资源描述:
《哈工程数字信号处理实验五》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、实验五谱分析一实验目的1研究不同类型的窗函数,研究一些不同的方法来测试窗的性能;2专注于有关窄带信号的几个不同的情形。二实验原理信号是无限长的,而在进行信号处理时只能采用有限长信号,所以需要将窄带信号“截断”。在信号处理中,“截断”被看成是用一个有限长的“窗口”看无限长的信倍号,或者从分析的角度是无限长的信号x(t)乘以有限长的窗函数w(t),由傅立叶变换性质可知:x(t)w(t)==1/2πX(jw)W(jw)如果x(t)是频带有限信号,而w(t)是频带无限函数,截断后的信号也必是频带无限信号,从而产生所谓的频谱泄露。频谱泄露是不可避免的,但是尽量减小,因此设计了不同的窗函数满足不
2、同的要求。从能量的角度,频谱泄露也是能量泄露,因为加窗后,是原来的信号集中在宅频带内的能量分散到无限的频带范围。三实验内容、程序及结果㈠用MATLAB编程绘制各种窗函数的形状及其幅度响应。(1)矩形窗;N=88;w=boxcar(N);n=[0:N-1];[H,W]=dtft(w,1024);subplot(211);stem(n,w),title('矩形窗形状');xlabel('n'),ylabel('w(n)');subplot(212);plot(W/2/pi,abs(H)),title('矩形窗幅度响应');xlabel('f'),ylabel('
3、W(w)
4、')⑵汉宁窗;
5、N=88;w=hanning(N);n=[0:N-1];[H,W]=dtft(w,1024);subplot(211);stem(n,w),title('汉宁窗形状');xlabel('n'),ylabel('w(n)');subplot(212);plot(W/2/pi,abs(H)),title('汉宁窗幅度响应');xlabel('f'),ylabel('
6、W(w)
7、')⑶汉明窗;N=88;w=hamming(N);n=[0:N-1];[H,W]=dtft(w,1024);subplot(211);stem(n,w),title('汉明窗形状');xlabel('n'),yla
8、bel('w(n)');subplot(212);plot(W/2/pi,abs(H)),title('汉明窗幅度响应');xlabel('f'),ylabel('
9、W(w)
10、')⑷巴特利特窗;N=88;w=bartlett(N);n=[0:N-1];[H,W]=dtft(w,1024);subplot(211);stem(n,w),title('巴特利特窗形状');xlabel('n'),ylabel('w(n)');subplot(212);plot(W/2/pi,abs(H)),title('巴特利特窗幅度响应');xlabel('f'),ylabel('
11、W(w)
12、')⑸布莱
13、克曼窗;N=88;w=blackman(N);n=[0:N-1];[H,W]=dtft(w,1024);subplot(211);stem(n,w),title('布莱克曼窗形状');xlabel('n'),ylabel('w(n)');subplot(212);plot(W/2/pi,abs(H)),title('布莱克曼窗幅度响应');xlabel('f'),ylabel('
14、W(w)
15、')⑹Triang窗;N=88;w=triang(N);n=[0:N-1];[H,W]=dtft(w,1024);subplot(211);stem(n,w),title('Triang窗形状')
16、;xlabel('n'),ylabel('w(n)');subplot(212);plot(W/2/pi,abs(H)),title('Triang窗幅度响应');xlabel('f'),ylabel('
17、W(w)
18、')⑺Kaiser窗;N=88;w=kaiser(N,2);n=[0:N-1];[H,W]=dtft(w,1024);subplot(211);stem(n,w),title(‘Kaiser窗形状');xlabel('n'),ylabel('w(n)');subplot(212);plot(W/2/pi,abs(H)),title('Kaiser窗幅度响应');xlabe
19、l('f'),ylabel('
20、W(w)
21、')⑻切比雪夫窗;N=88;w=chebwin(N,16);n=[0:N-1];[H,W]=dtft(w,1024);subplot(211);stem(n,w),title('切比雪夫窗形状');xlabel('n'),ylabel('w(n)');subplot(212);plot(W/2/pi,abs(H)),title('切比雪夫窗幅度响应');xlabel('f'),ylabel('
22、W(w)
23、')3,绘