学生实验报告
课程名称 数字信号处理与应用 实验日期 2012 年 12 月 26 日
注:学生做每个实验都必须填写实验报告。实验报告可以是纸质的,也可以是电子形式的。凡电子文本形式的实验报告须在学期末刻录进该学期教学资料光盘并存个人教学档案。
实验指导教师签名 年 月 日
实验四 用窗函数法设计 FIR数字滤波器
一、 实验目的
(1) 掌握用窗函数法设计FIR数字滤波器的原理和方法。
(2) 熟悉线性相位FIR数字滤波器特性。
(3) 了解各种窗函数对滤波特性的影响。
二、 实验原理、
滤波器的理想频率响应函数为Hd(ejω),则其对应的单位脉冲响应为
hd(n) =
窗函数设计法的基本原理是用有限长单位脉冲响应序列h(n)逼hd(n)。由于hd(n)往往是无限长序列,且是非因果的,所以用窗函数。w(n)将hd(n)截断,并进行加权处理:
h(n) = hd(n) w(n)
h(n)就作为实际设计的FIR数字滤波器的单位脉冲响应序列,其频率响应函数H(ejω)为
H(ejω) =
用窗函数法设计的滤波器性能取决于窗函数w(n)的类型及窗口长度N的取值。设计过程中,要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度N。
三、 实验内容
1.MATALAB程序
任务一
N1=15;
N2=33;
b1=fir1(14,1/4,hanning(15));
b2=fir1(32,1/4,hanning(33));
[H1,W]=freqz(b1,1);
H1_db=20*log10(abs(H1));
magH1=abs(H1);
phaH1=angle(H1);
[H2,W]=freqz(b2,1);
H2_db=20*log10(abs(H2));
magH2=abs(H2);
phaH2=angle(H2);
figure(1);
subplot(2,1,1);stem(b1);
title('N=15时,汉宁窗得到的FIR滤波器脉冲响应');
subplot(2,1,2);stem(b2);
title('N=33时,汉宁窗得到的FIR滤波器脉冲响应');
figure(2);
subplot(3,1,1);plot(W/pi,magH1);
title('N=15时汉宁窗得到的FIR滤波器的幅频特性')
subplot(3,1,2);plot(W/pi,H1_db);
title('N=15时汉宁窗得到的FIR滤波器的对数幅频特性')
subplot(3,1,3);plot(W/pi,phaH1);
title('N=15时汉宁窗得到的FIR滤波器的相频特性')
figure(3);
subplot(3,1,1);plot(W/pi,magH2);
title('N=33时汉宁窗得到的FIR滤波器的幅频特性')
subplot(3,1,2);plot(W/pi,H2_db);
title('N=33时汉宁窗得到的FIR滤波器的对数幅频特性')
subplot(3,1,3);plot(W/pi,phaH2);
title('N=33时汉宁窗得到的FIR滤波器的相频特性')
任务二
N=33;
b1=fir1(32,1/4,boxcar(33));
b2=fir1(32,1/4,hanning(33));
b3=fir1(32,1/4,bartlett(33));
b4=fir1(32,1/4,blackman(33));
[H1,W]=freqz(b1,1);
H1_db=20*log10(abs(H1));
magH1=abs(H1);
phaH1=angle(H1);
[H2,W]=freqz(b2,1);
H2_db=20*log10(abs(H2));
magH2=abs(H2);
phaH2=angle(H2);
[H3,W]=freqz(b3,1);
H3_db=20*log10(abs(H3));
magH3=abs(H3);
phaH3=angle(H3);
[H4,W]=freqz(b4,1);
H4_db=20*log10(abs(H4));
magH4=abs(H4);
phaH4=angle(H4);
figure(1);
subplot(2,2,1);stem(b1);
title('矩形窗得到的FIR滤波器脉冲响应')
subplot(2,2,2);stem(b2);
title('汉宁窗得到的FIR滤波器脉冲响应')
subplot(2,2,3);stem(b3);
title('三角窗得到的FIR滤波器脉冲响应')
subplot(2,2,4);stem(b4);
title('布拉克曼窗得到的FIR滤波器脉冲响应')
figure(2);
subplot(2,2,1);plot(W/pi,magH1);
title('矩形窗得到的FIR滤波器幅频特性')
subplot(2,2,2);plot(W/pi,magH2);
title('汉宁窗得到的FIR滤波器幅频特性')
subplot(2,2,3);plot(W/pi,magH3);
title('三角窗得到的FIR滤波器幅频特性')
subplot(2,2,4);plot(W/pi,magH4);
title('布拉克曼窗得到的FIR滤波器幅频特性')
figure(3);
subplot(2,2,1);plot(W/pi,H1_db);
title('矩形窗得到的FIR滤波器相频特性')
subplot(2,2,2);plot(W/pi,H2_db);
title('汉宁得到的FIR滤波器相频特性')
subplot(2,2,3);plot(W/pi,H3_db);
title('三角窗得到的FIR滤波器相频特性')
subplot(2,2,4);plot(W/pi,H4_db);
title('布拉克曼得到的FIR滤波器相频特性')
figure(4);
subplot(2,2,1);plot(W/pi,phaH1);
title('矩形窗得到的FIR滤波器相频特性')
subplot(2,2,2);plot(W/pi,phaH2);
title('汉宁窗得到的FIR滤波器相频特性')
subplot(2,2,3);plot(W/pi,phaH3);
title('三角窗得到的FIR滤波器相频特性')
subplot(2,2,4);plot(W/pi,phaH4);
title('布拉克曼窗得到的FIR滤波器相频特性')
2.实验波形图
任务一
任务二
四、 实验结论
1.N的大小决定了窗谱的主瓣宽度,N越大,窗谱的主瓣宽度越大
2.最小阻带衰减只有窗行决定,不受N的影响,过渡带宽度与N和窗形都有关,N越大,过渡带宽越小
3.由实验可知滤波特性:布拉克曼窗>汉宁窗>三角窗>矩形窗
五、 思考题
(1) 如果给定通带截止频率和阻带截止频率以及阻带最小衰减,如何用窗函数法设计线性相位低通滤波器?写出设计步骤。
答:①根据计算出所要设计的滤波器的过渡带
②查表找到满足的阻带最小衰减所需的窗函数
③利用已选的窗函数的过渡带宽除以得到低通滤波器单位脉冲响应的长度N
④根据求出h(n)
⑤再根据H(ejω) =求出响应函数H(ejω)
(2) 如果要求用窗函数法设计带通滤波器,且给定上、下边带截止频率为ω1和ω2,试求理想带通的单位脉冲响应hd(n)。
答:由算出,求出滤波器过渡带宽和滤波器长度N,最后由或求出,其中
实验四IIR数字滤波器设计及软件实现实验报告实验四IIR数字滤波器设计及软件实现实验报告一实验目的1熟悉用双线性变换法设计IIR数…
电子信息工程系学生实验报告实验课程名称数字信号处理实验项目名称IIR数字滤波器设计及软件实现实验时间20xx年06月01日班级座号…
西北工业大学matlab实验报告实验三IIR数字滤波器设计实验报告一实验目的1通过仿真冲激响应不变法和双线性变换法2掌握滤波器性能…
数字信号处理实验报告数字信号处理实验报告实验四IIR数字滤波器的设计学班学指导教生姓名级号师张志翔电子信息工程1203班12401…
实验报告姓名学号班级组号李鹏博20xx300704030411021实验名称课程名称实验室名称实验日期IIR数字滤波器设计数字信号…
重庆交通大学综合性设计性实验报告班级学号姓名实验项目名称滤波器设计综合实验应用实验项目性质综合性设计实验所属课程数字信号处理实验室…
数字信号处理实验报告姓名班级学号指导教师实验目的1熟悉用双线性变换法设计IIR数字滤波器原理和试验方法2掌握数字滤波器的计算机仿真…
实验六IIR数字滤波器的设计一实验目的从理论上讲任何的线性是不变LTI离散时间系统都可以看做一个数字滤波器因此设计数字滤波器实际就…
物理与电子信息工程学院实验报告实验课程名称数字信号处理实验名称FIR数字滤波器设计与软件实现班级姓名学号101234153成绩实验…
实验二数字图像的空间域滤波平滑滤波1实验目的1掌握图像滤波的基本定义及目的2理解空间域滤波的基本原理及方法3掌握进行图像的空域滤波…