实验二快速傅里叶变换(FFT)及其应用

2）求循环相关 >> n=0:1:15;

>> xan=exp(-(n-8).^2/2);

>> xbn=exp(-0.1*n).*sin(2*pi*0.0625*n); >> k=length(xbn); >> xak=fft(xan); >> xbk=fft(xbn);

>> rm=real(ifft(conj(xak).*xbk));

>> stem(n,rm);xlabel('n');ylabel('幅度');legend('循环相关');

（8）用FFT分别计算xa(n)（p?8,q?2）和xb(n)（a?0.1,f?0.0625）的自相关函数。

解答： >> n=0:1:15;

>> xan=exp(-(n-8).^2/2); >> k=length(xan); >> xak=fft(xan,2*k);

>> rm=real(ifft(conj(xak).*xak)); >> rm=[rm(k+2:2*k) rm(1:k)]; >> m=(-k+1):(k-1);

>> stem(m,rm);xlabel('m');ylabel('幅度');

（2） >> n=0:1:15;

>> xbn=exp(-0.1*n).*sin(2*pi*0.0625*n); >> k=length(xbn); >> xbk=fft(xbn,2*k);

>> rm=real(ifft(conj(xbk).*xbk)); >> rm=[rm(k+2:2*k) rm(1:k)]; >> m=(-k+1):(k-1);

>> stem(m,rm);xlabel('m');ylabel('幅度');