数字信号处理第7章下方案.ppt

频率响应为： 在线性相位FIR滤波器的设计中，N的奇、偶选择，用第一种还是第二种频率抽样，主要取决于待设计的滤波器，由使用者选择。 实际求 ，可采用离散傅立叶反变换 （IDFT）先求h（n），再求序列的傅立叶变换（离散时间傅立叶变换，DTFT） ，并可用FFT算法。 四种线性相位FIR DF特性 第一种情况，偶、奇，四种滤波器都可设计。 第二种情况，偶、偶，可设计低、带通滤波器，不能设计高通和带阻。 第三种情况，奇、奇，只能设计带通滤波器，其它滤波器都不能设计。 第四种情况，奇、偶，可设计高通、带通滤波器，不能设计低通和带阻。 I）根据截止频率对应的K及线性相位的要求（如对第一类线性相位FIR滤波器，第一种抽样），对给定的频率响应 等间隔抽样确定： 频率抽样法的设计步骤： III）求FIR滤波器的单位抽样响应： II） IV）求出 验证结果。 确定 共轭对称，h（n）为实数，且为 满足线性相位， h（n）也对称（偶或奇） MATLAB例程：当N为奇数时线性相位第一种频率抽样： alpha = -(N -1)/2 kwc=fix(wc*N /(2*pi))+1 Hrs=[ones(1,kwc),zeros(1,N-2*kwc+1),ones(1,kwc-1)] k1 = 0:floor((N-1)/2); k2 = floor((N-1)/2)+1:N-1; angH = [-alpha*(2*pi)/N*k1, alpha*(2*pi)/N*(N-k2)]; H = Hrs.*exp(j*angH); h = real(ifft(H,N)); [Hm,wm] = freqz(h,1,512); mag=abs(Hm); pha=angle(Hm); 过渡带加抽样点，可增大最小阻带衰减，但过渡带宽也变大。 过渡带宽等于频域抽样间隔。 要想同时减小过渡带宽，只有也增加抽样点数N。 增加抽样点数N使抽样点变密，减小过渡带宽度，但增加了计算量。 过渡带加抽样点的值，是一个优化的问题，可通过线性规划的方法得到。 3）过渡带抽样的优化设计 加过渡抽样点（过渡带） (a)一点过渡带； (b)二点过渡带； (c)三点过渡带 频率抽样法设计低通滤波器例子 (a) 和 （N＝33）； (b)设计结果 增加过渡带抽样点，可加大阻带衰减 不加过渡抽样点： 加一点： 加两点： 加三点： 优点： 可以在频域直接设计，适合于最优化设计； 缺点： 抽样频率只能等于 的整数倍， 或等于 的整数倍加上 ，不能 确保截止频率 的自由取值。 频率抽样法特点： 例：利用频率抽样法，设计一个FIR线性相位低通数字滤波器，截止频率为 抽样点数为N=33。 解：根据题意有： 频率抽样后的序列为 。 选 周期偶对称， N为奇数，则： h（n）偶对称。 由于 对称于 ，只取 即 区间，截止频率 对应于 按第一种频率抽样法， N=33，则有 必须周期偶对称 对称中心N/2=16.5 所以频率响应为： 过渡带宽为一个频率采样间隔 ： 最小阻带衰减为：-20dB 抽样点数取为N=33，过渡带无点频率抽样： 对称中心N/2=16.5 频率响应为： 过渡带抽样点增加一点：k=9， 过渡带宽为： 最小阻带衰减为：-40dB 抽样点数取为N=33，过渡带抽样点一点： 频率响应为： 继续增加过渡带抽样点，将导致过渡带宽增加，故需同时增加抽样点数。抽样点数取为N=65，过渡带抽样点二点： k=17， ，k=18， 过渡带宽为： 最小阻带衰减为：-60dB 抽样点数取为N=65，过渡带抽样点二点： * 将给定的理想频率响应 等间隔抽样，即： 根据 就可确定FIR滤波器的h（n）、系统函数 和FIR滤波器的特性： 2、频率抽样设计法 系统函数： 单位抽样响应： 频率响应： 代入内插函数： 频率响应为： 频率响应的抽样： I）在抽样点上，滤波器的实际频率响应与理想频率响应数值相等； II）在抽样点之间，滤波器的频率响应是各抽样点加权内插函数延伸叠加而成，因而有一定的逼近误差。 频率抽样特点： III）在靠近通带边缘的误差较大，而在通带和阻带误差较小。 1）线性相位的约束 如果我们设计的是线性相位的FIR数字滤波器，则其抽样值的幅度和相位一定要满足7.2 节所讨论的约束条件，这些条件已归纳在表7－1之中。 对于第一类线性相位FIR滤波器（ h（n）偶对称，N为奇数），有： 取 也是标量幅度函