MATLABIIR数字滤波器设计附DSP实现.docVIP

基于MATLAB的IIR数字滤波器的设计 （MATLAB在通信中的应用-----课程论文） 论文作者: 王紫艳 专业：： 基于MATLAB的IIR数字滤波器的设计(DSP)工具箱(signal processing toolbox)可以快速有效的实现数字滤波器的设计与仿真。MATLAB软件的频谱分析和滤波器的分析设计功能很强，从而使数字信号处理变得十分简单。本文比较了设计IIR数字滤波器的两种方法，介绍了利用双线性变换法设计IIR滤波器的原理、步骤，并在MATLAB下进行了仿真。 关键字：MATLAB IIR 数字滤波器 DSP 一二IIR数字滤波器的特征是具有无限持续时间的冲激响应，这种滤波器也称之为递归型滤波器。IIR数字滤波器的设计常采用两种方法：一、冲激响应不变法 二、双线性变换法。 2.1 冲激响应不变法 冲激响应不变法是从时域出发，要求数字滤波器的冲激响应h(n)对应于模拟滤波器的冲激响应ha(t)的抽样,即 h(n) = ha(nT) ，其中T是抽样周期。冲激响应不变法 Ha(s) → ha(t) → h(n)? →? H(z) ? 由于数字滤波器的冲激响应对应模拟滤波器冲激响应的抽样，由抽样定理可知数字滤波器的频率响应为模拟滤波器频率响应的周期延拓: 其中 为了避免混叠失真要求模拟滤波器的频谱限带于折叠频率以内，即要满足: 而实际的滤波器并非严格限带，所以用冲激响应不变法设计的数字滤波器不可避免地会产生混叠失真。所以此法只适于设计带限滤波器。其映射关系如下： 在MATLAB工具箱中，提供了专用函数impinvar来实现以上计算，其调用格式为： 格式一：[BZ,AZ]= impinvar(B,A,Fs) 功能：把具有[B,A]模拟滤波器传递函数模型转换为采样频率为Fs(HZ)的数字滤波器的传递函数模型[BZ,AZ]。采样频率默认值为Fs=1。 格式二：[BZ,AZ]= impinvar(B,A,Fs，TOL) 功能：利用指定的容错误差TOL来确定极点是否重复。如果设置的容差增大，则函数认为相邻的很近的极点为重复极点的可能性增大。默认的TOL=0.001，即0.1%。 2.2 双线性变化法原理 为了 (1)通过下面的变换公式把整个S平面压缩到中介S1平面的一条横带里（宽度从-π/T到π/T），其中C为常数可根据设计要求选取； (2) 通过 将此横带变换到整个Z平面，这样就使S平面和Z平面是一一对应的单值映射关系，消除了频谱混叠现象。 图2-2 双线性变换法的映射关系 可以证明，从SZ平面的映射关系为： ， 由于从S→Z和从Z→S的映射规则都是分式线性变换，因此称其为双线性变换。 MATLAB信号处理工具箱为实现双线型变换提供了函数bilinear,基本调用格式如下： 格式一：[Zd,Pd,Kd]= bilinear(Z,P,K,Fs) 功能：把模拟滤波器的零极点模型转换为数字滤波器的零极点模型。其中Fs为采样频率。 格式二：[numd,dend]= bilinear(num,den,Fs) 功能：把模拟滤波器的传递函数模型转换为数字滤波器的传递函数模型。 格式三：[Ad,Bd,Cd,Dd]= bilinear(A,B,C,D,Fs) 功能：把模拟滤波器的状态方程模型转换为数字滤波器的状态方程模型。 说明：以上三种调用格式中，可以在增设一个畸变频率Fp(HZ)输入参数。在进行双线性变换之前，对采样了频率进行畸变处理，以保证频率冲激响应在双线性变换前后，在Fp处具有良好的单值映射关系。 三IIR数字滤波器的设计步骤 (1)按一定规则将给出的数字滤波器的技术指标转换为模拟滤波器的技术指标； (2)根据转换后的技术指标设计模拟低通滤波器H(s)； (3)在按一定规则将H(s)转换为H(z)； 若所设计的数字滤波器是低通的，那么上述设计工作可以结束，若所设计的是高通、带通或者带阻滤波器，那么还有步骤： (4)将高通、带通或者带阻数字滤波器的技术指标先转化为低通滤波器的技术指标，然后按上述步骤(2)设计出模拟低通滤波器H(s)，再将H(s)转换为所需的H(z)。 四 IIR 数字滤波器的MATLAB 设计及实现 MATLAB信号工具箱提供了几个直接设计IIR数字滤波器的函数，直接调用这些函数就可以很方便地对滤波器进行设计。这里选取巴特沃斯法、切比雪夫I、切比雪夫Ⅱ、椭圆法四种方法进行比较。给出用上述方法设计数字滤波器的函数如下： Butterworth 滤波器： [N，Wc]=buttoed(Wp,Ws,Rp,Rs) [b,a]=butter(N,Wc) [b,a]=