实验三 窗函数的特性分析

1、 本科学生实验报告学号 * 姓名 * 学院 物电学院 专业、班级 *实验课程名称 数字信号分析与处理 教师及职称 * 开课学期 2015 至 2016学年 下 学期填报时间 2016 年 3 月 25 日云南师范大学教务处编印一、验设计方案实验序号实验三实验名称窗函数的特性分析实验时间2016/3/25实验室同析楼三栋313实验室1实验目的 分析各种窗函数的时域和频域特性，灵活应用窗函数分析信号频谱和设计FIR数字滤波器。2 实验原理、实验流程或装置示意图 在确定信号谱分析、随机信号功率谱估计以及FIR数字滤波器设计中，窗函数的选择对频谱分析和滤波器设计都起着重要的作用。在确定信号谱分析和随机

2、信号功率谱估计中，截短无穷长的序列会造成频率泄露，影响频率普分析的精确度和质量。合理选取窗函数的类型，可以改善泄露现象。在FIR数字滤波器设计中，截短无穷长的系统单位脉冲序列会造成FIR滤波器的幅度特性产生波动，且出出现过渡带。 【例1.3.1】 写出分析长度N=51点矩形窗的时域波行和频谱的MATLAB程序。 解 N=51;w=boxcar(N);W=fft(w,256);subplot(2,1,1);stem(0:N-1,w);subplot(2,1,2);plot(-128:127,abs(fftshift(W);运算结果如图1.3.1所示 图1.3.1 矩形窗的时域波形和频谱3实验设备

3、及材料 计算机，MATLAB软件4实验方法步骤及注意事项注意事项：（1） 在使用MATLAB时应注意中英输入法的切换，在中文输入法输入程序时得到的程序是错误的；（2） MATLAB中两个信号相乘表示为x.*u,中间有个.，同样两个信号相除也是如此；（3） 使用MATLAB编写程序时，应新建一个m文件，而不是直接在Comandante窗口下编写程序；（4） 在使用编程时，应该养成良好的编写习惯。 5实验数据处理方法 图像法6参考文献信号分析与处理MATLAB数值计算与方法二、报告1实验现象与结果实验内容 1.分析并绘出常用窗函数的时域波形特性。 N=51;w=boxcar(N);subplot(

4、3,1,1);stem(0:N-1,w);title('矩形窗的时域波形')w=hanning(N);subplot(3,1,2);stem(0:N-1,w);title('汉宁窗的时域波形')w=hamming(N);subplot(3,1,3);stem(0:N-1,w);title('汉明窗的时域波形') 实验运行结果如图 3.1 图3.1 矩形窗、汉宁窗及汉明窗的时域特性波形N=51;w=blackman(N);subplot(3,1,1);stem(0:N-1,w);title('布莱克曼窗的时域波形');w=bartl

5、ett(N);subplot(3,1,2);stem(0:N-1,w);title('Bartlett窗的时域波形');beta=2*N;w=Kaiser(N,beta);subplot(3,1,3);stem(0:N-1,w);title('凯泽窗的时域波形');实验结果如图3.2所示 图3.2 布莱克曼窗、Bartlett窗、凯泽窗时域特性波形3. 研究凯泽窗（Kaiser）的参数选择对其时域和频域的影响。 （1）固定beta=4，分别取N=20,60,110N=20;beta=4;w=Kaiser(N,beta);subplot(3,2,1);stem(0

6、:N-1,w);title('第3题凯泽窗N=20时域波形');W=fft(w,256);subplot(3,2,2);plot(-128:127,abs(fftshift(W);title('第3题凯泽窗N=20频域波形 ')N=60;w=Kaiser(N,beta);subplot(3,2,3);stem(0:N-1,w);title('第三题凯泽窗N=60波形');W=fft(w,256);subplot(3,2,4);plot(-128:127,abs(fftshift(W);title('第3题凯泽窗N=60频域波形 '

7、);N=110;w=Kaiser(N,beta);subplot(3,2,5);stem(0:N-1,w);title('第三题凯泽窗N=110波形');W=fft(w,256);subplot(3,2,6);plot(-128:127,abs(fftshift(W);title('第3题凯泽窗N=110频域波形 ')实验结果如图3.3 图3.3 N取不同值的时的时域和频域波形 （2）固定N=60，分别取beta=1,5,11.N=60beta=1;w=Kaiser(N,beta);subplot(3,2,1);stem(0:N-1,w);title('

8、第3题凯泽窗beta=1时域波形');W=fft(w,256);subplot(3,2,2);plot(-128:127,abs(fftshift(W);title('第3题凯泽窗beta=1频域波形 ')beta=5;w=Kaiser(N,beta);subplot(3,2,3);stem(0:N-1,w);title('第三题凯泽窗beta=5波形');W=fft(w,256);subplot(3,2,4);plot(-128:127,abs(fftshift(W);title('第3题凯泽窗beta=5频域波形 '); beta=1

9、1;w=Kaiser(N,beta);subplot(3,2,5);stem(0:N-1,w);title('第三题凯泽窗beta=11波形');W=fft(w,256);subplot(3,2,6);plot(-128:127,abs(fftshift(W);title('第3题凯泽窗beta=11频域波形 ');实验运行结果如图3.4 图3.4 beta取不同值的时的时域和频域波形4.某序列为，使用fft函数分析其频谱。 （1）利用不同宽度的N的矩形窗截短该序列，N分别为20,40,160，观察不同长度N的窗对谱分析结果的影响。实验matlab程序代码N=2

10、0;k=0:N-1;w=0.5*cos(11*pi*k)/N)+cos(9*pi*k)/N);W=fft(w,256);subplot(3,1,1);plot(-128:127,abs(fftshift(W);title('第4题 xkN=20频谱 ')subplot(3,1,2);N=40;k=0:N-1;w=0.5*cos(11*pi*k)/N)+cos(9*pi*k)/N);W=fft(w,256);plot(-128:127,abs(fftshift(W);title('第4题 xkN=40频谱 ')subplot(3,1,3);N=160;k=0:N-

11、1;w=0.5*cos(11*pi*k)/N)+cos(9*pi*k)/N);W=fft(w,256);plot(-128:127,abs(fftshift(W);title('第4题 xkN=160频谱 ')实验运行结果如图3.5 图3.5 不同宽度N的矩形窗对谱分析结果影响 （2）利用汉明窗重做（1）。 实验程序代码N=20;k=0:N-1;w=0.5*cos(11*pi*k)/N)+cos(9*pi*k)/N).*(0.54-0.46.*cos(2*pi.*k/(N-1);W=fft(w,256);subplot(3,1,1);plot(-128:127,abs(ffts

12、hift(W),'r');title('第4题 第二问xkN=20频谱 ')subplot(3,1,2);N=40;k=0:N-1;w=0.5*cos(11*pi*k)/N)+cos(9*pi*k)/N).*(0.54-0.46.*cos(2*pi.*k/(N-1);W=fft(w,256);plot(-128:127,abs(fftshift(W),'r');title('第4题 第二问xkN=40频谱 ')subplot(3,1,3);N=160;k=0:N-1;w=0.5*cos(11*pi*k)/N)+cos(9*pi*k

13、)/N).*(0.54-0.46.*cos(2*pi.*k/(N-1);W=fft(w,256);plot(-128:127,abs(fftshift(W),'r');title('第4题 第二问xkN=160频谱 ') 实验运行结果如图3.6 图3.6 不同宽度N的汉明窗对谱分析结果影响3 实验总结由实验结果可以看出矩形窗波形为方形，汉宁窗，汉明窗，布莱克曼窗，凯泽窗波形为正弦波形，Bartlett窗波形为三角形矩形窗，汉宁窗，汉明窗，布莱克曼窗，Bartlett窗的波形固定，一旦选择了某种窗函数，用它进行谱分析得到的频谱纹波或设计出的滤波器的阻带衰减是确定的。凯泽窗是一种可调窗，可以通过改变窗函数的形状来控制频谱纹波或阻带衰减指标，因而获