13数字滤波器设计技术分析.pptVIP

* * 常用窗函数的频谱 * * 用窗函数法设计FIR滤波器的步骤 给出希望设计的滤波器的频率响应函数Hd(ejω)。 根据允许的过渡带宽及阻带衰减，选定窗函数和N值。 计算hd(n) 如果Hd(ejω)不能用简单函数表示，可以用求和代替积分。 * * 将hd(n)与窗函数相乘得FIR数字滤波器的冲激响应h(n) 计算FIR数字滤波器的频率响应，并验证是否达到所要求的指标 * * 例 窗函数法设计FIR滤波器 例8.1 用窗函数法设计线性相位FIR低通滤波器，设N=11，ω=0.2?rad 解：理想数字低通滤波器 单位取样响应 * * 要求设计的FIR数字滤波器的单位取样响应 * * 用不同窗函数设计的FIR滤波器 用矩形窗时过渡带最窄，而阻带衰减最小，布莱克曼窗过渡带最宽，但阻带衰减加大。为保证有同样的过渡带，必须加大窗口长度N * * 8.5 IIR与FIR数字滤波器的比较 IIR滤波器存在着输出对输入的反馈，因此可以用比FIR滤波器少的阶数来满足技术指标，这样，IIR滤波器所用的存储单元和所需的运算次数都比FIR滤波器少。 例如用频率抽样法设计阻带衰减为－20dB的FIR滤波器，其阶数要33阶才能达到要求，而如果用双线性变换法设计一个切比雪夫IIR滤波器，只需4-5阶就可以达到同样的指标，所以FIR滤波器的阶数要高5－10倍。 FIR滤波器可得到严格的线性相位，而IIR滤波器则做不到这一点。IIR滤波器的选频特性越好，则相位的非线性就越严重。如果要求IIR滤波器具有线性相位，同时又要求它满足幅度要求，那么就必须用一个全通网络进行相位校正，这必然会大大增加滤波器的节数和复杂性。因此在需要严格线性相位的情况下应该选择FIR滤波器。 * * IIR滤波器必须采用递归结构实现，只有当所有极点都在单位圆内时滤波器才是稳定的。但实际中由于存在有限字长效应，滤波器有可能变得不稳定。而FIR滤波器主要采用非递归结构，因而从理论上以及从实际的有限精度的运算中，都是稳定的。另外，FIR滤波器可以采用快速傅立叶变换（FFT）来实现，在相同阶数下，运算速度可以快得多 * * IIR滤波器可利用模拟滤波器现成的设计公式、数据和表格，因而计算工作量较小，对计算工具要求不高。FIR滤波器没有现成的设计公式，窗函数法只给出窗函数的计算公式，但计算通带和阻带衰减仍无显式表达式。一般，FIR滤波器的设计只有计算机程序可以利用，因此对计算工具要求较高，要借助计算机来设计。另外，IIR滤波器主要是设计规格化的、频率特性为分段常数的标准低通、高通、带通、带阻和全通滤波器，而FIR滤波器可设计出理想正交变换器、理想微分器、线性调频器等各种网络，适应性较广。 * * 8.6 FIR数字滤波器的Matlab仿真实现 窗函数法设计FIR滤波器 FIR滤波器的优化设计 * * 8.6.1 窗函数法设计FIR滤波器 窗函数法是通过对理想滤波器的单位取样响应加窗来逼近理想滤波器的。 函数fir1用于设计标准的低通、带通、高通和带阻滤波器。 * * 函数fir1的调用格式 函数fir1的调用格式为 b= fir1(n, Wc, ftype, Windows) 其中，n为滤波器阶数，Wc为截止频率 ftype决定滤波器类型，ftype= high，设计高通FIR滤波器，ftype= stop，设计带阻FIR滤波器。 Windows指定窗函数类型，默认为Hamming窗；可选Hanning、Hamming、Blackman、triangle、bartlett和boxcar窗，每种窗都可以由Matlab的相应函数生成。 * * 例：设计butterworth低通滤波器 例8.3 设计一个15阶FIR低通滤波器，截止频率为0.2π 。 b= fir1(15, 0.2); freqz(b, 1, 512); 函数freqz(b, a, N)用于计算由a和b构成的数字滤波器的频率响应，并用图形方式分别表示其幅度响应和相位响应。 * * 例8.3 程序运行结果 程序运行结果如图8.13所示。 1) 信号消噪 滤波技术是信号消噪的基本方法。根据噪声频率分量的不同，可选用具有不同滤波特性的滤波器。当噪声的频率高于信号的频率时，应选用低通滤波器；反之，选用高通滤波器。当噪声的频率低于和高于信号的频率时，应选用带通滤波器。当噪声的频率处于信号的频率范围时，应选用带阻滤波器。 数字滤波应用 语音信号消噪 200 -200 0 2000 4000 6000 8000 200 -200 0 2000 4000 6000 8000 20 lg| H ( w ) |/dB 0 -50 -100 0 0.5 1 w / 2 p 20lg