滤波器设计方法.pptVIP

同样，取左半平面中的三对极点，得连续Butterworth低通滤波器 的系统函数： 同双线性变换公式代入，得离散连续Butterworth低通滤波器的系统函数： 第三十一页，共七十三页，2022年，8月28日 六阶Butterworth低通滤波器频率响应： ω=0.2π ? -0.56dB (0.937) ω=0.3π ? -15dB (0.1778) 阻带比连续滤波器下降快 ω = π ? Ω = ∞ 频率轴被压缩 第三十二页，共七十三页，2022年，8月28日 连续Butterworth滤波器幅度平方函数： 可用双线性变换的频率关系式： 直接代入得到离散Butterworth滤波器幅度平方函数： 式中 实际设计并不首先用，因为上式无法获得在z平面的极点，从而得到单位圆内的极点 ? 滤波器的系统函数H(z) 第三十三页，共七十三页，2022年，8月28日 双线性变换法的Butterworth，Chebyshev, elliptic滤波器逼近例子 所设计低通离散滤波器的技术指标： 即δ1=0.01， δ2=0.001， ωp=0.4π， ωs=0.6 例7.4 Butterworth逼近 省去具体计算过程，可得N=14，频率响应如图： 对数幅度图 第三十四页，共七十三页，2022年，8月28日 通带的幅度图： 群延迟图： 第三十五页，共七十三页，2022年，8月28日 例7.5 Chebyshev逼近 Chebyshev的类型： I型：频率响应 ---- 通带呈纹波特性，阻带单调 II型：频率响应 ----- 阻带呈纹波特性，通带单调 两种类型逼近阶数N均为八阶（比Butterworth要低） I型逼近的频率响应图： 第三十六页，共七十三页，2022年，8月28日 通带中幅度细节图： 群延迟图： 第三十七页，共七十三页，2022年，8月28日 II型逼近的频率响应图 通带中的细节图： 第三十八页，共七十三页，2022年，8月28日 群延迟图： 若通带和阻带均容许纹波 滤波器的阶数可以减少 两种类型逼近的零极点分布（I型零点对应于连续滤波器s = ∞） I型 II型 第三十九页，共七十三页，2022年，8月28日 例7.6 elliptic逼近 满足指标的有理函数最低阶逼近，N=6 通带中的细节图： 第四十页，共七十三页，2022年，8月28日 群延迟图： 零极点分布： 第四十一页，共七十三页，2022年，8月28日 讨论： （1）各种滤波器逼近的双线性变换法是IIR滤波器设计标准方法 （2）H(z)的极点在单位圆内（稳定性），零点在单位圆上 （3）非线性频率响应特性，非线性群延迟 （4）通带边缘处或过渡带出现最大群延迟（其它频带中群延迟变化较小） （5）II型Chebyshev具有最小的群延迟以及相应群延迟近似为常数的最宽通带区域 （6）elliptic逼近具有最低阶次的系统函数，因此计算量也最少。 第四十二页，共七十三页，2022年，8月28日 7.2 用窗函数法设计FIR滤波器 IIR滤波器设计 ----- 连续滤波器变换 （1）连续滤波器设计方法成熟 （2）IIR滤波器直接设计困难（在离散域） FIR滤波器设计 ------ 可以直接在离散域进行（频率域逼近） 线性相位条件 ? 避免了直接设计的复杂性 最简单的方法：窗函数法（window method） 基本思想：对理想滤波器的脉冲响应hd[n]（无限长，非因果序列） ? 截断，得到具有对称性的有限长因果序列h[n] ? FIR滤波器（线性相位、因果）------ 逼近 问题：时域截断 ?? 频域逼近（关系、影响） ？ 即， hd[n] ? h[n] Hd(ejω) ? H(ejω) 第四十三页，共七十三页，2022年，8月28日 理想滤波器的频率响应： 脉冲响应： ---- 非因果、无限长序列 ? 频率响应的不连续性 实际因果FIR滤波器的脉冲响应（因果、有限长）： 相当于： 其中窗函数： 矩形窗 第四十四页，共七十三页，2022年，8月28日 时域处理简单，讨论频域关系： 由傅立叶变换性质：时域相乘 ?? 频域卷积 得： 如图所示 （周期性） 若w[n]=1, 对于所有n, 即没有截断 则 W(ejω) 是周期为2π的δ函数（周期脉冲串） 并有 H(ejω) = Hd(ejω) 可见，W(ejω)（频域）能量集中（波形窄）----- 逼近性越好 但， W(ejω) 越窄 ? w[n]（时域）越长 ----- 矛盾性