《数字信号处理导论第7章》PPT.pptVIP

第7章 FIR 数字滤波器设计;IIR数字滤波器：;7.1 Fourier 级数法（窗函数法）;设理想低通滤波器的幅频为1，相频为零：;于是：; ;加窗处理后对频率响应的影响：; ; ;加窗函数的影响：;例1.设计低通 FIR DF, 令归一化截止频 率 ＝0.125, N＝10，20，40， 用矩形窗截短。;接上例：N＝10 分别用矩形窗 和Hamming 窗;高通：;令：;令：;例： 理想差分器及其设计;奇对称，纯虚函数;实际相频特性;例： 设计 Hilbert 变换器 ;优点：1. 无稳定性问题； 2. 容易做到线性相位； 3. 可以设计各种特殊类型的滤波器； 4. 方法特别简单。;7.2 窗函数;对窗函数的技术要求： 1. 3 dB 带宽 ：主瓣归一化幅度降到－ 3 dB 时的带宽；或直接用 。令 则 的单位为 ；;常用窗函数：;矩形窗;汉宁窗-布拉克曼窗比较;矩形窗-汉宁窗-布拉克曼窗比较;矩形窗-三角形窗比较;矩形窗-海明窗-凯泽窗比较; 窗函数;窗函数法的设计步骤;公式法：;解：;2）求hd(n);4）确定N 值;5）确定FIR滤波器的h(n);7.3 FIR DF设计的频率抽样法;;可指定：;转移函数、频率响应和给定的 的关系：;原抽样;用插值的方法得到所要的滤波器：; 为偶数：;窗函数法和频率抽样法 设计的 FIR DF 的频率响应都不理想，即通带不够平，阻带衰减不够大，过渡带过宽，频率边缘不能精确指定。因此我们要寻找新的设计方法。此方法即是Chebyshev 最佳一致逼近 法。该方法在数字信号处理中占有重要的定位，是设计 FIR DF 最理想的方法。由于课时限制，此方法不在课堂上讲解，有兴趣的同学可查阅参考书自行学习。;附：IIR和FIR数字滤波器的比较;与本章内容有关的MATLAB文件：;9．fir1.m 用“窗函数法”设计FIR DF。 调用格式： （1）b = fir1(N,Wn); (2) b = fir1(N,Wn,‘high’); (3) b = fir1(N,Wn, ‘stop’); N：阶次，滤波器长度为N＋1； Wn：通带截止频率，其值在0～1之间，1对应 Fs/2 b: 滤波器系数。;对格式（1），若Wn为标量，则设计低通滤波器，若 Wn是1×2的向量，则用来设计带通滤波器，若Wn是 1×L的向量，则可用来设计L带滤波器??这时，格式 （1）要改为: b = fir1(N,Wn, DC-1), 或 b = fir1(N,Wn, DC-0) 前者保证第一个带为通带，后者保证第一个带为阻带。 格式（2）用来设计高通滤波器， 格式（3）用来设计带阻滤波器。 在上述所有格式中，若不指定窗函数的类型，fir1自动选择Hamming窗。;10．fir2.m 本文件采用“窗函数法”设计具有任意幅 频相应的FIR 数字滤波器。其调用格式是： b = fir1(N, F, M); F是频率向量，其值在0～1之间，M是和F相对应 的所希望的幅频相应。如同fir1, 缺省时自动选用 Hamming窗。;N=30,90 时幅频响应响应及理想幅频响应；;11. remez.m 设计Chebyshev最佳一致逼近FIR滤波器、Hilbert变换器和差分器。调用格式是： (1) b=remez(N, F, A); (2) b=remez(N, F, A, W); （3）b=remez(N,F,A,W,‘Hilbert’); (4) b=remez(N, F, A,W, ‘differentiator) N是给定的滤波器的阶次，b是设计的滤波器的系数，其长度为N＋1；F是频率向量，A是对应F的各频段上的理想幅频响应，W是各频段上的加权向量。;F、A及W的指定方式和例7.4.1和7.4.2所讨论过