第六FIR数字滤波器的设计.ppt

表中 的系数 与原系数 和 的关系如下： （6-4-4） （6-4-5） （6-4-6） 将（6-4-3）式代入（6-4-1）式，得到 （6-4-7） （6-4-8） （6-4-9） 因而最优化的问题变成了选择 的系数使误差函数 的最大绝对值最小。 （6-4-10） 2、利用切比雪夫等波纹最佳逼近法设计  FIR滤波器的原理 交错定理：设 为 的线性组合 在A上能唯一的最佳逼近连续函数 的充分必要条件是加权误差函数 在A内至少 有r+1个极值频率点，即在A内必须存在r+1个频率点 ，其中 且 及 （6-4-11） 下面以情况1为例，讨论线性相位滤波器幅度响应极值数目的约束问题。 表示成 多项式 将余弦函数 将上式代入对情况1的 ，得 式中r=(N+1)/2，ak是与a (n)和amn有关的常数 求极值频率，将 对ω微分一次得 上式右边余弦函数多项式的最高阶次为r-2=(N-3)/2，应有r-2个根，sinω在ω=0和ω=π处有两个根(零点)，所以在区间[0,π]上，Hg(ω)应有r个极值点。 对逐段恒定的幅度特性，Hg(ω)的极值点也是E(ω)的极值点。 (6-4-12) 根据等波纹最佳一致逼近定理，误差函数至少应有 r+1个极值，因此E(ω)本身至少还应再增加一个独立极 值。对理想低通特性的Hd(ω)，Hg(ω)的边界频率ωp 和ωc也是E(ω)的极值点(不是Hg(ω)的极值点) ，这样 E(ω)就有r+2个极值频率(这种情况称为超波纹逼近)。 例如，N=13, (N+1)/2=7 时幅度特性最多有7个极值点，而误差函数E(ω)有9个极值点，如下图所示。 3、 求解方法及Remez算法 设已知在A上r+1个极值率为 将这些频率代入式 式中 按最佳逼近定理得 (6-4-13) 令an=B(n)，将(6-4-13)式写成矩阵形式 (6-4-14) 上式中有r+1个极值频率需要求解，实际求解起来很困难。 Remez算法的计算步骤 Remez算法：靠一次次迭代求得一组交错点即极值频率点，从而求出系数，而且在每一次求解极值频率的迭代过程中能够避免直接求解(6-4-14)式。 (1) 首先给出r+1个极值频率的初始估计值，通常在A上等间隔地取r+1个初始频率点, 即 由下式计算δ值 （6-4-15） 式中 一般不可能是最优值 (2) 将δ值及r+1个极值频率点代入(6-4-13)式，可求得 的离散值 利用拉格朗日（Lagrange）插值公式, 可由 的离散 (6-4-16) 值得到连续的 将 代入(6-4-10)式可求得误差函数E(ω)，如果在子集A的 所有频率上都有 说明δ是最佳的极值，初始估计值 恰好就是交错点组，工作可结束。 如果在某些频率点处， 说明初始估计的极值点不对，然后在这些点附近搜索局部极值点， 于是得到一组新的交错点，从而完成了一次迭代。 (3) 利用和步骤(2)相同的方法，把在各频率处使 的局部极值点作为新的极值频率的估计值，从而又得到一组 新的极值估计值。 重复以上步骤，每次迭代后的δ都是递增的, 最后收敛到 也就最佳一致逼近 . 迭代结束后，最后按新的最佳极值频率求出 B(n) （即 a(n)， ），从而求得FIR滤波器的单位 最佳极值，此时 的r个系数 脉冲响应h(n)。 在滤波器的ωp和ωs 及δp和δs等指标给定后，在利用Remez 交换算法开始设计前，需要能比较准确地确定N值。 对低通滤波器，当 进行估计： 时，N值可由下述经验公式 式中 (6-4-17) (6-4-18) 。 (6-4-19) 也可用下式进行粗略的估计 (6-4-20) 其中 对于带通滤波器，N值的估计可按以下经验公式进行：令ωp和ωs分别为带通滤波器的通带和阻带边频（取带通滤波器的两个过渡带中过渡带最小的边频），则有 (6-4-21) 式中 人有了知识，就会具备各种分析能力， 明辨是非的能力。 所以我们要勤恳读书，广泛阅读， 古人说“书中自有黄金屋。 ”通过阅读科技书籍，我们能丰富知识， 培养逻辑思维能力； 通过阅读文学作品，我们能提高文学鉴赏水平， 培养文学情趣； 通过阅读报刊，我们能增长见识，扩大自己的知识面。 有许多书籍还能