DSP试验 04008012
figure(2)
subplot(2,1,1); plot(n,yn_1);
xlabel('n');ylabel('y(n)'); subplot(2,1,2);
plot(n,abs(fft(yn_1))); xlabel('k');ylabel('Y(k)'); axis([0,600,0,300]); 运行结果：
11
DSP试验 04008012
（7）用FFT分别计算xa(n)(p=8,q=2)和 xb(n)（a＝0.1，f＝0.0625）的16点循环相关和线性相关，问一共有多少种结果，它们之间有何共同点？ n=0:1:15;
xan=exp(-(n-8).^2/2);
xbn=exp(-0.1*n).*sin(2*pi*0.0625*n); k=length(xbn);
xan1=[xan zeros(1,k-1)]; xbn1=[xbn zeros(1,k-1)]; xak=fft(xan1); xbk=fft(xbn1);
rm=real(ifft(conj(xak).*xbk)); rm1=[rm(k+1:2*k-1) rm(1:k)]; m=(-k+1):(k-1); subplot(2,1,1);
stem(m,rm1);xlabel('n'); subplot(2,1,2);
12
DSP试验 04008012
xak=fft(xan); xbk=fft(xbn);
rm=real(ifft(conj(xak).*xbk)); stem(n,rm);xlabel('n'); 运行结果：
（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); subplot(2,1,1);
13
DSP试验 04008012
stem(m,rm);xlabel('m'); 运行结果：
四，思考题
（1）实验中的信号序列xc(n)和xd(n)，在单位圆上的Z变换频谱|Xc(jω)|和|Xd(jω)|会相同吗？如果不同，你能说出哪一个低频分量更多一些吗？为什么？
答：不同。反三角波的低频分量更多一些。
（2）对一个有限长序列进行DFT等价于将该序列周期延拓后进行DFS展开，因为DFS也只是取其中一个周期来计算，所以FFT在一定条件下也可以用以分析周期信号序列。如果实正弦信号sin(2πfn),f=0.1用16点FFT来做DFS运算，得到的频谱时信号本身的真实谱吗？
14