FIR滤波器切比雪夫逼近设计的直接Remez算法研究.pdf

142 信 号 处 理 FIR滤波器切比雪夫逼近设计的直接Remez算法 ‘ 赖晓平 (山东大学威海分校控制工程系) !摘要】针对FIR数字涟波器的切比雪夫逼近设计问题，本文提出了一套递推计算公式，使用这套公式可 直接由一组极值频率确定滤波器的单位采样响应.从而使FIR识字滤波器设计的Remez算法变得直接、简单 和明了。分析和计算实例都表明，本文的算法和标准的Remez设计算法在计算性能上相差无几。 关键询:FIR数字滤波器 切比雪夫通近 Remez算法 一、引言 FIR数字滤波器的设计问题一直是语音及图象处理领域人们所关心的问题之一，人们也提出了各种各样的 设计算法，如切比雪夫逝近设计的Remez多重交换算法n.xl准等披纹设计的加权最小二乘算法131,随机抽样 递推最小_乘算法141及进化策略算法151等。在这些算法中，L.RRabiner等人提出的Remez多重交换算法以最 小的最人误差及等波纹特性成为标准的设计算法写进了许多教科书Is71中。但是，标准的Remez算法采用了插 值、傅}i.叶变换及多重交换等多种技术，使这一算法显得非常复杂。本文提出一种由一组极值频率直接确定 滤波器单位采样响应的算法，使Romez算法变得简单明了。 二、设计问题及标准Remez算法 FIR数字滤WN脉冲传递函数H(z)一艺hn():一’的频率特性为H(ei)二艺h(n)-lwP(eJ)exv()，其中 P(el)为共辐频特性，Oe)为其相频特性.简单考虑奇对称 (即N为奇数且h(n)=h(N-1-n))的情况. 这时滤波器的相频特性例e)=一夸。是线性特性(这一特性使得任一信号经过滤波器后不会产生畸变，这 在许多应用场合都是非常必要的.)a幅频特性为 P(e.)一艺acosnax， (1) 其中;二午十】，as一h(r一1)，a二2h(r-1一n)(n二1,2,...,;一1). FIR数字池波器的设计问题就是寻找h《(n)),使涟波器的频率特性在某一准则下，通近期望特性.在奇对 称情况下就是要寻找a《,卜使P(e)逼近期望幅频特性D(eJ)。因此，FIR数字滤波器的设计问题本质上就 是一个函数返近问题。切比雪夫福祈设计1【,,6,71着眼于在有盒义的频带A上，使加权误差函数 E(el)=W(e)ID(el)一P(e)1 (W(e灼为误差加权函数)较均匀一致，并且通过合理选择(a.),能使IE(e)I的最大值达到最小.若用s 来表示这个最小值，则切比雪夫退近设计问题可表示为: 5一黔瞥IE(e})I 求解切比雪夫逼近设计问题的依据是下面的交错点组定理/13.1,31 文.点组定理:设A为闭区间。‘。二‘上的一个闭子集，P(e)为如((I)式所示的二I阶三角多项式，则 P(el)成为A上的连续函数D(e})的最佳一致逼近的充要条件是.P(e)在A上至少存在(r+1)个频率点w6 口-W,,rq，使得 E(ell)=-E(el) (k=0,1，一，r-1), 并且 IE(el})1=黔IE(el)I。 这(}0个频率点即是 “交错点组.，’ 根据交错点组定理，若交错点(极值频率)组(wk,k二0,1...,弓是已知的，且设最大误差为8，则可以 写出下面的方程组: 中受国家自然科学基金资助，基金项目批准号:65;74002 信 号 处 理 143 W(ew)(D(el.)一P(e))=(-1)8 (k=0,1.、小 注意到P(e)的表达式(1)，上面的方程组可以写成矩阵形式 甘