R;(t,t +T) = Aj(t) Aj(t +T) / (w) eiwt dw i j = 1,2,3 Is that formulation can be written as matlab code by using FFT? t= time series tau= time lag w= angular frequency