phasicfft.m

利用ARMA、AR、MA模型

clear;clf;
tt = 0:1/2000:6;
yy = cos(2*pi*100.4*tt+pi/3)+cos(2*pi*150.6*tt+pi/2);
NFFT = 4096;
yy1 = yy(1:NFFT*2-1);
yy1 = yy1(:);
%vecter =  [1:NFFT,NFFT-1:-1:1];
%vecter = vecter/NFFT;
vecter =  conv(hanning(NFFT)',hanning(NFFT)');
vecter = vecter/NFFT/max(vecter);
for ii = 1:NFFT-1,
    yy2(ii) = vecter(ii)*yy1(ii) + vecter(NFFT+ii)*yy1(NFFT+ii);
end
yy2(NFFT) = vecter(NFFT)*yy1(NFFT);
yy2=[yy2(NFFT) yy2];
yy2_fft = fft(yy2,NFFT);
yy2_phase = phase(yy2_fft)*180/pi;
yy2_phase = mod(yy2_phase,360);
disp('中间点的测量值程序')
yy2_phase(206) 
yy2_phase(310)
disp('中间点的理论值')
tt=4095;
mod((2*pi*100.4*tt/2000+pi/3)*180/pi,360)
mod((2*pi*150.6*tt/2000+pi/2)*180/pi,360)