陷波器设计概述

FPGA教程目录

MATLAB教程目录

-----------------------------------------------------------------------

        陷波器是一种特殊的带阻滤波器,其阻带在理想情况下只有一个频率点，因此也被称为点阻滤波器。这种滤波器主要用于消除某个特定频率的干扰，由于陷波器频率特性的特殊性，它除了可采用双线性变换进行设计外，还可以采用所谓零极点配置的方法进行设计。

       一个理想的陷波滤波器的频率特性要在消除的信号频率点处，其值等于零；而在其他频率点处，其值等于1。由于数字滤波器的频率特性就是其单位冲激响应在单位圆上的Z变换，因此只需要在单位圆上相应于所需带阻滤波器阻带位置的频率处设置零点，就可以使滤波器的频率特性在所需阻带频率处为零。但是仅仅进行零点设置只考虑到了滤波器的阻带特性。为了得到非常陡峭的过渡带和常数幅度的通带特性，必须在Z平面上为每一个零点再配置一个相应的极点。Z平面单位圆附近的零点会在滤波器幅频特性的相应频率处产生陷落，零点离单位圆越近，陷落越深；而Z平面单位圆附近的极点会在滤波器幅频特性的相应频率处产生凸峰，极点离单位圆越近，凸峰越高。因此在完成了零点的配置后，为了抵消零点引起的陷落对滤波器通带范围内幅频特性的影响，还需要再配置相应的极点，由于滤波器稳定性的要求，极点必需配置在单位圆内，显然极点离单位圆越近则极点对零点的抵消作用越明显，得到的滤波器的阻带就越窄，过渡带就越陡峭。

      本函数的主要过程就是设计一个陷波滤波器，将叠加的60hzsin信号滤除。

       当带阻滤波器的阻带很窄时，又称为陷波滤波器。陷波滤波器的设计是以模拟滤波器为原型，通过一定变换转换为数字滤波器。陷波滤波器的设计有两种一是先利用模拟频域带阻变换法，再利用数字化法设计数字带阻型滤波器。方法二是直接从模拟低通原型滤波器通过s平面变换成数字带阻滤波器的z平面。

    这里假设模拟低通到模拟带阻的变换关系为：

图7 陷波滤波器频谱特性

从上面的仿真结果可以看到，陷波滤波器在60hz出凹陷，从而将60hz处的叠加信号滤除。通过陷波滤波器之后，其仿真波形如下所示：

图8 EMG信号滤除sin之后效果对比

    设计一个MATLAB函数，通过傅里叶变换得到EMG信号60hz附近的频谱图，这个函数的输入信号为EMG带干扰信号，输出为陷波滤波之后的信号，做傅里叶变换，得到其滤波前后的频谱图效果。

    这个问题的内容和函数二相似，其主要是将输入的信号数字化，并通过DFT变化通过陷波滤波之后的信号的频谱变化，这里对其理论不多做介绍了，通过仿真，其结果如下所示：

图8 陷波滤波之后的频谱图

    通过陷波滤波之后，信号的频谱图在60hz出叠加的sin分量基本被滤除。

    设计一个比对函数，将原始的EMG信号和通过陷波滤波之后的信号进行比较，分析器RMS值。

该模块为信号对比模块，就是将原始输入的EMG信号和叠加sin的EMG信号通过陷波滤波后的信号进行对比，通过分析RMS来判断其正确性，该系统的主要功能将两个输入的信号相减，并计算器均方根值。对比陷波滤波前后的波形图：

图9 陷波滤波前后的仿真效果对比图

从上图的仿真结果可以看到，滤波之后的仿真图和正常德EMG信号十分接近，在MATLAB中通过如下的代码计算RMS值：

RMS=sqrt(sum((EMG_vector-EMG_notch).^2)/length(EMG_vector));

运行代码，可以得到RMS值为：0.38。而将滤波之前的信号与标准的信号进行计算，其RMS值为：0.701，由此可见，通过陷波滤波之后，信号基本和正确的测试信号相似。