FIR数字滤波器的设计等波纹最佳逼近法.docxVIP

FIR 数字滤波器的设计--等波纹最佳逼近法一、等波最佳逼近的原理简介 等波纹最佳逼近法是一种优化设计法，即最大误差最小化准则，它克服了窗函数设计法和频率采样法的缺点，使最大误差（即波纹的峰值）最小化，并在整个逼近频段上均匀分布。用等波纹最佳逼近法设计的FIR 数字滤波器的幅频响应在通带和阻带都是等波纹的，而且可以分别控制通带和阻带波纹幅度，这就是等波纹的含义。最佳逼近是指在滤波器长度给定的条件下， 使加权误差波纹幅度最小化。与窗函数设计法和频率采样法比较，由于这种设计法使滤波器的最大逼近误差均匀分布，所以设计的滤波器性能价格比最高。阶数相同时，这种设计法使滤波器的最大逼近误差最小，即通带最大衰减最小，阻带最小衰减最大；指标相同时，这种设计法使滤波器阶数最低。 d等波纹最佳逼近法的设计思想 。用H (?) 表示希望逼近的幅度特性函数，要求设计线性相位 FIR 数字滤波器时, Hd (? ) 必须满足线性相位约束条件。用 d H ?? ?表示实际设计的滤波器的幅度特性函数。定义加权误差函数? ???为 ? ???? W ????Hd ???? H ???? 式中，W ?? ?为幅度误差加权函数，用来控制不同频带（一般指通带和阻带） 的幅度逼近精度。等波纹最佳逼近法的设计在于找到滤波器的系数向量h?n?， 使得在通带和阻带内的最大绝对值幅度误差? ??? 为最小，这也就是最大误差 最小化问题。 二、等波纹逼近法设计滤波器的步骤和函数介绍 根据滤波器的设计指标的要求：边界频率，通带最大衰减，阻带最大衰等 估计滤波器阶数n，确定幅度误差加权函数W ?? ? 采用Parks-McClellan 算法，获得所设计滤波器的单位脉冲响应h?n? 实现 FIR 数字滤波器的等波纹最佳逼近法的 MATLAB 信号处理工具函数为firpm 和firpmord。 firpm 函数采用数值分析中的多重交换迭代算法求解等 波纹最佳逼近问题，求的满足等波纹最佳逼近准则的FIR 数字滤波器的单位脉冲响应h?n?。firpmord 根据逼近指标,计算采用Parks-McClellan 算法等波纹最佳逼近滤波器的最低阶数，误差加权向量w,归一化边界频率向量f。 3 对 firpm 和firpmord 的说明firpm 函数功能：采用Parks-McClellan 算法设计FIR 滤波器函数格式：hn=firpm(n,f,m,w) n 是滤波器的阶 hn 是数字滤波器的单位脉冲响应，其长度为n+1 f 是希望滤波器的边界频率向量，要求 f 是单调增向量，并且从 0 开始，以 1 结束，1 对于数字频率? m 是与 f 对应的希望滤波器的幅度向量，m 和f 的长度相等，m?k ?表示希望滤波器在频率点 f ?k ?上的幅频响应，m 和 f 给出了希望滤波器的幅度特性。 w 是误差加权向量，其长度为 f 的一半。w(i) 表示对 m 中第 i 个频率段幅度逼近精度的加权值。w 缺省时，函数默认 w 全为 1，即每个频率段的逼近误差加权值相同。 firpmord 函数 函数功能：根据逼近指标，计算采用Parks-McClellan 算法等波纹最佳逼近滤波器的最低阶数n，误差加权向量 w 和归一化边界频率f。其返回参数作为firpm 函数的调用参数。 函数格式：[n,f,m,w]=firpmord(f,m,rip,fs) f 可以是归一化边界数字频率向量，也可以是模拟边界频率向量，但必须以 0 开始，以 1 结束或 fs/2 结束，并且其中省略了 0 和fs/2 两个频率点。 fs 是时域采样频率，单位Hz。fs 缺省时,函数默认fs=2Hz.但这是f 的长度 （包括省略 0 和 fs/2 两个频率点）是 m 的两倍，即 m 中的每个元素表示 f 给定的一个逼近频段上的希望逼近的幅度值。 rip 表示f 和 m 描述的各逼近段允许的波纹振幅（幅频响应最大偏差）。4.注意事项： 省略fs 时，f 必须是归一化的数字频率 有时计算的阶数n 略小，使设计结果达不到指标要求，这时要取n+1 或n+2 三、程序 低通滤波器设计 % Lowpass filter: wp=, ws=, N=26 % peak passband ripple is , peak stopband ripple is clc;clear;clear all; wp_l=*pi; ws_l=*pi; N1=26; deltal_1=; deltal_2=; %初始化参 数 f_l=[0 wp_l ws_l pi]/pi; %设定归一化的firpm函数参数deltal=[deltal_1,deltal_2]; a_l=[1 1 0 0]; %幅值 [n,f,m,w_l]=firpmord(f_l(2:3),