经典滤波器的设计和使用

1、*实践教学*兰州理工大学计算机与通信学院2015年春季学期信号处理课程设计报告题 目： 经典滤波器的设计和使用 班 级： 通信工程12级（ 1 ）班 姓名： 设计质量（30分）： 学号： 122501xx 说明书质量（10分）： 同组成员： 指导教师： 兰州理工大学课程设计报告目录摘要1一、基本原理21.1iir数字滤波器设计原理21.2fir数字滤波器设计原理3二、滤波器详细设计52.1 信号源52.2 低通滤波器62.3 带通滤波器92.3高通滤波器11总结13参考文献14附录15摘要滤波器广泛应用于电子电路中，具有平滑信号、选频、去噪和使交流信号变为直流信号的功能。可分为低通滤波器、高通

2、滤波器、带通滤波器和带阻滤波器，也可分为模拟滤波器和数字滤波器，数字滤波器分为iir数字滤波器和fir数字滤波器。设计滤波器之前先根据实际需要确定滤波器类型和指标，通过matlab分别设计低通、带通和高通三种滤波器，首先通过三种滤波器对时域信号滤波，然后对信号频域变换后滤波，最后对信号进行时域搬移后滤波。通过matlab得出的频谱图对信号的三种形式通过三种滤波器后的结果进行分析，通过分析比较各滤波器的优点和缺点。关键词：matlab iir滤波器 fir滤波器 频域变换 时域搬移一、基本原理数字滤波器是一个离散系统，其系统函数一般可表示为的有理多项式形式，即当ai；i=1，2，n都为0时，上式

3、描述的系统称为有限脉冲响应数字滤波器，简称fir(finite-impulse response)数字滤波器。当系数ai；i=1，2，n中至少有一个是非0时，上式描述的系统称为无限脉冲响应数字滤波器，简称iir(infinite-impulse response)数字滤波器。对于iir数字滤波器，一般满足mn，这时将系统称为n阶iir数字滤波器。对于fir数字滤波器，系统函数中的有理多项式的最高次幂m就是其阶数。在设计数字滤波器时，首先根据具体的滤波成分，确定待设计数字滤波器的技术指标，再求解出满足设计指标的离散系统的系统函数h(z)。数字滤波器的技术指标通常用通带截频，通带最大衰减ap，阻带

4、截频，阻带最小衰减as表示。1.1 iir数字滤波器设计原理iir滤波器的设计就是根据给定的数字滤波器技术指标，确定上式中滤波器的阶数n和系数ai，bj。在满足技术指标的条件下，滤波器的阶数应尽可能低。应为滤波器阶数越低，实现滤波器的成本就越低。由于模拟滤波器的设计技术已经非常成熟，且可得闭合形式的解，因此在设计iir滤波器时，一般是通过模拟滤波器来设计滤波器。设计方法是先将数字滤波器技术指标转换为对应的模拟滤波器技术指标，然后设计满足技术指标的模拟滤波器h(s)，再将模拟滤波器h(s)转换为对应的数字滤波器h(z)。因此，在iir数字滤波器设计中，模拟滤波器设计是基础，模拟滤波器到数字滤波器

5、的转换是核心。而模拟滤波器的设计都是通过设计模拟低通滤波器来实现，比较常见的模拟低通滤波器有巴特沃思和切比雪夫等。将模拟滤波器变换为数字滤波器的主要方法有脉冲响应不变法和双线性变换法。数字滤波器h(z)模拟滤波器h(s)待设计数字滤波器指标模拟滤波器指标频率转换设计模拟滤波器s到z域转换图1-1 iir数字滤波器的设计过程1.2 fir数字滤波器设计原理表1-1 四种线性相位fir滤波器的特性类型iiiiiiiv阶数m偶数奇数偶数奇数hk的对称性偶对称偶对称奇对称奇对称0.5m0.5m0.5m0.5m000.50.5a()关于=0的对称性偶对称偶对称奇对称奇对称a()关于=的对称性偶对称奇对称

6、奇对称偶对称a()的周期2424a(0)任意任意00a()任意00任意可适用的滤波器类型lp，hp，bp，bs等lp，bp微分器，hilbert(希尔伯特)变换器微分器，hilbert(希尔伯特)变换器，hp窗函数法设计fir滤波器的基本思想是在时域逼近理想滤波器的单位脉冲相应。首先根据待逼近的理想滤波器的频率响应，由idtft求出理想滤波器的单位脉冲响应，再将无限长的加窗截断得到有限长序列。为了获得线性相位fir滤波器，在窗函数法设计fir滤波器的过程中，需要将线性相位因子加入理想滤波器的频率响应，具体步骤如下：（1） 根据所需设计的滤波器，确定线性相位滤波器的类型（i型，ii型，iii型，

7、iv型）。（2） 确定理想滤波器的幅度函数。（3） 确定理想滤波器的相位。对i型和ii型线性相位fir滤波器，对iii型和iv型线性相位fir滤波器。（4） 计算。（5） 截断得其中wk是长度n=m+1的矩形窗。二、滤波器详细设计2.1 信号源产生一个连续信号，包含低频、中频、高频分量，对其进行采样，进行频谱分析分别设计三种高通、低通、带通滤波器对信号进行滤波处理，观察滤波后信号的频谱。合成信号频谱程序如下:t=(1:100)/2000;x=sin(2*pi*t*20)+sin(2*pi*t*200)+sin(2*pi*t*500);plot(t,x);n=1:100;t=n/2000;x=f

8、ft(x,512);w=(0:255)/256*1000;x=sin(2*pi*t*20)+sin(2*pi*t*200)+sin(2*pi*t*500);plot(w,abs(x(1:256);xlabel('hz');ylabel('频率响应幅度');title('合成信号频谱图');图2-1 合成信号波形图图2-2 合成信号频谱图2.2 低通滤波器根据数字滤波器的设计原理，首先将数字域的指标转化为模拟域的指标设计模拟低通滤波器，然后应用双线性变换法将模拟滤波器转化为数字滤波器将指标转换成归一化模拟低通滤波器的指标，通过归一化的模拟低通滤波器

9、阶数n和3db截止频率的计算，将模拟域频率变换成模拟低滤波器h(s)，并用双线性变换法将h(s)转换成数字低通滤波器h(z)，由此得到低通数字滤波器。iir低通滤波器的设计程序为:ft=8000;fp=1000;fs=1200;as=100;ap=1;wp=2*pi* fp/ft;ws=2*pi*fs/ft;fp=2*fp*tan(wp/2);fs=2*fs*tan(ws/2);n11,wn11=buttord(wp,ws,1,50,'s');b11,a11=butter(n11,wn11,'s');num11,den11=bilinear(b11,a11,0.

10、5);h,w=freqz(num11,den11);plot(w*8000*0.5/pi,abs(h);title('低通滤波器');图2-3 iir低通滤波器读入信号对信号傅里叶变换进行频谱分析，利用filter函数滤波，对得到的信号傅里叶变换进行频谱分析，与滤波前进行比较。图2-4 滤波前信号波形图2-5滤波后信号波形图2-6滤波后信号的频谱图2-7滤波前信号的频谱2.3 带通滤波器iir带通滤波器的设计程序为：ft=8000;fp1=1200;fp2=3000;fs1=1000;fs2=3200;as=100;ap=1;wp1=tan(pi*fp1/ft);wp2=tan

11、(pi*fp2/ft);ws1=tan(pi*fs1/ft);ws2=tan(pi*fs2/ft);w=wp1*wp2/ws2;bw=wp2-wp1;wp=1;ws=(wp1*wp2-w.2)/(bw*w);n12,wn12=buttord(wp,ws,1, 50,'s');b12,a12=butter(n12,wn12,'s');num2,den2=lp2bp(b12,a12,sqrt(wp1*wp2),bw);num12,den12=bilinear(num2,den2,0.5);h,w=freqz(num12,den12);plot(w*8000*0.5/

12、pi,abs(h);title('带通滤波器');图2-8 iir带通滤波器图2-9滤波后的信号波形图2-10滤波后信号的频谱读入信号对信号傅里叶变换进行频谱分析，利用filter函数滤波，对得到的信号傅里叶变换进行频谱分析，与滤波前进行比较.并在滤波前后用sound函数播放语音信号，观察滤波前后的变化。2.3 高通滤波器iir高通滤波器的设计程序为：ft=8000;fp=4000;fs=3500;wp1=tan(pi*fp/ft);ws1=tan(pi*fs/ft);wp=1;ws=wp1*wp/ws1;n13,wn13=cheb1ord(wp,ws,1,50,'s&

13、#39;);b13,a13=cheby1(n13,1,wn13,'s');num,den=lp2hp(b13,a13,wn13);num13,den13=bilinear(num,den,0.5);h,w=freqz(num13,den13);plot(w*21000*0.5/pi,abs(h);title('高通滤波器');图2-11 iir高通滤波器图2-12滤波后的信号波形图2-13滤波后的信号频谱读入信号对信号傅里叶变换进行频谱分析，利用filter函数滤波，对得到的信号傅里叶变换进行频谱分析，与滤波前进行比较。分析滤波前后信号的变化，波形图变窄，频谱变

14、化表明:信号的高频率段被保留，低频率段被滤除。这说明滤波器设计基本符合指标要求。总结通过几周的课程设计，我对matlab的使用有了很大提高，能够熟练掌握matlab编程原理和应用matlab处理各种问题，matlab中带有海量功能函数和命令，在使用其设计滤波器时，首先要熟悉滤波器函数的使用方法，matlab应用范围广使用人数多，涌现出各类相关书籍、查询手册和网络资料。同时通过对滤波器的设计，提高了对数字信号处理的能力，此次设计的滤波器有低通滤波器、带通滤波器和高通滤波器，三种滤波器为iir数字滤波器或fir数字滤波器，在设计iir滤波器之前需要根据实际分析得出数字滤波器的指标，在设计数字滤波器

15、之前要设计模拟滤波器，所以需要将数字滤波器指标转换为模拟滤波器指标，设计好模拟滤波器后再通过脉冲响应不变法和双线性变换法变为iir数字滤波器，fir线性相位滤波器分为四种类型：i型、ii型、iii型、iv型。可以使用窗函数法设计线性相位fir数字滤波器，设计的滤波器会有gibbs现象，严重影响滤波器的性能，在设计时应减少这种现象的影响，常用窗函数有矩形窗、hann窗、hamming窗、blackman窗、kaiser窗。各个窗函数的特性都不一样，设计滤波器之前应熟悉这些窗函数的性能，以使设计出的滤波器性能更加完美，也可以使用频率取样法设计线性相位fir数字滤波器。iir滤波器的重大缺点在于，由

16、于存在反馈其稳定性不能得到保证。另外，反馈还使iir滤波器的数字运算可能溢出。从性能上来说，iir滤波器传输函数的极点可位于单位圆内的任何地方，因此可用较低的阶数获得高的选择性，所用的存贮单元少，所以经济而效率高。但是这个高效率是以相位的非线性为代价的。选择性越好，则相位非线性越严重。相反，fir滤波器却可以得到严格的线性相位，然而由于fir滤波器传输函数的极点固定在原点，所以只能用较高的阶数达到高的选择性。在设计之初由于对matlab不熟悉和专业知识掌握不牢，设计相对艰难，进度缓慢，通过老师的指导、同学们的帮助和查阅资料，滤波器的设计才得以顺利进行，提高了独自处理问题的能力和团队合作能力，在

17、此对指导老师和同学们表示感谢。参考文献1 丁玉美. 数字信号处理. 西安电子科技大学出版社，20122 陈后金. 数字信号处理(第2版）. 高等教育出版社, 20123 a.v.oppenheim. 离散时间信号处理(第3版)英文版 ,电子工业出版社, 20114 v.k. ingle. 数字信号处理: 应用matlab(第3版)英文影印版,科学出版社, 20125 胡广书.数字信号处理.清华大学出版社(第3版), 2012.6 樊昌信. 通信原理. 北京:国防工业出版社,2010.7 刘敏. matlab 通信仿真与应用. 北京:国防工业出版社.附录iir低通滤波器滤波程序：f1=20;f2

18、=200;f3=500;t=(1:100)/2000;x1=sin(2*pi*t*f1);x1=sin(2*pi*t*f2);x1=sin(2*pi*t*f3);x=sin(2*pi*t*f1)+sin(2*pi*t*f2)+sin(2*pi*t*f3);n=1:100;t=n/2000;x=fft(x,512);w=(0:255)/256*1000;x=sin(2*pi*t*f1)+sin(2*pi*t*f2)+sin(2*pi*t*f3);ft=8000;fp=1000;fs=1200;as=100;ap=1;wp=2*pi*fp/ft;ws=2*pi*fs/ft;fp=2*fs*tan(

19、wp/2);fs=2*fs*tan(ws/2);n11,wn11=buttord(wp,ws,1,50,'s');b11,a11=butter(n11,wn11,'s');num11,den11=bilinear(b11,a11,0.5);h,w=freqz(num11,den11);z11=filter(num11,den11,x);m11=fft(z11);plot(abs(x),'g');title('滤波前信号的频谱');plot(abs(m11),'r');title('滤波后信号的频谱'

20、);plot(x);title('滤波前信号的波形');plot(z11);title('滤波后信号的波形');iir带通滤波器滤波程序：f1=20;f2=200;f3=500;t=(1:100)/2000;x1=sin(2*pi*t*f1);x2=sin(2*pi*t*f2);x3=sin(2*pi*t*f3);x=sin(2*pi*t*f1)+sin(2*pi*t*f2)+sin(2*pi*t*f3);n=1:100;t=n/2000;x=fft(x,512);w=(0:255)/256*1000;x=sin(2*pi*t*f1)+sin(2*pi*t*f2

21、)+sin(2*pi*t*f3);ft=8000;fp1=1200;fp2=3000;fs1=1000;fs2=3200;as=100;ap=1;wp1=tan(pi*fp1/ft);wp2=tan(pi*fp2/ft);ws1=tan(pi*fs1/ft);ws2=tan(pi*fs2/ft);w=wp1*wp2/ws2;bw=wp2-wp1;wp=1;ws=(wp1*wp2-w.2)/(bw*w);n12,wn12=buttord(wp,ws,1,50,'s');b12,a12=butter(n12,wn12,'s');num2,den2=lp2bp(b12,a12,sqrt(wp1*wp2),bw);num12,den12=bilinear(num2,den2,0.5);h,w=freqz(num12,den12);z12=filter(num12,den12,x);m12=fft(z12);plot(abs(x),'g');title('滤波前信号的频谱');plot(abs(m12),'r');title('滤波后信号的频谱');plot(x);title('滤波前信号的波形');plot(z12);title('