一种风机监测抗抖动校正方法与流程

技术特征：
1.一种风机监测抗抖动校正方法，其特征在于，包括如下步骤：步骤s1：获取风机监测时的数据流，利用长度值为l的窗口截取数据流得到截取数据，截取数据包括当前历元数据和历史历元数据；步骤s2：对截取数据进行修正，输出当前历元的校正数据；步骤s3：设定在计算下一历元校正数据时的窗口长度，设定规则如下：步骤s3.1：若当前历元的窗口长度值l＞2，则执行步骤s3.2，否则执行步骤s3.3；步骤s3.2：通过式6)得到敏感门限，通过式7)得到不同历元的第一数据差值和第二数据差值，式6)以及式7)如下：thres＝n*std(f')
ꢀꢀꢀꢀ
6)；其中，thres表示敏感门限；n表示倍数；std表示返回截取数据的标准差；f'表示截取数据；δf1'表示第一数据差值；δf
′2表示第二数据差值；f'(k)表示截取数据中当前历元为k的数据；若δf1'＜thres＜δf
′2，取当前历元的窗口长度值l＝1并返回步骤s1进行下一历元的计算，否则进入步骤s3.3；步骤s3.3：当l＜l
max
时，当前历元的窗口长度值l进行一个单位长度的增加后返回步骤s1进行下一历元的计算；l
max
表示设定的窗口长度最大值；当l≥l
max
时，当前历元的窗口长度值l不变并返回步骤s1进行下一历元的计算。2.根据权利要求1所述的风机监测抗抖动校正方法，其特征在于，所述步骤s2包括步骤s2.1、步骤s2.2以及步骤s2.3；步骤s2.1，对截取数据进行均值处理，得到均值；均值和截取数据进行差分得到残差；步骤s2.2，对残差进行处理得到傅里叶变换级数；步骤s2.3，利用均值和傅里叶变换级数进行叠加得到修正序列，根据修正序列输出校正数据。3.根据权利要求2所述的风机监测抗抖动校正方法，其特征在于，所述步骤s2.1中，通过式1)对截取数据进行均值处理，得到截取数据的均值，式1)如下：其中，mean表示截取数据的均值；l表示当前历元窗口的长度值；f'(k)表示截取数据中当前历元为k的数据。4.根据权利要求2所述的风机监测抗抖动校正方法，其特征在于，所述步骤s2.1中，通过式2)得到残差，式2)如下：f0＝f'-mean
ꢀꢀꢀꢀꢀ
2)；其中，f0表示残差；f'表示截取数据；mean表示截取数据的均值。5.根据权利要求2所述的风机监测抗抖动校正方法，其特征在于，所述步骤s2.2中，傅里叶变换级数通过式3)得到，式3)如下：
其中，ft(t)表示傅里叶变换级数；a0表示n＝0时的余弦分量振幅；level表示傅里叶分解层级；a
n
表示各级余弦分量振幅；b
n
表示正弦分量振幅；n表示级数；ω表示角频率；t表示正余弦函数时间轴。6.根据权利要求5所述的风机监测抗抖动校正方法，其特征在于，所述步骤s2.2中，level、ω以及t通过式4)计算得到，式4)如下：其中，ceil表示返回不小于目标值的整数；l表示当前历元窗口的长度值；p为比例系数；num表示离散数据点的数量；t表示傅里叶变换级数采用的拟合正余弦函数的周期长度。7.根据权利要求6所述的风机监测抗抖动校正方法，其特征在于，所述比例系数p∈[0.5，1]。8.根据权利要求2所述的风机监测抗抖动校正方法，其特征在于，所述步骤s2.3中，根据式5)得到修正序列，取修正序列的最后一位作为当前历元的校正数据输出，式5)如下：ft'(t)＝ft(t) mean
ꢀꢀꢀꢀꢀ
5)；其中，ft'(t)表示修正序列；mean表示截取数据的均值；ft(t)表示傅里叶变换级数。9.根据权利要求1所述的风机监测抗抖动校正方法，其特征在于，所述步骤s3.2中，所述倍数n为5-8。10.根据权利要求9所述的风机监测抗抖动校正方法，其特征在于，所述步骤s3.3中，l
max
为15-20。

技术总结
本发明提供了公开一种风机监测抗抖动校正方法，本发明中自适应调节长度值的窗口，能够对数据进行平稳化处理，得到风机监测的全部历元的校正数据，还能进行断点识别和数据断层的应对，即本发明增加了窗口长度以及敏感门限自适应功能和数据平稳化预处理功能，对数据进行预处理，使得被校正后的数据满足后续处理的条件，并且敏感门限跟随实时数据标准差而设定。定。定。


技术研发人员：匡宇龙 梁晓东 朱芳政 刘建峰 刘柯 雷孟飞 汤金毅 李荣学 孙永旭 吴勇生
受保护的技术使用者：湖南联智科技股份有限公司
技术研发日：2022.05.07
技术公布日：2022/7/29