FIR滤波器的设计培训讲义.ppt

6.3 FIR滤波器设计 IIR数字滤波器的设计方法利用了模拟滤波器的研究成果，设计方法简单而有效，能得到较好的幅度特性，特别是双线性变换法因没有频谱混叠而很受欢迎。但是FIR滤波器的相位特性是非线性的．相位特性不好控制是IIR波波器的缺点。 IIR滤波器若需线性相位，则要采用全通网络进行相位校正，使滤波器设计变得复杂，成本也高。而有限长脉冲响应FIR数字滤波器具有严格的线性相位，任意的幅度特性，工作稳定。另外，任何非因果有限长序列经过一定的延时都能变成因果的有限长序列，因而FIR滤波器总能用因果系统来实现，而且FIR滤波器由于单位脉冲响应是有限长的，可以用FFT算法来计算，从而大大提高运算效率。当然要取得很好的集成特性，FIR滤波器h(n)的阶次要高于IIR滤波器。 6.3.1 线性相位FIR数字滤波器的条件和特点 1. 线性相位条件 对于长度为N的h(n)，传输函数为 H(ejω)线性相位是指θ(ω)是ω的线性函数，即 θ(ω)=τω, τ为常数 第一类线性相位 ? 如果θ(ω)满足下式： θ(ω)=θ0-τω,θ0是起始相位 第二类线性相位 以上两种情况都满足群时延是一个常数，即 下面推导与证明满足第一类线性相位的条件是：h(n)是实序列且对(N-1)/2偶对称，即 h(n)=h(N-n-1) (6-44) ? 满足第二类线性相位的条件是： h(n)是实序列且对(N-1)/2奇对称，即 h(n)=-h(N-n-1) (6-45) (1) 第一类线性相位条件证明： 按照上式可以将H(z)表示为 (2) 第二类线性相位条件证明： 3. 线性相位FIR滤波器零点分布特点 第一类线性相位的系统函数满足 第二类满足 ，综合起来用下式表示： 6.3.2 利用窗函数法设计FIR滤波器 设理想滤波器传输函数为Hd(ejω),hd(n)是与其对应的单位脉冲响应，因此 相应的单位取样响应hd(n)为 我们实际实现的滤波器的单位取样响应为h(n),长度为N，其系统函数为H(z)， 以上就是用窗函数法设计FIR滤波器的思路。 将H(ejω)写成下式： 通过以上分析可知，对hd(n)加矩形窗处理后，H(ω)和原理想低通Hd(ω)差别有以下两点： (1)在理想特性不连续点ω=ωc附近形成过渡带。过渡带的宽度，近似等于RN(ω)主瓣宽度，即4π/N。 (2)通带内增加了波动，最大的峰值在ωc-2π/N处。阻带内产生了余振，最大的负峰在ωc+2π/N处。以上两点就是对hd(n)用矩形窗截断后，在频域的反映，称为吉布斯效应。这种效应直接影响滤波器的性能。通带内的波动影响滤波器通带中的平稳性，阻带内的波动影响阻带内的衰减，可能使最小衰减不满足技术要求。当然，一般滤波器都要求过渡带愈窄愈好。下面研究如何减少吉布斯效应的影响，设计一个满足要求的FIR滤波器。 直观上，增加矩形窗口的宽度，即加大N，可以减少吉布斯效应的影响。只要分析一下N加大时， RN(ω)的变化，就可以看到这一结论是否完全正确。我们讨论在主瓣附近的情况。 在主瓣附近(w接近于0)，由于 ，RN(ω)可近似为 以上分析说明，调整窗口长度N可以有效地控制过渡带的宽度，但不改变旁瓣的相对值。减少带内波动以及加大阻带的衰减只能从窗函数的形状上找解决方法。如果能找到的窗函数形状，使其谱函数的主瓣包含更多的能量，相应旁瓣幅度就减小了；旁瓣的减小可使通带、阻带波动减小，从而加大阻带衰减。但这样总是以加宽过渡带为代价的。 下面介绍几种常用的窗函数。设 h(n)=hd(n)w(n) ? 式中w(n)表示窗函数。 1.