测试系统技术 教学课件 作者 郭军 第1－8章 第6章.pptVIP

　　除了巴特沃思滤波器以外，在模拟滤波器中还有很多低通滤波器原型可以供数字滤波器设计时使用，例如切比雪夫(Chebyshev)滤波器、贝赛尔(Bessel)滤波器等。通过一些简单的变换，还可以将原型低通滤波器转化为高通或带通滤波器。但是由于在使用冲激响应不变变换法设计滤波器时，为了保证滤波器的幅频特性不会产生混叠，原型滤波器在ωs/2以上的幅频特性应该等于0(或者接近于0)，所以这种方法不适于实现高通滤波器的设计，这时可以用 下面介绍的双线形变换法进行设计。有关这些内容，这里就不深入介绍了，有兴趣的读者可以参阅有关滤波器综合以及数字信号处理方面的书籍。 　　2. 双线性变换法 　　双线性变换法(bilinear design method)的基本思路是将连续系统微分方程通过数值积分近似，导出与微分方程相近的差分方程，从而完成离散系统设计。它可以克服冲激响应不变变换法中的频率混叠问题。假设有一个一阶的连续时间系统，其微分方程为 ra′(t)+a0ra(t)=b0ea(t)　 　　　　(6.25) 其传递函数为 (6.26) 根据式(6.25)，可以将ra(t)表示成ra′(t)的积分: (6.27) 　　用t0=kT，t=(k+1)T代入上式，可以得到 (6.28) 　　当T很小的时候，式(6.28)中的积分可以用数值积分来代替，其中最简单的算法就是用直线逼近时间区间［kT，(k+1)T］内的函数ra′(τ)。这时，式(6.28)可以近似表示成 (6.29) 　　根据式(6.25)，又可以得到 　　　　ra′a(kT)=b0ea(kT)－a0ra(kT) 　　　　ra′[(k+1)T]=b0ea[(k+1)T]－a0ra［(k+1)T］ 将这两个等式代入式(6.29)，可以得到 将式子变换可得 上式实际上就是一个一阶差分方程。假设r(k)=ra(kT)，e(k)=ea(kT)，则有 对该等式两边同时求Z变换，可以得到 故 (6.30) 这就是与Ha(s)对应的数字滤波器的传输函数。 　　将式(6.30)表示的离散时间系统的传输函数与式(6.26)表示的连续时间系统的传输函数相比较，可以得到两者的对应关系: (6.31) 也就是说，通过映射关系 (6.32) 　　可以直接由连续时间系统的传输函数Ha(s)推导出离散时间系统的传输函数。这个映射关系属于双线性变换映射，所以这种滤波器设计方法被称为双线性变换法。虽然这里以一阶系统为例推导出这个映射关系，但是这种映射关系在高阶系统中依然存在。 　　双线性变换实际上是按照映射关系　　　　　将原来一个以复数s为自变量的复变函数Ha(s)映射成一个以z为自变量的复变函数，其反映射关系为 (6.33) 　　令s=σ+jΩ，代入上式，并对等式两边求绝对值，可以得到 (6.34) 　　由式(6.34)可以看出，如果σ＜0，则|z|＜1；如果σ=0，则|z|=1；如果σ＞0，则|z|＞1。 也就是说，这个映射将s平面中虚轴以左部分的点映射到了z平面的单位圆内部，将s平面中虚轴以右部分的点映射到了z平面的单位圆外部，将s平面中虚轴上的点映射到了z平面的单位圆上。由此可见，如果原来连续时间系统Ha(s)是一个稳定系统，其所有的极点都处于s平面的虚轴以左部分，则H(z)的所有极点一定处于z平面的单位圆内，离散时间系统一定也是稳定的。 　　下面再进一步研究两个系统频率响应之间的关系。 　　令s=jΩ，z=ejωT，带入式(6.32)得 故 或 　　Ω与ω之间的映射关系如图6.11所示。从图中可以看到，这种映射关系是一个一一对应的单调映射。当频率Ω或ω较小时，映射基本上呈线性关系，Ω≈ω。这时原来系统的频率特性可以得到保证。但是当频率较大时，映射出现了非线性关系，使得原来的连续系统在高频段的频率响应空 间产生了压缩，当Ω→∞时，ω→ωs/2。也就是说，映射到离散时间系统的频率分量不会超过ωs/2，所以，无论采样频率取什么值，都不会出现频谱的混叠现象，这是双线性变换法的优点。但是，由此带来的问题是离散时间系统的频率响应不再与连续系统一致，特别是在高频部分，频率特性会产生非线性失真。这一点在设计时必须加以考虑。一般的方法 图6.8 线性相位直接型结构 　　除了满足式(6.7)的FIR滤波器以外，还有一种FIR滤波器相位特性同样也是直线，只不过该直线不过原点。通常也可以将它纳入到线性相位FIR滤波器范畴内。这种FIR滤波器的单位冲击响应函数的波形关于　　　奇对称，即满足 h(k)=－h(N－1－k)　　　　　　　　(6.8) 　基于与上面的定理相同的分析过程，可以证明其相频特性为 (6.9) 　　由式(6.9)可见，这种系统的相频特性上多了π/2的相位移，所以这