哈工程数字信号处理实验五.docVIP

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