基于GN的IIR滤波器设计序列最小化方法.pdfVIP

第3l卷第4期 杭州电子科技大学学报 、roI．31．No．4 0f 2011年惦月 J,mmal HallgzIlou Aug．201I Dian五Umv∞fity 基于GN的IIR滤波器设计序列最小化方法 罗家浩，赖晓平，侯秀竹，杨吉祥 (杭州电子科技大学信息与控制研究所，浙江杭州310018) 摘要：无限冲击响应数字滤波器不具有内禀稳定性，因此在其设计中要考虑稳定性问题。无限冲 击响应滤波器的一种充分必要稳定条件就是稳定三角形条件。为了充分利用三角形稳定条件的 充分必要性和线性性质，该文针对无限冲击响应滤波器的最大误差最小化设计问题，提出将序列 最小化技术与高斯牛顿策略相结合的方法，把原问题转化成一系列的二阶锥规划子问题进行求 解。仿真实例说明了此方法的有效性。 关键词：无限冲击响应；数字滤波器；序列最小化技术；高斯牛顿策略 中图分类号：TN911．72文献标识码：A 文章编号：1001—9146(2011)04—0020—04 O 引 言 数字滤波器的设计问题是数字信号处理中的基本问题，IIR滤波器设计中通常采用首一多项式表 示滤波器传递函数的分母，针对这一分母多项式现存多种充分的稳定性约束条件【l。]。其中，正实性条 条件只是保证滤波器稳定的充分条件而非必要条件，因此那些具有更好性能的稳定滤波器可能被排除 在由上述条件确定的稳定域之外。另外一种表示滤波器传递函数分母的方法是将其由一连串二阶因子 级联而成，对于每个二阶因子均要求其零点位于单位圆内，即使该因子的两个多项式系数位于稳定三角 形H．51内。然而，采用这种级联二阶因子对分母参数化的方法也存在一些缺点，当传递函数的分母不止 含有一个二阶因子时，频率响应函数关于二阶因子的多项式系数是非线性的。为了克服这个困难，文献 a ofMinimizationthat theSee· 7中引入一种序列最小化方法(solves Sequential problemsahematelyupdates ond—Order Conic (Second—Order 此方法绝对收敛且能得到比现有一些方法更好的设计结果。 1设计问题的数学描述 设计滤波器频率响应函数分子和分母的阶数分别为M和N，用Hz)表示IIR滤波器的传递函数 如下： 。 收稿日期：2011—06一10 基金项目：国家自然科学基金资助项目；浙江省自然科学基金资助项目(Y1090109) 作者简介：罗家浩(1986一)，男，河南禹州人，在读研究生，数字滤波器设计． 第4期 罗家浩等：基于GN的IIR滤波器设计序列最小化方法 2l 一=而犏 1B(z)=bo+blz一”·+bMzd Iu P=1，2，…，P 【AP(z)=1+仅plz一1+ap2z—z 式中，P=N／2，bm J，a2 ， J b，…，bM]T，a=[o【…or,12’0【21，of．，…，’…，M Lo【1l，12，0【2l’0[,22，…，ff．pIO[,r2 坟∥‘mJ刀期望明狈翠州胜，开疋义：设D(m)为期望的频率响应，并定义： E(∞)=日(扩)一D(∞) (2) 为实际频率响应与期望频率响应