数字信号处理IIR系统设计ppt

1、第八章第八章 iiriir系统的设计系统的设计 iiriir系统的设计概述系统的设计概述 原型滤波器与模拟系统设计原型滤波器与模拟系统设计 模拟系统与数字系统的变换模拟系统与数字系统的变换 数字信号处理数字信号处理 iiriir系统的表达与计算系统的表达与计算 冲激响应长度无限冲激响应长度无限, , 采用系统函数方式表采用系统函数方式表 达，利用递归差分方程形式实现：达，利用递归差分方程形式实现： n n n n zaza zbzbb zh .1 . 1 1 1 10 .1.1 110 nyanxbnxbny iiriir系统的硬件实现系统的硬件实现 n n n n zaza zbzbb zh

2、 .1 . 1 1 1 10 iiriir滤波器的特点滤波器的特点 n n阶系统具有阶系统具有n n个不在原点的极点，在设计同个不在原点的极点，在设计同 样性能的滤波器时，样性能的滤波器时，iiriir滤波器的阶数通常低滤波器的阶数通常低 于于firfir滤波器，设计效率较高。滤波器，设计效率较高。 分母多项式系数不可能设置为对称，因此相分母多项式系数不可能设置为对称，因此相 频特性通常无法控制为线性，只能针对幅频特频特性通常无法控制为线性，只能针对幅频特 性进行设计。性进行设计。 设计目标幅频特性为：设计目标幅频特性为： 1 1 将设计目标转换为模拟系统的幅频特性：将设计目标转换为模拟系统的

3、幅频特性： 2 2 设计满足要求的模拟系统：设计满足要求的模拟系统： 3 3 将模拟系统转换为数字系统：将模拟系统转换为数字系统： 模拟变换法基本设计思想模拟变换法基本设计思想 j eh jh sh zhsh iiriir滤波器的模拟变换法设计滤波器的模拟变换法设计 问题：问题： 如何设计模拟系统？如何设计模拟系统？ 如何将数字幅频特性转换为模拟幅频特性？如何将数字幅频特性转换为模拟幅频特性？ 如何将模拟系统转换为数字系统？如何将模拟系统转换为数字系统？ 模拟滤波器基本设计步骤模拟滤波器基本设计步骤 将幅频特性的平方表达为将幅频特性的平方表达为2n2n阶多项式分式：阶多项式分式： n nn n

4、 nn ccc ddd c d mjh . . 22 1 2 0 22 1 2 0 2 2 2 2 利用各种优化逼近设计方法，可以得到利用各种优化逼近设计方法，可以得到 分式中的各优化系数。分式中的各优化系数。 根据幅频特性平方与系统函数的关系：根据幅频特性平方与系统函数的关系： )( 2 2 sgshshjhjhjh js 可将平方系统的系统函数表达为：可将平方系统的系统函数表达为： n nn n nn cscsc dsdsd shshsg . . 1 2 1 2 0 1 2 1 2 0 2 模拟滤波器基本设计步骤模拟滤波器基本设计步骤 和和 的零极点形成对称分布的零极点形成对称分布 shs

5、h 模拟滤波器基本设计步骤模拟滤波器基本设计步骤 得到优化的多项式系数后，可以求解出得到优化的多项式系数后，可以求解出 全部的零极点；全部的零极点； 选取虚轴左半的零极点就可以得到因果选取虚轴左半的零极点就可以得到因果 稳定的系统函数。稳定的系统函数。 模拟滤波器基本设计步骤模拟滤波器基本设计步骤 第八章第八章 iiriir系统的设计系统的设计 iir iir系统的设计概述系统的设计概述 原型滤波器与模拟系统设计原型滤波器与模拟系统设计 模拟系统与数字系统的变换模拟系统与数字系统的变换 数字信号处理数字信号处理 模拟原型滤波器模拟原型滤波器 当设计目标是对理想滤波器的逼近时，当设计目标是对理想

6、滤波器的逼近时， 通常可以采用原型滤波器进行变量代换设通常可以采用原型滤波器进行变量代换设 计。计。 原型滤波器定义：原型滤波器定义： 模拟低通滤波器，截止频率模拟低通滤波器，截止频率1 c 模拟原型滤波器：模拟原型滤波器： butterworth 最平坦滤波器最平坦滤波器 n m 2 2 1 1 n jh 2 1 1 对于任意阶数对于任意阶数n n， 随频率增加，系统增益单调下降；随频率增加，系统增益单调下降； 截止频率为截止频率为-3db-3db频率；频率； 随着随着n n值增大，可以逼近理想滤波器。值增大，可以逼近理想滤波器。 1 0 jh 2 1 1 2 jh 模拟原型滤波器：模拟原型

7、滤波器： butterworth 由幅频特性平方可以求出平方系统的零极点：由幅频特性平方可以求出平方系统的零极点： n m 2 2 1 1 .3, 11 2 mes jm n p .3, 1 2 1 2 mes n m j p 模拟原型滤波器：模拟原型滤波器： butterworth n s sg 2 2 1 1 系统函数为：系统函数为： pnpp ssssss sh . 1 21 模拟原型滤波器：模拟原型滤波器： butterworth z,p,kz,p,k=buttap(nbuttap(n) ) 得出得出n n阶阶butterworthbutterworth原型滤波器的零点原型滤波器的零点

8、z z （空向量），极点（空向量），极点p p向量和增益向量和增益k k（通常为（通常为 1 1）；）； butterworthbutterworth原型滤波器的设计仿真原型滤波器的设计仿真 b,ab,a=zp2tf(z,p,k) =zp2tf(z,p,k) 根据零极点分布给出滤波器的系统函数：根据零极点分布给出滤波器的系统函数： 对于对于butterworthbutterworth，b=1b=1； h,wh,w=freqs(b,afreqs(b,a) ) 根据模拟系统函数得出模拟频率响应。根据模拟系统函数得出模拟频率响应。 butterworthbutterworth原型滤波器的设计仿真原型

9、滤波器的设计仿真 仿真函数的使用举例仿真函数的使用举例 z,p,k=buttap(12); b,a=zp2tf(z,p,k); h,w=freqs(b,a);zplane(b,a); figure(2),plot(w,abs(h); axis(0,4,0,1.2); butterworthbutterworth原型滤波器的设计仿真原型滤波器的设计仿真 仿真函数的使用举例仿真函数的使用举例 butterworthbutterworth原型滤波器的设计仿真原型滤波器的设计仿真 butterworthbutterworth低通滤波器设计低通滤波器设计 对原型滤波器进行变量代换，可以得到指对原型滤波器

10、进行变量代换，可以得到指 定截止频率的低通模拟滤波器：定截止频率的低通模拟滤波器： n c n jh 22 /1 1 1 1 : 设计参数：设计参数： 滤波器阶数滤波器阶数 n n -3db -3db截止频率截止频率 c butterworthbutterworth低通滤波器设计低通滤波器设计 考虑滤波器参数与设计指标之间的关系：考虑滤波器参数与设计指标之间的关系： 通带波动通带波动 通带截止频率通带截止频率 阻带波动阻带波动 阻带截止频率阻带截止频率 在通带和阻带的边缘，可以得到：在通带和阻带的边缘，可以得到： p s p s p n cp 1 /1 1 2 s n cs 2 /1 1 利用

11、对数坐标，对纹波采用分贝为单位：利用对数坐标，对纹波采用分贝为单位： p n cp 1 . 02 10/1 s n cs 1 . 02 10/1 对上式联立求解，可以得到滤波器最低阶对上式联立求解，可以得到滤波器最低阶 数为：数为： sp s p n loglog 110log110log 2 1 1 . 0 1 . 0 butterworthbutterworth低通滤波器设计低通滤波器设计 相关的相关的matlabmatlab函数函数 n, wn = buttors(wp, ws, rp, rs, s) 给定滤波器的设计参数，上述函数可以求解出给定滤波器的设计参数，上述函数可以求解出 bu

12、tterworthbutterworth滤波器的最低阶数和滤波器的最低阶数和-3db-3db截止频率。截止频率。 22 2 1 1 n v m 22 1 1 n v jh n n阶阶chebyshevchebyshev多项式：多项式： 21 2 nnn vvv 1 v 12 2 2 v 模拟原型滤波器：模拟原型滤波器： chebyshev 1 没有零点，极点在单位圆内呈椭圆分布；没有零点，极点在单位圆内呈椭圆分布； 通带为等纹波，阻带为单调变化；通带为等纹波，阻带为单调变化； 模拟原型滤波器：模拟原型滤波器： chebyshev 1 z,p,kz,p,k=cheb1ap(n,rp) =cheb

13、1ap(n,rp) 给出通带波动为给出通带波动为rprp（dbdb）的）的n n阶阶1 1型型chebyshevchebyshev滤滤 波器的零点波器的零点z z，极点，极点p p和增益和增益k k； z,p,k=cheb1ap(10,1); b,a=zp2tf(z,p,k); h,w=freqs(b,a);zplane(b,a); figure(2),plot(w,abs(h); axis(0,4,0,1.2); chebyshevchebyshev 1 1型滤波器的设计仿真型滤波器的设计仿真 chebyshevchebyshev 1 1型低通滤波器的设计型低通滤波器的设计 滤波器设计参数：

14、滤波器设计参数： p pn v jh /1 1 22 滤波器阶数滤波器阶数 通带截止频率通带截止频率 通带纹波参数通带纹波参数 n 可以直接代入频率响应表达式；可以直接代入频率响应表达式； 利用通带波动指标确定纹波参数：利用通带波动指标确定纹波参数： 利用阻带衰减和阻带截止频率确定滤波器最小阶数：利用阻带衰减和阻带截止频率确定滤波器最小阶数： p p 2 1/1lg10 nv sns 22 1log10 chebyshevchebyshev 1 1型低通滤波器的设计型低通滤波器的设计 确定频率响应函数后，利用平方系统函确定频率响应函数后，利用平方系统函 数求系统极点：数求系统极点： cheby

15、shevchebyshev 1 1型低通滤波器的设计型低通滤波器的设计 2 22 2 /1 1 sg v jh pn 利用左半平面的极点构建因果稳定的系利用左半平面的极点构建因果稳定的系 统函数。统函数。 /11 1 1 22 2 n v m /11 /1 22 n n v v jh 将将chebyshevchebyshev 1 1型滤波器进行反转和变量型滤波器进行反转和变量 代换得到。代换得到。 模拟原型滤波器：模拟原型滤波器： chebyshev 2 阻带截止频率阻带截止频率 阻带为等纹波，通带单调下降；阻带为等纹波，通带单调下降； 1 s 2 2 1 s jh 模拟原型滤波器：模拟原型滤

16、波器： chebyshev 2 系统极点分布与系统极点分布与1 1型呈倒数关系；型呈倒数关系； 系统具有系统具有n n个零点分布在虚轴上；个零点分布在虚轴上； 模拟原型滤波器：模拟原型滤波器： chebyshev 2 z,p,kz,p,k=cheb2ap(n,rs) =cheb2ap(n,rs) 给出阻带波动为给出阻带波动为rsrs（dbdb）的）的n n阶阶2 2型型chebyshevchebyshev滤波器滤波器 的零点的零点z z，极点，极点p p和增益和增益k k； z,p,k=cheb2ap(10,25); b,a=zp2tf(z,p,k); h,w=freqs(b,a);zplan

17、e(b,a); figure(2),plot(w,abs(h),title(a); axis(0,4,0,1.2); chebyshevchebyshev 2 2型滤波器的仿真型滤波器的仿真 chebyshevchebyshev 2 2型滤波器的仿真型滤波器的仿真 为为chebyshevchebyshev型滤波器的综合形式，采用等纹波逼型滤波器的综合形式，采用等纹波逼 近设计；近设计； n n个极点分布在单位圆内，形成通带纹波；个极点分布在单位圆内，形成通带纹波； n n个零点分布在虚轴上，形成阻带纹波；个零点分布在虚轴上，形成阻带纹波； 22 2 1 1 n q m 22 1 1 n q j

18、h 模拟原型滤波器：模拟原型滤波器：elliptic z,p,kz,p,k=ellipap(n,rp,rsellipap(n,rp,rs) ) 给出通带波动给出通带波动rprp、阻带波动、阻带波动rsrs的的n n阶阶elliptic滤波器滤波器 的零点的零点z z，极点，极点p p和增益和增益k k； z,p,k=ellipap(6,0.5,25); b,a=zp2tf(z,p,k); h,w=freqs(b,a);zplane(b,a); figure(2),plot(w,abs(h); axis(0,4,0,1.2); elliptic滤波器的设计仿真滤波器的设计仿真 elliptic滤

19、波器的设计仿真滤波器的设计仿真 从原型到一般滤波器：变量代换从原型到一般滤波器：变量代换 原型滤波器：低通，截止频率为原型滤波器：低通，截止频率为1 1； 任何频率选择滤波器都可以通过频率变换由原型任何频率选择滤波器都可以通过频率变换由原型 滤波器得到：滤波器得到： 对频率响应函数，作变量代换：对频率响应函数，作变量代换： 对系统函数，作变量代换：对系统函数，作变量代换： f sfs 常用变换关系举例常用变换关系举例 原型到低通：原型到低通：p / 原型到高通：原型到高通： p ss/ / p ss p / 原型到带通：原型到带通： plpu pupl 2 原型到带阻：原型到带阻： 2 pup

20、l plpu 模拟滤波器的模拟滤波器的matlabmatlab设计仿真函数设计仿真函数 函数使用说明：函数使用说明： 用于频率选择滤波器的设计；用于频率选择滤波器的设计； 函数指定选用的原型；函数指定选用的原型； 需要设置滤波器类型和各截止频率；需要设置滤波器类型和各截止频率； 输出系统函数的分子输出系统函数的分子/ /分母多项式系数分母多项式系数b b和和a a。 设计仿真函数形式设计仿真函数形式 b,ab,a=butter(n,wp,sbutter(n,wp,s) ) n n阶阶butterworthbutterworth低通滤波器，低通滤波器，3db3db截止频率截止频率wpwp； ss代表模拟滤波器，若没有此标识则代表数字滤波代表模拟滤波器，若没有此标识则代表数字滤波 器；器； b,ab,a=butter(n,wpl,wpu,sbutter(n,wpl,wpu,s) 2n2n阶阶butterworthbutterworth带通滤波器带通滤波器 b,ab,a=butter(n,wp