Notchfilter.docVIP

Notch filter 程序说明 在信号处理中经常要把某些单频 窄带 干扰信号去除，例如求系统采集信号中的把工频信号滤除。实际上有一个很好的方法，便是使用陷波器。在附件中给出了陷波器的设计技术，并举了例子。 陷波器是无限冲击响应 IIR 数字滤波器，该滤波器可以用以下常系数线性差分方程表示： 1 式中： x n 和y n 分别为输人和输出信号序列；和为滤波器系数。 对式 1 两边进行z变换，得到数字滤波器的传递函数为： 2 式中：和分别为传递函数的零点和极点。 由传递函数的零点和极点可以大致绘出频率响应图。在零点处，频率响应出现极小值；在极点处，频率响应出现极大值。因此可以根据所需频率响应配置零点和极点，然后反向设计带陷数字滤波器。考虑一种特殊情况，若零点在第1象限单位圆上，极点在单位圆内靠近零点的径向上。为了防止滤波器系数出现复数，必须在z平面第4象限对称位置配置相应的共轭零点、共轭极点。 这样零点、极点配置的滤波器称为单一频率陷波器，在频率ωo处出现凹陷。而把极点设置在零的的径向上距圆点的距离为l-μ处，陷波器的传递函数为： 3 式 3 中μ越小，极点越靠近单位圆，则频率响应曲线凹陷越深，凹陷的宽度也越窄。当需要消除窄带干扰而不能对其他频率有衰减时，陷波器是一种去除窄带干扰的理想数字滤波器。 当要对几个频率同时进行带陷滤波时，可以按 2 式把几个单独频率的带陷滤波器 3 式串接在一起。 一个例子：设有一个输入，它由50Hz信号和100Hz信号组成。50Hz是一个干扰信号，要设计一个50 Hz的带陷滤波器，采样频率为400Hz。 因此z平面上的零极点可设置为 展开式为 它的传递函数为 因此分子系数是[1 1.414 1]；分母系数是[1 1.4126 0.999]。 差分方程有 程序清单有 见notch_filter.m 程序说明： freqz_m在很多书中都有介绍.如基于MATLAB的系统分析和设计--信号处理西安电子科技大学出版社维纳_恩格尔数字信号处理--使用MATLAB使用方法查阅一下书本,函数内容是:function [db, mag, pha, grd,w] freqz_m b,a ;%Modified version of freqz subroutine[H,w] freqz b,a,1000,whole ;H H 1:501 ; w w 1:501 ;mag abs H ;db 20*log10 mag+eps /max mag ;pha angle H ;grd grpdelay b,a,w ; 下面一个代码也是一个滤除50Hz及其整倍数的陷波器，供参考function data trapper data,f,sf,n ;for i 1:nfi i*f;omg 2*pi*fi/sf;num [1 -2*cos omg 1];den [1 -2*0.999*cos omg 0.999];data filter num,den,data ;end%%%%%%%%%%sf:采样频率f:要滤除的频率在主程序里调用即可.n大于100效果就差了.若要将49Hz也滤除,0.999是不是要减小?