哈工程数字信号处理第五次

《哈工程数字信号处理第五次》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、实验五谱分析一实验目的1研究不同类型的窗函数，研究一些不同的方法来测试窗的性能；2专注于有关窄带信号的几个不同的情形。二实验原理信号是无限长的，而在进行信号处理时只能采用有限长信号，所以需要将窄带信号“截断”。在信号处理中，“截断”被看成是用一个有限长的“窗口”看无限长的信倍号，或者从分析的角度是无限长的信号x(t)乘以有限长的窗函数w(t),由傅立叶变换性质可知：x(t)w(t)=1/2πX(jw)W(jw)如果x(t)是频带有限信号，而w(t)是频带无限函数，截断后的信号也必是频带无限信号，从而产生所谓的频谱泄露。频谱泄露是不可避免的，但是尽量减

2、小，因此设计了不同的窗函数满足不同的要求。从能量的角度，频谱泄露也是能量泄露，因为加窗后，是原来的信号集中在宅频带内的能量分散到无限的频带范围。三实验内容、程序及结果㈠用MATLAB编程绘制各种窗函数的形状及其幅度响应。w1=boxcar(30);subplot(421)stem(w1)title('juxingchang')xlabel('t');ylabel('w1');w2=hanning(30);subplot(422)stem(w2)title('hanningchuang')xlabel('t');ylabel('w2');w3=hamm

3、ing(30);subplot(423)stem(w3)title('hanmingchuang');xlabel('t');ylabel('w3');w4=bartlett(30);subplot(424)stem(w4)title('batelitechuang')xlabel('t');ylabel('w4');w5=blackman(30);subplot(425)stem(w5)title('balaikemanchuang')xlabel('t');ylabel('w5');w6=triang(30);subplot(426)stem(w

4、6)title('triang')xlabel('t');ylabel('w6');w7=kaiser(30,60);subplot(427)stem(w7)title('kaiser')xlabel('t');ylabel('w7');w8=chebwin(30,60);subplot(428)stem(w8)title('chebwin')xlabel('t');ylabel('w8');2,各种窗函数的幅频响应w1=boxcar(30);[X,W]=dtft(w1,300)subplot(421)plot(W/2/pi,abs(X));titl

5、e('juxingchuang');xlabel('w');ylabel('

6、W(jw)

7、');w2=hanning(30);[X,W]=dtft(w2,300)subplot(422)plot(W/2/pi,abs(X));title('hanningchuang');xlabel('w');ylabel('

8、W(jw)

9、');w3=hamming(30);[X,W]=dtft(w3,300)subplot(423)plot(W/2/pi,abs(X));title('hanmingchuang');xlabel('w');ylabel('

10、W(j

11、w)

12、');w4=bartlett(30);[X,W]=dtft(w4,300)subplot(424)plot(W/2/pi,abs(X));title('batelite');xlabel('w');ylabel('

13、W(jw)

14、');w5=blackman(30);[X,W]=dtft(w5,300)subplot(425)plot(W/2/pi,abs(X));title('bulaikemanchuang');xlabel('w');ylabel('

15、W(jw)

16、');w6=triang(30);[X,W]=dtft(w6,300)subp

17、lot(426)plot(W/2/pi,abs(X));title('triang');xlabel('w');ylabel('

18、W(jw)

19、');w7=kaiser(30,60);[X,W]=dtft(w7,300)subplot(427)plot(W/2/pi,abs(X));title('kaiser');xlabel('w');ylabel('

20、W(jw)

21、');w8=chebwin(30,60);[X,W]=dtft(w8,300)subplot(428)plot(W/2/pi,abs(X));title('chebwin');xlabel

22、('w');ylabel('

23、W(jw)

24、');3矩形窗的幅频响应w1=boxcar(10);[X,W]=d