《信号与系统课件--第七章有限长单位脉冲响应(FIR)数字滤波器的设计方法》课件.ppt

* * * * * * * * * * * * * * * * * * * * * * * * * * 四，窗函数法设计步骤 1，给定所要求的 2，先求得 或者使用计算机：在给定 当M很大时，有 频域抽样导致时域序列周期(M)延拓: 3，由过渡区宽度??以及阻带最小衰减的要求，选定窗函数的形状，并估计N值，一般N要试探几次最后确定。 一般 4，求得所设计的FIR的 5，求得 （还应验算是否满足设计要求，不行时则改变N或变换窗函数） 其中A 与窗函数有关，如矩形窗的 A=4? 例 1 用矩形窗设计一个线性相位带通滤波器 -ωc≤ω-ω0≤ωc ０≤ωω0-ωc, ω0+ωcω≤π （1） 设计N为奇数时的h(n)。  （2） 设计N为偶数时的h(n)。  （3） 若改用海明窗设计，求以上两种形式的h(n)表达式。 五 举例 解 根据该线性相位带通滤波器的相位 可知该滤波器只能是h(n)=h(N-1-n)即h(n)偶对称的情况，h(n)偶对称时，可为第一类和第二类滤波器，其频响 （1）当N为奇数时，h(n)=h(N-1-n)，可知H(ejω)为第一类线性相位滤波器，H(ω)关于ω=0, π, 2π有偶对称结构。题目中仅给出了Hd(ejω)在 0～π上的取值，但用傅里叶反变换求hd(n)时， 需要Hd(ejω)在一个周期［-π,π］或［0, 2π］上的值，因此, Hd(ejω)需根据第一类线性相位滤波器的要求进行扩展，扩展结果为 ω0-ωc≤ω≤ω0+ωc, -ω0-ωc≤ω≤-ω0+ωc -ω0+ωcωω0-ωc, -π≤ω-ω0-ωc, ω0+ωcω≤π 则 h(n)=hd(n)RＮ(n) (2) N 为偶数时，H(ejω)为第二类线性相位滤波器，H(ω)关于ω=0 呈偶对称。所以, Hd(ejω)在［-π,π］之间的扩展同上， 则hd(n)也同上，即: (3) 若改用海明窗 则 N为奇数时 N为偶数时 上面两个表达式形式虽然完全一样，但由于N为奇数时， 对称中心点α=(N-1)/2 为整数，N为偶数时，α为非整数，因此N在奇数和偶数情况下，滤波器的单位脉冲响应的对称中心不同， 在 0≤n≤N-1 上的取值也完全不同。 例 2 根据下列技术指标，设计一个FIR低通滤波器。  通带截止频率ωp=0.2π，通带允许波动Ap=0.25dB；  阻带截止频率ωs=0.3π，阻带衰减As=50dB。  解 查表可知，海明窗和布拉克曼窗均可提供大于 50 dB的衰减。但海明窗具有较小的过渡带从而具有较小的长度N。 根据题意，所要设计的滤波器的过渡带为 由表可知，利用海明窗设计的滤波器的过渡带宽Δω=8π/N，所以低通滤波器单位脉冲响应的长度为 3 dB通带截止频率为 理想低通滤波器的单位脉冲响应为 海明窗为 则所设计的滤波器的单位脉冲响应为 N=80 所设计的滤波器的频率响应为 利用计算机编程实现，结果如图所示。图（a）是理想低通滤波器的单位脉冲响应hd(n); 图 (b)是海明窗函数; 图 (c) 是实际低通滤波器的单位脉冲响应h(n); 图 (d)是实际低通滤波器的幅频特性|H(ejω)|，以dB为单位。滤波器长N=80，实际阻带衰减为As=53dB，通带波动为Ap=0.0316 dB，均满足设计要求。 图 例2中低通滤波器设计结果 §7.3 频率抽样设计法 一、基本思想方法 ? 窗函数法------从时域开始去设计 ?频率抽样法---直接从频域去设计 ? 有限长 代替 无限长 代替 ?有限多点 无限多点 逼近 力求 插值恢复 F变换 重要的是窗函数 重要的是插值函数 二、频率抽样设计法 1，由 抽样得到N个H(k)值 2，由N个H(k)唯一确定有限长的h(n) 3，由H(k)值按插值公式恢复 内插结果 平坦段：误差不大 突变段：误差很大 过渡采样点 例 3 利用频率采样法，设计一个线性相位低通FIR数字滤波器，其理想频率特性是矩形的 0≤ω≤ωc 其他 已知ωc=0.5π，采样点数为奇数N=33。试求各采样点的幅值Hk及相位θk，也即求采样值H(k)。 解 N=33, 且低通滤波