第4章 节 5-7 数字滤波器的原理和设计方法 数字信号处理 .pptVIP

第四章 数字滤波器的原理和设计方法;4.5 IIR数字滤波器的频率变换;设计这样一些选频滤波器的传统方法有两种。;我们曾经用双线性变换将模拟系统函数变换成数字系统函数，对于低通、高通、带通和带阻型选频滤波器来说，也可以利用一种和双线性变换非常相象的有理变换从低通数字滤波器变换出来。;频率变换中的变量代换公式必须满足下列条件：;2)低通→高通的z平面变换;4)低通→带阻的z平面变换;4.6 FIR数字滤波器的设计方法;FIR滤波器的主要缺点是，必须用很长冲激响应的FIR滤波器才能很好地逼近锐截止滤波器，这意味着需要很大的运算量。 另一个缺点是，线性相位FIR滤波器的时延不一定总是样本间隔的整数倍，在某些信号处理应用中，这种非整数时延会带来一些不希望有的问题。 从以上讨论看出，我们最感兴趣的是具有线性相位的FIR滤波器。对非线性相位的FIR滤波器，一般可以用IIR滤波器来代替，因而这里不去讨论它。 FIR数字滤波器的设计方法与IIR数字滤波器的设计方法很不一样，它不能利用模拟滤波器的设计方法。 FIR数字滤波器的设计方法主要有窗函数法、频率取样法和等波纹逼近法等3种，本章主要介绍窗函数法，也简要地介绍频率取样法。;4.6.1 窗函数法;我们以一个截止频率为ωc的线性相位的理想矩形幅度特性的低通滤波器为例来加以讨论，设时延为α，即;现在需要寻找一个有限长序列h(n)来逼近hd(n)，h(n)应满足FIR滤波器的基本条件，即它是偶对称或奇对称的，以满足线性相位的要求，它还应当是因果的。这样，有;即FIR数字滤波器的频谱函数是理想低通滤波器的频谱函数与窗函数频谱的卷积。采用不同的窗函数，对应的H(ejω)有不同的形状。;理想低通滤波器的频率响应可表示为：;只要看几个特殊的频率点，就可以看出H(ω)的一般情况。特别要注意卷积过程给H(ω)造成的起伏现象。;1、先来看ω=0时零频率处的响应值H(0)。根据式(4.81)，H(0)等于图4.49中(a)与(b)两个函数乘积的积分，即WR(θ)在θ=-ωc到θ=+ωc这一段的面积，当ωc2π/N时(这个条件一般能满足)，H(0)实际上就很近似于WR(θ)全部θ从-π到+π的面积。 2、再看ω=ωc时的卷积值，这时Hd(θ)正好与WR(ω-θ)的一半重叠，如图4.49(d)所示，因此卷积值正好是零频响应H(0)的一半，即H(ωc)/H(0)＝0.5，如图4.49 (f)所示。 3、当ω在通带截止频率ωc以内，即ω=ωc-2π/N时，WR(ω-θ)的整个主瓣都在Hd(θ) 的通带内，如图4.49(d)所示，因此卷积结果有最大值，这时频率响应出现正肩峰。 4、对于ω=ωc-2π/N、WR(ω-θ)的主瓣全部在Hd(θ)的通带外，如图4.49(e)所示，在通带内旁瓣负的面积大于正的面积，因此卷积值达到最大负值，H(ω)在这里出现负肩峰。;5、当ω进一步增大时，卷积值也将随着WR(ω-θ)的旁瓣在通带内的面积的变化而变化，这样就造成H(ω)以零值为中心的上下起伏波动。;从以上分析及图4.49(f)可以看出，理想低通滤波器经加窗处理后，主要受到加窗处理两方面的影响。;以矩形窗为例，它的频谱为;为了减小波纹幅度，一方面可以加大窗的长度N，但效果并不显著；另一方面可采用不同的窗函数来改善不均匀收敛性。图4.50所示的是几种常用的窗函数：;它们的定义式和频谱函数分述如下：;3、汉宁(Hanning)窗(升余弦窗);当N1时，N-1≈N，得到窗谱的幅度函数为;4、哈明(Hamming)窗，又称改进的升余弦窗;其频率响应的幅度函数为;其频谱的幅度函数为;;图4.53所示的是用这5种窗函数设计的低通FIR数字滤波器的频率响应特性。窗函数的长度N＝51，理想低通滤波器的截止频ωc=π/2。 从图中可看出，用矩形窗设计的滤波器的过渡带最窄，但阻带衰减指标最差，仅有-21dB左右。而用布莱克曼窗设计的阻带衰减指标最好，可达-74dB，但过渡带最宽，约为矩形窗的3倍。;;6、凯泽(Kaiser)窗;当n=(N-1)/2(即中点)时，ω(n)=1，当n从中点向两边变化时，ω(n)逐渐减小，β越大，ω(n)变化愈快。当n=0和n=N-1时，ω(0)= ω(N-1)=1/I0(β)。在上图中，β=5.44的曲线接近于哈明窗，β=8.5的曲线与Blackman窗相近，而β=0的曲线就是矩形窗。参数β选得越大，则ω(n)窗越窄，而频谱的旁瓣越小，但主瓣宽度也相应增加。因而改变β值就可对主瓣宽度与旁瓣衰减进行选择。;可以看出，阻带最小衰减只由窗形状决定，不受N的影响，而过渡带的宽度则随窗宽N的增加而减小。;窗函数设计FIR数字滤波器的步骤：;窗函数法计算中的主要问题;;第二个问题的解决办法是用求和来代替积分，以便在计算机上计算，也就是要计