《現代数字信号处理》仿真题(3-4章).docVIP

3.17 (1) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %FFT法与直接法求自相关函数的比较 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% clc; f1=0.15; f2=0.17; f3=0.26; n=1:32; d=1; %白噪声的方差取1 N=length(n); A1=sqrt((10^3)*2); A2=A1; A3=sqrt((10^2.7)*2); s=A1*cos(2*pi*f1*n+unifrnd(0,2*pi))+A2*cos(2*pi*f2*n+unifrnd(0,2*pi))+A3*cos(2*pi*f3*n+unifrnd(0,2*pi))+normrnd(0,1,1,N); %直接法求自相关函数 r1=zeros(1,N); k=1; for k1=0:N-1 for k2=k1+1:N r1(k)=r1(k)+s(k2)*(s(k2-k1)); end r1(k)=r1(k)/N; k=k+1; end %FFT法求自相关函数 s=[s,zeros(1,N)]; r2=fft(s,2*N); r2=abs(r2).^2/N; r3=ifft(r