数字信号处理教案第16次课.docVIP

4.4从模拟滤波器低通原型到各种数字滤波器的频率变换 ?? 对于模拟滤波器，已经形成了许多成熟的设计方案，如巴特沃兹滤波器，切比雪夫滤波器等，每种滤波器都有自己的一套准确的计算公式，同时，也已制备了大量归一化的设计表格和曲线，为滤波器的设计和计算提供了许多方便，因此在模拟滤波器的设计中，只要掌握原型变换，就可以通过归一化低通原型的参数，去设计各种实际的低通、高通、带通或带阻滤波器。这一套成熟、有效的设计方法，也可通过前面所讨论的各种变换应用于数字滤波器的设计，具体过程如下： ? ?? 也可把前两步合并成一步，直接从模拟低通归一化原型通过一定的频率变换关系，完成各类数字滤波器的设计，设计过程如下： ????????????????????? 频率变换 模拟原型? ? 数字低通、高通、带通、带阻? ??下面举例讨论应用模拟滤波器低通原型，设计各种数字滤波器的基本方法，着重讨论双线性变换法。 ????双线性变换法由模拟原型设计数字滤波器的四个步骤: ??? 1）确定数字滤波器的性能要求，确定各临界频率{ωk}。 2）由变换关系将{ωk}映射到模拟域，得出模拟滤波器的临界频率值{Ωk}。 ??? 3）根据{Ωk}设计模拟滤波器的Ha(s) 。 ??? 4) 把 Ha(s)变换成数字滤波器传递函数H(z) 。 ? 低通变换 ??例1.设采样周期 T=250μs(fs=4khz) ,设计一个三阶巴特沃兹LP滤波器,其3dB截止频率fc=1khz。分别用脉冲响应不变法和双线性变换法求解。 解:a.脉冲响应不变法 由于脉冲响不变法的频率关系是线性的，所以可直接按Ωc =2πfc设计Ha(s)。根据上节的讨论，以截止频率Ωc 归一化的三阶巴特沃兹 滤波器的传递函数为： 以 s/Ω 代替其归一化频率，得： ???????????? ? 也可以查表得到巴特沃兹多项式的系数，之后以 s/Ω代替归一化频率，即得Ha(s)。 将 Ωc=2πfc 代入，就完成了模拟滤波器的设计，但为简化运算，减小误差积累，fc数值放到数字滤波变换后代入。 ? ?? 为进行脉冲响应不变法变换，计算Ha(S)分母多项式的根，将上式写成部分分式结构： 对照前面学过的脉冲响应不变法中的部分分式形式有 将上式系数代入数字滤波器的传递函数： , Si ：极点 并将 Ωc=ωc/T代入，计算得：?? 合并上式后两项，并将代入，计算得： ?????????? ??? 可见，H（Z）与采样周期T有关，T越小，H（Z）的相对增益越大，这是不希望的。为此，实际应用脉冲响应不变法时稍作一点修改，即求出H（Z）后，再乘以因子T,使H（Z）只与ωc 有关，即只与fc和fs的相对值 fc/fs 有关，而与采样频率fs无直接关系。例如，fs=4KHZ,fc=1KHZ 与 fs=4KHZ,fc=10KHZ 的数字滤波器具有相同的传递函数，这一结论适合于所有的数字滤波器设计。 最后得： ??????????? ? b. 双线性变换法 ?(一)首先确定数字域临界频率 ?(二)根据频率的非线性关系，确定预畸的模拟滤波器临界频率 ?(三) 以S/Ωc代入归一化的三阶巴特沃模拟器传递函数 　???? ??并将Ωc=2/T代入上式。 (四)将双线性变换关系代入，求H(Z)。 ??? ???????? ????? ??????????????????????????????? ??? 图3.11为两种设计方法所得到的频响，对于双线性变换法，由于频率的非线性变换，使截止区的衰减越来越快，最后在折 叠频率处 (Z=-1,ω=π)形成一个三阶传输零点，这个三阶零点正是模拟滤波器在处的三阶传输零点通过映射形成的。因此，双线性变换法使过渡带变窄，对频率的选择性改善，而脉冲响应不变法存在混淆，且没有传输零点。 高通变换???????????????????????????????????? HYPERLINK /zq/signal/new/course/three/3_23.htm \l a变换方法的选用 设计高通、带通、带阻等数字滤波器时，有两种方法： ① 先设计一个相应的高通、带通或带阻模拟滤波器，然后通过脉冲响应不变法或双线性变换法转换为数字滤波器。 ??? 模拟原型模拟高通、带通、带阻数字高通、带通、带阻 设计方法同上面讨论的低通滤波器的设计。 即确定 {ωk} 转换为相应的{Ωk} ???????????????????????????????????? 高通、带通、带阻 模拟滤波器的