中学课件数字信号处理5.pptVIP

式中，Ω0 为带阻滤波器阻带的几何对称中心角频率; B为带阻滤波器的阻带宽度，它与低通原型中的截止频率Ωc成反比。  设数字带阻的中心频率为ω0，数字带阻滤波器的上、下边带的截止频率分别为ω2和ω1， 则将式（5-77）代入式（5-78）、 式（5-79），可得: （5-80） （5-81） 考虑到模拟带阻到数字带阻是阻带中心频率相对应的映射关系，则有 （5-82） 将式（5-80）、 式（5-81）和式（5-82）代入式（5-75）及式（5-76），并应用一些标准三角恒等式运算后可得: 所以，在设计时，要给定中心频率和带宽或者是中心频率和边带频率， 利用式（5-83）和式（5-84）来确定D1和E1两常数， 然后利用式（5-74）的变换，把模拟低通系统函数一步变成数字带阻系统函数 （5-85） 式中，Ha(s)为模拟低通原型传递函数。  可以看出，数字带阻滤波器的极点数（或阶数）将是模拟低通滤波器极点数的两倍。 下面来讨论模拟低通滤波器与数字带阻滤波器频率之间的关系。令s=jΩ，z=ejω代入式（5-74），经推导后可得 (5-86) 其变换关系曲线如图5-24所示，其映射关系为: Ω=0→ω=0, ω=π Ω=±∞→ω=ω0 也就是说，低通滤波器的通带（Ω=0 附近）映射到带阻滤波器的阻带范围之外（ω=0, π），低通滤波器的阻带（Ω=±∞）映射到带阻滤波器的阻带上（ω=ω0 附近）。 表 5-4 利用双线性变换法从截止频率为Ωc的低通原型模拟 滤波器到实际数字滤波器的变换 式中，Ωc=2πfc。上式也可由巴特沃思滤波器的幅度平方函数求得。  为了进行脉冲响应不变法变换，将上式进行因式分解并表示成如下的部分分式形式： 将此部分分式系数代入(5-40)式就得到 式中，ωc=ΩcT=2πfcT=0.5π是数字滤波器数字频域的截止频率。 将上式两项共轭复根合并，得 从这个结果我们看到，H(z)只与数字频域参数ωc有关，也即只与临界频率fc与采样频率fs的相对值有关，而与它们的绝对大小无关。 例如fs=4kHz，fc=1 kHz与fs=40 kHz，fc=10kHz的数字滤波器将具有同一个系统函数。这个结论适合于所有的数字滤波器设计。 将ωc=ΩcT=2πfcT=0.5π代入上式，得 这个形式正好适合用一个一阶节及一个二阶节并联起来实现。 脉冲响应不变法由于需要通过部分分式来实现变换，因而对采用并联型的运算结构来说是比较方便的。  图5-18给出了脉冲响应不变法得到的三阶巴特沃思数字低通滤波器的频响幅度特性，同时给出例5-5双线性变换法设计的结果。由图可看出，脉冲响应不变法存在微小的混淆现象，因而选择性将受到一定损失，并且没有传输零点。 图 5-18 三阶巴特沃思数字低通滤波器的频响 5.6.2 模拟低通滤波器变换成数字高通滤波器 由表5-3可知，由低通模拟原型到模拟高通的变换关系为 （5-52） 式中，Ωc为模拟低通滤波器的截止频率，Ωc′为实际高通滤波器的截止频率。  根据双线性变换原理，模拟高通与数字高通之间的S平面与Z平面的关系仍为 （5-53） 把变换式（5-52）和变换式（5-53）结合起来，可得到直接从模拟低通原型变换成数字高通滤波器的表达式，也就是直接联系s与z之间的变换公式 （5-54） 式中， 。由此得到数字高通系统函数为 式中，Ha(s)为模拟低通滤波器传递函数。 可以看出，数字高通滤波器和模拟低通滤波器的极点数目（或阶次）是相同的。  根据双线性变换，模拟高通频率与数字高通频率之间的关系仍为 则 又因 故 （5-55） 下面讨论模拟低通滤波器与数字高通滤波器频率之间的关系。令s=jΩ，z=ejω， 代入式（5-54），可得 或 （5-56） 其变换关系曲线如图5-19所示。由图可看出，Ω=0 映射到ω=π, 即z=-1上; Ω=∞映射到ω=0, 即z=1上。通过这样的频率变换后就可以直接将模拟低通变换为数字高通，如图5-20所示。 图5-19 从模拟低通变换到数字高通时频率间关系的曲线 还应当明确一点，所谓高通数字滤波器，并不是ω高到∞都通过， 由于数字域存在折叠频率ω=π，对于实数响应的数字滤波器， ω由π到2π的部分只是ω由π到0的镜像部分。因此，有效数字域仅只是从ω=0 到ω=π，高通也仅指这一端的高端，即到ω=π为止的部分。 图 5-20 模拟低通变换到数