发布时间 : 星期三 文章实验二快速傅里叶变换(FFT)及其应用更新完毕开始阅读8bedea61cf2f0066f5335a8102d276a2012960f8
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('幅度');