FIR滤波器设计课件.pptVIP

理想数字带通滤波器 理想的线性相位带通DF的频率特性为： 1 其冲激响应hd(n) 加矩形窗处理后，得到 分析： 若选择相位有相移的理想带通DF频率特性为： 此时的脉冲响应 加矩形窗处理后 分析 此时的h(n)一定为反对称序列 当N为奇数时，对应第三种线性相位， 当N为偶数时，对应第四种线性相位， 7.4 加窗对系统频率响应的影响 根据频域卷积定理，加窗后，滤波器的频率响应 现在我们以低通滤波器为例来讨论： 加窗后，频率响应发生了什么变化 加什么样的窗，可以使变化减至最小 7.4.1 矩形窗 矩形窗口的频率特性为 用幅度响应和相位响应的乘积表示为 矩形窗（2） 当w很小时， 当w很大时，WR(w)为周期函数 主瓣 矩形窗处理后的频率响应 根据频域卷积定理可得 Wr(w-wc) Wr(w-wc+2 π/N) Wr(w-wc-2 π/N) 加窗后的低通滤波器频谱 几个特殊频率点 w=0处，响应值 为窗函数频谱Wr(w-θ)和理想低通滤波器频率特性Hd(θ)的乘积的积分，可近似看作Wr(θ)在- π到π的全部积分面积 w=wc处， Hd(θ)刚好与Wr(w-θ)的一半重叠，因此H(wc)=0.5H(0) w=wc-2 π/N处， Wr(w-θ)的全部主瓣在Hd(θ)的通带之内，因此卷积结果有最大值，频率响应出现正肩峰 w=wc+2 π/N处， Wr(w-θ)的全部主瓣在Hd(θ)的通带之外，通带内的旁瓣负的面积大于正的面积，因此卷积结果有最负值，频率响应出现负肩峰 几个特殊频率点（2） 当wwc+2 π/N后， Wr(w-θ)的左边旁瓣的起伏部分扫过通带，卷积值围绕零值而波动 当wwc+2 π/N时， Wr(w-θ)的右旁瓣进入通带，卷积值围绕H(0)而波动 加矩形窗对理想低通滤波器的影响 使理想频率特性不连续点处边沿加宽，形成过渡带，过渡带的宽度等于窗的频谱主瓣宽度4π/N 在截止频率wc的两边 处，H(w)出现肩峰，肩峰的两侧形成起伏振荡，振荡幅度取决于旁瓣相对幅度，振荡多少，取决于旁瓣的多少 增加截取长度，则主瓣附近的窗的频率响应 可见改变N，只能改变窗的主瓣宽度，w坐标的比例和Wr(w)的绝对值大小，而不能改变主瓣与旁瓣的相对比例 矩形窗的频谱示意图 各种窗函数 矩形窗阶段造成府肩峰为8.95%，阻带最小衰减为21dB，不符合工程需要 为了加大阻带衰减，只能改善窗函数形状，使窗谱尽量逼近冲击函数，即绝大部分能量集中在频谱中点 一般窗函数满足两项要求： 窗谱主瓣尽可能的窄，以得到较陡的过渡带 尽量减少最大旁瓣的相对幅度 一般而言，上面两项要求不能同时满足 矩形窗 三角形(BARTLETT)窗 升余弦窗（汉宁Hanning窗）-1 左移 右移 倒余弦 升余弦窗（汉宁Hanning窗）-2 升余弦窗（汉宁Hanning窗）-3 由于频谱是由三个互有频移的不同幅值的矩形窗函数相加而成，这样使旁瓣大大抵消，从而能量相当有效地集中在主瓣内。 其代价:主瓣加宽一倍，可达到减少肩峰，余振，提高阻带衰减。缺点：过滤带加大 改进的升余弦窗(汉明Hanning窗)-1 其频谱函数为 其幅度函数为 III类线性相位：反对称脉冲响应，M为奇数 这种情况下，beta=pi/2，alpha=(M-1)/2是整数 h(n)=-h(M-1-n)， 0≤n≤M-1 Hr(0)=Hr(pi)=0, 因此这种滤波器不适合设计低通或高通滤波器 exp(jpi/2)=j，这种特性非常适合设计希尔伯特变换器和微分器 IV类线性相位：反对称脉冲响应，M为偶数 这种情况和II类似，有 Hr(0)=0 and exp(jpi/2)=j. 因此这种类型适合用于设计数字希尔伯特变换器和微分器 MATLAB实现 Hr_type1:求I类线性相位的Hr(w) 调用格式：[Hr,w,a,L]=Hr_type1(h) Hr_type2:求II类线性相位的Hr(w) 调用格式：[Hr,w,b,L]=Hr_type2(h) Hr_type3:求III类线性相位的Hr(w) 调用格式：[Hr,w,c,L]=Hr_type3(h) Hr_type4:求IV类线性相位的Hr(w) 调用格式：[Hr,w,d,L]=Hr_type4(h) 小结 了解了线性相位FIR滤波器的各种特性，便可根据实际需要选择合适的FIR滤波器，同时设计时要遵循有关约束条件。 如：第3、4种情况，对于任何频率都有固定的π/2相移，一般微分器及90°相移器采用这两种情况，而选频性滤波器则用第1、2种情况。 (1)设计线性相位的低通Digtal Filter 从幅度特性考虑，只能选择第1种或第2种 第一种： 第二种 (2)设计线性相位的