FIR滤波器约束Minimax设计算法.pdfVIP

　 第 24卷 第 2 期 系统工程与电子技术 Systems Engineering and Electronics 　 Vol.24 , No.2 2002 　 　 收稿日期:2000-10-24　　修订日期:2001-01-31 基金项目:国家自然科学基金;山东省自然科学基金资助课题 作者简介:赖晓平(1965-),男 ,副教授 ,博士 ,主要研究方向为控制理论与应用 ,数字滤波器设计 ,神经网络及电力系统负载预测。 　文章编号:1001-506X(2002)02-0084-05 FIR滤波器约束 Minimax设计算法 赖晓平 (山东大学威海分校控制工程系 ,山东 威海 264209) 摘　要:考虑了一类具有频域等式约束的 FIR滤波器Minimax 设计问题 ,得到了最优滤波器应满足的一个充分 必要条件。根据该条件 ,并应用参数估计中的正交投影算法以及数值逼近理论中的 Remez 多重交换算法 , 提出了 FIR滤波器约束Minimax 设计的OPRemez 算法。设计实例表明设计算法是十分有效的 , 并具有很好的数值特性 , 可 设计阶数很大的最优约束 FIR滤波器。 关键词:FIR滤波器;约束Minimax 设计;OPRemez算法 中图分类号:TN713.7　　文献标识码:A Constrained Minimax Design Algorithm for FIR Filters LAI Xiao-ping (Department of Control Engineering , Shandong University at Weihai , Weihai 264209 , China) Abstract:This paper studies the Minimax design of a class of FIR filters with frequency domain equation constraints.A suffi- cient and necessary condition for the optimal filter is obtained.According to this condition , anOPRemez algorithm for the design of optimal constrained filters is proposed , which is a combination of the orthogonalized projection(OP)algorithm and the Remezmultiple exchange algorithm.Simulation results demonstrate the effectiveness of the proposed designmethod. Keywords:FIR filters;constrainedMinimax design;OPRemez algorithm 1　设计问题描述 FIR滤波器的显著特点之一是容易实现线性相位 , 因而 在要求滤波后波形不产生畸变的许多重要场合得到了广泛 的应用。线性相位FIR滤波器的设计问题就是选择滤波器 的系数使其幅频特性在某一准则下逼近期望特性。一个阶 数 N 为奇数 、脉冲响应对称的 FIR滤波器 , 其幅频特性可表 示为 H(ω,α)= ∑ r-1 n=0 αncosnω (1) 式中　 r=(N-1)/ 2+1;α=(α0 , α1 , … , αr-1)———一组与 滤波器脉冲响应相联系的系数。 FIR滤波器的 Minimax 设计 就是选择α使H(ω,α)与期望特性在某频带上的最大加权 误差最小。 在许多应用场合 ,为了保证某些特定频率下滤波器有特 定的幅值 , 滤波器的设计除了要使最大加权误差最小外 , 还 要加上一些约束条件 , 使设计问题变成如下的约束 Minimax 逼近问题 ,即 min α f(α)=max ω∈A |W(ω)[ H(ω, α)-D(ω)] | (2) s.t.　H( ωj , α)= H j　(j =1 , 2 , … , M , M r) (3) 式中　A———区间[ 0 , π] 的一个闭子集; ωj ∈[ 0 , π] (1≤ j ≤ M r)———约束频率;D(ω)———定义在 A 上的期望 幅频特性;W(ω) 0 ———定义在 A 上的加权函数。E(ω, α)≡ W(ω)[ H(ω, α)(D(ω)] ———加权误差函数。约束频 率构成的集合记为 Ψ={ ω1 , ω2 , … , ωM}。 例如 , 在图像处理等应用中 , 低通滤波器的设计要求信 号的直流分量经过滤波器后保持不变。另