数字信号处理-第6章I

1第6章无限长单位脉冲响应IIR滤波器设计2概述:许多信息处理过程，如信号的过滤，检测、预测等都要用到滤波器，数字滤波器是数字信号处理中使用得最广泛的一种线性系统，是数字信号处理的重要基础。实现方法主要有两种：数字信号处理硬件和计算机软件。3数字滤波器(DF)的定义输入和输出均是数字信号，通过一定运算关系(数值运算)，改变输入数字信号所含频率成份的相对比例或滤除某些频率成份的器件。数字滤波器的特点(相对模拟滤波器)：精度高、稳定、体积小、重量轻、不要求阻抗匹配。数字滤波器处理模拟信号通过A/DC和D/AC，使用数字滤波器对模拟信号的处理。4(1)一般分类经典滤波器：输入信号中的有用的频率成分和希望滤除的频率成分占用不同的频带，通过选频滤波器达到滤波的目的。现代滤波器：信号和干扰的频带相互重叠，要利用信号的统计分布规律，从干扰中最佳提取信号，如：维纳滤波器、卡尔曼滤波器和自适应滤波器等。(2)从滤波器的结构分类IIR/FIR(3)从滤波器的功能分类分为低通、高通、带通、带阻滤波器1、数字滤波器的分类5DF的传输函数是以2为周期，低通的中心频带处于2的整数倍处，高通的中心频带处于的奇数倍附近。0π67确定传输函数H(z)的过程称为数字滤波器设计.在大多数应用中，关键的问题是用一个可实现的传输函数去逼近给定的滤波器幅度响应指标，而滤波器的相位响应可以通过级联全通滤波器来校正。一种广泛应用的IIR滤波器设计方法是将一个模拟的原型传输函数转换为一个数字的传输函数，FIR滤波器的设计则是基于对指定幅度响应的直接逼近在设计数字传输函数H(z)之前，有两个关键的问题需要考虑：1.分析使用数字滤波器的整个系统的需求，确定合理的滤波器频率响应指标.2.确定所设计的滤波器是FIR还是IIR数字滤波器。8NiNiiiinybinxany01)()()(数字滤波器的数学描述：1)差分方程NMzdzcAzbzaZHNiiMiiNiiiMiii一般111110)1()1(1)(2)系统函数91110.7072pcs)(jeH通带过渡带阻带低通滤波器的技术要求10阻带截止频率通带截止频率::sp2s1p)(1)(10jjeHeH在阻带中，阻带频率范围：）－在通带中（，通带频率范围：11阻带内允许最小衰减通带内允许最大衰减::spdBeHeHdBeHeHspjjjj)()(lg20)()(lg200s0p12，则归一化为将1)(H0jedBeHdBeHjjp)(lg20)(lg20ssp－－cpc,322称＝，＝时，当幅度下降到dB为3dB通带截止频率。sc,,p通称为边界频率。13)(jeHc11+p1-psps低通数字滤波器的典型幅度指标14dBs=-20log10(s)dBppp11lg201511/A211cps)(jeH归一化的数字低通滤波器幅度响应指标162101log20pdBAs10log20dB17数字滤波器的设计步骤1）按照实际需要确定滤波器的性能要求。2）用一个因果稳定系统的H(z)或h(n)去逼近这个性能要求，即求h(n)的表达式确定系数、或零极点、，以使滤波器满足给定的性能要求3）用一个有限精度的运算去实现这个系统函数。包括选择运算结构：如级联型、并联型、卷积型、频率采样型以及快速卷积（FFT）型等;选择合适的字长和有效数字的处理方法等。iaibicid18IIR数字滤波器设计方法：1）先设计一个合适的模拟滤波器，然后变换成满足预定指标的数字滤波器。由于模拟的网络综合理论已经发展得很成熟，模拟滤波器有简单而严格的设计公式，设计起来方便、准确、可将这些理论推广到数字域，作为设计数字滤波器的工具。DF指标转换相应的AF指标)(sHaH（z）spsp,,,spsp,,,（）（）方法：脉冲响应不变法、双线性变化法注：仅适合IIR－DF，而不适合FIR－DF。192）最优化设计方法分两步：a)确定一种最优准则，如最小均方误差准则，即使设计出的实际频率响应的幅度特性（与所要求的理想频率响应的均方误差最小，此外还有其他多种误差最小准则，b)在此最佳准则下，求滤波的系数和通过不断地迭代运算，改变、，直到满足要求为止。|)(|jeH|)(|jdeHmin)()(21MijdjiieHeHiaibiaib零极点累试法频域逼近法、时域逼近法20以上两种设计方法中，第一种比较重要，因为数字滤波器在很多场合所要完成的任务与模拟滤波器相同，如作低通、高通、带通及带阻网络等，这时数字滤波也可看作是“模仿”模拟滤波器。在IIR滤波器设计中，采用这种设计方法目前最普遍。由于计算机技术的发展，最优化设计方法的使用也逐渐增多。21§6.2常用模拟低通滤波器特性为了方便学习数字滤波器，先讨论几种常用的模拟低通滤波器设计方法，高通、带通、带阻等模拟滤波器可利用变量变换方法，由低通滤波器变换得到。模拟滤波器的设计指标有sspp,,,221110.7072pcs)(jHa通带过渡带阻带图：低通滤波器的技术要求23阻带截止频率通带截止频率::sp，阻带频率范围：，通带频率范围：sp024阻带内允许最小衰减通带内允许最大衰减::spdBjHjHdBjHjHsaapaa)()0(lg20)()0(lg20sp25，则归一化为将1)0(HajdBjHdBjHsapa)(lg20)(lg20sp－－cpc,322称＝，＝时，当幅度下降到dB为3dB通带截止频率。26几种典型模拟滤波器（AF）的特点（对相同的阶数N）（1）巴特沃斯（butterworth）AF：幅频特性单调下降，但选择性最差；27（2）切比雪夫I（ChebyshevI）AF：幅频响应在通带内等波纹，阻带单调下降，选择性比切比雪夫II好。2829（3）切比雪夫II（ChebyshevII）AF：幅频响应在通带内单调下降，阻带等波纹；3031（4）椭圆（ellip）AF：选择性最好，通带阻带均为等波纹。32（5）贝塞尔（Bessel）AF：通带内有较好的线性相位特性。33模拟滤波器的设计相应的AF指标)(sHa（）spsp,,,？34模拟滤波器的设计就是根据一组设计规范设计模拟系统函数Ha(s)，使其逼近某个理想滤波器特性。因果系统中式中ha(t)为系统的冲激响应，是实函数。∴不难看出0)()(dtethjHtjaa0sincos)()(dttjtthjHaa)()(jHjHaa35定义振幅平方函数式中Ha(s)—模拟滤波器系统函数Ha(jΩ)—滤波器的频率响应|Ha(jΩ)|—滤波器的幅频响应又S=jΩ，Ω2=－S2∴A(Ω2)=A(-S2)|S=jΩ)1()()()()()()()()()(222jsaaaaaaasHsHjHjHAjHjHjHA36问题：由A(-S2)→Ha(S)对于给定的A(-S2)，先在S复平面上标出A(-S2)的极点和零点，由(1)式知，A(-S2)的极点和零点总是“成对出现”，且对称于S平面的实轴和虚轴，选用A(-S2)的对称极、零点的任一半作为Ha(s)的极、零点，则可得到Ha(s)。为了保证Ha(s)的稳定性，应选用A(-S2)在S左半平面的极点作为Ha(s)的极点。37零点的分布则无此限制，只和滤波器的相位特性有关，如果要求是最小相位延迟特性，则Ha(s)应取左半面零点，若无特殊要求，则可将对称零点的任一半（应为共轭对）取为Ha(s)的零点。38由此看出，由确定的方法如下：（1）由得到象限对称的s平面函数；（2）将因式分解，得到各零极点，将左半平面极点归于。轴上的零点或者极点都为偶次，应取一半（应为共轭对）作为的零点或极点。（3）按照与的低频或高频特性的对比就可以确定出增益常数。（4）由求出的零点，极点及增益常数，则可完全确定系统函数2)()()(jHsHsHajsaa2)(jHa)(sHa2)(jHa)(sHa)(sHa)(jHa)(sHaj)(sHa39例根据以下幅度平方函数确定系统函数2)(jHa)(sHa)36)(49()25(16)(22222jHa解：)36)(49()25(16)()()(22222sssjHsHsHajsaa)6)(7()25()(Hj5s67s20asssks＝，零点为，＝其极点为4)()(000kjHsHasa可以确定Step1:Step2:Step3:4042131004)6)(7()25(4)(H222asssssss411）巴特沃兹滤波器(Butterworth滤波器)(巴特沃兹逼近)特点：具有通带内最大平坦的振幅特性，且随f↗，幅频特性单调↘。其幅度平方函数：NcajjjHA22211)()(三种模拟低通滤波器的设计42图1巴特沃兹滤波器振幅平方函数43通带:使信号通过的频带阻带：抑制噪声通过的频带过渡带：通带到阻带间过渡的频率范围Ωc：3dB截止频率。过渡带为零，阻带|H(jΩ)|=0通带内幅度|H(jΩ)|=cons.，H(jΩ)的相位是线性的。理想滤波器44图1中，N增加，通带和阻带的近似性越好，过渡带越陡。通带内，分母Ω/Ωc1,(Ω/Ωc)2N1，A(Ω2)→1。过渡带和阻带，Ω/Ωc1，(Ω/Ωc)2N1,Ω增加，A(Ω2)快速减小。Ω=Ωc,，,幅度衰减，相当于3db衰减点。21)0()(2AAc2121)(2A45振幅平方函数的极点：令分母为零，得式中，k＝0，1，2，，，（2N－1）可见，Butterworth滤波器的振幅平方函数有2N个极点，它们均匀对称地分布在|S|=Ωc的圆周上。NcaajSSHSH2)(11)()()(2121221)()1(NkjccNkejS46例：N=3阶BF振幅平方函数的极点分布47考虑到系统的稳定性，知DF的系统函数是由S平面左半部分的极点（SP3，SP4，SP5）组成的，它们分别为：系统函数为：令，得归一化的三阶BF：如果要还原的话，则有3254323,,jcpcpjcpeSSeS))()(()(5433pppcaSSSSSSsH1c1221)(23SSSsHa1)/(2)/(2)/(1)(23cccassssH48由于各滤波器的幅度特性不同，为使设计统一，将所有频率归一化，采用对3dB截止频率归一化，归一化后的表示为c)(sHaccNkckcassssHj,)(1)(10＝其中，＝令c称为归一化频率；jp令，p称为归一化复变量，这样归一化的巴特沃斯的传输函数为4910)(1)(NkkapppH式中，为归一化极点，1,,2,1,0,)(21221NkeppNkjkkNNNappbpbpbbpH1122101)(则只要知道阶数N，由表6.2.1查表可得到)(pHa及各极点。(6.2.11)50问题：阶数N如何确定？阶数N的大小主要影响幅度特性下降速度，应由技术指标确定。sspp,,,51NcajjjH2211)(dBjHjHdBjHjHsaapaa)()0(lg20)()0(lg20sp52102102101101spNcsNcp解上面两式，得到1101101010spNsp53令110110,1010