无限长单位冲激晌应 (IIR) 数字滤波器的设计方法.pptVIP

[例 6-7] 给定模拟滤波器指标与例 6-6 相同 , 即 求用巴特沃思滤波器实现时，所需的阶次、滤波器的截止频率 Ωc 以及滤波器的系统函数Ha(s)的表达式。 解 滤波器的 1dB 截止频率为 ,15dB 的截止频率为 ，而Ωc为待求的巴特沃思滤波器 3dB 的截止频率。 (1) 求 N。由题意可写出 由此得出 将它们代入巴特沃思滤波器幅度平方函数的表达式 (6-73) 式中 , 得 将此两式对 N 求解 , 可得 取 N=6。 2) 求 3dB 处的截止频率Ωc 。把 N=6 代入到阻带指标要求的公式中 , 阻带指标满足要求 , 但是由于 N=65.885, 因而通带指标一定超过性能要求。将 N=6 代人阻带条件 , 可得 则有 （3）由于 N 和Ωc都已求出, 可以用查表法来求解。查巴特沃思滤波器的表6-4（分母为 N 阶多项式）或表6-6（分母为因是分解形式），可得到归一化 Ωc=1 的系统函数为 利用(6-110)式 可得到截止频率为 时的系统函数 代入Ωc 的值 , 得 与例 6-6 比较 , 在同样性能指标下 , 巴特沃思滤波器所用的阶次 (N=6) 比切贝雪夫滤波器的阶次 (N=4) 要多。 * dengsong 设为阻带截止频率，即当 时，将上面 的 的表达式代入（6-80）式，可 得 再将（6-87）式代入，可得 由此得出 由于 ，所以由（6-81）式的第二式有 再 由此，并考虑（6-86）式，可得 如果要求阻带边界频率上衰减越大（即 A越大），也就是过渡带内幅度特性越陡，则所需的阶数N越高。 或者对 求解，可得 这里 是切贝雪夫滤波器的通带宽度，但不是3dB带宽， 可以求出3bB带宽为（A= ） 注意，只有当时 ，才能采用（6-90）式来求解 （因为满足 ） 给定后，就可求得滤波器 当然,已知 后，也可先求 的极点分布，即可求得 我们由幅度平方函数来求解，将 代入 ，再考虑 （6-80）式可得 令其分母多项式为零，即 即 令 Ch(.)和sh(.)分别为双曲余弦函数与双曲正弦函数。 由（6-92）式可得 为了导出 与 的关系，把（6-92）式代入（6-91）式， 且考虑到 的定义（见（6-81）式） 由此得出 由（6-94a）式,且考虑 可解出 则有 把a代入(6-94b)式,可得 因而 把(6-95)式 的及(6-96)式的 代入(6-93)式可得 的极点为 即 其中 将(6-98a)式.(6-98b)式分别取平方再化简,可得 在S平面的极点分布满足的关系式 这是一个椭圆方程,由于双曲余弦CH（.）大于双曲正弦CH(.) 故长轴为 (在虚轴上),短轴为 (在实轴上). 取 在s左半平面的极点，即(6-97)式的左半平面的 就是 极点。 实际上,用(6-99)式求解a,b并不方便,还可以再加简化,考虑到 其中 查数学手册知,若 则有 将此式取反变换有 由此得出 将 代入(6-102a)式、(6-102b)式,可得 这样,由(6-103)式可求得 ，由 可求得a,b,从而利用(6-101) 式求得 的极点. 求出幅度平方函数的极点后, 的极点就是s不面左半平面的 诸极点 ,从而得到切贝雪夫滤波器的系数函数为 常数K可由 和 的低频或高频特性对比求得. 也可用下面的方法直接求出(6-105a)式的系数K: 将(6-80)式开平方,并代入 再考虑(6-105a)式,则有 可以看出,第一个等号后的公式,其分母多项式的首项 的系数不为1,这是因为 的首项 系数为 , 因而其 项的系数为 ，整个分母多项式 项的系数 则为 ， 而第二个等号后的分式,其分母的首项 的