FIR数字滤波器的设计教学讲义.pptVIP

FIR数字滤波器的设计;群延时为;故有 ;式7.1是使FIR滤波器具有q(w)=－tw 线性相位的必要且充分条件。;h(n)以(N-1)/2为轴呈偶对称;式7.2是使FIR滤波器具有q(w)=b－tw（ b=±p/2）线性相位 的必要且充分条件。 ;要求h(n)序列以n=(N-1)/2为奇对称中心，时延t=(N-1)/2个抽样周期。 当n=(N-1)/2时代入式7.2;直接画网络结构，有N次乘法与N-1次加法;令n=N-1-m;图7.3 线性相位FIR滤波器网络结构（直接型） h(n)为偶对称，N为奇数;将z=ejω 代入，并利用欧拉公式;其中;2π;直接画出网络结构，有N次乘法与N次加法;图7.4 线性相位FIR滤波器网络结构 h(n)为偶对称，N为偶数;提出因子;注意：这种滤波器不能用于高通与带阻， 因为H(p)=0 ，而以上二者在w=p处不为0。;此时;图7.5 线性相位FIR滤波器网络结构 h(n)为奇对称，N为奇数;进行化简;注意：这种滤波器同样不能用于高通与带阻滤波器。;图7.6 线性相位FIR滤波器网络结构 h(n)为奇对称，N为偶数;幅度函数;任何一种线性相位FIR滤波器的群时延恒定。;则;或;⑴ 零点 zi 既不在实轴上，也不在单位圆上，如图ri≠1，wi≠0;其中;jIm(z);这时零点为实数，共轭值就是本身。 因而他们的基本因子为：;jIm(z);线性相位FIR滤波器的零点只能有以上四种情况;§7.2 窗口法（傅氏级数法）;此冲激响应是无限长，但要求的是有限长冲激响应滤波器，要由hd(n)得到FIR滤波器的冲激响应h(n)，最直接的方法就是将hd(n)截短，即令（假设N为奇数）;这样便得到了FIR数字滤波器的冲激响应h(n)，但此滤波 器的频率响应H(ejω)肯定与理想滤波器的频率响应Hd(ejω) 有差异。（因为h(n)与hd(n)的差异）;级数求和后利用欧拉公式，得到;图中阴影所示面积，即为积分的值，当w变化时， 此曲线左右移动，此面积也就发生变化。;当w逐渐增大，随着图中不同正负，不同大小的旁瓣 移出和移入积分区间，使得H(ejw)的大小产生波动。;w继续增大，主瓣开始移出积分区间，因此H(ejω) 迅速下降， 进入过渡带。;即主瓣的中心移到了wc处，此时区间内曲线下的面积近似等于w=0时面积的一半，因此H(ejwc)≈ H(ej0)/2 ;⑸ w继续增大到p，H(ejw)随着区间内旁瓣的移动而在阻带内 波动 ;由图可知，加了矩形窗后所到的FIR数字滤波器的频率响应H(ejω)与理想的频响Hd(ejω)之间产生了差异，表现在H(ejω)出现了肩峰、过渡带以及通带和阻带内的波动，这就是所谓的吉布斯现象。我们当然希望肩峰和波动尽可能小，过渡带尽可能窄，这样才能更接近理想特件。;出现的这些差异与哪些出素有关？;因此，窗口法设计FIR滤波器;因此，我们只需考察w(n)和W(w)的表示式即可。;以上为对称中心在n＝0处的非因果矩形窗;1;1;1;频谱特性如图，由于这三部分频谱的相加，使总频谱的旁瓣大大抵销，从而使能量有效地集中在主瓣内，但其代价是使主瓣与矩形窗主瓣相比加宽了一倍，为8p/N;其幅度频谱为：;汉宁窗a＝0.5；对于哈明窗a＝0.54。;布莱克曼;矩形窗函数的傅氏变换 （N=51） 旁瓣峰值衰减 -13dB 主瓣宽度4p/N （书表达不准）;三角形窗函数的傅氏变换 （N=51） 旁瓣峰值衰减 -25dB 主瓣宽度8p/N;汉宁窗函数的傅氏变换 （N=51） 旁瓣峰值衰减 -31dB 主瓣宽度8p/N;哈明窗函数的傅氏变换 （N=51） 旁瓣峰值衰减 -41dB 主瓣宽度8p/N;布莱克曼窗函数的傅氏 变换(N=51) 旁瓣峰值衰减 -57dB 主瓣宽度12p/N;以上几种窗函数都是以一定的主瓣加宽为代价来换取某种程度的旁瓣抑制。;四、设计方法;n;n;例7-2-1：设计一个线性相位FIR低通滤波器，给定抽样频率为Ws=2p×1.5×104(rad/s)，通带截止频率为Wp=2p×1.5×103(rad/s)，阻带截止频率Wst=2p×3×103(rad/s)，阻带衰减d2不小于-50dB。幅度特性如图：;理想滤波器截止频率Wc为两个肩峰值处频率的中点， 所以可近似有 Wc≈ (Wp + Wst)/2= 2p×2.25×103(rad/s);哈明窗过渡带宽为 Dw=wst -wp = 6.6p/N = 0.4p-0.2p = 0.2p，求得N = 33 t=(N-1)/2=16;(6)由h(n)求H(ejω)，检验各项指标是否满足要求，如不满足要改变N或改变窗形状重新计算。 H(ejω)图形如下，满足要求。;是全通滤波器，因此一个高通滤波器相当于一个全通滤波器减去一个低