一种基于奇异值分解的轴承故障信号降噪方法

1.本发明属于机械故障诊断领域，特别涉及了一种轴承故障信号降噪方法。


背景技术：

2.轴承是机械使用中最广泛且最易损坏的机械零件之一，其状态检测与故障诊断是机械设备故障诊断技术的重要内容。常用的轴承故障诊断方法有振动、噪声分析法、油液分析法、温度测定法、声发射分析法、轴承间隙测定法、油膜电阻测定法、光学纤维测定法、奇异值分解法等。
3.近年来，奇异值分解(singular value decomposition，svd)方法在轴承的故障诊断中表现出独特的优势。许多基于svd的轴承故障诊断方法也在实际应用中表现出良好的诊断特性，例如2005年杨宇、于德介等人结合经验模态分解(empirical mode decomposition，emd)和奇异值分解实现了滚动轴承故障诊断；2014年王超、孔凡让等人提出改进的奇异值分解算法，丰富了奇异值分解在轴承故障诊断领域的应用。而针对求解奇异值的奇异向量矩阵的应用并未展开过多的探索，使得现阶段的故障诊断方法中定位奇异值的突变位置存在不足。


技术实现要素：

4.为了解决上述背景技术提到的技术问题，本发明提出了一种基于奇异值分解的轴承故障信号降噪方法。
5.为了实现上述技术目的，本发明的技术方案为：
6.一种基于奇异值分解的轴承故障信号降噪方法，包括以下步骤：
7.(1)对轴承振动信号进行采样；
8.(2)对采样信号构造hankel矩阵a：
[0009][0010]
其中，x
a
为第a个采样点，a＝1,2,
…
,n，n为采样点数，1＜n＜n；
[0011]
(3)对矩阵a进行奇异值分解a＝u
t
σv，其中，σ为奇异值矩阵，u和v为奇异向量矩阵，t表示转置，u＝(u
ij
),1≤i,j≤m，u
ij
为矩阵u的第i行第j列的元素，m＝n
‑
n+1，v＝(v
ls
),1≤l,s≤n，v
ls
为矩阵v的第l行第s列的元素；且有：
[0012]
[0013][0014]
其中，表示矩阵a删除第j行后的矩阵，表示矩阵a删除第s列后的矩阵，σ
i
(a)、σ
k
(a)、σ
l
(a)、σ
t
(a)为矩阵a的第i个、第k个、第l个、第t个奇异值，为矩阵的第k个奇异值，为矩阵的第t个奇异值；
[0015]
则将a＝u
t
σv改写为
[0016][0017]
其中，u
c
和v
c
为矩阵a的第c个奇异向量，σ
c
为矩阵a的第c个奇异值，c＝1,2,
…
,q；
[0018]
(4)令矩阵a经奇异值分解后的第c个子矩阵a
c
为
[0019][0020]
其中，x
c,a
为子矩阵a
c
相应位置处的元素，a＝1,2,
…
,n；
[0021]
令p
c,1
＝[x
c,1 x
c,2
ꢀ…ꢀ
x
c,n
]，q
c,n
＝[x
c,n+1 x
c,n+2
ꢀ…ꢀ
x
c,n
]，则重构后的第c个子矩阵p
c
为
[0022][0023]
(5)将步骤(3)中得到的q个奇异值σ
c
按照从大到小的顺序排序，令
[0024]
b
d
＝σ
d
‑
σ
d+1
,d＝1,2,
…
,q
‑1[0025]
其中，σ
d
为排序后的第d个奇异值；
[0026]
则奇异值突变位置为
[0027]
(6)利用步骤(4)得到的重构后的子矩阵p
c
和步骤(5)得到的奇异值突变位置d，重构出降噪后的信号p：
[0028][0029]
采用上述技术方案带来的有益效果：
[0030]
本发明求解奇异值的奇异向量矩阵更加简单、快捷，该部分的模块化可直接从实验平台作为产品插入其他系统，适用范围广。通过本发明得到的降噪后的信号进行故障诊断的精确度更好。
附图说明
[0031]
图1是本发明的流程图；
[0032]
图2是实施例中的奇异值差分谱图。
具体实施方式
[0033]
以下将结合附图，对本发明的技术方案进行详细说明。
[0034]
本发明设计了一种基于奇异值分解的轴承故障信号降噪方法，如图1所示。以小波信号来进行说明，构造响应函数如下：
[0035][0036][0037]
n(t)为高斯白噪声，ξ为粘滞阻尼比，k为冲击响应次数，f0和f
n
为载波频率，τ
n
和τ0为时间参数，t0为冲击响应周期，在实际信号中表现为故障频率。此处取ξ＝0.37，k＝5，f
n
＝0.385hz，f0＝200hz，τ
n
＝0.5s，τ0＝1s，t0＝0.05s，以采样频率f
s
＝8000hz对信号采样2048点，进而构造故障信号。
[0038]
对于采样点信号，构造1024
×
1025的hankel矩阵：
[0039][0040]
对矩阵a进行奇异值分解a＝u
t
σv，其分解出的两个奇异向量矩阵u＝(u
ij
),1≤i,j≤1024和v＝(v
ls
),1≤l,s≤1025的元素如下：
[0041][0042][0043]
其中，表示矩阵a删除第j行后的矩阵，表示矩阵a删除第s列后的矩阵，将将a＝u
t
σv改写为
[0044][0045]
其中，u
c
和v
c
为矩阵a的第c个奇异向量，σ
c
为矩阵a的第c个奇异值，c＝1,2,
…
,q，q＝1024。
[0046]
令矩阵a经奇异值分解后的第c个子矩阵a
c
为
[0047][0048]
令p
c,1
＝[x
c,1 x
c,2
ꢀ…ꢀ
x
c,1024
]，q
c,1024
＝[x
c,1025 x
c,1026
ꢀ…ꢀ
x
c,2048
]，则重构后的第c个子矩阵
[0049]
将q个奇异值σ
c
按照从大到小的顺序排序，令b
d
＝σ
d
‑
σ
d+1
,d＝1,2,
…
,q
‑
1则奇异值突变位置为然后重构出降噪后的信号由于绘图比例的关系，只显示前100个较重要的奇异值，且奇异值大小为原值的并给出相应的差分谱，如图2所示，从图中可知应取前2个信号进行重构。
[0050]
对于本领域技术人员而言，显然本发明不限于上述示范性实施例的细节，而且在不背离本发明的精神或基本特征的情况下，能够以其他的具体形式实现本发明。因此，无论从哪一点来看，均应将实施例看作是示范性的，而且是非限制性的，本发明的范围由所附权利要求而不是上述说明限定，因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。
[0051]
此外，应当理解，虽然本说明书按照实施方式加以描述，但并非每个实施方式仅包含一个独立的技术方案，说明书的这种叙述方式仅仅是为清晰起见，本领域技术人员应当将说明书作为一个整体，各实施例中的技术方案也可以经适当组合，形成本领域技术人员可以理解的其他实施方式。