第六章无限脉冲响应数字滤波器的设计课件

2024/7/17第第1页页一、数字滤波器一、数字滤波器 是指输入、输出均为数字信号，通过一定运算关系改变输入信是指输入、输出均为数字信号，通过一定运算关系改变输入信号所含频率成分的相对比列或者虑除某些频率成分的器件。号所含频率成分的相对比列或者虑除某些频率成分的器件。二、数字滤波器的分类二、数字滤波器的分类1.1.经典滤波器经典滤波器n从功能上分：低通、高通、带通、带限从功能上分：低通、高通、带通、带限n从实现的网络结构或单位脉冲响应：从实现的网络结构或单位脉冲响应：IIRIIR、FIRFIR2.2.现代滤波器现代滤波器维纳滤波器、卡尔曼滤波器、自适应滤波器等最佳滤波器。维纳滤波器、卡尔曼滤波器、自适应滤波器等最佳滤波器。第一节第一节 数字滤波器的基本概念数字滤波器的基本概念2024/7/17第第2页页三、标准理想滤波器的特性三、标准理想滤波器的特性低通低通高通高通带通带通带阻带阻01-01-01012024/7/17第第3页页通带截止频率通带截止频率阻带截止频率阻带截止频率四、滤波器的技术指标四、滤波器的技术指标过渡带过渡带00.7073dB通带截止频率通带截止频率2024/7/17第第4页页如将|H(ej0)|归一化为1，上式则表示成：2024/7/17第第5页页第二节第二节 模拟滤波器的设计模拟滤波器的设计 模拟滤波器的理论和设计方法已发展得相当成熟，且有若干典型的模拟滤波器供我们选择，如巴特沃斯(Butterworth)滤波器、切 比 雪 夫(Chebyshev)滤 波 器、椭 圆(Cauer)滤波器、贝塞尔(Bessel)滤波器等。2024/7/17第第6页页模拟低通滤波器的设计指标有 ，其中 分别称为通带截止频率和阻带截止频率，是通带 中的最大衰减系数，是阻带 的最小衰减系数一、模拟低通滤波器设计指标及逼近方法一、模拟低通滤波器设计指标及逼近方法2024/7/17第第7页页滤波器的技术指标给定后，需要设计一个传输函数，希望其幅度平方函数满足给定的指标 和 ，一般滤波器的单位冲激响应为实数，因此幅度平方函数可表示成：2024/7/17第第8页页 巴巴特特沃沃斯斯滤滤波波器器是是指指在在通通带带内内的的幅幅度度特特性性具具有有最最大大平平坦坦特特性性，是是一一全全极极点点型型滤滤波波器器，巴巴特特沃沃斯斯低低通通滤滤波波器器的的幅幅度度平平方方函函数数用下式表示：用下式表示：巴特沃斯幅度特性和N的关系二、巴特沃斯低通滤波器的设计二、巴特沃斯低通滤波器的设计2024/7/17第第9页页将幅度平方函数写成s的函数：此式表明幅度平方函数有2N个极点，极点sk用下式表示：极点分布：2N个极点在S平面上是象限对称分布在半径为c的圆上；极点间的角度间隔为/N rad；极点决不会落在虚轴上；保证了系统的稳定性 N为奇数则实轴上有极点，N为偶数实轴上无极点。三阶巴特沃斯滤波器极点分布2024/7/17第第10页页 为形成稳定的滤波器，2N个极点中只取s平面左半平面的N个极点构成Ha(s)，而右半平面的N个极点构成Ha(s)。Ha(s)的表示式为：设N=3，极点有6个，它们分别为取s平面左半平面的极点s0,s1,s2组成Ha(s)：2024/7/17第第11页页 由于各滤波器的幅频特性不同，为使设计统一，将所有的频率归一化。这里采用对3dB截止频率c归一化，归一化后的Ha(s)表示为 式中，s/c=j/c。令=/c，称为归一化频率；令p=j，p称为归一化复变量，这样归一化巴特沃斯的传输函数为式中pk为归一化极点，用下式表示：2024/7/17第第12页页2024/7/17第第13页页令,则N可表示为：上式确定的N可能有小数部分，应取大于等于N的最小整数。关于3dB截止频率c，如果技术指标中没有给出，可以由下面推导得到：2024/7/17第第14页页(1)根据技术指标p,p,s和s，求出滤波器的阶数N。(2)求出归一化极点pk，得到归一化传输函数Ha(p)。(3)将Ha(p)去归一化。将p=s/c代入Ha(p)，得到实际的滤波器传输函数Ha(s)。低通巴特沃斯滤波器的设计步骤如下：低通巴特沃斯滤波器的设计步骤如下：2024/7/17第第15页页巴特沃斯归一化低通滤波器参数 2024/7/17第第16页页2024/7/17第第17页页2024/7/17第第18页页 例 已知通带截止频率fp=5kHz，通带最大衰减p=2dB，阻带截止频率fs=12kHz，阻带最小衰减s=30dB，按以上技术指标设计巴特沃斯低通滤波器。解(1)确定阶数N。2024/7/17第第19页页 (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.23612024/7/17第第20页页 (3)为将Ha(p)去归一化，先求3dB截止频率c：将p=s/c代入Ha(p)中得到：2024/7/17第第21页页切比雪夫型滤波器幅度平方函数可表示为：三、切比雪夫型滤波器的设计切比雪夫型滤波器幅频特性 2024/7/17第第22页页式中为小于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)2024/7/17第第23页页 下图为N=0,4,5时切比雪夫多项式特性，由图可见：(1)切比雪夫多项式的过零点在|x|1的范围内；(2)当|x|1时，|CN(x)|1,在|x|1时，CN(x)是双曲线函数，随x单调上升。2024/7/17第第24页页 平方幅度函数与三个参数即,p和N有关。其中与通带内允许的波动大小有关，定义允许的通带波纹用下式表示：切比雪夫型与巴特沃斯低通曲线2024/7/17第第25页页 设阻带的起始点频率(阻带截止频率)用s表示，在s处的A2(s)为：令s=s/p，由s1，有可以解出 2024/7/17第第26页页 3dB截止频率用c表示通常取c1，因此上式中仅取正号，得到3dB截止频率计算公式：以上p,和N确定后，可以求出滤波器的极点，并确定Ha(p)，p=s/p。2024/7/17第第27页页 设Ha(s)的极点为si=i+ji，可以证明：上式是一个椭圆方程，长半轴为pch(在虚轴上)，短半轴为psh(在实轴上)。令bp和ap分别表示长半轴和短半轴，可推导出：2024/7/17第第28页页设N=3，平方幅度函数的极点分布如图所示。为稳定，用左半平面的极点构成Ha(p)，即由上式可导出：c=2 N-1，从而得归一化的传输函数为：去归一化后的传输函数为2024/7/17第第29页页1)确定技术要求p,p,s和p是=p时的衰减系数，s是=s时的衰减系数：这里p就是前面定义的通带波纹。归一化频率切比雪夫切比雪夫型滤波器设计步骤：型滤波器设计步骤：2024/7/17第第30页页2)求滤波器阶数N和参数 由幅度平方函数可得：2024/7/17第第31页页3)求归一化传输函数Ha(p)为求Ha(p)，先求出归一化极点pk,k=1,2,:,N。根据定义可得：将极点pk代入即得：4)将Ha(p)去归一化，得到实际的Ha(s)，即 其中p=2024/7/17第第32页页 例 设计低通切比雪夫滤波器，要求通带截止频率fp=3kHz，通带最大衰减p=0.1dB，阻带截止频率fs=12kHz，阻带最小衰减s=60dB。(2)求阶数N和：解：(1)滤波器的技术要求：2024/7/17第第33页页 (3)求Ha(p):求出N=5时的极点pi，代入上式，得到：(4)将Ha(p)去归一化，得到：2024/7/17第第34页页1)低通到高通的频率变换 和之间的关系为上式即是低通到高通的频率变换公式，如果已知低通G(j)，高通H(j)则用下式转换：四、模拟滤波器的频率变换 模拟高通、带通、带阻滤波器的设计2024/7/17第第35页页 模拟高通滤波器的设计步骤：(1)由高通滤波器的技术指标（通带下限频率p，阻带上限频率s，通带最大衰减p，阻带最小衰减s）确定相应低通滤波器的设计指标：低通滤波器通带截止频率p=1/p；低通滤波器阻带截止频率s=1/s；通带最大衰减仍为p，阻带最小衰减仍为s。(2)设计归一化低通滤波器G(p)。(3)求模拟高通的H(s)。将G(p)转换成归一化高通H(q)，为去归一化，将q=s/c代入H(q)中，得 2024/7/17第第36页页 解 高通技术要求：fp=200Hz,p=3dB;fs=100Hz,s=15dB 归一化频率：低通技术要求：例 设计高通滤波器：fp=200Hz,fs=100Hz，幅度特性单调下降，fp处最大衰减为3dB，阻带最小衰减s=15dB。2024/7/17第第37页页 设计归一化低通巴特沃斯滤波器G(p)求模拟高通H(s)：2024/7/17第第38页页2)低通到带通的频率变换 低通与带通滤波器的幅度特性如图所示带通与低通滤波器的幅度特性 2024/7/17第第39页页由对应关系可得与的关系为：由上表及上式有：上式称为低通到带通的频率变换公式。利用该式将带通的边界频率转换成低通的边界频率。2024/7/17第第40页页 1)确定模拟带通滤波器的技术指标：带通上限频率u，带通下限频率l 下阻带上限频率 s1,上阻带下限频率 s2 通带中心频率20=lu，通带宽度B=ul 与以上边界频率对应的归一化边界频率如下：模拟带通的设计步骤模拟带通的设计步骤2024/7/17第第41页页 (2)确定归一化低通技术要求：s与-s的绝对值可能不相等，一般取绝对值小的s，这样保证在较大的s处更能满足要求。通带最大衰减仍为p，阻带最小衰减亦为s。(3)设计归一化低通G(p)。(4)将G(p)转换成带通H(s)：2024/7/17第第42页页 例 设计模拟带通滤波器，通带带宽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=62024/7/17第第43页页 2)模拟归一化低通技术要求：取s=1.833,p=3dB,s=15dB。3)设计模拟归一化巴特沃斯低通滤波器G(p)：2024/7/17第第44页页 取N=3，查表得：(4)求模拟带通H(s)：2024/7/17第第45页页 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 2024/7/17第第46页页 根据与的对应关系，可得到：2024/7/17第第47页页1)确定模拟带阻滤波器的技术要求，即：下通带截止频率l,上通带截止频率u阻带下限频率s1，阻带上限频率s2阻带中心频率 ，阻带宽度B=ul它们相应的归一化边界频率为 l=l/B,u=u/B,s1=s1/B;s2=s2/B,20=ul以及通带最大衰减p和阻带最小衰减s。带阻滤波器设计步骤2024/7/17第第48页页 (2)确定归一化模拟低通技术要求，即：取s和s的绝对值较小的s；通带最大衰减为p，阻带最小衰减为s。(3)设计归一化模拟低通G(p)。(4)将G(p)转换成带阻滤波器H(s)。2024/7/17第第49页页 例 设计模拟带阻滤波器，其技术要求为：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=252024/7/17第第50页页 (2)归一化低通的技术要求：(3)设计归一化低通滤波器G(p):(4)带阻滤波器的H(s)为2024/7/17第第51页页第三节第三节 脉冲响应不变法设计脉冲响应不变法设计IIR数字低通滤波器数字低通滤波器 数字滤波器的设计可先按设计要求设计一个模拟滤波器，再按一定转换关系将模拟滤波器转换成数字滤波器。为了保证转换后的H(z)稳定且满足技术要求，对转换关系提出两点要求：1)因果稳定的模拟滤波器转换成数字滤波器，仍是因果稳定的。2)数字滤波器的频率响应模仿模拟滤波器，s平面的虚轴映射z平面的单位圆，相应的频率之间成线性关系。2024/7/17第第52页页 设模拟滤波器的传输函数为Ha(s),相应的单位冲激响应是ha(t)设模拟滤波器Ha(s)只有单阶极点，且分母多项式的阶次高于分子多项式的阶次，将Ha(s)用部分分式表示：式中si为Ha(s)的单阶极点。将Ha(s)进行逆拉氏变换得：对ha(t)进行等间隔采样，采样间隔为T，得到：2024/7/17第第53页页对h(n)进行Z变换，得到数字滤波器的系统函数H(z)：设ha(t)的采样信号用 表示：对 进行拉氏变换，得到:上式表示采样信号的拉氏变换与相应的序列的Z变换之间的映射关系可用下式表示：2024/7/17第第54页页 模拟信号ha(t)的傅里叶变换Ha(j)和其采样信号的傅里叶变换之间的关系如下：将s=j代入上式，得上式表明将模拟信号ha(t)的拉氏变换在s平面上沿虚轴按照周期s=2/T延拓后，再按 映射关系映射到z平面上，就得到H(z)。可称为标准映射关系。2024/7/17第第55页页下面进一步分析这种映射关系。设得到:因此得到:那么 =0,r=1 0,r0,r1 另外，注意到z=esT是一个周期函数，可写成为任意整数2024/7/17第第56页页z=esT,s平面与z平面之间的映射关系频率混叠现象2024/7/17第第57页页 假设 没有频率混叠现象，即满足由 ，并将关系式s=j代入，=T，代入得到：为避免采样间隔较小而导致 H(z)过大，令：则：2024/7/17第第58页页 如果模拟滤波器的二阶基本节的形式为极点为 可以推导出相应的数字滤波器二阶基本节为：如果模拟滤波器二阶基本节的形式为 极点为 相应的数字滤波器二阶基本节为：一般Ha(s)的极点是复数且共轭成对的形式出现，将一对复数共轭极点放在一起，形成一个二阶基本节。2024/7/17第第59页页 例 已知模拟滤波器的传输函数Ha(s)为用脉冲响应不变法将其转换成数字滤波器的系统函数。解 首先将Ha(s)写成部分分式：极点为：那么H(z)的极点为：由上式并经整理得：2024/7/17第第60页页也可以首先将Ha(s)写成如下形式，如极点s1，2=1j1,则再由上分析得H(z)为：T=1sT=0.1s2024/7/17第第61页页脉冲响应不变法的特点：优点：优点：T保证了模拟角频率与数字角频率之间的线性关系保证了模拟角频率与数字角频率之间的线性关系缺点：缺点：在在 附近产生混叠失真附近产生混叠失真所以，一般只适用于带限滤波器，高通、带阻滤波器不能直接所以，一般只适用于带限滤波器，高通、带阻滤波器不能直接进行变换进行变换所以，脉冲响应不变法适用情况：所以，脉冲响应不变法适用情况：带限滤波器带限滤波器H(s)要能用部分分式展开成单阶极点要能用部分分式展开成单阶极点2024/7/17第第62页页双线性变换的含义双线性变换的含义n第一次映射第一次映射:将整个将整个S平面压缩到平面压缩到S1平面的一条横带内平面的一条横带内（-/T,+/T）n第二次映射第二次映射:将将S1平面映射到平面映射到Z平面平面双线性映射的变换式双线性映射的变换式n第一次映射：第一次映射：正切变换实现频率压缩n第二次映射第二次映射:双线性变换法的公式：双线性变换法的公式：第四节第四节 双线性变换法设计双线性变换法设计IIR数字低通滤波器数字低通滤波器2024/7/17第第63页页S平面S1平面Z平面2024/7/17第第64页页双线性变换法的频率变换关系双线性变换法幅度和相位特性非线性映射2024/7/17第第65页页系数关系表 2024/7/17第第66页页2024/7/17第第67页页 例 试分别用脉冲响应不变法和双线性不变法将下图所示的RC低通滤波器转换成数字滤波器。解 首先按照图写出该滤波器的传输函数Ha(s)为 利用脉冲响应不变法转换，数字滤波器的系统函数为 利用双线性变换法转换，数字滤波器的系统函数H2(z)为2024/7/17第第68页页H1(z)和H2(z)的网络结构分别如图所示数字滤波器H1(z)和H2(z)的幅频特性 2024/7/17第第69页页 利用模拟滤波器设计IIR数字低通滤波器的步骤1)确定数字低通滤波器的技术指标：通带截止频率p、通带衰减p、阻带截止频率s、阻带衰减s。2)将数字低通滤波器的技术指标转换成模拟低通滤波器的技术指标。如果采用双线性变换法，边界频率的转换关系为 3)按照模拟低通滤波器的技术指标设计模拟低通滤波器。4)将模拟滤波器Ha(s)，从s平面转换到z平面，得到数字低通滤波器系统函数H(z)。2024/7/17第第70页页 解 (1)用脉冲响应不变法设计数字低通滤波器。数字低通的技术指标为 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。指定采用巴特沃斯低通滤波器。试分别用脉冲响应不变法和双线性变换法设计滤波器。2024/7/17第第71页页 设计巴特沃斯低通滤波器，先计算阶数N及3dB截止频率c。取N=6。由p和p求得3dB截止频率c=0.7032rad/s。根据阶数N=6，查表得归一化传输函数为2024/7/17第第72页页用脉冲响应不变法将Ha(s)转换成H(z)。首先将Ha(s)进行部分分式，并按脉冲响应不变法得到：去归一化：将p=s/c代入Ha(p)中，得实际传输函数：2024/7/17第第73页页 (2)用双线性变换法设计数字低通滤波器。数字低通技术指标仍为 p=0.2rad,p=1dB;s=0.3rad,s=15dB 模拟低通的技术指标为2024/7/17第第74页页 设计巴特沃斯低通滤波器。阶数N计算如下：取N=6。由s和s求c=0.7662rad/s。查表得归一化传输函数Ha(p)，去归一化，将p=s/c代入Ha(p)，得实际的Ha(s)2024/7/17第第75页页 用双线性变换法将Ha(s)转换成数字滤波器H(z)：双线性变换法的幅度特性脉冲响应不变法的幅度特性2024/7/17第第76页页第五节第五节 数字高通、带通和带阻滤波器设计数字高通、带通和带阻滤波器设计具体设计步骤如下：1)确定所需类型数字滤波器的技术指标。2)将数字滤波器技术指标转换成相应的模拟滤波器技术指标，转换公式为3)将所需模拟滤波器指标转换成模拟低通滤波器技术指标4)设计模拟低通滤波器。5)将模拟低通转换成所需类型的模拟滤波器。6)采用双线性变换法，将所需类型的模拟滤波器转换成所需类型的数字滤波器。2024/7/17第第77页页 例 设计一个数字高通滤波器，要求通带截止频率p=0.8rad，通 带 衰 减 不 大 于 3dB，阻 带 截 止 频 率s=0.44rad,阻带衰减不小于15dB。希望采用巴特沃斯型滤波器。解 (1)数字高通的技术指标为 p=0.8rad,p=3dB;s=0.44rad,s=15dB (2)模拟高通的技术指标计算如下：令T=1，则有2024/7/17第第78页页(3)模拟低通滤波器的技术指标计算如下：将p和s对3dB截止频率c归一化，这里c=p,(4)设计归一化模拟低通滤波器G(p)。模拟低通滤波器的阶数N计算如下：2024/7/17第第79页页 查表得到归一化模拟低通传输函数G(p)为 去归一化，将p=s/c代入上式得到：(5)将模拟低通转换成模拟高通。将上式中G(s)的变量换成1/s，得到模拟高通Ha(s)：(6)用双线性变换法将模拟高通转换成数字高通H(z)：实际上(5)、(6)两步可合并成一步，即2024/7/17第第80页页 例 设计一个数字带通滤波器，通带范围为0.3rad到0.4rad，通带内最大衰减为3dB，0.2rad以下和0.5rad以上为阻带，阻带内最小衰减为18dB。采用巴特沃斯型模拟低通滤波器。解 (1)数字带通滤波器技术指标为 通带上截止频率 u=0.4rad 通带下截止频率 l=0.3rad 阻带上截止频率 s2=0.5rad 阻带下截止频率 s1=0.2rad 通带内最大衰减p=3dB，阻带内最小衰减s=18dB。2024/7/17第第81页页 (2)模拟带通滤波器技术指标如下：设T=1，则有(通带中心频率)(带宽)将以上边界频率对带宽B归一化，得到 u=3.348,l=2.348;s2=4.608,s1=1.498;0=2.8042024/7/17第第82页页 (3)模拟归一化低通滤波器技术指标：归一化阻带截止频率归一化通带截止频率p=1p=3dB,s=18dB (4)设计模拟低通滤波器：查表得归一化低通传输函数G(p)2024/7/17第第83页页 (5)将归一化模拟低通转换成模拟带通：(6)通过双线性变换法将Ha(s)转换成数字带通滤波器H(z)。下面将(5)、(6)两步合成一步计算：将上式代入(5)中的转换公式，得将上面的p等式代入G(p)中，得 2024/7/17第第84页页 例 设计一个数字带阻滤波器，通带下限频率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=1，则有2024/7/17第第85页页阻带中心频率平方为 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=3dB2024/7/17第第86页页 (4)设计模拟低通滤波器：(5)将G(p)转换成模拟阻带滤波器Ha(s)：(6)将Ha(s)通过双线性变换，得到数字阻带滤波器H(z)。2024/7/17第第87页页第六节第六节 IIR 数字滤波器的直接设计法数字滤波器的直接设计法 1.零极点累试法称为零极点累试法。在确定零极点位置时要注意：1)极点必须位于单位圆内，保证数字滤波器因果稳定；2)复数零极点必须共轭成对，保证系统函数有理式的系数是实的。a)零极点分布；b)幅度特性 2024/7/17第第88页页 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)个待定的系数，求它们的原则是使E最小。下面研究采用第一式网络结构，如何求出(4K+1)系数。2024/7/17第第89页页由上式可知E是(4K+1)个未知数的函数，用下式表示：上式表示4K个系数组成的系数向量。为推导公式方便，令 为选择A使E最小，令2024/7/17第第90页页 设k是 的第k个分量(ak或bk或ck或dk)，因为 ,式中H*i表示对Hi函数共轭。将上式具体写成对ak,bk,ck,dk的偏导，得到：式中，k=1,2,3,:,K；i=1,2,3,:,N。2024/7/17第第91页页同理求得 由于系统函数是一个有理函数，极、零点均以共轭成对的形式存在，对于极点z1,一定有下面关系：2024/7/17第第92页页 例 设计低通数字滤波器，其幅度特性如图所示。截止频率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,:,2024/7/17第第93页页 待求的参数是A，a1,b1,c1,d1。设初始值=(0000.25)T经过90次迭代，求得E=1.2611，系统函数零、极点位置为 零点0.67834430j0.73474418;极点0.75677793j1.3213916 为使滤波器因果稳定，将极点按其倒数搬入单位圆内，再进行62次优化迭代，求得结果为 零点0.82191163j0.56961501;极点0.89176390j0.19181084;Ag=0.11733978,E=0.56731 N=29，取k=1，系统函数为2024/7/17第第94页页 3.在时域直接设计IIR数字滤波器 设我们希望设计的IIR数字滤波器的单位脉冲响应为hd(n)，要求设计一个单位脉冲响应h(n)充分逼近hd(n)。下面我们介绍这种设计方法。设滤波器是因果性的，系统函数为式中a0=1，未知系数ai和bi共有N+M+1个，取h(n)的一段，0np-1，使其充分逼近hd(n)，用此原则求解M+N+1个系数。2024/7/17第第95页页令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的非线性函数2024/7/17第第96页页考虑到iM时，bi=0，一般表达式为：设x(n)为给定的输入信号，yd(n)是相应的希望的输出信号，x(n)和yd(n)长度分别为M和N，实际滤波器的输出为y(n)，下面按照y(n)和yd(n)的最小均方误差求解滤波器的最佳解，设均方误差用E表示：上式中x(n),0nM1；yd(n)，0nN-1 为选择h(n)使E最小，令2024/7/17第第97页页 由上式得到2024/7/17第第98页页 例 设计数字滤波器，要求在给定输入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.31 3h(2)+9h(3)=0.03 2024/7/17第第99页页 解联立方程，得 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)在五项以后幅度值很小。