基于随机矩阵理论与累积和的悬浮系统早期故障检测方法技术方案

本发明专利技术公开了基于随机矩阵理论与累积和的悬浮系统早期故障检测方法，包括以下步骤：S1、基于随机矩阵理论和大数据构建初始矩阵；S2、通过预设的移动窗口将初始矩阵转换为批量随机矩阵；S3、根据随机矩阵理论单环定理将批量随机矩阵转换为非Hermite矩阵；S4、计算非Hermite矩阵的特征值，并根据所述特征值得到该特征值的MSR，并将MSR作为悬浮系统的健康状态值；S5、依据得到的悬浮系统的健康状态值，通过累积和函数实现悬浮系统早期故障检测。本发明专利技术通过随机矩阵理论获得悬浮系统的健康状态值，并通过悬浮系统健康状态的历史值，利用累积和函数实现早期故障检测，具有检测简单、可靠性高的优点。

Early fault detection method of suspension system based on random matrix theory and cumulative sum

【技术实现步骤摘要】
基于随机矩阵理论与累积和的悬浮系统早期故障检测方法
本专利技术涉及磁悬浮列车
，特别是涉及一种基于随机矩阵理论与累积和的悬浮系统早期故障检测方法。
技术介绍
随着磁悬浮列车的普及，悬浮系统的安全性和可靠性越来越受到关注。在磁浮列车运行过程中，一旦悬浮系统发生故障，列车将无法运行。如果能在悬浮系统故障发生前检测出系统早期故障，将很大程度上避免这类事情发生。但悬浮系统健康程度在其工作寿命期间是逐渐降低其性能直至发生故障。其退化过程通常受到两种变化的影响，即个体变异性和时间变异性。因此，如何准确的检测悬浮系统的早期故障是目前急需解决的问题。
技术实现思路
有鉴于此，本专利技术提出基于随机矩阵理论与累积和的悬浮系统早期故障检测方法，利用随机矩阵理论和累积和函数实现早期故障检测并获得早期故障时刻，具有检测简单、可靠性高的优点。一方面，本专利技术提供了基于随机矩阵理论与累积和的悬浮系统早期故障检测方法，包括以下步骤：S1、基于随机矩阵理论和大数据构建初始矩阵；S2、通过预设的移动窗口将初始矩阵转换为批量随机矩阵；S3、根据随机矩阵理论单环定理将批量随机矩阵转换为非Hermite矩阵；S4、计算非Hermite矩阵的特征值，并根据所述特征值得到该特征值的MSR，并将MSR作为悬浮系统的健康状态值；S5、依据得到的悬浮系统的健康状态值，通过累积和函数实现悬浮系统早期故障检测。进一步地，所述初始矩阵通过如下表达方式描述：式中，A表示初始矩阵，T表示矩阵转置，R表示变量的数量，Nm表示第m个时间序列中的样本数，，表示R行Nm列矩阵集合，X1、X2、Xk、XR分别为初始矩阵A的元素，其表达式为：式中，Xk表示第k个变量在第m个时间序列中的一维随机矩阵，xk1、xk2、分别为一维随机矩阵Xk的元素。进一步地，所述批量随机矩阵的表达式为:Aj＝(Yj1,Yj2,...,Yjk,...,YjR)T∈CR×d(3)式中，Yj1、Yj2、Yjk、YjR分别为批量随机矩阵Aj中长度为d的向量，d表示移动窗口的宽度，CR×d表示R行d列矩阵集合，j∈[1,Nm-d+1]是批量随机矩阵Aj的序列号，其中：Yjk＝(yk1,yk2,...,ykd)∈C1×d(4)式中，yk1、yk2、ykd为Yjk的元素，C1×d表示l行d列矩阵集合。进一步地，所述移动窗口的宽度大小通过如下公式获取：d＝R*c(6)式中，c表示行列式的比值，并通过如下途径获取：1)设定c的最小值为：式中，ε∈(0，1]为防止d过大的比例因素；2)将选定的时间序列作为对象，令c从c的最小值到1之间以0.1为步长，获得不同的MSR曲线；3)比较获得的不同MSR曲线，取最佳c使曲线单调递减；4)使用其他时间序列来确定c是否合适。进一步地，所述步骤S3中非Hermite矩阵具体通过如下途径获得：当得到c后，对批量随机矩阵Aj中的元素进行归一化处理以获得标准的非Hermite矩阵的表达式为：式中，Bj1、Bj2、Bjk、BjR分别为非Hermite矩阵中长度为d的向量。进一步地，所述非Hermite矩阵的特征值通过如下步骤获得：设的特征值为其中i＝1,...,n，它们每一个都按照递增次序排列；依据如下公式，计算获得非Hermite矩阵的特征值式中，为的奇异值等价矩阵，为的转置矩阵，为的转置矩阵，其中，通过如下公式求取：式中，U为哈雷矩阵。进一步地，定义非Hermite矩阵的特征值的MSR为MSRj,q，用于反映矩阵的特征值分布，其通过如下公式求取：式中，MSRj,q表示第q时刻的悬浮系统健康状态值，是特征值的半径。进一步地，所述累计和函数被定义为：式中，是第q时刻向上的累计和，是第q时刻向下的累积和，是第q-1时刻向上的累计和，是第q-1时刻向下的累计和，u0是求得的不同MSR的平均值，s是求得的不同MSR的3倍标准差。进一步地，用于悬浮系统早期故障检测，并通过如下过程获得悬浮系统早期故障发生时刻：设定阈值θ，当且时，第q时刻就是悬浮系统早期故障发生时刻。本专利技术提供的基于随机矩阵理论与累积和的悬浮系统早期故障检测方法，主要通过两步实现早期故障检测：第一步，通过随机矩阵理论获得悬浮系统的健康状态值；第二步，通过悬浮系统健康状态的历史值，利用累积和函数实现早期故障检测，具有检测简单、可靠性高的优点。附图说明构成本专利技术的一部分的附图用来提供对本专利技术的进一步理解，本专利技术的示意性实施例及其说明用于解释本专利技术，并不构成对本专利技术的不当限定。在附图中：图1为本专利技术实施例提供的基于随机矩阵理论与累积和的悬浮系统早期故障检测方法的流程图。具体实施方式需要说明的是，在不冲突的情况下，本专利技术中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本专利技术。为更好地理解本专利技术，预先做出如下解释：图1本专利技术实施例提供的基于随机矩阵理论与累积和的悬浮系统早期故障检测方法的流程图。如图1所示，基于随机矩阵理论与累积和的悬浮系统早期故障检测方法，包括以下步骤：S1、基于随机矩阵理论和大数据构建初始矩阵；S2、通过预设的移动窗口将初始矩阵转换为批量随机矩阵；S3、根据随机矩阵理论单环定理将批量随机矩阵转换为非Hermite矩阵(非厄米特矩阵)；S4、计算非Hermite矩阵的特征值，并根据所述特征值得到该特征值的MSR(meanspecialradius，平均谱半径)，并将MSR作为悬浮系统的健康状态值；S5、依据得到的悬浮系统的健康状态值，通过累积和函数实现悬浮系统早期故障检测。从其上论述可知，本专利技术主要通过两步实现早期故障检测：第一步，通过随机矩阵理论(RandomMatrixTheory，RMT，通过数理统计方法分析复杂系统的特征值和能谱，以获得复杂系统的内在属性，可用于分析系统的相关状态，以反映系统的运行状态)获得悬浮系统的健康状态值；第二步，通过悬浮系统健康状态的历史值，利用累积和函数实现早期故障检测，具有检测简单、可靠性高的优点。具体地，上述步骤通过如下过程实现：假设变量的数量是R，第k个变量在第m个时间序列中的一维随机矩阵的表达如下：式中，Xk表示第k个变量在第m个时间序列中的一维随机矩阵，xk1、xk2、分别为一维随机矩阵Xk的元素，Nm表示第m个时间序列中的样本数。初始矩阵描述为：式中，A表示初始矩阵，T表示矩阵转置，R表示变量的数量，表示R行Nm列矩阵集合，X1、X2、Xk、XR分别为初始矩阵A的元素。设定移动窗口的宽度大小为d，通过宽度为d的移动窗口将初始矩阵A转换为批量随机矩阵Aj：本文档来自技高网...

【技术保护点】
1.基于随机矩阵理论与累积和的悬浮系统早期故障检测方法，其特征在于，包括以下步骤：/nS1、基于随机矩阵理论和大数据构建初始矩阵；/nS2、通过预设的移动窗口将初始矩阵转换为批量随机矩阵；/nS3、根据随机矩阵理论单环定理将批量随机矩阵转换为非Hermite矩阵；/nS4、计算非Hermite矩阵的特征值，并根据所述特征值得到该特征值的MSR，并将MSR作为悬浮系统的健康状态值；/nS5、依据得到的悬浮系统的健康状态值，通过累积和函数实现悬浮系统早期故障检测。/n

【技术特征摘要】
1.基于随机矩阵理论与累积和的悬浮系统早期故障检测方法，其特征在于，包括以下步骤：
S1、基于随机矩阵理论和大数据构建初始矩阵；
S2、通过预设的移动窗口将初始矩阵转换为批量随机矩阵；
S3、根据随机矩阵理论单环定理将批量随机矩阵转换为非Hermite矩阵；
S4、计算非Hermite矩阵的特征值，并根据所述特征值得到该特征值的MSR，并将MSR作为悬浮系统的健康状态值；
S5、依据得到的悬浮系统的健康状态值，通过累积和函数实现悬浮系统早期故障检测。


2.根据权利要求1所述的基于随机矩阵理论与累积和的悬浮系统早期故障检测方法，其特征在于，所述初始矩阵通过如下表达方式描述：



式中，A表示初始矩阵，T表示矩阵转置，R表示变量的数量，Nm表示第m个时间序列中的样本数，表示R行Nm列矩阵集合，X1、X2、Xk、XR分别为初始矩阵A的元素，其表达式为：



式中，Xk表示第k个变量在第m个时间序列中的一维随机矩阵，xk1、xk2、分别为一维随机矩阵Xk的元素。


3.根据权利要求2所述的基于随机矩阵理论与累积和的悬浮系统早期故障检测方法，其特征在于，所述批量随机矩阵的表达式为:
Aj＝(Yj1,Yj2,...,Yjk,...,YjR)T∈CR×d(3)
式中，Yj1、Yj2、Yjk、YjR分别为批量随机矩阵Aj中长度为d的向量，d表示移动窗口的宽度，CR×d表示R行d列矩阵集合，j∈[1,Nm-d+1]是批量随机矩阵Aj的序列号，其中：
Yjk＝(yk1,yk2,...,ykd)∈Cl×d(4)
式中，yk1、yk2、ykd为Yjk的元素，C1×d表示l行d列矩阵集合。


4.根据权利要求3所述的基于随机矩阵理论与累积和的悬浮系统早期故障检测方法，其特征在于，所述移动窗口的宽度大小通过如下公式获取：
d＝R*c(6)
式中，c表示行列式的比值，并通过如下途径获取：
1)设定c的最小值为：



式中，ε∈(0，1]为防止d过大的比例因素；
2)将选定的时间序列作为对象，令...

【专利技术属性】
技术研发人员：龙志强，王平，李晓龙，苗欣，江守亮，刘纪龙，
申请(专利权)人：中国人民解放军国防科技大学，中车青岛四方机车车辆股份有限公司，
类型：发明
国别省市：湖南;43