§3．3 从模拟滤波器低通原型到各种数字滤波器的频率变换（原型变换）.pptVIP

§3．3 从模拟滤波器低通原型到各种数字滤波器的频率变换（原型变换） 也可以查表得到。由手册中查出巴特沃兹多项式的系数，之后以 代替归一化频率，即得 。 将 代入，就完成了模拟滤波器的设计，但为简化运算，减小误差积累，fc数值放到数字滤波变换后代入。 可见，H（Z）与采样周期T有关，T越小，H（Z）的相对增益越大，这是不希望的。为此，实际应用脉冲响应不变法时稍作一点修改，即求出H（Z）后，再乘以因子T,使H（Z）只与 有关，即只与fc和fs的相对值 有关，而与采样频率fs无直接关系。 例如， 与 的数字滤波器具有相同的传递函数，这一结论适合于所有的数字滤波器设计。 最后得： （1）频率变换关系 由于倒数关系不改变模拟滤波器的稳定性，因此，也不会影响双线变换后的稳定条件，而且 轴仍映射在单位圆上，只是方向颠倒了。 例 通带损耗 时， N 3时, 系统函数为： wc 2*1000*tan 2*pi*400/ 2*1000 ; wt 2*1000*tan 2*pi*317/ 2*1000 ; [N,wn] cheb1ord wc,wt,0.5,19,s ; [B,A] cheby1 N,0.5,wn,high,s ; [num,den] bilinear B,A,1000 ; [h,w] freqz num,den ; f w/pi*500; plot f,20*log10 abs h ; axis [0,500,-80,10] ; grid; xlabel ylabel 幅度/dB 解:首先确定上下边界频率 代入上式,求中心频率 再求模拟低通滤波器的边界频率 反归一化的系统函数: 则带阻数字滤波器为 P.88例7 w1 95/500; w2 105/500; [B,A] butter 1,[w1, w2],stop ; [h,w] freqz B,A ; f w/pi*500; plot f,20*log10 abs h ; axis [50,150,-30,10] ; grid; xlabel 频率/Hz ylabel 幅度/dB 3、稳定性证明 同时，这一变换也满足稳定性要求，设 由于上式完全是实数，所以是映射在S平面 轴上。 其中分子永远非负的 ， 因此 的正负决定于分母 由此证明了，S左半平面映射在单位圆内，而右半平面映射在单位圆外，这种变换关系是稳定的变换关系，可用它来完成带通的变换，如图1。 4、带通滤波器设计 设计带通时，一般只给出上、下边带的截止频率 作为设计要求。 为了应用以上变换，首先要将上下边带参数 换算成中心频率 及模拟低通截止频率 。 为此将 代入变换关系式： 由于 在模拟低通中是一对镜象频率， 代入上面两等式，求出 又 同时也就是模拟低通的截止频率 ， 有了这两个参数就可完成全部计算。 5、例子 采样 fs 400kHz，设计一巴特沃兹带通滤波器，其3dB边界频率分别为f2 90kHz,f1 110kHz,在阻带f3 120kHz处最小衰减大于10dB。 解：确定数字频域的上下边带的角频率 求中心频率： 求模拟低通的通带截止频率 与阻带边界频率 ： 从 频率增加了约1.05倍，衰减增加了（10-3）dB，故选用二阶巴特沃兹滤波器可满足指标（查表） 归一化的系统函数： 代入 ， 代入变换公式 P88 例6 带通滤波器设计 3dB 10dB f2 f1 f3 H ejω f w1 2*400*tan 2*pi*90/ 2*400 ; w2 2*400*tan 2*pi*110/ 2*400 ; wr 2*400*tan 2*pi*120/ 2*400 ; [N,wn] buttord [w1 w2],[0 wr],3,10,s ; [B,A] butter N,wn,s ; [num,den] bilinear B,A,400 ; [h,w] freqz num,den ; f w/pi*200; plot f,20*log10 abs h ; axis [40,160,-30,10] ; grid; xlabel 频率/kHz ylabel 幅度/dB 巴特沃兹带通滤波器 频率/kHz 幅度/ dB 四．带阻变换 1、变换原理 把带通的频率关系倒置就得到带阻变换。 给定 2、举例 例 一数字滤波器采样频率fs 1kHz，要求滤除100Hz的干扰，其３dB的边界频率为95Hz和105Hz，原型归一化低通滤波器为 * * 对于模拟滤波器,已经形成了许多成熟的设计方案，如巴特沃兹滤波器,切比雪夫滤波器,考尔滤波器,每种滤波