资源描述:
《哈工程数字信号处理第五次》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
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