一种基于交叉熵的水电机组振动信号降噪方法及系统技术方案

本发明专利技术公开了一种基于交叉熵的水电机组振动信号降噪方法及系统，涉及信号处理技术领域，包括以下步骤：利用经验模态分解法对原始振动信号进行分解处理，获得若干模态分量；计算若干模态分量与所述原始振动信号的交叉熵并计算阈值，筛选有效分量；利用奇异值分解法对所述有效分量进行二次降噪，得到二次降噪后的分量；将所述二次降噪后的分量进行重构，得到降噪后的振动信号。本发明专利技术结合经验模态分解法、奇异值分解法对振动信号进行降噪处理，具有一定的鲁棒性，可避免信号中有效信息的损失。失。失。

【技术实现步骤摘要】
一种基于交叉熵的水电机组振动信号降噪方法及系统


[0001]本专利技术涉及信号处理
，更具体的说是涉及一种基于交叉熵的水电机组振动信号降噪方法及系统。

技术介绍

[0002]水电机组是水电的心脏，对其进行定期的故障诊断工作是保障水力水电事业健康运转的前提。振动信号之中隐藏了水电机组的大部分故障信息，但是它却经常被淹没在水电机组运行中的众多噪声中，从而使得大部分的信号失真无法获取确切的故障信息。因此在进行后续的故障诊断工作前，对水电机组振动信号进行降噪的预处理工作是必要的步骤。
[0003]因此，如何对振动信号进行很好地降噪处理，简化计算过程，避免有效信息的损失是本领域技术人员亟需解决的技术问题。

技术实现思路

[0004]有鉴于此，本专利技术提供了一种基于交叉熵的水电机组振动信号降噪方法及系统，结合经验模态分解法、奇异值分解法对振动信号进行降噪处理，具有一定的鲁棒性，减少信号中有效信息的损失。
[0005]为了实现上述目的，本专利技术提供如下技术方案：
[0006]一种基于交叉熵的水电机组振动信号降噪方法，包括以下步骤：
[0007]利用经验模态分解法对原始振动信号进行分解处理，获得若干模态分量；
[0008]计算若干模态分量与所述原始振动信号的交叉熵并计算阈值，筛选有效分量；
[0009]利用奇异值分解法对所述有效分量进行二次降噪，得到二次降噪后的分量；
[0010]将所述二次降噪后的分量进行重构，得到降噪后的振动信号。
[0011]上述技术方案达到的技术效果为：结合经验模态分解法、奇异值分解法能够很好地对信号进行降噪处理，且具有一定的鲁棒性，解决降噪过程中有效信息的损失问题。
[0012]可选的，所述获得若干模态分量，具体包括以下步骤：
[0013]确定原始振动信号在分析时段内所有的极大值点、极小值点，利用三次样条线将所有的极大值点、极小值点分别用平滑的线段连接，形成上包络线和下包络线；所述上包络线和下包络线囊括所有的极大值点、极小值点；
[0014]计算所述上包络线和下包络线的均值，得到第一信号h1：
[0015]h1＝x(t)
‑
m1ꢀꢀꢀ
(1)；
[0016]判断所述第一信号h1是否满足预设条件，若满足则将第一信号h1作为第一个本征模态函数，记为c1(t)；若不满足则重复以上操作直至满足预设条件；
[0017]将c1(t)从原始振动信号中分离，得到第二信号r1：
[0018]r1＝x(t)
‑
c1ꢀꢀꢀ
(2)；
[0019]不断重复以上操作，得到若干模态分量：
[0020][0021]直至得到一个单调函数，无法从信号中提取本征模态函数时，结束循环：
[0022][0023]式中，r
n
(t)为残余函数，表示信号的平稳趋势；x(t)为原始振动信号；m1表示上、下包络线的均值。
[0024]可选的，所述预设条件为：
[0025]在所述分析时段内，极值点的个数和过零点的个数相等或最多相差1；
[0026]在任意时刻，由所述分析时段内的极大值点形成的上包络线和极小值点形成的下包络线之间的平均值为零，即所述分析时段的信号基于时间轴局部上下对称。
[0027]可选的，所述筛选有效分量，具体为：
[0028]若随机变量有两个独立的概率分布P(x)与Q(x)，则交叉熵表示为：
[0029][0030]当两个概率分布均服从高斯分布时，交叉熵可表示为：
[0031][0032]式中，μ1、μ2、σ1、σ2分别表示两个概率密度分布函数的均值与标准差；
[0033]假设处理信号中的有用信号与背景噪声的子带对数能量分布均符合高斯分布，将有用信号的对数能量概率分布记为P
S
，背景噪声信号的对数能量概率分布记为P
N
，均值分别为μ
S
、μ
N
，标准差分别为σ
S
、σ
N
，规定测度表达式为：
[0034][0035]阈值表达式为：
[0036][0037]式中，E(
·
)表示均值；ρ
S,N
表示最终的对称交叉熵测度，n表示模态分量的个数，H
i
表示第i个模态分量与原信号的相关系数；
[0038]选取高于阈值的分量作为可疑分量，进行二次降噪。
[0039]可选的，所述得到二次降噪后的分量，具体包括以下步骤：
[0040]将长度为N的时间序列，重构为m
×
n阶的Hankle矩阵，如下：
[0041][0042]根据正交化方法有：
[0043]K＝UDV
T
ꢀꢀꢀ
(10)；
[0044]式中，N＝m+n+1，m、n表示矩阵的阶数，U∈R
m
×
m
，V∈R
n
×
n
，D∈R
m
×
n
，D＝(diag(σ1,σ2,...,σ
q
),0)，0代表零矩阵；q＝min(m,n)；σ1≥σ2≥...≥σ
q
＞0代表矩阵的奇异值；
[0045]基于奇异值差分谱，确定最大突变点作为有效信息与噪声信息的分界点，选择分界点前的部分进行叠加重构，得到二次降噪后的分量。
[0046]上述技术方案达到的技术效果为：利用奇异值分解法(SVD)进行二次降噪，减少信号有效信息的损失。
[0047]本专利技术还公开了一种基于交叉熵的水电机组振动信号降噪系统，包括：第一分解模块、筛选模块、第二分解模块、重构模块，且各结构依次相连；
[0048]所述第一分解模块，利用经验模态分解法对原始振动信号进行分解处理，获得若干模态分量；
[0049]所述筛选模块，用于计算若干模态分量与所述原始振动信号的交叉熵并计算阈值，筛选有效分量；
[0050]所述第二分解模块，利用奇异值分解法对所述有效分量进行二次降噪，得到二次降噪后的分量；
[0051]所述重构模块，用于将所述二次降噪后的分量进行重构，得到降噪后的振动信号。
[0052]经由上述的技术方案可知，与现有技术相比，本专利技术公开提供了一种基于交叉熵的水电机组振动信号降噪方法及系统，基于经验模态分解法与交叉熵对振动信号进行分解、选取有效分量，并利用奇异值分解法进行二次降噪，减少信号有效信息的损失；引入交叉熵概念并对计算过程进行简化，将其作为有效分量筛选的判定标准。本技术方案结合了经验模态分解法、奇异值分解法，能够很好地对振动信号进行降噪处理，且具有一定的鲁棒性，解决降噪过程中有效信息的损失问题。
附图说明
[0053]为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于交叉熵的水电机组振动信号降噪方法，其特征在于，包括以下步骤：利用经验模态分解法对原始振动信号进行分解处理，获得若干模态分量；计算若干模态分量与所述原始振动信号的交叉熵并计算阈值，筛选有效分量；利用奇异值分解法对所述有效分量进行二次降噪，得到二次降噪后的分量；将所述二次降噪后的分量进行重构，得到降噪后的振动信号。2.根据权利要求1所述的一种基于交叉熵的水电机组振动信号降噪方法，其特征在于，所述获得若干模态分量，具体包括以下步骤：确定原始振动信号在分析时段内所有的极大值点、极小值点，利用三次样条线将所有的极大值点、极小值点分别用平滑的线段连接，形成上包络线和下包络线；所述上包络线和下包络线囊括所有的极大值点、极小值点；计算所述上包络线和下包络线的均值，得到第一信号h1：h1＝x(t)
‑
m1ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)；判断所述第一信号h1是否满足预设条件，若满足则将第一信号h1作为第一个本征模态函数，记为c1(t)；若不满足则重复以上操作直至满足预设条件；将c1(t)从原始振动信号中分离，得到第二信号r1：r1＝x(t)
‑
c1ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)；不断重复以上操作，得到若干模态分量：直至得到一个单调函数，无法从信号中提取本征模态函数时，结束循环：式中，r
n
(t)为残余函数，表示信号的平稳趋势；x(t)为原始振动信号；m1表示上、下包络线的均值。3.根据权利要求2所述的一种基于交叉熵的水电机组振动信号降噪方法，其特征在于，所述预设条件为：在所述分析时段内，极值点的个数和过零点的个数相等或最多相差1；在任意时刻，由所述分析时段内的极大值点形成的上包络线和极小值点形成的下包络线之间的平均值为零，即所述分析时段的信号基于时间轴局部上下对称。4.根据权利要求1所述的一种基于交叉熵的水电机组振动信号降噪方法，其特征在于，所述筛选有效分量，具体为：若随机变量有两个独立的概率分布P(x)与Q(x)，则交叉熵表示为：
当两个概率分布均服从高斯分布时，交叉熵可表示为：式中，μ1、μ2、σ1、σ2分别表示两个概率密度分布函数的均值与标准差；假设处理信号中的有用信号与背...

【专利技术属性】
技术研发人员：王斌，张文静，吴凤娇，陈帝伊，任晟民，徐哲熙，
申请(专利权)人：西北农林科技大学，
类型：发明
国别省市：