实验五:FIR 数字滤波器设计与软件实现
一、实验目的
(1)掌握用窗函数法设计FIR 数字滤波器的原理和方法。
(2)掌握用等波纹最佳逼近法设计FIR 数字滤波器的原理和方法。
(3)掌握FIR 滤波器的快速卷积实现原理。
(4)学会调用MA TLAB 函数设计与实现FIR 滤波器。
二、实验内容及步骤
(1)调用信号产生函数xtg 产生具有加性噪声的信号xt ,并显示xt 及其频谱。观察xt 的频谱,确定滤波器指标参数,设计低通滤波器,从高频噪声中提取xt 中的单频调幅信号,要求信号幅频失真小于0.1dB ,将噪声频谱衰减60dB 。
(2)根据滤波器指标选择合适的窗函数,计算窗函数的长度N ,调用MATLAB 函数fir1设计一个FIR 低通滤波器,调用函数fftfilt 实现对xt 的滤波。绘图显示滤波器的频响特性曲线、滤波器输出信号的幅频特性图和时域波形图。
(3)改用等波纹最佳逼近法重新设计滤波器。
(4)采样频率Fs=1000Hz ,采样周期T=1/Fs , p 20.24p f ωπ=T =π,αp=0.1dB , s 20.3s f ωπ=T =π,αs=60dB 。
三、实验程序框图:如图一所示
四、实验程序及实验结果
图一 实验程序流程图
图二 stg 信号与噪声