loading
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('幅度');