利用FDATool进行数字滤波器的设计以及仿真实验.docVIP

利用FDATool进行IIR数字滤波器的设计以及仿真实验 姓名：许锦旋 专业：物理电子学 IIR数字滤波器具有无限宽的冲激响应 ，与模拟滤波器相匹配 ，所以IIR滤波器的设计可以采取在模拟滤波器设计基础上进一步变换的方法。其设计方法主要有经典设计法、直接设计法和最大平滑滤波器设计法，现在主要介绍利用FDATool的设计方法。 1　FDATool界面设计 1.1　FDATool的介绍 FDATool (Filter DesignAnalysis Tool)是MATLAB信号处理工具箱里专用的滤波器设计分析工具 ，MATLAB6.0以上的版本还专门增加了滤波器设计工具箱 (Filter Design Toolbox)。 FDATool可以设计几乎所有的常规滤波器，包括FIR和IIR的各种设计方法 。它操作简单 ,方便灵活。 FDATool界面分为两大部分，一部分是Design Filter,在界面的下半部，用来设置滤波器的设计参数；另一部分则是特性区，在界面的上半部分，用来显示滤波器的各种特性。 Design Filter部分主要分为：FilterType(滤波器类型 )选项，包括Lowpass(低通 )、Highpass(高通 )、Bandpass(带通 )、Bandstop(带阻 )和特殊的FIR滤波器。 Design Method(设计方法 )选项，包括IIR滤波器的Butterowrth(巴特沃思 )法、Chebyshev Type I(切比雪夫Ｉ型 )法、Chebyshev Type Ⅱ (切比雪夫Ⅱ型 )法、Elliptic(椭圆滤波器 )法和FIR滤敞器的Equiripple法、Least-Squares(最小乘方 )法、Window(窗函数 )法。 Filter Order(滤波器阶数 )选项，定义滤波器的阶数，包括Specify Order(指定阶数 )和Minimum Order(最小阶数 )。在Specify Order中填入所要设计的滤波器的阶数(N阶滤波器， Specify Order=N -1)。如果选择Minimum Order，则MATLAB根据所选择的滤波器类型自动使用最小阶数。 Frequency Specifications选项，可以详细定义频带的各参数，包括采样频率和频带的截止频率。它的具体选项由Filter Type选项和Design Method选项决定。例如Bandpass(带通 )滤波器需要定义Fstop1(下阻带截止频率 )、Fpass1(通带下限截止频率 )、Fpass2(通带上限截止频率 )、Fstop2(上阻带截止频率 )，而Lowpass(低通 )滤波器只需要定义Fstop1、Fpass1。采用窗函数设计滤波器时，由于过渡带是由窗函数的类型和阶数所决定，所以只需定义通带截止频率，而不必定义阻带参数。 Magnitude specifications选项，可以定义幅值衰减的情况。例如设计带通滤波器时，可以定义Wstop1(频率Fstop1处的幅值衰减 )、Wpass(通带范围内的幅值衰减 )、Wstop2(频率Fstop2处的幅值衰减 )。当采用窗函数设计时 ,通带截止频率处的幅值衰减固定为6db，所以不必定义。 Window Specifications选项 ，当选取采用窗函数设计时，该选项可定义，它包含了各种窗函数。 1.2　IIR滤波器设计实例 下面以一个IIR滤波器的设计为例来说明如何使用MATLAB设计数字滤波器。要求设计一个3阶的低通ChebyshelⅠ滤波器，要求通带是0到 20Ｈｚ，采样频率为 1000Hz，Rp =0.5。 本例中，首先在Filter Type中选择Lowpass(低通滤波器 )；在Design Method选项中选择IIR，接着在相邻的右则选项中选择ChebyshevI(切比雪夫I型 )；指定Filter Order项中的Specify Order=3；由于采用的是切比雪夫设计，不必在Options中选择；然后在Frequency Specifications中选择Unit为Hz，给出采样频率Fs =1000，Fpass=20；最后在Magnitude specifications中选择Unit为db，Apass =0.5。设置完成后点击Design Filter即可得到所设计的IIR滤波器。通过菜单选项Analysis可以在特性区看到所设计的幅频响应、相频响应、冲激响应和零极点配置等特性，如图 1、2、3所示。设计完成后将结果保存为lowpassfilter.fda文件。 图1：低通滤波器的幅频特性 图2：低通滤波器的相频特性 图3：低通滤波器的冲击响应 1.3　设计实例分析 由图 （1）可知，所设计的低通滤波器在0-20Hz的通带