利用FFT成功实现拓扑识别（六）--为自己设计的滤波器做仿真验证

为自己设计的滤波器做仿真验证

• 前文中我们已经设计好了2、4和6阶的数字带通滤波器，实际效果如何昵？可以通过mathlab工具进行仿真和验证

第一步：导出滤波器

• 在mathlab的fdatool中设计好滤波器后，首先要对其幅频特性，带宽等基本特性做到心中有数，无误后，选择 文件–>导出———>选择对象，默认为Hd的名字,如下图：

这样，我们设计的滤波器就医Hd的变量名称导出了，

第二步：mathlab引入滤波器

• 名称为Hd的滤波器变量可以在mathlab的工作区使用，首先通过tf函数从Hd变量中得到滤波器系数：

   							[b,a] = tf(Hd)
  

再通过filters函数就可以实现滤波了。

								Sf=filter(b,a,S)

如下图：

滤波后，可以plot绘出滤波前后的波形对比

可以看出，滤波器，833hz的信号叠加着50Hz工频信号上，而滤波后，工频信号基本不存在了。

第三步：FFT验证滤波效果

我们利用滤波完成的信号再次进行FFT运算，结果如下：

50hz的频谱已经从结果中消失，效果还是特别明显的。