matlab滤波器设计

1、matlab滤波器设计 在现代通信系统中，因为信号中常常混有各种复杂成分，所以无数信号的处理和分析都是基于而举行的。但是，传统的数字滤波器的设计用法繁琐的公式计算，转变参数后需要重新计算，从而在设计滤波器尤其是高阶滤波器时工作量很大。利用matlab信号处理箱(signal processing toolbox)可以迅速有效地实现数字滤波器的设计与。 2 数字滤波器及传统设计办法数字滤波器可以理解为是一个计算程序或算法，将代表输入信号的数字时光序列转化为代表输出信号的数字时光序列，并在转化过程中，使信号按预定的形式变幻。数字滤波器有多种分类，按照数字滤波器冲激响应的时域特征，可将数字滤波器分为

2、两种，即无限长冲激响应(iir)滤波器和有限长冲激响应(fir)滤波器。iir数字滤波器具有无限宽的冲激响应，与模拟滤波器相匹配，所以iir滤波器的设计可以实行在模拟滤波器设计的基础上进一步变换的办法。其设计办法主要有经典设计法、挺直设计法和最大平滑滤波器设计法。fir数字滤波器的单位脉冲响应是有限长序列。它的设计问题实质上是确定能满足所要求的转移序列或脉冲响应的常数问题，设计办法主要有窗函数法、频率采样法和等波纹最佳靠近法等。在对滤波器实际设计时，囫囵过程的运算量是很大的。设计阶数较高的iir滤波器时，计算量更大，设计过程中转变参数或滤波器类型时都要重新计算。设计完成后对已设计的滤波器的频率

3、响应要举行校核。要得到幅频、相频响应特性，运算量也是很大的。平常所要设计的数字滤波器，阶数和类型并不一定是彻低给定的，无数时候要按照设计要求和滤波效果不断地调节，以达到设计的最优化。在这种状况下，滤波器设计就要举行大量复杂的运算，单纯的靠公式计算和编制容易的程序很难在短时光内完成。利用matlab强大的计算功能举行计算机辅助设计，可以迅速有效地设计数字滤波器，大大地简化了计算量。3 iir滤波器的matlab设计3.1 fdatool界面设计3.1.1 fdatool的介绍fdatool(filter design&analysis tool)是matlab信号处理工具箱里专用的滤波器

4、设计分析工具，matlab 6.o以上的版本还特地增强了滤波器设计工具箱(filter design toolbox)。fdatool可以设计几乎全部的常规滤波器，包括fir和iir的各种设计办法。它操作容易，便利灵便。fdatool。界面总共分两大部分，一部分是design filter。在界面的下半部，用来设置滤波器的设计参数；另一部分则是特性区，在界面的上半部分，用来显示滤波器的各种特性。design filter部分主要分为：filter type(滤波器类型)选项，包括lowpass(低通)、highpass(高通)、bandpass(带通)、bandstop(带阻)和特别的fir滤

5、波器。design method(设计办法)选项，包括iir滤波器的butterwotth(巴特沃思)法、chebyshev type i(切比雪夫i型)法、chebyshev type ii(切比雪夫ii型)法、elliptic(椭圆滤波器)法和fir滤波器的equiripple法、least-squares(最小乘方)法、window(窗函数)法。filter order(滤波器阶数)选项，定义滤波器的阶数，包括specify order(指定阶数)和minimum order(最小阶数)。在specify order中填入所要设计的滤波器的阶数(n阶滤波器，specify order=n

6、-1)。假如挑选minimum order，则matlab按照所挑选的滤波器类型自动用法最小阶数。frequency specifications选项，可以具体定义频带的各参数，包括采样频率和频带的截止频率。它的详细选项由filter type选项和design method选项打算。例如bandpass(带通)滤波器需要定义fstop1(下阻带截止频率)、fpass1(通带下限截止频率)、fpass2(通带上限截止频率)、fstop2(上阻带截止频率)，而lowpass(低通)滤波器只需要定义fstop1、fpass1。采纳窗函数设计滤波器时，因为过渡带是由窗函数的类型和阶数所打算，所以只需

7、定义通带截止频率，而不必然义阻带参数。magnitude specifications选项，可以定义幅值衰减的状况。例如设计带通滤波器时，可以定义wstop1(频率fstop1处的幅值衰减)、wpass(通带范围内的幅值衰减)、wstop2(频率fstop2处的幅值衰减)。当采纳窗函数设计时，通带截止频率处的幅值衰减固定为6db，所以不必然义。window specifications选项，当选取采纳窗函数设计，该选项可定义，它包含了各种窗函数。3.1.2 iir滤波器设计实例本文以一个iir滤波器的设计为例解释如何用法matlab设计数字滤波器。要求设计一个10阶的带通chebyshev i

8、滤波器，它的通带范围是100到200hz，采样频率为1000hz，rp=0.5。本例中，首先在filter type中挑选bandpass(带通滤波器)；在design method选项中挑选iir，接着在相邻的右则选项中挑选chebyshev i(切比雪夫i型)；指定filter order项中的specify order=10；因为采纳的是切比雪夫设计，不必在options中挑选；然后在frequency specifications中挑选unit为hz，给出采样频率fs=1000，通带fpass1=100和fpass2=200；最后在magnitude specifications中挑选

9、unit为db，apass=0.5。设置完成后点击design filter即可得到所设计的iir滤波器。通过菜单选项analysis可以在特性区看到所设计的幅频响应、相频响应、冲击响应和零极点配置等特性，1所示。设计完成后将结果保存为filterl.fda文件。图1滤波器的幅频、相频和冲激响应(特性区)3.2 程序设计法在matlab中，对各种滤波器的设计都有相应的计算振幅、相位和冲激响应的函数，可以用来做滤波器程序设计。上例的iir滤波器的冲激响应可用程序设计如下：%ehebyshevl bpfn=10; %阶数为10rp=0.5; %幅值衰减为0.5wn100 200/500；b,a=c

10、hebyl(n,rp,wn);y,timpz(b,a,101);stem(t,y,'.');在matlab环境下运行该程序即可得到滤波器的时域冲激响应。因为篇幅所限，这里不再详述源程序。3.3 设计实例分析由图1(a)可知，这种滤波器在100-200hz的通带范围内是等波纹的，而在阻带中是单调的，这是chebyshev i滤波器的幅频特性。由图1(b)可知，在100-200hz的范围内相移较小，其曲线近似一条直线，失真较小；当频率超过这一范围时，相移较大，而且其曲线是非直线的，所以失真也较大。图1(c)得到了chebyshev i滤波器的时域冲激响应，在5ms-100ms有冲激

11、响应，超过这一范围的冲激响应近似为零，进而实现了chebyshev i带通滤波器的设计。4 simulink仿真4.1 simulink仿真切例通过调用simulink中的功能模块，可以构成数字滤波器的仿真框图。在仿真过程中，双击各功能模块，随时转变参数，获得不同状态下的仿真结果。例如原始信号x(t)=2sin(0.05wt)+w(t)，w(t)为随机信号，幅值为0.2，通过一传递函数为h(z)=0.15/(1-0.8z-1)的滤波器可得到2的仿真结果。其中仿真过程中可导入fdatool所设计的滤波器文件。4.2 仿真分析比较图2中(b)和(c)的波形可知，输入的原始信号经过滤波器滤波后，(b)中波形的毛刺部分(即干扰噪声)被滤除，输出的信号更临近正弦波，如(c)中所示波形。由此解释，传递函数为h(z)=0.15/(1-0.8z-1)的滤波器的设计是恰当的。5 结论利用matlab的强大运算功能，基于matlab