[发明专利]一种基于SES-LOF的旋转机械监测数据噪点检测方法

权利要求书

1.一种基于SES-LOF的旋转机械监测数据噪点检测方法，其特征在于，包括以下步骤：

1)获得一段旋转机械监测信号作为原始信号x(t)，其中t＝1,…,N，N为该段信号采样点个数；根据幅值大小，依次计算原始信号x(t)数据段中每个点的局部异常因子值；具体为，当t＝j时，点x(j)的局部异常因子值计算表达式如下：

其中，k为LOF算法的参数，其值取为5；

N_k(x(j))是点x(j)的邻域，即x(j)的第k距离以内的所有点，包括第k距离上的点；

lrd_k(o)与lrd_k(x(j))分别为点o与点x(j)的局部可达密度，根据式(2)、式(3)分别计算：

其中，o'为o的邻域，k_distance(o')表示点o'的第k距离，k_distance(o)表示点o的第k距离，d(o,o')表示点o与o'的欧氏距离，d(x(j),o)表示点x(j)与点o的欧氏距离，N_k(o)是点o的邻域，max{·}表示取所有值中的最大值；

2)将局部异常因子值大于1的所有点标记为可疑噪点，设有M个可疑噪点，则得到可疑噪点集合C_p；

3)从原始信号x(t)中剔除一个可疑噪点得到信号y_i(n)：任取一可疑噪点x_i(j)∈C_p,i＝1,…,M，其中i表示可疑噪点集合C_p的第i个点，j表示该点在原始信号x(t)中的位置，则原始信号x(t)剔除点x_i(j)后，得到信号y_i(n),n＝1,…,N-1；

4)对信号y_i(n)进行希尔伯特变换得到

然后对进行傅里叶变换，得到zf_i：

5)计算信号y_i(n)的Shannon熵：由S＝-∑plnp得：

其中，S_i为剔除可疑噪点x_i(j)后，得到信号y_i(n)的Shannon熵，表示求zf_i²的均值；

6)重复步骤3)至步骤5)，直至遍历完集合C_p中所有可疑噪点，得到Shannon熵的一个集合

7)从所有可疑噪点中筛选出真正噪点：根据3σ准则，设置阈值其中std(·)表示计算标准差；若S_i＞T，则S_i对应的可疑噪点x_i(j)为真实噪点，最终根据阈值T检测出所有真正噪点。