有限长单位抽样响应数字滤波器的设计演示文件修改版.pptVIP

增加两点过渡带抽样点， 且增加抽样点数为N=65， 令H(17)=0.5886 H(18)=0.1065 过渡带宽： 阻带衰减：-60dB 窗函数设计法和频率抽样设计法的优点是简单直观、易于实现。但同时二者存在以下不足：所设计滤波器的通带和阻带的截止频率不易精确控制；通带和阻带内的波纹幅度不能分开控制；由于所设计的滤波器一般在阻带的截止频率附近刚好满足衰减要求，而在大于该截止频率的频段，衰减加大，因而存在很大的衰减裕量，即这两种设计方法存在较大的资源浪费，所设计的滤波器性价比不高。 7.5 等波纹最佳逼近设计法 针对窗函数设计法和频率抽样设计法存在的缺点，我们希望有一种方法，在滤波器长度 给定的情况下，设计出的滤波器是所能得到的所有滤波器中性能最佳的，即寻求一种最优化的设计方法，而等波纹最佳逼近设计法就是这样一种方法。 7.5.1 设计思想 设计滤波器的过程，其实就是用一个实际的频率响应函数H(ejω)（或对应的h(n) ）对所希望的频率响应函数Hd(ejω) （或对应的hd (n) ）的一个逼近过程。而最优化设计过程则是依据某一优化准则，改变一组变量值，进行迭代运算，寻求一种最佳逼近的过程。 函数逼近准则 （1）均方误差最小准则：窗函数设计法 （2）插值逼近准则：频率抽样设计法 这两种逼近有一个共同的缺点：在整个要求逼近的区间（0~π）内误差分布是不均匀的，特别是在幅度特性具有跃变点的过渡区附近，误差最大。 （3）最大误差最小化准则：等波纹最近逼近设计法，使得在要求逼近的整个范围内，逼近误差分布是均匀的，所以也称为最佳一致意义下的逼近。 在设计滤波器时，对通带和阻带的性能要求是不同的，为了统一使用最大误差最小化准则，则采用误差函数加权的办法，使得不同频段（如通带和阻带）的加权误差最大值是相等的。 设所希望的幅度频率响应为Hd(ω) ，实际设计得到的频率响应为H(ω) ，误差的加权函数为W(ω) ，则定义加权逼近误差函数为 （7.5.1） 使式(7.5.1)所定义的加权逼近误差函数 的最大绝对值在各个要求逼近的频带内达到最小，即 式中，||E(ω)||为得到的最小值，A为所要求逼近的通带和阻带。 7.5.2 交错点组定理 若H(ω)是r个余弦函数的线性组合，即 A是[0,π]内的一个闭区间（包括通带、阻带，但不包括过渡带），Hd(ω)是A上的一个连续函数，则Hd(ω)的唯一和最佳一致逼近的充要条件是：加权逼近误差函数E(ω)在A中至少存在r+2个交错频率点ω0ω1…ωr+1,使得 图7.5.1低通滤波器的最佳逼近波形 r=(N-1)/2=7(N=15) 交错点组： 0、ω1、ω2、ωp、ωst、ω3、ω4、ω5、ω6、π （1）交错点数最多可能是r+3个 （2）交错点比出现在ω=ωp、ωst处 （3）除ω=0或ω=π外，滤波器在通带和阻带内必是等波纹的 7.5.3 设计方法 若希望设计的滤波器为线性相位理想低通滤波器，其幅度响应函数为 现在的任务就是寻求一个H(ω) ，使其在通带和阻带内最佳一致逼近Hd(ω) ，当然为了保证线性相位，频率响应函数 要满足表7.1.1中的相应条件。图7.5.1中，δ1为通带内波纹峰值、 δ2为阻带内波纹峰值。因此在设计滤波器时，共有五个指标参数，即ωp 、 ωst 、δ1、δ2和单位抽样响应h(n)长度N。 对于低通滤波器，误差的加权函数为 定义为 用这种加权方式，逼近误差的最大绝对值在通带和阻带上均为δ=δ1 。 依据交错点组定理，若知道A上的r+2 （此时r=(N-1)/2）个交错点频率，则依据式(7.5.3)，得 （7.5.4） 将上式写成矩阵形式 瑞米兹（Remez）算法可以设计任何最优的线性相位FIR滤波器，是一种非常实用的算法，将该算法的步骤归纳如下： （1）在频率子集A上等间隔地取r+2个频率点ωi（i=0,1, …,r+1 ）作为交错点组（极值点频率）的初始猜测值。将这r+2个 代入下式求出 （7.5.5） 式中 此时的δ是相对于第一次猜测的交错点组产生的偏差，一般不是最优值。接下来依据重心形式的拉格朗日插值公式求得H(ω) ，即 （7.5.6） （7.5.7） 式中 把求得的H(ω)代入式(7.5.1) 中，求出加权逼近误差函数E(ω)。若对这组交错点组中的所有频率都满足 | E(ω )||δ|，则说明δ是纹波的极值，且初始猜测ωi（i=0,1, …,r+1）恰好是交错点组频率，计算即可结束。 （7.5.8） （2）对上次确定的交错点组 ωi（i=0,1, …,r+1）中的每一点，都检查在其附近是否存在使| E(ω )||δ|的频率点，如存在，则