无限脉冲响应数字滤波器的设计.pptx

《无限脉冲响应数字滤波器的设计.pptx》由会员分享，可在线阅读，更多相关《无限脉冲响应数字滤波器的设计.pptx（99页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、2023/2/12第1页三、标准理想滤波器的特性低通高通带通带阻01-01-0101第1页/共99页2023/2/12第2页通带截止频率阻带截止频率四、滤波器的技术指标过渡带00.7073dB通带截止频率第2页/共99页2023/2/12第3页如将|H(ej0)|归一化为1，上式则表示成：第3页/共99页2023/2/12第4页第二节 模拟滤波器的设计 模拟滤波器的理论和设计方法已发展得相当成熟，且有若干典型的模拟滤波器供我们选择，如巴特沃斯(Butterworth)滤波器、切比雪夫(Chebyshev)滤波器、椭圆(Cauer)滤波器、贝塞尔(Bessel)滤波器等。第4页/共99页2023

2、/2/12第5页模拟低通滤波器的设计指标有 ，其中 分别称为通带截止频率和阻带截止频率，是通带 中的最大衰减系数，是阻带 的最小衰减系数一、模拟低通滤波器设计指标及逼近方法第5页/共99页2023/2/12第6页滤波器的技术指标给定后，需要设计一个传输函数，希望其幅度平方函数满足给定的指标 和 ，一般滤波器的单位冲激响应为实数，因此幅度平方函数可表示成：第6页/共99页2023/2/12第7页 巴特沃斯滤波器是指在通带内的幅度特性具有最大平坦特性，是一全极点型滤波器，巴特沃斯低通滤波器的幅度平方函数用下式表示：巴特沃斯幅度特性和N的关系二、巴特沃斯低通滤波器的设计第7页/共99页2023/2/

3、12第8页将幅度平方函数写成s的函数：此式表明幅度平方函数有2N个极点，极点sk用下式表示：极点分布：2N个极点在S平面上是象限对称分布在半径为c的圆上；极点间的角度间隔为/N rad；极点决不会落在虚轴上；保证了系统的稳定性 N为奇数则实轴上有极点，N为偶数实轴上无极点。三阶巴特沃斯滤波器极点分布第8页/共99页2023/2/12第9页 为形成稳定的滤波器，2N个极点中只取s平面左半平面的N个极点构成Ha(s)，而右半平面的N个极点构成Ha(s)。Ha(s)的表示式为：设N=3，极点有6个，它们分别为取s平面左半平面的极点s0,s1,s2组成Ha(s)：第9页/共99页2023/2/12第1

4、0页 由于各滤波器的幅频特性不同，为使设计统一，将所有的频率归一化。这里采用对3dB截止频率c归一化，归一化后的Ha(s)表示为 式中，s/c=j/c。令=/c，称为归一化频率；令p=j，p称为归一化复变量，这样归一化巴特沃斯的传输函数为式中pk为归一化极点，用下式表示：第10页/共99页2023/2/12第11页第11页/共99页2023/2/12第12页令,则N可表示为：上式确定的N可能有小数部分，应取大于等于N的最小整数。关于3dB截止频率c，如果技术指标中没有给出，可以由下面推导得到：第12页/共99页2023/2/12第13页(1)根据技术指标p,p,s和s，求出滤波器的阶数N。(2

5、)求出归一化极点pk，得到归一化传输函数Ha(p)。(3)将Ha(p)去归一化。将p=s/c代入Ha(p)，得到实际的滤波器传输函数Ha(s)。低通巴特沃斯滤波器的设计步骤如下：第13页/共99页2023/2/12第14页巴特沃斯归一化低通滤波器参数 第14页/共99页2023/2/12第15页第15页/共99页2023/2/12第16页第16页/共99页2023/2/12第17页 例 已知通带截止频率fp=5kHz，通带最大衰减p=2dB，阻带截止频率fs=12kHz，阻带最小衰减s=30dB，按以上技术指标设计巴特沃斯低通滤波器。解(1)确定阶数N。第17页/共99页2023/2/12第1

6、8页 (2)其极点为归一化传输函数为 上式分母可以展开成为五阶多项式，或者将共轭极点放在一起，形成因式分解形式。直接查表更简单，由N=5，直 接 查 表 得 到 极 点：-0.3090j0.9511,-0.8090j0.5878;-1.0000式 b0=1.0000,b1=3.2361,b2=5.2361,b3=5.2361,b4=3.2361第18页/共99页2023/2/12第19页 (3)为将Ha(p)去归一化，先求3dB截止频率c：将p=s/c代入Ha(p)中得到：第19页/共99页2023/2/12第20页切比雪夫型滤波器幅度平方函数可表示为：三、切比雪夫型滤波器的设计切比雪夫型滤波

7、器幅频特性 第20页/共99页2023/2/12第21页式中为小于1的正数，表示通带内幅度波动的程度，愈大，波动幅度也愈大。p称为通带截止频率。CN(x)称为N阶切比雪夫多项式，定义为当N=0时，C0(x)=1；当N=1时，C1(x)=x；当N=2时，C2(x)=2x 21；当N=3时，C3(x)=4x 3 3x。由此可归纳出高阶切比雪夫多项式的递推公式为C N+1(x)=2xCN(x)C N-1(x)第21页/共99页2023/2/12第22页 下图为N=0,4,5时切比雪夫多项式特性，由图可见：(1)切比雪夫多项式的过零点在|x|1的范围内；(2)当|x|1时，|CN(x)|1,在|x|1

8、时，CN(x)是双曲线函数，随x单调上升。第22页/共99页2023/2/12第23页 平方幅度函数与三个参数即,p和N有关。其中与通带内允许的波动大小有关，定义允许的通带波纹用下式表示：切比雪夫型与巴特沃斯低通曲线第23页/共99页2023/2/12第24页 设阻带的起始点频率(阻带截止频率)用s表示，在s处的A2(s)为：令s=s/p，由s1，有可以解出 第24页/共99页2023/2/12第25页 3dB截止频率用c表示通常取c1，因此上式中仅取正号，得到3dB截止频率计算公式：以上p,和N确定后，可以求出滤波器的极点，并确定Ha(p)，p=s/p。第25页/共99页2023/2/12第

9、26页 设Ha(s)的极点为si=i+ji，可以证明：上式是一个椭圆方程，长半轴为pch(在虚轴上)，短半轴为psh(在实轴上)。令bp和ap分别表示长半轴和短半轴，可推导出：第26页/共99页2023/2/12第27页设N=3，平方幅度函数的极点分布如图所示。为稳定，用左半平面的极点构成Ha(p)，即由上式可导出：c=2 N-1，从而得归一化的传输函数为：去归一化后的传输函数为第27页/共99页2023/2/12第28页1)确定技术要求p,p,s和p是=p时的衰减系数，s是=s时的衰减系数：这里p就是前面定义的通带波纹。归一化频率切比雪夫切比雪夫型滤波器设计步骤：型滤波器设计步骤：第28页/

10、共99页2023/2/12第29页2)求滤波器阶数N和参数 由幅度平方函数可得：第29页/共99页2023/2/12第30页3)求归一化传输函数Ha(p)为求Ha(p)，先求出归一化极点pk,k=1,2,:,N。根据定义可得：将极点pk代入即得：4)将Ha(p)去归一化，得到实际的Ha(s)，即 其中p=第30页/共99页2023/2/12第31页 例 设计低通切比雪夫滤波器，要求通带截止频率fp=3kHz，通带最大衰减p=0.1dB，阻带截止频率fs=12kHz，阻带最小衰减s=60dB。(2)求阶数N和：解：(1)滤波器的技术要求：第31页/共99页2023/2/12第32页 (3)求Ha

11、(p):求出N=5时的极点pi，代入上式，得到：(4)将Ha(p)去归一化，得到：第32页/共99页2023/2/12第33页1)低通到高通的频率变换 和之间的关系为上式即是低通到高通的频率变换公式，如果已知低通G(j)，高通H(j)则用下式转换：四、模拟滤波器的频率变换 模拟高通、带通、带阻滤波器的设计第33页/共99页2023/2/12第34页 模拟高通滤波器的设计步骤：(1)由高通滤波器的技术指标（通带下限频率p，阻带上限频率s，通带最大衰减p，阻带最小衰减s）确定相应低通滤波器的设计指标：低通滤波器通带截止频率p=1/p；低通滤波器阻带截止频率s=1/s；通带最大衰减仍为p，阻带最小衰

12、减仍为s。(2)设计归一化低通滤波器G(p)。(3)求模拟高通的H(s)。将G(p)转换成归一化高通H(q)，为去归一化，将q=s/c代入H(q)中，得 第34页/共99页2023/2/12第35页 解 高通技术要求：fp=200Hz,p=3dB;fs=100Hz,s=15dB 归一化频率：低通技术要求：例 设计高通滤波器：fp=200Hz,fs=100Hz，幅度特性单调下降，fp处最大衰减为3dB，阻带最小衰减s=15dB。第35页/共99页2023/2/12第36页 设计归一化低通巴特沃斯滤波器G(p)求模拟高通H(s)：第36页/共99页2023/2/12第37页2)低通到带通的频率变换

13、 低通与带通滤波器的幅度特性如图所示带通与低通滤波器的幅度特性 第37页/共99页2023/2/12第38页由对应关系可得与的关系为：由上表及上式有：上式称为低通到带通的频率变换公式。利用该式将带通的边界频率转换成低通的边界频率。第38页/共99页2023/2/12第39页 1)确定模拟带通滤波器的技术指标：带通上限频率u，带通下限频率l 下阻带上限频率 s1,上阻带下限频率 s2 通带中心频率20=lu，通带宽度B=ul 与以上边界频率对应的归一化边界频率如下：模拟带通的设计步骤模拟带通的设计步骤第39页/共99页2023/2/12第40页 (2)确定归一化低通技术要求：s与-s的绝对值可能

14、不相等，一般取绝对值小的s，这样保证在较大的s处更能满足要求。通带最大衰减仍为p，阻带最小衰减亦为s。(3)设计归一化低通G(p)。(4)将G(p)转换成带通H(s)：第40页/共99页2023/2/12第41页 例 设计模拟带通滤波器，通带带宽B=2200rad/s，中心频率0=21000rad/s，通带内最大衰减p=3dB，阻带s1=2830rad/s,s2=21200rad/s，阻带最小衰减s=15dB。解：1)模拟带通的技术要求：0=21000rad/s,p=3dB s1=2830rad/s,s2=21200rad/s,s=15dB B=2200rad/s;0=5,s1=4.15,s2

15、=6第41页/共99页2023/2/12第42页 2)模拟归一化低通技术要求：取s=1.833,p=3dB,s=15dB。3)设计模拟归一化巴特沃斯低通滤波器G(p)：第42页/共99页2023/2/12第43页 取N=3，查表得：(4)求模拟带通H(s)：第43页/共99页2023/2/12第44页 3)低通到带阻的变换 低通与带阻滤波器的幅频特性如图所示。l和u分别是下通带和上通带截止频率，s1和s2分别为阻带的下限频率和上限频率，0为阻带中心频率，20=ul，阻带带宽B=ul，B作为归一化参考频率。相应的归一化边界频率为：u=u/B,l=l/B,s1=s1/B,s2=s2/B;20=ul

16、 第44页/共99页2023/2/12第45页 根据与的对应关系，可得到：第45页/共99页2023/2/12第46页1)确定模拟带阻滤波器的技术要求，即：下通带截止频率l,上通带截止频率u阻带下限频率s1，阻带上限频率s2阻带中心频率 ，阻带宽度B=ul它们相应的归一化边界频率为 l=l/B,u=u/B,s1=s1/B;s2=s2/B,20=ul以及通带最大衰减p和阻带最小衰减s。带阻滤波器设计步骤第46页/共99页2023/2/12第47页 (2)确定归一化模拟低通技术要求，即：取s和s的绝对值较小的s；通带最大衰减为p，阻带最小衰减为s。(3)设计归一化模拟低通G(p)。(4)将G(p)

17、转换成带阻滤波器H(s)。第47页/共99页2023/2/12第48页 例 设计模拟带阻滤波器，其技术要求为：l=2905rad/s,s1=2980rad/s,s2=21020rad/s,u=21105rad/s,p=3dB,s=25dB。试设计巴特沃斯带阻滤波器。解 (1)模拟带阻滤波器的技术要求：l=2905,u=21105;s1=2980,s2=21020;20=lu=4+21000025,B=ul=2200;l=l/B=4.525,u=u/B=5.525;s1=s1/B=4.9,s2=5.1;20=lu=25第48页/共99页2023/2/12第49页 (2)归一化低通的技术要求：(3

18、)设计归一化低通滤波器G(p):(4)带阻滤波器的H(s)为第49页/共99页2023/2/12第50页第三节 脉冲响应不变法设计IIR数字低通滤波器 数字滤波器的设计可先按设计要求设计一个模拟滤波器，再按一定转换关系将模拟滤波器转换成数字滤波器。为了保证转换后的H(z)稳定且满足技术要求，对转换关系提出两点要求：1)因果稳定的模拟滤波器转换成数字滤波器，仍是因果稳定的。2)数字滤波器的频率响应模仿模拟滤波器，s平面的虚轴映射z平面的单位圆，相应的频率之间成线性关系。第50页/共99页2023/2/12第51页 设模拟滤波器的传输函数为Ha(s),相应的单位冲激响应是ha(t)设模拟滤波器Ha

19、(s)只有单阶极点，且分母多项式的阶次高于分子多项式的阶次，将Ha(s)用部分分式表示：式中si为Ha(s)的单阶极点。将Ha(s)进行逆拉氏变换得：对ha(t)进行等间隔采样，采样间隔为T，得到：第51页/共99页2023/2/12第52页对h(n)进行Z变换，得到数字滤波器的系统函数H(z)：设ha(t)的采样信号用 表示：对 进行拉氏变换，得到:上式表示采样信号的拉氏变换与相应的序列的Z变换之间的映射关系可用下式表示：第52页/共99页2023/2/12第53页 模拟信号ha(t)的傅里叶变换Ha(j)和其采样信号的傅里叶变换之间的关系如下：将s=j代入上式，得上式表明将模拟信号ha(t

20、)的拉氏变换在s平面上沿虚轴按照周期s=2/T延拓后，再按 映射关系映射到z平面上，就得到H(z)。可称为标准映射关系。第53页/共99页2023/2/12第54页下面进一步分析这种映射关系。设得到:因此得到:那么 =0,r=1 0,r0,r1 另外，注意到z=esT是一个周期函数，可写成为任意整数第54页/共99页2023/2/12第55页z=esT,s平面与z平面之间的映射关系频率混叠现象第55页/共99页2023/2/12第56页 假设 没有频率混叠现象，即满足由 ，并将关系式s=j代入，=T，代入得到：为避免采样间隔较小而导致 H(z)过大，令：则：第56页/共99页2023/2/12

21、第57页 如果模拟滤波器的二阶基本节的形式为极点为 可以推导出相应的数字滤波器二阶基本节为：如果模拟滤波器二阶基本节的形式为 极点为 相应的数字滤波器二阶基本节为：一般Ha(s)的极点是复数且共轭成对的形式出现，将一对复数共轭极点放在一起，形成一个二阶基本节。第57页/共99页2023/2/12第58页 例 已知模拟滤波器的传输函数Ha(s)为用脉冲响应不变法将其转换成数字滤波器的系统函数。解 首先将Ha(s)写成部分分式：极点为：那么H(z)的极点为：由上式并经整理得：第58页/共99页2023/2/12第59页也可以首先将Ha(s)写成如下形式，如极点s1，2=1j1,则再由上分析得H(z

22、)为：T=1sT=0.1s第59页/共99页2023/2/12第60页脉冲响应不变法的特点：优点：T保证了模拟角频率与数字角频率之间的线性关系缺点：在 附近产生混叠失真所以，一般只适用于带限滤波器，高通、带阻滤波器不能直接进行变换所以，脉冲响应不变法适用情况：带限滤波器H(s)要能用部分分式展开成单阶极点第60页/共99页2023/2/12第61页双线性变换的含义n第一次映射:将整个S平面压缩到S1平面的一条横带内（-/T,+/T）n第二次映射:将S1平面映射到Z平面双线性映射的变换式n第一次映射：正切变换实现频率压缩n第二次映射:双线性变换法的公式：第四节 双线性变换法设计IIR数字低通滤波

23、器第61页/共99页2023/2/12第62页S平面S1平面Z平面第62页/共99页2023/2/12第63页双线性变换法的频率变换关系双线性变换法幅度和相位特性非线性映射第63页/共99页2023/2/12第64页系数关系表 第64页/共99页2023/2/12第65页第65页/共99页2023/2/12第66页 例 试分别用脉冲响应不变法和双线性不变法将下图所示的RC低通滤波器转换成数字滤波器。解 首先按照图写出该滤波器的传输函数Ha(s)为 利用脉冲响应不变法转换，数字滤波器的系统函数为 利用双线性变换法转换，数字滤波器的系统函数H2(z)为第66页/共99页2023/2/12第67页H

24、1(z)和H2(z)的网络结构分别如图所示数字滤波器H1(z)和H2(z)的幅频特性 第67页/共99页2023/2/12第68页 利用模拟滤波器设计IIR数字低通滤波器的步骤1)确定数字低通滤波器的技术指标：通带截止频率p、通带衰减p、阻带截止频率s、阻带衰减s。2)将数字低通滤波器的技术指标转换成模拟低通滤波器的技术指标。如果采用双线性变换法，边界频率的转换关系为 3)按照模拟低通滤波器的技术指标设计模拟低通滤波器。4)将模拟滤波器Ha(s)，从s平面转换到z平面，得到数字低通滤波器系统函数H(z)。第68页/共99页2023/2/12第69页 解 (1)用脉冲响应不变法设计数字低通滤波器

25、。数字低通的技术指标为 p=0.2rad,p=1dB;s=0.3rad,s=15dB 模拟低通的技术指标为 T=1s,p=0.2rad/s,p=1dB;s=0.3rad/s,s=15dB例 设计低通数字滤波器，要求在通带内频率低于0.2时，幅度误差不大于1dB；在频率0.3到之间的阻带衰减大于15dB。指定采用巴特沃斯低通滤波器。试分别用脉冲响应不变法和双线性变换法设计滤波器。第69页/共99页2023/2/12第70页 设计巴特沃斯低通滤波器，先计算阶数N及3dB截止频率c。取N=6。由p和p求得3dB截止频率c=0.7032rad/s。根据阶数N=6，查表得归一化传输函数为第70页/共99

26、页2023/2/12第71页用脉冲响应不变法将Ha(s)转换成H(z)。首先将Ha(s)进行部分分式，并按脉冲响应不变法得到：去归一化：将p=s/c代入Ha(p)中，得实际传输函数：第71页/共99页2023/2/12第72页 (2)用双线性变换法设计数字低通滤波器。数字低通技术指标仍为 p=0.2rad,p=1dB;s=0.3rad,s=15dB 模拟低通的技术指标为第72页/共99页2023/2/12第73页 设计巴特沃斯低通滤波器。阶数N计算如下：取N=6。由s和s求c=0.7662rad/s。查表得归一化传输函数Ha(p)，去归一化，将p=s/c代入Ha(p)，得实际的Ha(s)第73

27、页/共99页2023/2/12第74页 用双线性变换法将Ha(s)转换成数字滤波器H(z)：双线性变换法的幅度特性脉冲响应不变法的幅度特性第74页/共99页2023/2/12第75页第五节 数字高通、带通和带阻滤波器设计具体设计步骤如下：1)确定所需类型数字滤波器的技术指标。2)将数字滤波器技术指标转换成相应的模拟滤波器技术指标，转换公式为3)将所需模拟滤波器指标转换成模拟低通滤波器技术指标4)设计模拟低通滤波器。5)将模拟低通转换成所需类型的模拟滤波器。6)采用双线性变换法，将所需类型的模拟滤波器转换成所需类型的数字滤波器。第75页/共99页2023/2/12第76页 例 设计一个数字高通滤

28、波器，要求通带截止频率p=0.8rad，通带衰减不大于3dB，阻带截止频率s=0.44rad,阻带衰减不小于15dB。希望采用巴特沃斯型滤波器。解 (1)数字高通的技术指标为 p=0.8rad,p=3dB;s=0.44rad,s=15dB (2)模拟高通的技术指标计算如下：令T=1，则有第76页/共99页2023/2/12第77页(3)模拟低通滤波器的技术指标计算如下：将p和s对3dB截止频率c归一化，这里c=p,(4)设计归一化模拟低通滤波器G(p)。模拟低通滤波器的阶数N计算如下：第77页/共99页2023/2/12第78页 查表得到归一化模拟低通传输函数G(p)为 去归一化，将p=s/c

29、代入上式得到：(5)将模拟低通转换成模拟高通。将上式中G(s)的变量换成1/s，得到模拟高通Ha(s)：(6)用双线性变换法将模拟高通转换成数字高通H(z)：实际上(5)、(6)两步可合并成一步，即第78页/共99页2023/2/12第79页 例 设计一个数字带通滤波器，通带范围为0.3rad到0.4rad，通带内最大衰减为3dB，0.2rad以下和0.5rad以上为阻带，阻带内最小衰减为18dB。采用巴特沃斯型模拟低通滤波器。解 (1)数字带通滤波器技术指标为 通带上截止频率 u=0.4rad 通带下截止频率 l=0.3rad 阻带上截止频率 s2=0.5rad 阻带下截止频率 s1=0.2

30、rad 通带内最大衰减p=3dB，阻带内最小衰减s=18dB。第79页/共99页2023/2/12第80页 (2)模拟带通滤波器技术指标如下：设T=1，则有(通带中心频率)(带宽)将以上边界频率对带宽B归一化，得到 u=3.348,l=2.348;s2=4.608,s1=1.498;0=2.804第80页/共99页2023/2/12第81页 (3)模拟归一化低通滤波器技术指标：归一化阻带截止频率归一化通带截止频率p=1p=3dB,s=18dB (4)设计模拟低通滤波器：查表得归一化低通传输函数G(p)第81页/共99页2023/2/12第82页 (5)将归一化模拟低通转换成模拟带通：(6)通过

31、双线性变换法将Ha(s)转换成数字带通滤波器H(z)。下面将(5)、(6)两步合成一步计算：将上式代入(5)中的转换公式，得将上面的p等式代入G(p)中，得 第82页/共99页2023/2/12第83页 例 设计一个数字带阻滤波器，通带下限频率l=0.19,阻带下截止频率s1=0.198，阻带上截止频率s2=0.202，通带上限频率u=0.21，阻带最小衰减s=13dB，l和u处衰减p=3dB。采用巴特沃斯型。解 (1)数字带阻滤波器技术指标：l=0.19rad,u=0.21rad,p=3dB;s1=0.198rad,s2=0.202rad,s=13dB (2)模拟带阻滤波器的技术指标：设T=

32、1，则有第83页/共99页2023/2/12第84页阻带中心频率平方为 20=lu=0.421阻带带宽为 B=u-l=0.07rad/s 将以上边界频率对B归一化：l=8.786,u=9.786,s1=9.186,s2=9.386;20=lu=85.98 (3)模拟归一化低通滤波器的技术指标：p=1,p=3dB第84页/共99页2023/2/12第85页 (4)设计模拟低通滤波器：(5)将G(p)转换成模拟阻带滤波器Ha(s)：(6)将Ha(s)通过双线性变换，得到数字阻带滤波器H(z)。第85页/共99页2023/2/12第86页第六节 IIR 数字滤波器的直接设计法 1.零极点累试法称为零

33、极点累试法。在确定零极点位置时要注意：1)极点必须位于单位圆内，保证数字滤波器因果稳定；2)复数零极点必须共轭成对，保证系统函数有理式的系数是实的。a)零极点分布；b)幅度特性 第86页/共99页2023/2/12第87页 2.频域幅度平方误差最小法设计IIR数字滤波器 设IIR滤波器由K个二阶网络级联而成，系统函数用H(z)表示，式中A是常数；ai,bi,ci,di是待求的系数；Hd(e j)是希望设计的滤波器频响。如果在(0，)区间取N点数字频率i,i=1,2,:,N，在这N点频率上，比较|Hd(e j)|和|H(e j)|，写出两者的幅度平方误差E为第一式中共有(4K+1)个待定的系数，

34、求它们的原则是使E最小。下面研究采用第一式网络结构，如何求出(4K+1)系数。第87页/共99页2023/2/12第88页由上式可知E是(4K+1)个未知数的函数，用下式表示：上式表示4K个系数组成的系数向量。为推导公式方便，令 为选择A使E最小，令第88页/共99页2023/2/12第89页 设k是 的第k个分量(ak或bk或ck或dk)，因为 ,式中H*i表示对Hi函数共轭。将上式具体写成对ak,bk,ck,dk的偏导，得到：式中，k=1,2,3,:,K；i=1,2,3,:,N。第89页/共99页2023/2/12第90页同理求得 由于系统函数是一个有理函数，极、零点均以共轭成对的形式存在

35、，对于极点z1,一定有下面关系：第90页/共99页2023/2/12第91页 例 设计低通数字滤波器，其幅度特性如图所示。截止频率s=0.1rad。解 考虑到通带和过渡带的重要，在00.2区间，每隔0.01取一点i值，在0.2区间每隔0.1取一点i值，并增加一点过渡带，在=0.1处|Hd(e j)|=0.5。1.0,=0,0.01,0.02,:,0.09 0.5,=0.1 0.0,=0.11,0.12,:,0.19 0.0,=0.2,0.3,:,第91页/共99页2023/2/12第92页 待求的参数是A，a1,b1,c1,d1。设初始值=(0000.25)T经过90次迭代，求得E=1.261

36、1，系统函数零、极点位置为 零点0.67834430j0.73474418;极点0.75677793j1.3213916 为使滤波器因果稳定，将极点按其倒数搬入单位圆内，再进行62次优化迭代，求得结果为 零点0.82191163j0.56961501;极点0.89176390j0.19181084;Ag=0.11733978,E=0.56731 N=29，取k=1，系统函数为第92页/共99页2023/2/12第93页 3.在时域直接设计IIR数字滤波器 设我们希望设计的IIR数字滤波器的单位脉冲响应为hd(n)，要求设计一个单位脉冲响应h(n)充分逼近hd(n)。下面我们介绍这种设计方法。设

37、滤波器是因果性的，系统函数为式中a0=1，未知系数ai和bi共有N+M+1个，取h(n)的一段，0np-1，使其充分逼近hd(n)，用此原则求解M+N+1个系数。第93页/共99页2023/2/12第94页令p=M+N+1,则 令上面等式两边z的同幂次项系数相等，可得N+M+1个方程：h(0)=b0 h(0)a1+h(1)=b1 h(0)a2+h(1)a1+h(2)=b2 上式表明h(n)是系数ai,bi的非线性函数第94页/共99页2023/2/12第95页考虑到iM时，bi=0，一般表达式为：设x(n)为给定的输入信号，yd(n)是相应的希望的输出信号，x(n)和yd(n)长度分别为M和N

38、，实际滤波器的输出为y(n)，下面按照y(n)和yd(n)的最小均方误差求解滤波器的最佳解，设均方误差用E表示：上式中x(n),0nM1；yd(n)，0nN-1 为选择h(n)使E最小，令第95页/共99页2023/2/12第96页 由上式得到第96页/共99页2023/2/12第97页 例 设 计 数 字 滤 波 器，要 求 在 给 定 输 入 x(n)=3,1的 情 况 下，输 出yd(n)=1,0.25,0.1,0.01,0。解 设h(n)长度为p=4列出方程：10h(0)+3h(1)=3.25 3h(0)+10h(1)+3h(2)=0.85 3h(1)+10h(2)+3h(3)=0.3

39、1 3h(2)+9h(3)=0.03 第97页/共99页2023/2/12第98页 解联立方程，得 h(n)=0.3333,0.0278,0.0426,0.0109将h(n)以及M=1,N=2代入(6.6.15),(6.6.16)式中，得 a1=0.1824,a2=0.1126 b0=0.3333,b1=0.0330 滤波器的系统函数为相应的差分方程为 y(n)=0.3333x(n)+0.0330 x(n1)0.1824y(n1)+0.1126y(n2)当x(n)=3,1时，输出y(n)为 y(n)=0.9999,0.2499,0.1,0.0099,0.0095,0.0006,0.0012:将y(n)与给定yd(n)比较，y(n)的前五项与yd(n)的前五项很相近，y(n)在五项以后幅度值很小。第98页/共99页2023/2/12第99页感谢您的观看！第99页/共99页