数字信号实验5谱分析.docxVIP

实验五一、实验名称：谱分析二、实验目的：1 研究不同类型的窗函数，研究一些不同的方法来测试窗的性能； 2 专注于有关窄带信号的几个不同的情形。三、实验原理：信号是无限长的，而在进行信号处理时只能采用有限长信号，所以需要将窄带信号“截断”。在信号处理中，“截断”被看成是用一个有限长的“窗口”看无限长的信倍号，或者从分析的角度是无限长的信号x(t)乘以有限长的窗函数w(t),由傅立叶变换性质可知： x(t)w(t)==1/2πX(jw)W(jw)如果x(t)是频带有限信号，而w(t)是频带无限函数，截断后的信号也必是频带无限信号，从而产生所谓的频谱泄露。频谱泄露是不可避免的，但是尽量减小，因此设计了不同的窗函数满足不同的要求。从能量的角度，频谱泄露也是能量泄露，因为加窗后，是原来的信号集中在宅频带内的能量分散到无限的频带范围。四、实验内容及结果分析：1和2、用matlab编程绘制各种窗函数的形状(1)矩形窗和汉宁窗形状及其幅频响应：N=60;n=0:N-1;w1=boxcar(N);[H1,W1]=dtft(w1,600);subplot(221),stem(n,w1);grid,title(矩形窗形状),xlabel(n),ylabel(w(n));subplot(222),plot(W1/2/pi,abs(H1));grid,title(矩形窗幅度响应),xlabel(f),ylabel(|H(w)|)w2=hanning(N);[H2,W2]=dtft(w2,600);subplot(223),stem(n,w2);grid,title(汉宁窗形状),xlabel(n),ylabel(w(n));subplot(224),plot(W2/2/pi,abs(H2));grid,title(汉宁窗幅度响应),xlabel(f),ylabel(|H(w)|)汉明窗和巴特利特窗形状及其幅频响应：N=60;n=0:N-1;w1=hamming(N);[H1,W1]=dtft(w1,600);subplot(221),stem(n,w1);grid,title(汉明窗形状),xlabel(n),ylabel(w(n));subplot(222),plot(W1/2/pi,abs(H1));grid,title(汉明窗幅度响应),xlabel(f),ylabel(|H(w)|)w2=bartlett(N);[H2,W2]=dtft(w2,600);subplot(223),stem(n,w2);grid,title(巴特利特窗形状),xlabel(n),ylabel(w(n));subplot(224),plot(W2/2/pi,abs(H2));grid,title(巴特利特窗幅度响应),xlabel(f),ylabel(|H(w)|)布莱克曼窗和triang窗形状及其幅频响应：N=60;n=0:N-1;w1=blackman(N);[H1,W1]=dtft(w1,600);subplot(221),stem(n,w1);grid,title(布莱克曼窗形状),xlabel(n),ylabel(w(n));subplot(222),plot(W1/2/pi,abs(H1));grid,title(布莱克曼窗幅度响应),xlabel(f),ylabel(|H(w)|)w2=triang(N);[H2,W2]=dtft(w2,600);subplot(223),stem(n,w2);grid,title(Triang窗形状),xlabel(n),ylabel(w(n));subplot(224),plot(W2/2/pi,abs(H2));grid,title(Triang窗幅度响应),xlabel(f),ylabel(|H(w)|)kaiser窗和切比雪夫窗形状及其幅频响应：N=60;n=0:N-1;w1=kaiser(N);[H1,W1]=dtft(w1,600);subplot(221),stem(n,w1);grid,title( Kaiser窗形状),xlabel(n),ylabel(w(n));subplot(222),plot(W1/2/pi,abs(H1));grid,title( Kaiser窗幅度响应),xlabel(f),ylabel(|H(w)|)w2=chebwin(N);[H2,W2]=dtft(w2,600);subplot(223),stem(n,w2);grid,title(切比雪夫窗形状),xlabel(n),ylabel(w(n));subplot(224),plot(W2/2/pi,abs(H2));grid,title(切比雪夫窗幅度响应),xlabel(f),y