基于灰狼优化算法的故障检测方法、系统、装置及介质与流程

1.本发明涉及故障检测技术领域，尤其是一种基于灰狼优化算法的故障检测方法、系统、装置及介质。


背景技术：

2.现如今，对机械设备的故障诊断分析大多基于振动信号，而振动信号具有非线性、非平稳性等特点，利用它可获取充分表达信号特征的信息。现有技术中，对于振动信号的处理，一般采用傅里叶变换等进行时频转化方法，例如基于daubechies小波变换的滚动轴承故障特征提取方法。这种方法在处理过程中会出现混叠现象从而产生误差，且该类误差由傅里叶算法原理产生，不可避免。由于对振动信号的特征提取存在较大误差，因此，现有的关于振动信号的故障检测方法往往并不准确。此外，对于特征提取后的图像化结果，单纯采用神经网络进行特征识别时，依赖于较多的网络层和神经元才能完整囊括所有信息，并且也无法判定各神经元之间与各网络层之间的耦合关系，同时还可能存在过拟合和欠拟合问题，同样影响了故障检测的准确度。


技术实现要素：

3.本发明的目的在于至少一定程度上解决现有技术中存在的技术问题之一。
4.为此，本发明实施例的一个目的在于提供一种基于灰狼优化算法的故障检测方法，该方法一方面通过对振动信号进行分解，获取重复特征波形和重复时间序列两个不同尺度的信息，避免了傅里叶算法的混叠现象产生的误差，提高了信号特征提取的准确性；另一方面通过灰狼优化算法对时频特征图进行优化，从整体搜索的视角出发，能够更好地完成故障特征信息的提取，从而提高了设备故障检测的准确度。
5.本发明实施例的另一个目的在于提供一种基于灰狼优化算法的故障检测系统。
6.为了达到上述技术目的，本发明实施例所采取的技术方案包括：
7.第一方面，本发明实施例提供了一种基于灰狼优化算法的故障检测方法，包括以下步骤：
8.获取待检测设备的第一振动信号，并通过后缀树算法对所述第一振动信号进行分解，得到第一重复特征波形和第一重复时间序列；
9.根据所述第一重复特征波形和所述第一重复时间序列确定所述第一振动信号的第一时频特征图；
10.通过灰狼优化算法对所述第一时频特征图进行优化处理，提取出所述第一时频特征图中的故障特征信息；
11.根据所述故障特征信息确定待检测设备的故障类型。
12.进一步地，在本发明的一个实施例中，所述通过后缀树算法对所述第一振动信号进行分解，得到第一重复特征波形和第一重复时间序列这一步骤，其具体包括：
13.通过平均分布或高斯分布对所述第一振动信号进行编码得到第一时域信号；
14.通过后缀树算法对所述第一时域信号进行分解，得到多个故障波形信息和对应的时刻信息，并根据所述故障波形信息和所述时刻信息构建第一后缀树；
15.对所述第一后缀树的各个节点进行遍历，获取重复出现的故障波形信息作为第一重复特征波形，并确定所述第一重复特征波形的第一重复时间序列。
16.进一步地，在本发明的一个实施例中，所述对所述第一后缀树的各个节点进行遍历，获取重复出现的故障波形信息作为第一重复特征波形，并确定所述第一重复特征波形的第一重复时间序列这一步骤，其具体包括：
17.从所述第一后缀树的根节点开始，通过深度优先的嵌套遍历算法对所述第一后缀树的各个节点进行遍历；
18.获取重复出现的故障波形信息作为第一重复特征波形，并确定所述第一重复特征波形对应的多个时刻信息；
19.根据所述第一重复特征波形对应的多个时刻信息确定所述第一重复特征波形的第一重复时间序列。
20.进一步地，在本发明的一个实施例中，所述根据所述第一重复特征波形和所述第一重复时间序列确定所述第一振动信号的第一时频特征图这一步骤，其具体为：
21.根据预设的重复长度范围，对所述第一重复特征波形和所述第一重复时间序列进行归一化处理，得到所述第一振动信号的第一时频特征图。
22.进一步地，在本发明的一个实施例中，所述通过灰狼优化算法对所述第一时频特征图进行优化处理，提取出所述第一时频特征图中的故障特征信息这一步骤，其具体包括：
23.根据所述第一时频特征图确定初始灰狼种群；
24.计算所述初始灰狼种群中各个灰狼个体的适应度，进而根据所述适应度确定第一灰狼个体、第二灰狼个体以及第三灰狼个体；
25.根据所述第一灰狼个体、所述第二灰狼个体以及所述第三灰狼个体对所述初始灰狼种群进行迭代更新，当达到预设的迭代次数，根据当前各个灰狼个体的位置确定狩猎目标位置；
26.根据所述狩猎目标位置在所述第一时频特征图中确定所述故障特征信息。
27.进一步地，在本发明的一个实施例中，所述根据所述第一灰狼个体、所述第二灰狼个体以及所述第三灰狼个体对所述初始灰狼种群进行迭代更新这一步骤，其具体包括：
28.根据所述第一灰狼个体、所述第二灰狼个体以及所述第三灰狼个体的位置更新各个灰狼个体的位置；
29.根据更新后各个灰狼个体的位置重新计算各个灰狼个体的适应度，并重新确定第一灰狼个体、第二灰狼个体以及第三灰狼个体。
30.进一步地，在本发明的一个实施例中，所述根据所述故障特征信息确定待检测设备的故障类型这一步骤，其具体为：
31.将所述故障特征信息输入到预先训练好的故障特征识别模型中，得到待检测设备的故障类型。
32.第二方面，本发明实施例提出了一种基于灰狼优化算法的故障检测系统，包括：
33.信号分解模块，用于获取待检测设备的第一振动信号，并通过后缀树算法对所述第一振动信号进行分解，得到第一重复特征波形和第一重复时间序列；
34.时频特征图确定模块，用于根据所述第一重复特征波形和所述第一重复时间序列确定所述第一振动信号的第一时频特征图；
35.故障特征信息提取模块，用于通过灰狼优化算法对所述第一时频特征图进行优化处理，提取出所述第一时频特征图中的故障特征信息；
36.故障类型确定模块，用于根据所述故障特征信息确定待检测设备的故障类型。
37.第三方面，本发明实施例提供了一种基于灰狼优化算法的故障检测装置，包括：
38.至少一个处理器；
39.至少一个存储器，用于存储至少一个程序；
40.当所述至少一个程序被所述至少一个处理器执行时，使得所述至少一个处理器实现上述的一种基于灰狼优化算法的故障检测方法。
41.第四方面，本发明实施例还提供了一种计算机可读存储介质，其中存储有处理器可执行的程序，所述处理器可执行的程序在由处理器执行时用于执行上述的一种基于灰狼优化算法的故障检测方法。
42.本发明的优点和有益效果将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到：
43.本发明实施例获取待检测设备的第一振动信号，通过后缀树算法对第一振动信号进行分解得到第一重复特征波形和第一重复时间序列，然后根据第一重复特征波形和第一重复时间序列确定第一时频特征图，进而通过灰狼优化算法提取出第一时频特征图中的故障特征信息，从而可以根据该故障特征信息确定待检测设备的故障类型。本发明实施例一方面通过对振动信号进行分解，获取重复特征波形和重复时间序列两个不同尺度的信息，避免了傅里叶算法的混叠现象产生的误差，提高了信号特征提取的准确性；另一方面通过灰狼优化算法对时频特征图进行优化，从整体搜索的视角出发，能够更好地完成故障特征信息的提取，从而提高了设备故障检测的准确度。
附图说明
44.为了更清楚地说明本发明实施例中的技术方案，下面对本发明实施例中所需要使用的附图作以下介绍，应当理解的是，下面介绍中的附图仅仅为了方便清晰表述本发明的技术方案中的部分实施例，对于本领域的技术人员来说，在无需付出创造性劳动的前提下，还可以根据这些附图获取到其他附图。
45.图1为本发明实施例提供的一种基于灰狼优化算法的故障检测方法的步骤流程图；
46.图2为本发明实施例提供的后缀树算法的分解示意图；
47.图3为本发明实施例提供的叶节点递归调用顺序示意图；
48.图4为本发明实施例提供的时频特征示意图；
49.图5为本发明实施例提供的灰狼优化算法中灰狼个体和猎物位置关系示意图；
50.图6为本发明实施例提供的灰狼优化算法的狩猎过程示意图；
51.图7为本发明实施例提供的一种基于灰狼优化算法的故障检测系统的结构框图；
52.图8为本发明实施例提供的一种基于灰狼优化算法的故障检测装置的结构框图。
具体实施方式
53.下面详细描述本发明的实施例，实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能理解为对本发明的限制。对于以下实施例中的步骤编号，其仅为了便于阐述说明而设置，对步骤之间的顺序不做任何限定，实施例中的各步骤的执行顺序均可根据本领域技术人员的理解来进行适应性调整。
54.在本发明的描述中，多个的含义是两个或两个以上，如果有描述到第一、第二只是用于区分技术特征为目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量或者隐含指明所指示的技术特征的先后关系。此外，除非另有定义，本文所使用的所有的技术和科学术语与本技术领域的技术人员通常理解的含义相同。
55.参照图1，本发明实施例提供了一种基于灰狼优化算法的故障检测方法，具体包括以下步骤：
56.s101、获取待检测设备的第一振动信号，并通过后缀树算法对第一振动信号进行分解，得到第一重复特征波形和第一重复时间序列；
57.具体地，运用传感器采集待检测设备的第一振动信号f(n)，n为采样点数。单次处理的数据量要求至少包含2个关注信号的完整周期，一般应达到关注特征周期的3～5倍以上，以获取更多的故障重复特征。
58.进一步作为可选的实施方式，通过后缀树算法对第一振动信号进行分解，得到第一重复特征波形和第一重复时间序列这一步骤，其具体包括：
59.a1、通过平均分布或高斯分布对第一振动信号进行编码得到第一时域信号；
60.a2、通过后缀树算法对第一时域信号进行分解，得到多个故障波形信息和对应的时刻信息，并根据故障波形信息和时刻信息构建第一后缀树；
61.a3、对第一后缀树的各个节点进行遍历，获取重复出现的故障波形信息作为第一重复特征波形，并确定第一重复特征波形的第一重复时间序列。
62.具体地，按照第一振动信号的幅值区间，依照平均分布或者高斯分布按预设编码位数对第一振动信号进行编码。
63.按照平均分布对第一振动信号进行编码的公式如下：
64.c1(n)＝int[(f(n)
‑
f
1，min
)*l
code
/(f
1，max
‑
f
1，min
)]
[0065]
按照高斯分布对第一振动信号进行编码的公式如下：
[0066]
c1(n)＝int[igd((f(n)
‑
μ1)/σ1)*l
code
]
[0067]
其中，c1(n)为n次采样数据对应的第一码带的编码值，也即第一时域信号，int()为取整函数，l
code
为预设的编码位数，f
1,max
为第一振动信号的值域最大值，也即第一码带编码时的值域最大值，f
1,min
为第一振动信号的值域最小值，也即第一码带编码时的值域最小值；μ1为第一码带高斯分布均值，σ1为第一码带的高斯分布方差，igd为标准正态分布积分概率查表返回函数。
[0068]
同时获得残差信号r1(n)：r1(n)＝f(n)
‑
c1(n)。
[0069]
将编码后的c1(n)进行后缀树重建。后缀树构建过程中遵循三个传递规则：
[0070]
relu1：当向根节点root插入新后缀时使用。active_node保持为root，active_edge被设置为即将被插入的新后缀的首字符，active_length减1。
[0071]
relu2：分裂(split)一条边并且插入(insert)一个新的节点时，如果该新节点不是当前步骤中创建的第一个节点，则将先前插入的节点与该新节点通过一个特殊的指针连接，称为后缀连接(suffix link)，在图示中通常以虚线绘制。
[0072]
relu3：当从active_node不为根节点root的节点分裂边时，沿着后缀连接(suffix link)的方向寻找节点，如果存在一个节点，则设置该节点为active_node；如果不存在，则设置active_node为root。active_edge和active_length保持不变。
[0073]
下面对本发明实施例采用的后缀树算法进行介绍。存在一个原始数据列t＝t1t2...t
n
，其中t
i
(1≤i≤n 1)，n为数据长度，欲从t1到t
n
将原始数据依次分解为n 1个不重复子序列，第n 1个子序列为规定的终止符，用
‘
#’表示。为便于表达，相关符号说明如下：
[0074]
o：(root)根节点，序列起点，无具体意义；
[0075]
p：(acitve_piont)活动点，指定活动起始点；
[0076]
n：(active_node)活动节点，指定某一子节点；
[0077]
e：(active_edge)活动边，指定序列连接方向；
[0078]
l：(active_edge)活动长度，指定序列移动的数据量；
[0079]
r：(remainder)剩余后缀数，指明未连接的后缀的数目；
[0080]
#：终止符；
[0081]
street(t)：最终分解结果。
[0082]
从根节点o开始，从左到右依次分解原始数据序列，直到产生第n 1个序列。如下式：
[0083]
stree(t)＝(f
i
,f
i
,g
i
),i∈[1,n]
[0084]
其中，f
i
表示主边序列，f
i
表示子边序列，g
i
表示数据i的连接方式，包括活动点p和剩余后缀数r的取值。当数据处于t
i
(1≤i≤n)位置时，通过如下传递方式，完成各条边的连接。
[0085]
1)i＝1时
[0086]
p1＝(o,'f1',1)，r＝1。起始位置选定根节点o；活动边e设为
‘
f1’
；活动长度l和剩余后缀数r设置为1，表明只需传入一个数据量。stree(t1)＝(f1,g1)。
[0087]
2)i>1时
[0088]
①
即t
i
为stree(t
i
‑1)后新出现的数据，设置p
i
＝(o,'f1',i)，r＝i。直接连接在全部主边之后，从f1开始迭代更新p
ik
＝(o,'f
k
',i
‑
k 1)，其中k表示主边数目，k∈(0,i]，k＝k 1,r＝r
‑
1，直至i＝k，停止更新p
ik
。在对已有边进行延展后，创建一个新的主序列f
k 1
，以根节点o作为起始点，以t
i
作为其第一条边。stree(t
i
)＝(f
i
,g
i
)。
[0089]
②
t
i
∈t
i
‑1，即t
i
已经在前缀中出现，被视为重复数据。从主链f1的第一个数据t0开始，寻找与t
i
重复的数据出现的位置，取j表示边长l。设置p
i
＝(o,'f1',j)，r＝1。由于未知t
i 1
的数据特点，暂时只能取定活动点p
i
和剩余后缀数r，并不能给出具体的数据序列的方向。由j可断定所有前j个主链中均存在与t
i
重复的数据，因此，j从1开始，迭代更新p
ij
＝(o,'f
j
',i
‑
j)，l∈(0,j]，l＝l
‑
1，r＝1，每次更新中，主边续接t
i
，r值不变，恒为1，只对t
i
这一个数据量进行操作。此步骤可参考规则1，活动点为根节点，活动边设置为新后缀的初始数据，操作一次后，活动边长减1。该过程未创建新边。活动点的更新用于连接后缀数据，不能作为创建新边的标准。
[0090]
③
在第
②
步的前提下，创建子序列f
i
，f
i
均以活动点p
ij
为起点，从主边f
j
中分离，沿用主边f
j
的前缀数据。子序列f
i
中导入数据同主边f
j
一致，将数据t
i
续接在序列后即可。
[0091]
stree(t
i
)＝(f
i
,f
i
,g
i
),i＜n。
[0092]
④
建立后缀连接(suffix link)。重复数据t
i
出现后，每一主边均会产生一个或多个子边节点p
fi
，如果该节点不是当前数据插入过程中创建的第一个节点p
f1
，则将先前插入的节点与该新节点通过一个特殊的指针连接，称为后缀连接。
[0093]
⑤
当从n不为根节点o的节点分裂边时，沿着后缀连接(suffix link)的方向寻找节点，如果存在一个节点，则设置该节点为n；如果不存在，则设置n为o。e和l保持不变。
[0094]
⑥
循环上述
②‑⑤
步，直到i＝n时，完成了所有数据的分解。
[0095]
以
‘
abcabxabcd’字符串为例，完成字符串的后缀树分解，如图2所示，序号为主从边编号，
‘
#’为结束符，空心圆表示父节点，三角形编号表示叶节点，箭头表示搜索顺序。
[0096]
加入终止符后，总共分成11个子串，如下表1所示：
[0097][0098][0099]
表1
[0100]
由表1可以清晰看出，原字符串可分解为若干个不重复子串，子串所能携带的信息量也不尽相同。后缀树分解，仅仅是对原数据的处理，完成了数据编码的再编排过程，以牺牲存储量来换取特征数据段的读取时间，完整保留了所有数据信号。
[0101]
进一步作为可选的实施方式，对第一后缀树的各个节点进行遍历，获取重复出现的故障波形信息作为第一重复特征波形，并确定第一重复特征波形的第一重复时间序列这一步骤，其具体包括：
[0102]
a31、从第一后缀树的根节点开始，通过深度优先的嵌套遍历算法对第一后缀树的各个节点进行遍历；
[0103]
a32、获取重复出现的故障波形信息作为第一重复特征波形，并确定第一重复特征
波形对应的多个时刻信息；
[0104]
a33、根据第一重复特征波形对应的多个时刻信息确定第一重复特征波形的第一重复时间序列。
[0105]
进一步作为可选的实施方式，通过深度优先的嵌套遍历算法对第一后缀树的各个节点进行遍历这一步骤，其具体包括：
[0106]
b1、创建与第一时域信号等长的重复时刻存储数组；
[0107]
b2、创建与第一后缀树的非叶节点数量一致的重复特征记录数组；
[0108]
b3、运行depth_first嵌套函数，输入节点编号和父节点重复波形长度，输出重复时刻起始位置、重复波形长度、波形重复次数、节点重复字符串长度以及多个重复时刻；
[0109]
其中，重复特征记录数组用于存储重复时刻起始位置、重复波形长度、波形重复次数以及节点重复字符串长度，重复时刻存储数组用于存储重复时刻。
[0110]
具体地，算法伪代码如下：
[0111]
depth_first嵌套函数(简写为df)：
[0112]
输入：节点编号nnodeid，父节点重复波形长度nfathernoderepeatlength；
[0113]
输出：重复时刻起始位置，重复波形长度nnoderepeatlength，波形重复次数nwrindex，重复时刻存储数组(全局变量)；
[0114]
1：重复波形长度＝当前节点字符长度 父节点重复波形长度
[0115]
2：记录节点重复时刻起始位置
[0116]
3：按照节点编号访问该节点的第一个子节点
[0117]
4：波形重复次数＝0
[0118]
5：repeat
[0119]
6：if子节点是非叶节点
[0120]
7：递归调用子节点df函数，获取子节点波形重复次数
[0121]
8：波形重复次数＝波形重复次数 子节点波形重复次数
[0122]
9：else//子节点是叶节点
[0123]
10：波形重复次数＝波形重复次数 1
[0124]
11：将叶节点编号填入重复时刻存储数组
[0125]
12：end if
[0126]
13：until所有子节点均已访问
[0127]
14：return重复波形长度
[0128]
深度优先算法：
[0129]
1：创建与待处理信号等长的重复时刻存储数组；
[0130]
2：创建与非叶节点数量一致的重复特征记录数组，包含如下信息：{重复波形字符串终止地址，节点重复字符串长度，重复时刻起始位置，波形重复次数}；
[0131]
3：运行df(根节点编号，0)//父节点重复波形长度为0；
[0132]
4：获得故障的特征波形与对应重复时间序列的特征。
[0133]
时频特征提取算法将后缀树算法中创建的节点仅遍历了一次。由于后缀树算法的复杂度为o(n)，所以时频特征提取算法的复杂度同样是o(n)。
[0134]
下面以abcabxabcd字符串的后缀树结构为例，进一步说明本发明实施例时频特征
提取算法的基本原理。
[0135]
df函数递归调用顺序如图3所示，依次为df(0,0)，df(1,0)，df(3,2)，df(4,1)，df(2,0)，df(5,0)，df(0,0)，每条主边中，访问完叶节点均需回到父节点处，等待下一次节点搜索。
[0136]
根据图2中非叶节点次序完成子序列特征数组，分解结果见表2，表3。
[0137]
1234567891011abcabxabcd#1742853961011
[0138]
表2
[0139]
表2第1行为原字符串
‘
abcabxabcd’各字符编号，第2行为原字符串，第3行为提取后的重复性特征分解序列，其值对应图2中叶节点编号，也对应值所示采样点的时刻。该数组的存储量与原始数据长度一致，其最后一位始终对应字符串终止符，该位可以省略。
[0140]
numw_endrelenstart_rtvreptimesnode_str12213ab22143b33312c43242c53172c
[0141]
表3
[0142]
表3中记录了每一个非叶节点的重复特征数组信息，包括重复特征波形终止点，重复特征波形长度relen，重复时间特征向量起始点start_rtv，重复次数reptimes。为了便于理解，表3中添加了节点编号num与节点字符串node_str。表3中节点编号与节点字符串均与图2中标注一致，均可由后缀树算法直接获取。
[0143]
进一步作为可选的实施方式，故障检测方法还包括以下步骤：
[0144]
获取超过预设的重复次数阈值的重复特征波形的最大重复长度，当最大重复长度大于等于预设的长度阈值且残差编码的次数小于等于预设的编码次数阈值，对残差编码的参数进行更新后再进行残差编码和信号分解。
[0145]
具体地，检验超过预设重复次数波形的最大重复长度，如该长度小于预设值或者残差分析次数i大于预设值则终止分析。如该长度大于等于预设值且残差分析次数i小于等于预设值，则对第i次残差信号r
i
(n)继续执行分解。
[0146]
其中残差编码时参数修正的规则如下：
[0147]
如果依照平均分布进行编码，则采用下式更新上下限：
[0148][0149]
如果依照高斯分布进行编码，则：
[0150]
a)可重新依照平均分布进行残差编码，同样采用上式更新上下限。
[0151]
b)依然依照高斯分布编码，采用下式更新参数：
[0152]
μ
i
＝0
[0153][0154]
再利用下式进行残差编码：
[0155]
c
i
(n)＝int[(igd((r
i
‑1(n)
‑
μ
i
)/σ
i
)
‑
0.5)*2*l
code
]
[0156]
s102、根据第一重复特征波形和第一重复时间序列确定第一振动信号的第一时频特征图。
[0157]
具体地，输出故障特征包含重复特征波形与重复时间序列两个信息。如表3中重复特征波形长度relen最大为3，对应重复特征波形终止点w_end为3，因此该特征的波形为表2第2行前3位所示的“abc”，而对应的重复时间序列特征为表2第3行从重复时间特征向量起始点start_rtv即第1位开始的共重复次数reptimes的2个序列，即“1,7”。因此重复波形长度为3采样点的故障特征为：
[0158]
{“abc”，“1,7”}
[0159]
同理，重复长度为2的特征有：
[0160]
{“ab”，“1,7,4”}、{“bc”，“2,8”}
[0161]
重复特征波形与重复时间序列可单独或共同作为故障特征用于后续故障诊断。结合动力学模型，也可将多个重复时间序列的交、并集及其补集等集合运算结果作为故障特征。其中某时刻最大波形重复长度越长意味着该时刻附近更长时间内的数据结构稳定性越好。当干扰近似白噪声时，该时刻的瞬时频率也相对越低。反之某时刻最大波形重复长度越短，意味该时刻数据结构越不稳定，该时刻的瞬时信号也越接近于突变或冲击信号。
[0162]
进一步作为可选的实时方式，根据第一重复特征波形和第一重复时间序列确定第一振动信号的第一时频特征图，并根据第一时频特征图构建训练样本集这一步骤，其具体为：
[0163]
根据预设的重复长度范围，对第一重复特征波形和第一重复时间序列进行归一化处理，得到第一振动信号的第一时频特征图。
[0164]
具体地，第一时频特征图的横坐标表示时刻信息，第一时频特征图的纵坐标表示重复波形长度，第一时频特征图的像素点的色彩值表示对应时刻的数据在对应重复波形长度的重复特征波形中的参与程度，从而可以实现第一振动信号的故障特征可视化。
[0165]
如图4所示为本发明实施例提供的时频特征示意图。将字符串数据类比于振动信号，图4中颜色深浅代表着不同时刻数据在当前重复长度中的活跃程度。第2行“b”字符的对应的色彩值较亮，表示其处很可能存在频率较低的周期性信号。第1行“x”、“d”字符处数值存在低值，表明该时刻信号极不稳定，可能存在冲击等突变信号。第1行2处“a”字符数值也存在低值，但第二行对应位置值反而更大，表示“a”字符对应信号隶属于后续的“ab”字符，在处理数据中不会脱离“ab”单独出现。
[0166]
后缀树分解后，得到具有时频特征的图像结果。图像中出现的低频凹坑位置，可视为故障件受到冲击的时刻；图像中出现的高频高亮的区域，可视为系统的固有振动信号。对于一个旋转系统而言，旋转周期小，振动频率高，很难直接判断出图像各色块之间的差异(色块表征为信号时频特征)，因此，本发明实施例引入灰狼优化算法提取出时频特征图中的故障特征信息所在的区域，从而可以准确判断具体的故障类型。
[0167]
s103、通过灰狼优化算法对第一时频特征图进行优化处理，提取出第一时频特征图中的故障特征信息。
[0168]
具体地，灰狼优化算法是(grey wolf optimizer，gwo)由澳大利亚格里菲斯大学学者mirjalili等人于2014年提出来的一种群智能优化算法，其通过模拟狼群捕食猎物的狩猎活动来进行相关的特征提取及优化，具有较强的收敛性能、参数少、易实现等特点。步骤s103具体包括以下步骤：
[0169]
s1031、根据第一时频特征图确定初始灰狼种群；
[0170]
s1032、计算初始灰狼种群中各个灰狼个体的适应度，进而根据适应度确定第一灰狼个体、第二灰狼个体以及第三灰狼个体；
[0171]
s1033、根据第一灰狼个体、第二灰狼个体以及第三灰狼个体对初始灰狼种群进行迭代更新，当达到预设的迭代次数，根据当前各个灰狼个体的位置确定狩猎目标位置；
[0172]
s1034、根据狩猎目标位置在第一时频特征图中确定故障特征信息。
[0173]
进一步作为可选的实施方式，根据第一灰狼个体、第二灰狼个体以及第三灰狼个体对初始灰狼种群进行迭代更新这一步骤，其具体包括：
[0174]
c1、根据第一灰狼个体、第二灰狼个体以及第三灰狼个体的位置更新各个灰狼个体的位置；
[0175]
c2、根据更新后各个灰狼个体的位置重新计算各个灰狼个体的适应度，并重新确定第一灰狼个体、第二灰狼个体以及第三灰狼个体。
[0176]
具体地，gwo的优化过程是模拟狼群狩猎过程，拟定狼群各层阶级成员，并协调其具体分工，优化过程包含了灰狼的社会等级分层、跟踪、包围和攻击猎物等过程，具体如下：
[0177]
构建gwo优化模型时，首先需构建灰狼社会等级层次模型。计算种群每个个体的适应度，将狼群中适应度最好的三匹灰狼依次标记为(即第一灰狼个体)、β(即第二灰狼个体)、δ(即第三灰狼个体)，而剩下的灰狼标记为ω。也就是说，灰狼群体中的社会等级从高往低排列依次为：β、δ、ω。gwo的优化过程主要由每代种群中的最优的三个解(即β、δ)来指导完成。
[0178]
灰狼个体捜索猎物时会逐渐地接近猎物并包围它，该行为的数学模型如下：
[0179][0180][0181]
其中，t为当前迭代次数，向量乘积采用hadamard乘积操作；和是协同系数向量；xp表示猎物(即狩猎目标)的位置向量；x(t)表示当前灰狼的位置向量；协同系数向量和的计算公式如下：
[0182][0183][0184]
在整个迭代过程中a由2线性降到0；r1和r2是[0，1]中的随机向量。
[0185]
和向量的改变会在二维坐标中决定灰狼与猎物的位置关系，如图5所示，灰狼的初始位置位于(x,y)，猎物位于可能点的中心位置。通过决策者调整和向量改变灰狼的位置来更新其与猎物的相对位置关系。例如设定和来更改灰狼位置点为(x*
‑
x,y*)。
[0186]
灰狼具有识别潜在猎物(即狩猎目标，也即最优解)位置的能力，搜索过程主要依靠β、δ灰狼的指引来完成。但是很多问题的解空间特征是未知的，灰狼是无法确定猎物(最优解)的精确位置。为了模拟灰狼的搜索行为，假设β、δ具有较强识别潜在猎物位置的能力。因此，在每次迭代过程中，保留当前种群中的最好三只灰狼(β、δ)，然后根据它们的位置信息来更新其它灰狼个体(包括ω)的位置。该行为的数学模型可表示如下：
[0187][0188][0189][0190]
式中：x
β
,x
δ
分别表示当前种群中α,β,δ的位置向量；x表示灰狼的位置向量；d
β
,d
δ
分别表示当前候选灰狼与最优三条狼之间的距离；当|a|>1时，灰狼之间尽量分散在各区域并搜寻猎物。当|a|＜1时，灰狼将集中捜索某个或某些区域的猎物。
[0191]
如图6所示为灰狼优化算法的狩猎过程示意图，从图6中可看出，候选解的位置最终落在被α,β,δ定义的随机圆位置内。总的来说，α,β,δ需首先预测出猎物(潜在最优解)的大致位置，然后其它候选狼在当前最优的那只狼的指引下在猎物附近随机地更新它们的位置。
[0192]
构建攻击猎物模型的过程中，根据包围猎物的数学模型可知，a值的减少会引起a的值也随之波动。换句话说，a是一个在区间[
‑
a，a]上的随机向量，其中a在迭代过程中呈线性下降。当a在[
‑
1，1]区间上时，则灰狼个体(search agent)的下一时刻位置可以在当前灰狼与猎物之间的任何位置上。
[0193]
灰狼主要依赖α,β,δ的信息来寻找猎物。它们开始分散地去搜索猎物位置信息，然后集中起来攻击猎物。对于分散模型的建立，通过|a|>1使其捜索代理远离猎物，这种搜索方式使gwo能进行全局搜索。gwo算法中的另一个搜索系数是c。根据包围猎物的数学模型可知，c向量是在区间范围[0,2]上的随机值构成的向量，此系数为猎物提供了随机权重，以便増加(|c|>1)或减少(|c|＜1)。这有助于gwo在优化过程中展示出随机搜索行为，以避免算法陷入局部最优。值得注意的是，c并不是线性下降的，c在迭代过程中是随机值，该系数有利于算法跳出局部，特别是算法在迭代的后期显得尤为重要。
[0194]
总的来说，在狼群搜索过程中，首先给定灰狼一个随机位置，随着迭代反馈，α,β,δ监测最优解可能出现的位置，每一个监测者更新它们与最优解之间的距离，通过参数a从2减为0的过程中不断迭代优化，同时，监测者判断狼群的行动是远离猎物还是靠近猎物，最后，灰狼优化算法输出最终的狩猎目标位置。根据该狩猎目标位置即可在第一时频特征图
中确定特定故障特征的区域，从而可以提取出故障特征信号。
[0195]
s104、根据故障特征信息确定待检测设备的故障类型。
[0196]
进一步作为可选的实施方式，根据故障特征信息确定待检测设备的故障类型这一步骤，其具体为：
[0197]
将故障特征信息输入到预先训练好的故障特征识别模型中，得到待检测设备的故障类型。
[0198]
具体地，从第一时频特征图中提取出故障特征信息后，可直接将该故障特征信息输入到预先训练好的故障特征识别模型中，实现故障类型的识别。故障特征识别模型可以利用神经网络对已知故障类型的振动信号的故障特征信息进行训练得到，训练过程非本发明重点，在此不做赘述。
[0199]
可以理解的是，由于本发明实施例通过灰狼优化算法从第一时频特征图中提取出了故障特征信息，相比直接对第一时频特征图进行模型识别而言，能够降低识别模型的复杂程度和训练时长，提高故障识别的准确度。
[0200]
以上对本发明的方法步骤进行了说明。可以认识到，本发明实施例一方面通过对振动信号进行分解，获取重复特征波形和重复时间序列两个不同尺度的信息，避免了傅里叶算法的混叠现象产生的误差，提高了信号特征提取的准确性；另一方面通过灰狼优化算法对时频特征图进行优化，从整体搜索的视角出发，能够更好地完成故障特征信息的提取，从而提高了设备故障检测的准确度。
[0201]
与现有技术相比，本发明实施例还具有以下优点：
[0202]
1)后缀树算法提取故障信号具有实时性。适度增加信号的存储量来达到快速提取重复故障信号的目的，在硬件条件满足情况下，可达到毫秒级。
[0203]
2)信号的完整性。对采集的时域信号进行特征提取的整个过程均无信号的损失，保留了所关注信号的所有信息。
[0204]
3)自动化程度高。信号采集、处理、特征提取、类别识别等过程可在线依次完成，完成降低人工劳动强度和劳动成本的目标。
[0205]
4)通用性强。以整体搜寻的思想判断各个决策者和各个特征间的相互关系，更大程度保留了故障特征之间的耦合关系，故障类别的整体特征可迁移至别的机械设备或测试系统之中，仍有较好的作用效果。
[0206]
5)与复杂智能识别方法相比，基于灰狼优化算法的故障识别方法相对简易，参数值少，训练速度快，针对性强。
[0207]
参照图7，本发明实施例提供了一种基于灰狼优化算法的故障检测系统，包括：
[0208]
信号分解模块，用于获取待检测设备的第一振动信号，并通过后缀树算法对第一振动信号进行分解，得到第一重复特征波形和第一重复时间序列；
[0209]
时频特征图确定模块，用于根据第一重复特征波形和第一重复时间序列确定第一振动信号的第一时频特征图；
[0210]
故障特征信息提取模块，用于通过灰狼优化算法对第一时频特征图进行优化处理，提取出第一时频特征图中的故障特征信息；
[0211]
故障类型确定模块，用于根据故障特征信息确定待检测设备的故障类型。
[0212]
上述方法实施例中的内容均适用于本系统实施例中，本系统实施例所具体实现的
功能与上述方法实施例相同，并且达到的有益效果与上述方法实施例所达到的有益效果也相同。
[0213]
参照图8，本发明实施例提供了一种基于灰狼优化算法的故障检测装置，包括：
[0214]
至少一个处理器；
[0215]
至少一个存储器，用于存储至少一个程序；
[0216]
当上述至少一个程序被上述至少一个处理器执行时，使得上述至少一个处理器实现上述的一种基于灰狼优化算法的故障检测方法。
[0217]
上述方法实施例中的内容均适用于本装置实施例中，本装置实施例所具体实现的功能与上述方法实施例相同，并且达到的有益效果与上述方法实施例所达到的有益效果也相同。
[0218]
本发明实施例还提供了一种计算机可读存储介质，其中存储有处理器可执行的程序，该处理器可执行的程序在由处理器执行时用于执行上述一种基于灰狼优化算法的故障检测方法。
[0219]
本发明实施例的一种计算机可读存储介质，可执行本发明方法实施例所提供的一种基于灰狼优化算法的故障检测方法，可执行方法实施例的任意组合实施步骤，具备该方法相应的功能和有益效果。
[0220]
本发明实施例还公开了一种计算机程序产品或计算机程序，该计算机程序产品或计算机程序包括计算机指令，该计算机指令存储在计算机可读存介质中。计算机设备的处理器可以从计算机可读存储介质读取该计算机指令，处理器执行该计算机指令，使得该计算机设备执行图1所示的方法。
[0221]
在一些可选择的实施例中，在方框图中提到的功能/操作可以不按照操作示图提到的顺序发生。例如，取决于所涉及的功能/操作，连续示出的两个方框实际上可以被大体上同时地执行或上述方框有时能以相反顺序被执行。此外，在本发明的流程图中所呈现和描述的实施例以示例的方式被提供，目的在于提供对技术更全面的理解。所公开的方法不限于本文所呈现的操作和逻辑流程。可选择的实施例是可预期的，其中各种操作的顺序被改变以及其中被描述为较大操作的一部分的子操作被独立地执行。
[0222]
此外，虽然在功能性模块的背景下描述了本发明，但应当理解的是，除非另有相反说明，上述的功能和/或特征中的一个或多个可以被集成在单个物理装置和/或软件模块中，或者一个或多个功能和/或特征可以在单独的物理装置或软件模块中被实现。还可以理解的是，有关每个模块的实际实现的详细讨论对于理解本发明是不必要的。更确切地说，考虑到在本文中公开的装置中各种功能模块的属性、功能和内部关系的情况下，在工程师的常规技术内将会了解该模块的实际实现。因此，本领域技术人员运用普通技术就能够在无需过度试验的情况下实现在权利要求书中所阐明的本发明。还可以理解的是，所公开的特定概念仅仅是说明性的，并不意在限制本发明的范围，本发明的范围由所附权利要求书及其等同方案的全部范围来决定。
[0223]
上述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个
人计算机，服务器，或者网络设备等)执行本发明各个实施例上述方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器(rom，read
‑
only memory)、随机存取存储器(r am，random access memory)、磁碟或者光盘等各种可以存储程序代码的介质。
[0224]
在流程图中表示或在此以其他方式描述的逻辑和/或步骤，例如，可以被认为是用于实现逻辑功能的可执行指令的定序列表，可以具体实现在任何计算机可读介质中，以供指令执行系统、装置或设备(如基于计算机的系统、包括处理器的系统或其他可以从指令执行系统、装置或设备取指令并执行指令的系统)使用，或结合这些指令执行系统、装置或设备而使用。就本说明书而言，“计算机可读介质”可以是任何可以包含、存储、通信、传播或传输程序以供指令执行系统、装置或设备或结合这些指令执行系统、装置或设备而使用的装置。
[0225]
计算机可读介质的更具体的示例(非穷尽性列表)包括以下：具有一个或多个布线的电连接部(电子装置)，便携式计算机盘盒(磁装置)，随机存取存储器(ram)，只读存储器(rom)，可擦除可编辑只读存储器(eprom或闪速存储器)，光纤装置，以及便携式光盘只读存储器(cdrom)。另外，计算机可读介质甚至可以是可在其上打印上述程序的纸或其他合适的介质，因为可以例如通过对纸或其他介质进行光学扫描，接着进行编辑、解译或必要时以其他合适方式进行处理来以电子方式获得上述程序，然后将其存储在计算机存储器中。
[0226]
应当理解，本发明的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实施方式中，多个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或固件来实现。例如，如果用硬件来实现，和在另一实施方式中一样，可用本领域公知的下列技术中的任一项或他们的组合来实现：具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻辑电路，具有合适的组合逻辑门电路的专用集成电路，可编程门阵列(pga)，现场可编程门阵列(fpga)等。
[0227]
在本说明书的上述描述中，参考术语“一个实施方式/实施例”、“另一实施方式/实施例”或“某些实施方式/实施例”等的描述意指结合实施方式或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施方式或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施方式或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施方式或示例中以合适的方式结合。
[0228]
尽管已经示出和描述了本发明的实施方式，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施方式进行多种变化、修改、替换和变型，本发明的范围由权利要求及其等同物限定。
[0229]
以上是对本发明的较佳实施进行了具体说明，但本发明并不限于上述实施例，熟悉本领域的技术人员在不违背本发明精神的前提下还可做作出种种的等同变形或替换，这些等同的变形或替换均包含在本技术权利要求所限定的范围内。