2022年matlab滤波器设计方案与仿真.docx

《2022年matlab滤波器设计方案与仿真.docx》由会员分享，可在线阅读，更多相关《2022年matlab滤波器设计方案与仿真.docx（19页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精品学习资源封面欢迎下载精品学习资源作者： PanHongliang仅供个人学习目录欢迎下载精品学习资源1 引言.52 关于 MATLAB53 数字滤波的基本概念 64 设计方案 . .74.1 数字滤波器设计的基本步骤 .7欢迎下载精品学习资源4.1.1确定指标.74.1.2模型靠近.74.1.3实现性能分析和运算机仿真 .74.2 基于 MATLAB的 FIR 数字滤波器的设计与仿真 . 74.3 基于 MATLAB的 IIR 数字滤波器的设计. 104.3.1 IIR数字滤波器的设计原理 .104.3.2 IIR数字滤波器的传统设计方法 .104.3.3 IIR数字滤波器的设计 .114

2、.3.4 IIR数字滤波器的程序设计 .114.3.5 IIR数字滤波器的仿真 .125 FIR数字滤波器与 IIR 数字滤波器的比较 . .126 结论 13参考文献 错误！未定义书签；摘 要数字滤波器是一个离散系统；该系统能对输入的离散信号进行处理，从而获得所需的有用信息；现代数字滤波器的设计大体可以分为IIR 和 FIR 两大类，可以用软件和硬件两种方法来实现，而选用MATAB 信号处理工具箱为设计通用滤波器带来了极大的便利；本文按设计指标要求设计了滤波器，其中IIR 接受巴特沃什， FIR 接受布莱克曼窗进行设计，得出了与之对应的幅度响应曲线和相位响应曲线，最终对IIR 和 FIR 的

3、实现形式和性能等方面进行比较；关键词：MATLAB ；IIR ；FIRAbstractThedigitalfilteris a discretesystem.Thesystemcanbe abletohandle discretesignals.So itcanachieverequiredimportantinformation.Thereare欢迎下载精品学习资源twomajorkindsofdesignprincipleofdigitalfilter,whicharefiniteimpulse responseFIRand infiniteimpulseresponseIIR.Themo

4、derndigitalfilter can be receivedbytwokindsofmethodofsoftwareand hardware.Butusing MATLABsignaldisposingtool case to design the digital filter is moreconvenient and universally applied.The main body of thepaper is demanded to design adigital filter according to the designing index. IIR adopts Butter

5、worth and FIR adopts the Blackman window to design the digital filter.Finally,carryoutcomparisononIIRandtheFIRrealizationandfunction aspect.Keywords : MATLAB； IIR ； FIR欢迎下载精品学习资源1 引言理想滤波器就是一个让输入信号中的某些有用频谱重量无任何变化的通过，同时又能完全抑制另外那些不需要的成分的具有某种选择性的器件、网络或运算机硬件支持的运算程序；依据对不同信号的处理可将滤波器分为模拟滤波器和数字滤波器两种；模拟滤波器和数字

6、滤波器的概念相同，只是信号的形式和实现滤波的方法不同；数字滤波器是指输入输出都是数字信号的滤波器；滤波器的滤波原理就是依据信号与噪声占据不同的频带，将噪声的频率放在滤波器的阻带中而由于阻带的响应为零，这样就滤去了噪声；数字滤波器可以分为两大类：一类是经典滤波器，即一般的滤波器，特点是输入信号中有用的频率成分和期望滤去的频率成分各占不同的频率带，通过一个合适的选频滤波器达到滤波的目的，这种滤波器主要是无限冲激响应滤波器和有限冲激响应滤波器；另外一类滤波器是现代滤波器，当信号和干扰的频带相互重叠，经典滤波器不能完成对干扰的有效去除时，可以接受现代滤波 器，这些滤波器可以依据立刻信号内部的一些统计分

7、布规律，从干扰中正确的提取信号；这种滤波器主要有维纳滤波器、卡尔曼滤波器、自适应滤波器等10 ；与模拟滤波器相同，数字滤波器从功能上可以分为四类，即低通（LP）、高通（ HP）、带通（ BP）、带阻（ BS）滤波器等；数字滤波器在数字信号处理的各种应用中发挥着特殊重要的作用，任何数字系统都可看成一个数字滤波器，因此数字滤波器的应用相当广泛；设计任务：设计一数字带通滤波器，分别用IIR 和 FIR 两种类型来实现；2 关于 MATLABMATLAB是 Matrix Laboratory（矩阵试验室）的缩写；它是美国MathWorks公司推出的用于数值运算和图形处理的数学运算环境，现已进展成为一个

8、具有高性能数值运算和可视化功能的科学运算环境，问题的提出和解答只需以数学方式表达和描述，不需要大量原始而传统的编程过程；该软件包括：数值分析、矩阵运算、数字信号处理、建模和系统把握等应用程序；MATLAB语言简洁、紧凑、使用便利，有丰富的运算符号和函数库，仍具有结构化的把握语句和面对对象编程的特点；随着其版本的不断提高，各种工具箱的扩充和完善，欢迎下载精品学习资源MATLAB 的功能越来越强大，从而被广泛应用于仿真技术、自动把握和数字信号处理等领域；MATLAB 是一个交互式系统，特殊适用于争论、解决工程和数学问题，典型应用包括：一般的数值运算、算法原型以及通过矩阵公式解决一些特殊问题，促进了

9、统计、数字信号处理等科学的进展； MATLAB 功能强大、简洁易学、编程效率高，深受广大科技工作者的欢迎；特殊是 MATLAB 仍具有信号分析工具箱，不需具备很强的编程才能，就可以很便利地进行信号分析、处理和设计；MATLAB 的信号处理工具箱供应了丰富而简便的设计、实现IIR 和 FIR 的方法，使原先烦琐的程序设计简化成函数的调用，特殊是滤波器的表达方式和滤波器形式之间的相互转换显得特殊简便，为滤波器的设计和实现开创了一片宽敞的天地；数字滤波器接受传统的设计方法要进行大量复杂的运算，而利用MATLAB 强大的运算功能进行运算机帮忙设计，就可以快速要效的设计数字滤波器，大大的简化运算量，直观

10、简便； MATLAB的信号处理工具箱包含了各种经典的和现代的数字信号处理技术，是一个特殊优秀的算法争论和帮忙设计的工具；在设 计数字滤波器时，通常接受信号处理工具箱供应的设计模拟和数字滤波器的MATLAB函数接受编程的方法和仿真实现，亦可以利用信号处理工具箱供应的滤波器设计和分析工具（ Filter Design & Analysis Tool）实现；MATLAB的一种重要的工作方式就是M 文件的编程工作方式； M 文件有两种形式，一种是脚本文件，另一种是函数文件；M 文件的扩展名为“ .m”；M 文件可以通过任何纯文本编辑器进行编辑，MATLAB也自带有文本编辑器，使用“ edit”命令即可

11、开启；3 数字滤波的基本概念滤波器的种类很多，从功能上可以分为低通、高通、带通和带阻滤波器，每一种又有模拟滤波器和数字滤波器两种形式；假如滤波器的输入和输出都是离散时间信号，就该滤波器的冲击响应也必定是离散的，这种滤波器称之为数字滤波器；该滤波器通过对时域中离散的采样数据作差分运算实现滤波；一个线性时不变数字滤波器可以用常系数线性差分方程表示为欢迎下载精品学习资源式中 xn,yn 辨论是输入和输出信号序列；ai，bk 分别是滤波系数；当上式 bk 中系数全部为零时，就有这种形式的滤波器为有限冲激响应滤波器，简称FIR 型；此事系统的输出仅与输入有关；它的实现一般接受非递归算法；4 设计方案4.

12、1 数字滤波器设计的基本步骤4.1.1 确定指标在设计一个滤波器之前，必需第一依据工程实际的需要来确定滤波器 的技术指标；在很多实际应用中，数字滤波器常被用来实现选频操作；因此指标的形式一般在频域中给出幅度和相位响应；4.1.2 模型靠近一旦确定了技术指标，就可以利用已学习过的基本原理和关系式，提出一个滤波器模型来靠近给定的技术指标；这是滤波器设计所要争论的主要问题；4.1.3 实现性能分析和运算机仿真以上两步的结果得到的滤波器，通常是以差分方程、系统函数或脉冲响应描述的；依据这个描述就可以分析其频率特性和相位特性，以验证设计结果是否中意指标要求；或者利用运算机仿真实现设计的滤波器，再分析滤波

13、结果来判断；4.2 基于 MATLAB的 FIR 数字滤波器的设计与仿真利用 MATLAB 设计滤波器，可以随时对比设计要求和滤波器特性调整参数， 直观简便，极大地减轻了工作量，有利于滤波器设计的优化；FIR 滤波器的设计方法有窗函数法、最优化设计法及约束最小二乘靠近法；MATLAB信号处理工具箱中供应了滤波器设计方法的工具函数，编程中可以依据设计要求直接调用相应函数，便利快捷；本设计数字滤波器的性能指标为：，要求用窗函数设计一个高通滤波器；本例接受凯泽窗函数设计该高通滤波器，原程序节选如下：As=50；ws=0.2*Pi；欢迎下载精品学习资源wp=0.3*Pi ；tr_width=wp-ws

14、 ；M=ceilAs-7.95*2*pi/14.36*tr_width+l+l；disp滤波器的长度为 ,nu2mstrM ； beta=0.1102*AS-8.7；n=【0:l:M-l:disp 线性相位斜率为 ,nu2mstrbeta；w_kai=kaiserM,beta ；wc=ws+wp/2:hd=ideal_lppi,M-ideal_lpwc,M ；h=hd.*w_kai ；db,mag,pha,grd,w=freqz_mh,1；delta_w=2*pi/l000 ；Rp=-mindbwp/delta_w+l:l:501 ；disp 实际通带波动为 ,nu2mstrRp；As=-ro

15、undmaxdbl:l:ws/delta_w+l ；disp 最小阻带衰减为 ,nu2mstrAs ；subplotl,l,l ；subplot2,2,l；stemn,hd；title 理想脉冲响应 ；axis0,M-l -0,40.8 ；ylabel hdn； subplot2,2,2: stemn,w_kai:title 凯泽窗 :axis0,M-101.1 ；ylabel w；n subplot2,2,3；欢迎下载精品学习资源stemn,h；title 实际脉冲响应 ；axis【0,M-l -0.40.8 ；xlabel n； ylabel h；n subplot2,2,4；plotw/

16、pi,db ；title 幅度响应 /dB ；axis01 -10010；grid；欢迎下载精品学习资源xlabel ylabel以为单位的频率 :分贝数 /dB ；欢迎下载精品学习资源程序运行后得到的结果如图1 所示；由图 1 可知滤波器的数据：线性相位斜率为 4.5513，滤波器的长度为61，参数beta 为 4.55126，实际通带波动为0.040436，最小阻带衰减为 50，由此可知设计的滤波器达到了要求；程序执行的结果如下图所示：滤波器幅频特性相频特性滤波前后的波形FIR 的冲激响应从滤波成效图看 ,所设计的滤波器基本排除了噪声和干扰的影响,客观的得出真实信号； 在设计中 假如该滤波

17、器的特性不中意要求，原有的参数必需作相应的调整，在程序中只需对参数进行重新设定，就可以得到所需要的滤波器；在实际应用中，假如需要对某一信号源进行特定的滤波，并要检测滤波效果，应用传统方法比较繁琐；在MATLAB 环境下，可先用软件模拟产生信号源，再设计滤波器对其进行滤波；欢迎下载精品学习资源在 MATLAB 环境下，数字滤波器的设计已经变得特殊简洁和高效率，而且FIR 滤波器的冲激响应是有限长序列，其系统函数为一个多多项式，它所含的极点多为原点，所以 FIR 滤波器是稳固的；4.3 基于 MATLAB 的 IIR 数字滤波器的设计4.3.1 IIR数字滤波器的设计原理IIR 一般接受递归型的结

18、构，系统的输入与输出听从N阶差分方程：相应的传输函数为：设计 IIR 数字滤波器就是要确定传输函数中的系数aj 、bi ，或者零极点增益ci 、dj 、A，使滤波器的频率特性中意给定的性能指标要求；设计原理主要包 括两个方面：一是依据设计指标，先设计出相应的模拟滤波器, 再通过脉冲响应不变法或者双线性变化法转换成对应的数字滤波器；二是先选择一种最优准就，如最小均方误差准就，再在此准就下求出滤波传输函数的系数；4.3.2 IIR数字滤波器的传统设计方法IIR数字滤波器具有无线宽的冲击响应，与模拟滤波器相匹配；所以IIR滤波器的设计可以实行在模拟滤波器设计的基础上进一步变换的方法；例如利 用脉冲响

19、应不变法设计IIR数字滤波器，其基本设计思路是直接设计一个数字滤波器并让它的时间特性靠近一个模拟滤波器；为了达到时间特性的正确逼近，把模拟滤波器的冲击响应均匀取值，作为数字滤波器的单位脉冲响应， 即：然后将 hn 通过 z 变换即糗得 Hz ；因此，脉冲响应不变法实际是时域取样法，整个过程先是依据给定的指标设计一个模拟滤波器，进而按以下变换的 次序，最终糗得数字滤波器的系统函数Hz, 即；在对滤波器进行设计时，整个过程运算量是很大的；设计完成后要对已设计的滤波器的频率响应进行校核，得到幅频相频响应特性，运算量也是很大 的；通常，待设计的数字滤波器，阶数和类型并不愿定是完全给定的，很多时欢迎下载

20、精品学习资源候都是要依据设计要求和滤波成效不断进行调整，以达到设计的最优化；在这种情形下，滤波器的设计就要进行大量复杂的运算，单纯的靠公式运算和编制简洁的程序很难在短时间使问题得到解决；4.3.3 IIR数字滤波器的设计FDAToolFilter Design & Analysis Tool是 MATLAB信号处理工具箱里专用的滤波器设计分析工具， MATLAB 6.0 以上的版本仍特地增加了滤波器设计工具箱Filter Design Toolbox；FDATool 可以设计几乎全部的常规滤波器，包括 FIR 和 IIR 的各种设计方法；它操作简洁，便利灵敏；下面就利用本文要求设计一个 10

21、阶的带通 Chebyshev I滤波器，它的通带范畴是 100到 200Hz，采样频率为 1000Hz，Rp=0.5；本例中，第一在Filter Type中选择Bandpass带通滤波器 ；在 Design Method选项中选择 IIR ，接着在相邻的右就选项中选择Chebyshev. 切比雪夫 I型 ；指定 Filter Order项中的Specify Order=10；由于接受的是切比雪夫设计，不必在Options中选择；然后在 Frequency Specifications中选择 Unit 为 Hz，给出采样频率 Fs= 1000， 通带 Fpass1 = 100和 Fpass2=2

22、00；最终在 Magnitude Specifications中选择Unit 为 db， Apass=0.5；设置完成后点击 DesignFilter即可得到所设计的 IIR 滤波器；通过菜单项项 Analysis可以在特性区看到所设计的幅频响应、相频响 应、冲激响应和零极点配置等特性，如图1 所示；设计完成后将结果储存为filter1.fda文件；图 1 滤波器的幅频、相频和冲激响应 特性区4.3.4 IIR数字滤波器的程序设计上例 IIR 滤波器的冲激响应可用程序设计如下：%chebyshev1 bpf n=10； % 阶数为 10Rp=0.5； %幅值衰减为 0.5Wn=100 200/

23、500 ；b,a=cheby1n,Rp,Wn；y,t=impzb,a,101；欢迎下载精品学习资源stemt,y, . ；由图 1a 可知，这种滤波器在 100-200Hz 的通带范畴内是等波纹的 , 而在阻带中是单调的，这是 Chebyshev I 滤波器的幅频特性；由图 1b 可知，在100-200Hz 的范畴内相移较小，其曲线近似一条直线，失真较小；当频率超过这一范畴时，相移较大，而且其曲线是非直线的，所以失真也较大；图 1c 得到了 Chebyshev I 滤波器的时域冲激响应，在 5100ms有冲激响应，超过这一范畴的冲激响应近似为零，进而实现了 Chebyshev I 带通滤波器的

24、设计；4.3.5 IIR数字滤波器的仿真通过调用 Simulink中的功能模块，可以构成数字滤波器的仿真框图；在仿真过程中，双击各功能模块，随时转变参数，获得不同状态下的仿真结果；例如原始信号 xt=2sin0.05Pt+wt ，wt 为随机信号，幅值为 0.2 ，通过一传-1递函数为 H z=0.15/1-0.8z 的滤波器可得到如图 2 的仿真结果；其中仿真过程中可导入 FDATool所设计的滤波器文件；图 2 Simulink仿真结果-1比较图 2 中b 和c 的波形可知，输入的原始信号经过滤波器滤波后，b 中波形的毛刺部分 即干扰噪声 被滤除，输出的信号更接近正弦波，如c 中所示波形；

25、由此说明，传递函数为H z=0.15/1-0.8z 的滤波器的设计是恰当的；5 FIR数字滤波器与 IIR 数字滤波器的比较从性能上来说， IIR滤波器传递函数包括零点和极点两组可调因素，对极点的惟一限制是在单位圆内；因此可用较低的阶数获得高的选择性，所用的储备单 元少，运算量小，效率高；但是这个高效率是以相位的非线性为代价的；选择 性越好，就相位非线性越严肃；FIR 滤波器传递函数的极点固定在原点，是不能动的，它只能靠转变零点位置来转变它的性能；所以要达到高的选择性，必须用较高的阶数；对于同样的滤波器设计指标，FIR 滤波器所要求的阶数可能比 IIR 滤波器高 5-10 倍，结果，成本较高，

26、信号延时也较大；假如按线性相位要求来说，就 IIR 滤波器就必需加全通网络进行相位校正，同样要大大增加滤波器的阶数和复杂性；而 FIR 滤波器却可以得到严格的线性相位；欢迎下载精品学习资源从结构上看， IIR滤波器必需接受递归结构来配置极点，并保证极点位置在单位圆内；由于有限字长效应，运算过程中将对系数进行舍入处理，引起极点的偏移；这种情形有时会造成稳固性问题，甚至产生寄生振荡；相反，FIR滤波器只要接受非递归结构，不论在理论上仍是在实际的有限精度运算中都不存在稳固性问题，因此造成的频率特性误差也较小；此外FIR 滤波器可以接受快速傅里叶变换算法，在相同阶数的条件下，运算速度可以快得多；另外，

27、也应看到， IIR滤波器虽然设计简洁，但主要是用于设计具有分段常数特性的滤波器，如低通、高通、带通及带阻等，往往脱离不了模拟滤波器 的格局；而 FIR 滤波器就要灵敏得多，特殊是他易于适应某些特殊应用，如构成数字微分器或希尔波特变换器等，因而有更大的适应性和宽敞的应用领域；从上面的简洁比较可以看到IIR与 FIR 滤波器各有所长，所以在实际应用时应当从多方面考虑来加以选择；从使用要求上来看，在对相位要求不敏捷的 场合，如语言通信等，选用IIR较为合适，这样可以充分发挥其经济高效的特点；对于图像信号处理，数据传输等以波形携带信息的系统，就对线性相位要 求较高；假如有条件，接受FIR 滤波器较好；

28、当然，在实际应用中可能仍要考虑更多方面的因素；6 结论利用 MATLAB 的信号处理工具箱强大的信号处理功能，接受编程的方法设计数字滤波器，使设计达到了最优化；并且它可以快速的实现数字滤波器的仿真，使设计达到了最简化；本文设计的数字滤波器接受IIR 和 FIR 两种方法来实现，其中 IIR 接受巴特沃什， FIR 接受布莱克曼窗函数法，按设计指标要求进行设计；利用 MATLAB 接受编程的方法实现；通过对IIR 数字滤波器和 FIR 数字滤波器的比较可知IIR 数字滤波器幅频特性较好，而相频特性曲线呈非线性，会使信号产生失真，阶数比FIR 数字滤波器低，系统不稳固； FIR 数字滤波器相位线性

29、，信号在传输过程中不会产生失真，系统稳固，阶数高，而幅频特性难以中意要求；参考文献 1蒋志凯数字滤波与卡尔曼滤波 M ，北京：中国科学技术出版社，1993：100-120 2赵健，李勇数字信号处理 M ，北京：清华高校出版社，2006 ：103-155 欢迎下载精品学习资源 3陈怀琛数字信号处理教程MA TLAB释义与实现 M ，北京：电子工业出版社，2004 ： 106-123 4邓华 MATLAB通信仿真及应用实例解读 M ，北京：人民邮电出版社，1997： 123- 156 5胡光书数字信号处理：理论、算法与实现 M ，北京：清华高校出版社，1997 ：136-155 6楼天顺基于MA

30、TLAB的系统分析与设计 - 信号处理 M ，西安：西安电子科技高校出版社， 1998： 133-142 7徐明远 MATLAB仿真在通信与电子工程中的应用 M ，西安：西安电子科技高校出版社， 2005： 154-187 8陈后金数字信号处理 M ，北京：高等训练出版社，2004： 124-177 9 美Sanjit K.Mitra著，孙烘，余翔宇等译数字信号处理试验指导书（MA TLAB版） M ，北京：电子工业出版社，2005： 71-85 10陈桂明，张明照，戚红雨应用MA TLAB语言处理数字信号与数字图像 M ，北京：科学出版社， 2001： 165-200. 11赵健，李毅，杨菊

31、平数字信号处理学习辅导 M ，北京：清华高校出版社，欢迎下载精品学习资源2006： 93-95版权申明欢迎下载精品学习资源本文部分内容，包括文字、图片、以及设计等在网上搜集整理；版权为潘宏亮个人全部This article includes some parts, including text,pictures, and design. Copyright is Pan Hongliangs personalownership.用户可将本文的内容或服务用于个人学习、争论或观看，以及其他非商业性或非盈利性用途，但同时应遵守著作权法及其他相关法律的规定，不得侵害本网站及相关权益人的合法权益；除此以

32、欢迎下载精品学习资源外，将本文任何内容或服务用于其他用途时，须征得本人及相关权益人的书面许可，并支付酬劳；Users may use the contents or services of this article for personal study, research or appreciation, and other non-commercial or non-profit purposes, but at the same time, they shall abide by the provisions of copyright law and other relevant laws

33、, and shall not infringe upon the legitimate rights of this website and its relevant obligees. In addition, when any content or service of this article is used for other purposes, written permission and remuneration shall be obtained from the person concerned and the relevant obligee.转载或引用本文内容必需是以新闻

34、性或资料性公共免费信息为使用目的的合理、善意引用，不得对本文内容原意进行曲解、修 改，并自负版权等法律责任；Reproduction or quotation of the content of this article must be reasonable and good-faith citation for the use of news or informative public free information. It shall not misinterpret or modify the original intention of the content of this article, and shall bear legal liability such as copyright.欢迎下载精品学习资源欢迎下载