《数字信号处理—理论与实践》课件第5章.pptVIP

即N≥100，因此取N=100。程序清单如下：N=100;n=0:9;xn=0.5.^n;Xk=fft(xn,N);k=0:N－1;wk=2*k/N;subplot(1,2,1);plot(wk,abs(Xk));xlabel(′ω/π′);ylabel(′幅度′);

title(′100点DFT的幅频特性′);subplot(1,2,2);plot(wk,angle(Xk));xlabel(′ω/π′);

ylabel(′相位′);

title(′100点DFT的相频特性′);

运行结果如图5-18所示。图5-18例5-12的频谱特性从图5-18可见，x(n)的频谱变化很慢，所以在用DFT对

该信号作谱分析时，可以适当降低频率分辨率的要求，从而使N减小一些。比如当N取16时，其幅频特性和相频特性如图5-19所示，可以看出其与N=100的情况相差很小。图5-19例5-12的16点DFT的频率特性2.对无限长时域离散信号进行频谱分析

如果是对无限长序列作谱分析，则需要将序列截取为有限长的一段，再用DFT作频谱分析。但这种截取必然会带来失真，因此得到的是近似频谱，其性能也与截取时所用的窗函数有关。

例5-13假设x(n)=cos(πn/10)+sin(πn/6)+cos(2πn/5)，试用DFT对其进行频谱分析。要求分别使用矩形窗和海明窗进行截取，并假设有三种不同的截取长度，即N=16,64和128。解程序清单如下：

N1=16;N2=64;N3=128;

%矩形窗截取

n=0:N1－1;

xn1=cos(pi*n/10)+sin(pi*n/6)+cos(2*pi*n/5);

Xk1=fft(xn1,2048);%2048点DFT作为连续频谱的近似

wk=2*(0:2047)/2048;

subplot(3,2,1);plot(wk,abs(Xk1));xlabel(′ω/π′);

ylabel(′幅度′);

title(′矩形窗截取（长度为16）′);n=0:N2－1;

xn2=cos(pi*n/10)+sin(pi*n/6)+cos(2*pi*n/5);

Xk2=fft(xn2,2048);%2048点DFT作为连续频谱的近似

subplot(3,2,3);plot(wk,abs(Xk2));xlabel(′ω/π′);ylabel(′幅度′);

title(′矩形窗截取（长度为64）′);

n=0:N3-1;

xn3=cos(pi*n/10)+sin(pi*n/6)+cos(2*pi*n/5);

Xk3=fft(xn3,2048);%2048点DFT作为连续频谱的近似subplot(3,2,5);plot(wk,abs(Xk3));xlabel(′ω/π′);

ylabel(′幅度′);

title(′矩形窗截取（长度为128）′);

%海明窗截取

wn=hamming(N1);

xn1=xn1.*wn′;

Xk1=fft(xn1,2048);%2048点DFT作为连续频谱的近似

wk=2*(0:2047)/2048;

subplot(3,2,2);plot(wk,abs(Xk1));xlabel(′ω/π′);

ylabel(′幅度′);

title(′海明窗截取（长度为16）′);wn=hamming(N2);

xn2=xn2.*wn′;

Xk2=fft(xn2,2048);%2048点DFT作为连续频谱的近似

subplot(3,2,4);plot(wk,abs(Xk2));xlabel(′ω/π′);ylabel(′幅度′);

title(′海明窗截取（长度为64）′);

wn=hamming(N3);

xn3=xn3.*wn′;

Xk3=fft(xn3,2048);%2048点DFT作为连续频谱的近似

subplot(3,2,6);plot(wk,abs(Xk3));xlabel(′ω/π′