利用FFT成功实现拓扑识别（四）--增加数字滤波器

增加数字滤波器

• 在电力的拓扑识别应用中，传感器采集数据不可避免伴随着各种工频干扰信号，如何精确的识别目标信号是考研系统能否实现精确识别的基础，然而，无论采用何种传感器，都会存在不同特性的噪声，用户用电活动会随着应用环境的变化而变化。因此，我们需要通过设计滤波器、针对传感器原始测量数据进行“优化”。本篇基于mathlab的感觉对数字滤波器的设计方法进行阐述，着重针对MATLAB
  FDAtool工具箱的使用方法进行说明

**

mathlab滤波器设计

• mathlab的滤波器设计功能强大，支持FIR、IIR、加窗、等波纹、最小二乘、Butterworth、Chebyshev、椭圆、脉冲整形等多种规格需求的设计，主要通过子模块fdatool工具（后续版本可能升级为filterdisigner）实现。

fdatool的使用

在mathlab中命令行窗口输入fdatool（新建议为filterDesigner）
选择带通滤波器，选择IIR滤波器，指定阶数为6,选择分析—幅值响应：即可看到该滤波器的幅频特性如图。

可是设计的滤波器在实际系统中的应用效果如何？这就需要在实际系统中进行验证。

选择分析，滤波器系数，再编辑—》转换为单节，输出系数

系数的个数为阶数有关，2阶为3个一组，4阶为5个，6阶为7个，分子对应b系数，分母对应a系数，选择目标—》生成c头文件
输出系数头文件：如下为6阶的输出系数。

利用a、b系数写c函数

对应方程式为：

实测效果

：利用滤波后的数据进行FFT运算，结果如下：

可以看出，信号中的50hz工频信号及其多次谐波信号基本滤除了