教案iis数字滤波器.pptVIP

如果模拟滤波器二阶基本节的形式为  (6.3.15)  则对应的数字滤波器二阶基本节的形式为   利用以上这些变换关系，可以简化设计，使实现结构中无复数乘法器。 极点为 (6.3.16) * 脉冲响应不变法的优点是频率变换关系是线性的，即ω=ΩT，如果不存在频谱混叠现象，用这种方法设计的数字滤波器会很好地重现原模拟滤波器的频响特性。另外一个优点是数字滤波器的单位脉冲响应完全模仿模拟滤波器的单位冲激响应波形，时域特性逼近好。但是，有限阶的模拟滤波器不可能是理想带限的，所以，脉冲响应不变法的最大缺点是会产生不同程度的频率混叠失真，其适合用于低通、带通滤波器的设计，不适合用于高通、带阻滤波器的设计。 * 用脉冲响应不变法设计数字滤波器方法是： 如果已知模拟滤波器的系统函数Ha(s) 1. 将Ha(s)部分分式展开， 2.再按公式得到数字滤波器系统函数H(z)。（T要选合适） 如果已知数字滤波器的设计指标：ωp, αp , ωs, αs 1、将数字滤波器性能指标变换为中间模拟滤波器的性能指标。采用: 2、设计出符合要求的中间模拟滤波器的系统函数Ha(s)。 3、将Ha(s)展成部分分式形式，利用公式求H(z) 。 （T可以随便选取） * 　　【例6.3.1】 已知模拟滤波器的系统函数Ha(s)为   用脉冲响应不变法将Ha(s)转换成数字滤波器的系统函数H(z)。 　　解 首先将Ha(s)写成部分分式：   极点为  * 那么H(z)的极点为   按照(6.3.4)式，并经过整理，得到：   式中, T是采样间隔，若T选取过大，则会使ω=π附近频谱混叠现象严重。这里选取T=1 s和T=0.1 s两种情况，以便进行比较。设T=1 s时用H1(z)表示，T=0.1 s时用H2(z)表示，则 * 转换时，也可以直接按照(6.3.14)、(6.3.16)式进行转换。首先将Ha(s)写成（6.3.15）式的形式，令极点　s1, 2= －σ1±jΩ1，则 * 再按照（6.3.16）式，H(z)为   将T=1 s、 T=0.1 s分别代入H(z)中，得到H1(z)和H2(z)，其结果和前面得到的H1(z)、 H2(z)完全一样。将Ha(jΩ)、H1(ejω)和H2(ejω)的幅度特性用它们的最大值归一化后，如图6.3.3所示。由图6.3.3（a）可见，模拟滤波器Ha(s)通带很窄，但阻带衰减慢，拖了很长的尾巴，不是带限滤波器。图6.3.3（b）表示的是两种采样频率（T=0.1 s, 1 s），转换成数字滤波器的损耗函数，它的横坐标是对π归一化的数字频率。 * 图6.3.3 例6.3.1中不同采样频率转换结果对比 * 　　【例6.3.2】 用脉冲响应不变法设计数字低通滤波器，要求通带和阻带具有单调下降特性，指标参数如下: ωp=0.2π rad， ? p=1 dB, ωs=0.35π rad, ? s=10 dB。 　　解 例6.3.1仅是将给定的模拟滤波器转换成数字滤波器，本例才是用脉冲响应不变法设计数字滤波器。根据间接设计法的基本步骤求解。 　　（1） 将数字滤波器设计指标转换为相应的模拟滤波器指标。设采样周期为T，由（6.3.7）式得到：  * 　　(2） 设计相应的模拟滤波器，得到模拟系统函数Ha(s)。根据单调下降要求，选择巴特沃斯滤波器。设计过程与例6.2.1完全相同，求出阶数N=4。求解计算留做读者练习。 　　（3） 按照（6.3.1）和（6.3.4）式，将模拟滤波器系统函数Ha(s)转换成数字滤波器系统函数H(z)：  如上求解计算相当复杂。本例调用MATLAB信号处理工具箱函数进行设计。设计程序ep632.m如下。读者可以改变程序中的T值，观察T的大小与频谱混叠失真的关系。 * 图6.3.4 例6.3.2设计的模拟和数字滤波器的损耗函数 * 6.4 用双线性变换法设计IIR数字滤波器 为了消除脉冲响应不变法带来的混叠影响，我们希望模拟频率与数字频率的映射关系满足如下要求：s平面整个虚轴—→z平面单位圆一周； Ω=0 —→ ω=0 Ω: 0→+∞ —→ ω: 0→π Ω: 0→ -∞ —→ ω: 0→-π 当Ω1从π/T经过0变化到 π/T 时，Ω则由－∞经过0变到+∞，实现了s平面上整个虚轴完全压缩到s1平面上虚轴的±π/T之间的转换。 * 代入s=jΩ, s1=jΩ1, 得到: 再通过　　　　从s1