用于机械故障诊断的三阶张量的奇异值分解方法及系统技术方案

本发明专利技术涉及故障诊断技术领域，公开了一种用于机械故障诊断的三阶张量的奇异值分解方法及系统，该方法包括以下步骤：1）以设定的采用频率和采样点数从多个通道对机械工作时所产生的信号进行采集；2）将待诊断的多个通道所采集得到的信号构建为一个三阶张量信号；3）对所构建的三阶张量信号采用张量奇异值分解公式进行自适应分解以获得分量信号；4）对比所获得的分量信号的频率与机械故障特征频率是否一致，以诊断机械是否存在故障。本发明专利技术能够为机械设备故障诊断和状态监测提供基础支撑，实现多通道信号的有效分解和故障特征识别，避免机械设备漏诊问题，从而防止因设备故障引起重大事故发生。大事故发生。大事故发生。

【技术实现步骤摘要】
用于机械故障诊断的三阶张量的奇异值分解方法及系统


[0001]本专利技术涉及故障诊断领域，尤其是涉及一种用于机械故障诊断的三阶张量的奇异值分解方法及系统。

技术介绍

[0002]随着科技的快速发展，自动化程度日益提高，机械设备在工业生产中产生的作用也逐渐增大。机械在运转过程中产生的任何故障都有可能导致人员伤亡和巨大的经济损失。因此，对其运行状态的监测和故障诊断具有重要的意义。机械的故障振动信号是非平稳信号且含有大量噪声等干扰信号，这给基于信号处理的故障诊断技术带来极大的困难，并且工程实际中的信号就是此类信号，所以研究信号处理方法具有很大的工程应用价值。
[0003]近年来，研究学者们将矩阵奇异值分解算法应用于机械故障诊断领域并取得了一定的进展，但是现有的矩阵奇异值分解方法存在着一定的不足，其中最主要的体现便是它无法处理多通道信号，无法提取多通道信号之间的固有耦合关系。而在工程实际中，往往测取的是多通道信号，这就难以满足工程的实际需求。

技术实现思路

[0004]本专利技术一方面所要解决的技术问题是提供了一种用于机械故障诊断的三阶张量的奇异值分解方法，能够为机械设备故障诊断和状态监测提供基础支撑，实现多通道信号的有效分解和故障特征识别，避免机械设备漏诊问题，从而防止因设备故障引起重大事故发生。
[0005]为了解决上述技术问题，本专利技术提供了一种用于机械故障诊断的三阶张量的奇异值分解方法，该方法包括以下步骤：1）以设定的采用频率和采样点数从多个通道对机械工作时所产生的信号进行采集；2）将待诊断的多个通道所采集得到的信号构建为一个三阶张量信号；3）对所构建的三阶张量信号采用张量奇异值分解公式进行自适应分解以获得分量信号；所述张量奇异值分解公式包括以下步骤：31）构建任意阶张量与任意阶张量的o模式积乘法规则；32）根据任意阶张量与任意阶张量的o模式积乘法规则，得出三阶张量奇异值分解公式；33）根据三阶张量奇异值分解公式确定分解后张量的阶数和维数；4）对比所获得的分量信号的频率与机械故障特征频率是否一致，以诊断机械是否存在故障。
[0006]进一步优选地，步骤31）中，首先设定m阶张量中有o阶的维数与n阶张量中o阶的维数对应相等，它们之间的乘法定义为o模式积，则m阶
张量与n阶张量之间的o模式积为：
[0007]其中，中的p和q分别表示张量和张量中参与张量o模式积的阶数顺序为正序还是逆序，若为
‘1’
则表示正序，若为
‘2’
则表示逆序，且o=1, 2,
ꢀ…
, min(m, n)－1；I1、I2、
…
、I
m
分别表示张量的第一阶维数、第二阶维数、
…
第m阶维数；同理，J1、J2、
…
、J
n
分别表示张量的第一阶维数、第二阶维数、
…
第n阶维数；k1、 k2、
…
、k
m+n
‑
2o
分别表示o模式积后结果张量的第一阶维数、第二阶维数、
…
第m+n
‑
2o阶维数；为张量的第二阶至第m阶的元素形式表达式；为张量的第一阶至第n
‑
1阶的元素形式表达式；s1表示参与o模式积运算的第1个阶，S1为相应阶的维数；同理，s
o
表示参与o模式积运算的第o个阶，S
o
为相应阶的维数。
[0008]优选地，所述m阶张量与与n阶张量之间的o模式积的阶数为正序或逆序时的参数计算公式如下所示：
[0009]其中，和表示张量的第二阶至第m阶参数，如中的s1对应于张量的第二阶i2，以此类推，对应于张量的第m阶i
m
；表示的值等于张量的第o+1阶维数，表示的值等于张量的第m
‑
o+1阶维数；和表示张量的第一阶至第n阶参数；表示的值等于张量的第o阶维数，表示的值等于张量的第n
‑
o阶维数。
[0010]进一步优选地，步骤31）中，张量的乘法规则至少包括以下四种：、、和。
[0011]其中，表示张量和参与o模式积运算的阶数顺序都为正序；表示张量和参与o模式积运算的阶数顺序分别为正序和逆序；表示张量和参与o模式积运算的阶数顺序分别
为逆序和正序；表示张量和参与o模式积运算的阶数顺序都为逆序。
[0012]优选地，所述三阶张量的奇异值分解公式将三阶张量分解成、和三个张量，具体表达式如下：
[0013]进一步优选地，根据所得到三阶张量的三阶维数I1、I2和I3以及选定的o、p、q参数，进一步确定、和三个张量的阶数及各个阶的维数，具体表达式如下：
[0014]其中，表示张量为二阶张量，张量的第一阶维数为I1，第二阶维数为I2；表示张量为四阶张量，第一阶维数为I1，第二阶维数为I2，第三阶维数为I2，第四阶维数为I1；表示张量为三阶张量，且第一阶维数为I1，第二阶维数为I2，第三阶维数为I3；表示第一阶维数为I2，第二阶维数为I1，第三阶维数为I3；表示四阶张量的第一阶维数为I2，第二阶维数为I3，第三阶维数为I2，第四阶维数为I3；
表示四阶张量的第一阶维数为I3，第二阶维数为I2，第三阶维数为I2，第四阶维数为I3；表示二阶张量的第一阶维数和第二阶维数都为I3。
[0015]本专利技术第二方面公开了一种张量数据分解系统，所述系统包括：数据获取模块，用于获取机械运作时所产生的张量信号；张量分解模块，用于对所述张量信号进行处理，通过构建任意阶张量与任意阶张量的o模式积乘法规则；根据任意阶张量与任意阶张量的o模式积乘法规则，得出三阶张量奇异值分解公式；根据三阶张量奇异值分解公式确定分解后张量的阶数和维数。
[0016]本专利技术还公开了一种可读存储介质，该可读存储介质上存储有可执行指令，该可执行指令被执行时采用本专利技术第一方面所述的方法步骤。
[0017]本专利技术还公开了一种电子设备，包括存储器、处理器以及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述程序时适于采用本专利技术第一方面所述方法的步骤。
[0018]通过上述技术方案，本专利技术的用于机械故障诊断的三阶张量的奇异值分解方法通过以设定的采用频率和采样点数从多个通道对机械工作时所产生的信号进行采集；将待诊断的多个通道所采集得到的信号构建为一个三阶张量信号；对所构建的三阶张量信号采用张量奇异值分解公式进行自适应分解以获得分量信号；对比所获得的分量信号的频率与机械故障特征频率是否一致，以诊断机械是否存在故障，采用该诊断方法能够处理现有的矩阵奇异值分解方法中所无法处理多通道信号的缺陷，从而能够为机械设备的故障诊断和状态监测提供基础支撑，并且实现了多通道信号的有效分解和故障特征识别，避免机械设备漏诊问题，从而防止因设备故障引起重大事故的发生。
[0019]本专利技术的其他特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
...

【技术保护点】

【技术特征摘要】
1.一种用于机械故障诊断的三阶张量的奇异值分解方法，其特征在于，该方法包括以下步骤：1）以设定的采用频率和采样点数从多个通道对机械工作时所产生的信号进行采集；2）将待诊断的多个通道所采集得到的信号构建为一个三阶张量信号；3）对所构建的三阶张量信号采用张量奇异值分解公式进行自适应分解以获得分量信号；所述张量奇异值分解公式包括以下步骤：31）构建任意阶张量与任意阶张量的o模式积乘法规则；32）根据任意阶张量与任意阶张量的o模式积乘法规则，得出三阶张量奇异值分解公式；33）根据三阶张量奇异值分解公式确定分解后张量的阶数和维数；4）对比所获得的分量信号的频率与机械故障特征频率是否一致，以诊断机械是否存在故障。2.根据权利要求1所述的用于机械故障诊断的三阶张量的奇异值分解方法，其特征在于，步骤31）中，首先设定m阶张量中有o阶的维数与n阶张量中o阶的维数对应相等，它们之间的乘法定义为o模式积，则m阶张量与n阶张量之间的o模式积为：；其中，中的p和q分别表示张量和张量中参与张量o模式积的阶数顺序为正序还是逆序，若为
‘1’
则表示正序，若为
‘2’
则表示逆序，且o=1, 2,
ꢀ…
, min(m, n)－1；I1、I2、
…
、I
m
分别表示张量的第一阶维数、第二阶维数、
…
第m阶维数；同理，J1、J2、
…
、J
n
分别表示张量的第一阶维数、第二阶维数、
…
第n阶维数；k1、 k2、
…
、k
m+n
‑
2o
分别表示o模式积后结果张量的第一阶维数、第二阶维数、
…
第m+n
‑
2o阶维数；为张量的第二阶至第m阶的元素形式表达式；为张量的第一阶至第n
‑
1阶的元素形式表达式；s1表示参与o模式积运算的第1个阶，S1为相应阶的维数；同理，s
o
表示参与o模式积运算的第o个阶，S
o
为相应阶的维数。3.根据权利要求2所述的用于机械故障诊断的三阶张量的奇异值分解方法，其特征在于，所述m阶张量与n阶张量之间的o模式积的阶数为正序或逆序时的参数计算公式如下所示：
；其中，和表示张量的第二阶至第m阶参数，如中的s1对应于张量的第二阶i2，以此类推，对应于张量的...

【专利技术属性】
技术研发人员：孟力，杨博淙，杨康定，王飞彪，
申请(专利权)人：赛腾机电科技常州有限公司，
类型：发明
国别省市：