信号处理四大算法教程第7讲FIR滤波器设计.pptVIP

有限冲击响应滤波器 FIR的设计 数字滤波器 数字滤波器是对数字信号实现滤波的线性时不变系统。数字滤波实质上是一种运算过程，实现对信号的运算处理。输入数字信号（数字序列）通过特定的运算转变为输出的数字序列，因此，数字滤波器本质上是一个完成特定运算的数字计算过程，也可以理解为是一台计算机。描述离散系统输出与输入关系的卷积和差分方程只是给数字信号滤波器提供运算规则，使其按照这个规则完成对输入数据的处理。时域离散系统的频域特性: 其中 、 分别是数字滤波器的输出序列和输入序列的频域特性（或称为频谱特性）, 是 数字滤波器的单位取样响应的频谱，又称为数字滤波器的频域响应。输入序列的频谱 经过滤波后 ,因此，只要按照输入信号频谱的特点和处理信号的目的， 适当选择 ，使得滤波后的 满足设计的要求，这就是数字滤波器的滤波原理。 数字滤波器根据其冲激响应函数的时域特性，可分为两种，即无限长冲激响应(IIR)数字滤波器和有限长冲激响应(FIR)数字滤波器。 FIR数字滤波器的理论和设计 FIR滤波器在保证幅度特性满足技术要求的同时，还有严格的线性相位特性。设FIR滤波器单位脉冲响应h(n)长度为N，其系统函数 H(z)是z-1的N-1次多项式，它在z平面上有N-1个零点，原点z=0是N-1阶重极点，因此，H(z)永远稳定。稳定和线性相位特性是FIR滤波器的特出优点。 1、FIR如何实现线性相位？ 对于长度为N的h(n)，传输函数为： 式中，Hg(ω)称为幅度特性，θ(ω)称为相位特性。 注意： Hg(ω) ≠|H(ejω)|,因为Hg(ω) 是ω的实函数，其值可正可负，而|H(ejω)|总为正值。 H(ejω)线性相位是指θ(ω)是ω的线性函数，即 θ(ω) = -τω， τ为常数 第一类线性相位 θ(ω) = θ0- τω， θ0为起始相位 第二类线性相位 可证明，满足第一类线性相位的条件是：h(n)是实序列且对(N-1)/2偶对称，即h(n)=h(N-n-1) 满足第二类线性相位的条件是：h(n)是实序列且对(N-1)/2偶对称，即h(n)=-h(N-n-1)。 第一类线性相位条件证明： θ(ω) = -τω 将h(n)=h(N-n-1)代入上式 = 令m=N-n-1 = = 将z=ejω代入上式 = = 群时延为 ,满足条件。 第二类线性相位条件θ(ω) = θ0- τω，采用类似方法证明可得： 2.FIR中Hg(ω)的特点 由于h(n)的长度N取奇数还是取偶数，对Hg(ω)的特性有影响，因此对于两类线性相位，可以分为四种情况，如下图： 在ω =π处为0，不能用于高通、带阻 所有都可以 奇对称单位冲击响应 h(n)=-h(N-1-n) 在ω =0,2π处为0，不能用于低通、带阻 在ω =0, π,2π处为0，只能用于带通 FIR滤波器的直接型结构 y(n) h(N-1) h(N-2) h(1) h(0) x(n) Z-1 Z-1 y(n) h(N-1) h(N-2) h(1) h(0) x(n) Z-1 Z-1 Z-1 h(2) 3.利用窗函数法设计FIR滤波器 步骤： 1）如果给出通带、阻带衰减和边界频率的要求，可选用理想滤波器作为逼近函数，从而用理想滤波器的特性作傅立叶变换，求出Hd(n)。 ，|ω|≤ ωc ωc为3dB截止频率 ，ωc ≤ |ω|≤ л 单位取样响应为： 为保证线性相位，取 2)根据对过渡带即阻带衰减的要求，选择合适的窗函数，并估计窗口长度N。 设待求滤波器的过渡带用△ω表示，它近似等于窗函数主瓣宽度。同时因为过渡带△ω近似与窗口长度N成反比(N ≈ A/△ω)，A取决于窗口形式。如，矩形窗A=4л，哈明窗A=8л等，A参数选择参见后面的表格。 按照过渡带及阻带衰减情况，选择窗函数形式的原则是：在保证阻带衰减满足要求的情况下，尽量选择主瓣窄的窗函数。 3）计算h(n) = hd(n) w(n) 其中w(n)即为窗函数。 4）计算 其中x(n-k)为输入序列x(