设计IIR滤波器,实现对存在加性干扰的时域离散信号进行滤波。
已知带加性干扰的信号用x(n)表示,x(n)=xs(n)+η(n),式中xs(n)是有用的信号,是一个0~0.2πrad的带限信号。η(n)是一个干扰信号,其频谱分布在0.3πrad以上。要求设计一个巴特沃斯IIR数字滤波器对信号x(n)进行滤波,将干扰η(n)滤除。要求在xs(n) 所在的通带内滤波器幅度平坦,在0.2πrad处幅度衰减不大于1dB,在噪声所在的0.3πrad 以上的频带内滤波器幅度衰减大于等于40 dB。
步骤:1.根据题目要求确定要设计的数字滤波器的技术指标(低通滤波器指标:wp=0.2πrad,ws=0.3πrad,αp=1dB,αs=40dB);
2. 用双线性变换法频率转换公式,将DF技术指标转换为巴特沃斯AF的技术指标;
3. 调用MATLAB函数buttord和butter,设计该模拟滤波器;
4. 根据1所确定的技术指标,调用MATLAB函数buttord和butter,直接设计数字滤波器,观察设计结果与上面用双线性变换法的设计结果是否相同。
5. 滤波仿真:调用MATLAB工具箱函数filter对下面给出的带加性干扰的信号x(n)进行滤波,观察滤波效果(绘制滤波前后信号的时域和频域波形)。a. 滤波仿真:调用MATLAB工具箱函数filter对下面给出的带加性干扰的信号进行滤波,观察滤波效果(绘制滤波前后信号的时域和频域波形)。
实验信号x(n)的128个样值:
-0.0289 0.3943 0.9965 1.1266 0.9995 1.0891 1.2262 1.0699
0.8990 0.7685 0.7844 0.9471 1.4317 1.6765 1.7629 -3.2903
3.4122 4.5403 -2.1667 -2.0584 4.6694 2.0368 -0.4864 0.1427
2.5652 -1.8980 -0.0527 -1.4730 2.7884 -6.4092 2.9084 -5.1428
1.5929 0.0014 -0.6010 -4.3059