实验三：FIR数字滤波器设计与软件实现.pptVIP

2012 12 11 FIR数字滤波器设计与软件实现 主讲 ：王国林 1．实验目的 （1）掌握用窗函数法设计FIR数字滤波器的原理和方法。 （2）掌握用等波纹最佳逼近法设计FIR数字滤波器的原理和方法。 （3）掌握FIR滤波器的快速卷积实现原理。 （4）学会调用MATLAB函数设计与实现FIR滤波器。 FS=1000 T=1/FS 用窗函数法或者等波纹最佳逼近法设计FIR滤波器hn 对信号xt滤波：yt=fftfilt(hn,xt) Xt=xtg 产生信号xt，并显示xt及其频谱 1、计算并绘图显示滤波器损耗函数 2、绘图显示滤波器的输出信号yt及其频谱 实验内容及步骤 （1）认真复习第七章中用窗函数法和等波纹最佳逼近法设计FIR数字滤波器的原理； （2）调用信号产生函数xtg产生具有加性噪声的信号xt，并自动显示xt及其频谱，如图10.5.1所示； 图10.5.1 具有加性噪声的信号x(t)及其频谱如图 （3）请设计低通滤波器，从高频噪声中提取xt中的单频调幅信号，要求信号幅频失真小于0.1dB，将噪声频谱衰减60dB。先观察xt的频谱，确定滤波器指标参数。 （4）根据滤波器指标选择合适的窗函数，计算窗函数的长度N，调用MATLAB函数fir1设计一个FIR低通滤波器。并编写程序，调用MATLAB快速卷积函数fftfilt实现对xt的滤波。绘图显示滤波器的频响特性曲线、滤波器输出信号的幅频特性图和时域波形图。 （5）重复（3），滤波器指标不变，但改用等波纹最佳逼近法，调用MATLAB函数remezord和remez设计FIR数字滤波器。并比较两种设计方法设计的滤波器阶数。 窗函数法设计FIR滤波器 窗函数法是通过对理想滤波器的单位取样响应加窗来逼近理想滤波器的。 函数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的相应函数生成。 前面介绍了FIR滤波器的优化设计方法，通过迭代的方法求解FIR滤波器，过程十分复杂。 在Matlab中，可以调用函数remez实现滤波器的设计。 b= remez(n, f, m) 函数remez采用Parks-McClellan算法设计线性相位FIR滤波器，n为滤波器阶数，其幅频特性由f和m指定。 实验结果与分析 具有加性噪声的信号x(t)及其频谱如图 分析： 用窗函数法设计滤波器，滤波器长度 Nb=184。滤波器损耗函数和滤波器输出yw(nT)分别如图所示。 用等波纹最佳逼近法设计滤波器，滤波器长度 Ne=91。滤波器损耗函数和滤波器输出ye(nT)分别如图所示。 两种方法设计的滤波器都能有效地从噪声中提取信号，但等波纹最佳逼近法设计的滤波器阶数低得多，当然滤波实现的运算量以及时延也小得多，从图10.5.3(b)和（d）可以直观地看出时延差别。