[发明专利]陷波滤波器的设计方法

说明书

技术领域

本发明涉及一种陷波滤波器的设计方法，特别是涉及一种能够设计出对风机测试中的50Hz及其二到九次谐波进行滤除的陷波滤波器的设计方法。

背景技术

风机是风能转换为机械能的动力机械，风机运行是否正常直接影响到风力发电的产量，风机故障可能会导致机组本身的损坏，还有可能造成其它一些更严重的后果。因此，为了保障风机的安全运行，对风机的运行状况进行振动状态检测和故障诊断是非常重要的。

风机机组振动状态检测是通过在风机机组中设置振动传感器以对机组振动信号定期进行监测，从而掌握机组的运行状态。而在振动传感器采集信号时的主要噪声源是50Hz的工作频率及其高次谐波，这些噪声常以电磁波的辐射形式对电子设备造成比较严重的干扰，而这种干扰就被称之为工频干扰，严重的会导致风机设备运行异常。

现有技术中虽然存在对工频干扰进行滤除的算法，但是基本都是仅仅对50Hz的基波进行滤除而不能滤除高次谐波，并且都不能保证具有极窄的阻带，从而使得信号中的其它有用的频谱成分也会被滤除、受到损伤，以及输出容易受到初值影响还有容易漏检噪声源等弊端。因此，现有技术中还不能对风机机组中的工频及其高次谐波的干扰进行有效地滤除，严重影响到风机设备的正常运行。

发明内容

本发明要解决的技术问题是为了克服现有技术中对工频干扰进行滤除的算法基本都不能保证具有极窄的阻带，从而使得信号中的其它有用的频谱成分也会被滤除、受到损伤，以及输出容易受到初值影响还有容易漏检噪声源的缺陷，提供一种能够设计出对风机测试中的50Hz及其二到九次谐波进行滤除的陷波滤波器的设计方法。

本发明是通过下述技术方案来解决上述技术问题的：

本发明提供了一种陷波滤波器的设计方法，其特点在于，该设计方法包括以下步骤：

S₁、对K对中的第i对具有相同的相位延时a、截止频率分别为w_i1和w_i2的离散时间理想低通滤波器的时域单位脉冲响应函数进行变换以生成一对长度均为N的非递归低通滤波器的函数h_i1（n）和h_i2（n），其中i、K均为正整数且i≤K，a=（N-1）/2，在每一对离散时间理想低通滤波器中w_i1均大于w_i2；

S₂、计算h_i3（n）=h_i1（n）-h_i2（n），其中h_i3（n）为上限截止频率为w_i1、下限截止频率为w_i2的带通滤波器的函数；

S₃、对全通滤波器的时域单位脉冲响应函数进行变换以生成长度为N的函数h₄（n），计算其中h₅（n）为该陷波滤波器的时域单位脉冲响应函数；

S₄、对h₅（n）进行变换以得到该陷波滤波器。