基于稀疏退化建模的设备性能退化评估方法及系统

1.本发明涉及性能退化评估的技术领域，具体地，涉及基于稀疏退化建模的设备性能退化评估方法及系统，尤其涉及基于稀疏退化建模和平方包络谱融合的设备性能退化评估方法。


背景技术：

2.设备服役性能退化评估一般的策略是通过构造全寿命周期健康指数来监测设备的早期故障以及实现单调的退化评估。如何构造一个合适的全寿命周期健康指数来准确评估设备退化水平一直是国内外的研究热点。但是由于早期故障特征微弱以及强环境噪声的存在，构造具有早期故障退化趋势和单调退化趋势的健康指数有一定的挑战性。
3.在公开号为cn107356431a的专利文献中公开了一种基于admm与稀疏组合学习的滚动轴承性能退化评估方法，所述方法包括以下步骤：(1)对滚动轴承的正常模式进行学习，利用学到的知识建立正常模式的知识库；(2)建立基于各知识的量化指标求解器，求解被检测的滚动轴承的量化指标φ；(3)对量化指标φ进行分析，实现滚动轴承性能退化程度的评估。但所构造的健康指数不具有解释性，需要复杂的特征提取和组合学习过程。
4.在公开号为cn108398268a的专利文献中公开了一种基于堆叠去噪自编码器和自组织映射的轴承性能退化评估方法，它用于评估轴承退化技术领域。本发明解决了传统hi曲线构建中存在的提取退化特征需依赖大量专家经验和有监督训练，标签选择需依赖人工参与的问题。本发明的6个去噪自编码机构建堆叠去噪自编码器来对原始振动信号数据进行多层的特征提取，训练集数据对网络进行预训练后利用bp算法微调参数，将输出的100维特征输入som网络训练得到各时间点对应的hi，并构建训练集的hi曲线；将测试集数据输入训练好的堆叠去噪自编码器和som网络中得到各时间点处的hi，并构建hi曲线；分别对训练集和测试集的hi曲线进行平滑处理，获得平滑处理后的hi曲线。本发明可以应用于评估轴承性能退化领域用。但所构造的健康指数不具有解释性以及需要样本标签信息。
5.因此，需要提出一种新的技术方案以改善上述技术问题。


技术实现要素：

6.针对现有技术中的缺陷，本发明的目的是提供一种基于稀疏退化建模的设备性能退化评估方法及系统。
7.根据本发明提供的一种基于稀疏退化建模的设备性能退化评估方法，所述方法包括如下步骤：
8.步骤s1：以固定的时间间隔采集振动信号的方式，收集设备全寿命周期的时域振动数据矩阵x∈rn×m，表示如下：
[0009][0010]
其中，n表示设备全寿命周期过程中采样次数；m表示每次采样的样本数；x
j，
.∈r1×m表示第j次从设备采集到的振动样本序列，j＝1,2
…
,n；x
j,k
表示第j次采样中的第k个样本点，k＝1,2,
…
,m；
[0011]
步骤s2：基于希尔伯特变换和快速傅里叶变换将n次振动数据序列转化为平方包络谱幅值ses∈r
n∈m
，表示如下：
[0012][0013]
其中，ses
j,k
表示第j次采样中的第k个样本点所对应的平方包络谱幅值，由x
j,k
由基于希尔伯特变换和快速傅里叶变换得到；
[0014]
步骤s3：基于平方包络谱幅值融合的全寿命周期健康指数hi∈rn×1定义如下：
[0015][0016]
其中，hij(j＝1,2
…
,n)为第j次采样时间点的健康指数；w∈rm×1＝[w1,w2,
…
,wm]
′
为平方包络谱幅值的权重；wk(k＝1,2,
…
,m)为每次采样中的第k个平方包络谱幅值的权重；
[0017]
步骤s4：基于全寿命周期振动数据的平方包络谱幅值矩阵ses，构建一种稀疏退化凸优化模型自动确定平方包络谱幅值的权重w；
[0018]
步骤s5：根据步骤s4构建的稀疏退化凸优化模型求解平方包络谱幅值的权重w；
[0019]
步骤s6：根据步骤s3的全寿命周期健康指数定义计算hi，对设备全寿命周期进行早期故障检测和单调退化评估。
[0020]
优选地，所述步骤s1中的设备全寿命周期实验过程中共采集n次，每次采集m个振动数据序列，设备的全寿命周期振动数据矩阵x∈rn×m表示如下：
[0021][0022]
其中，n表示设备全寿命周期过程中采样次数；m表示每次采样的样本数；x
j，
.∈r1×m表示第j次从设备采集到的振动样本序列，j＝1,2
…
,n；x
j,k
表示第j次采样中的第k个样本点，k＝1,2,
…
,m。
[0023]
优选地，所述步骤s2包括如下步骤：
[0024]
步骤s2.1：首先利用希尔伯特变换运算，将n次的振动数据序列转化为解析信号矩阵e∈rn×m表示如下：
[0025][0026]
其中，|
·
|表示复数取模符号；x∈rn×m表示设备全寿命周期的时域振动数据矩阵；x
j，
.∈r1×m(j＝1,2
…
,n)表示第j次从设备采集到的振动样本序列；e
j,k
表示第j次采样中的第k个样本点所对应的解析信号，由x
j,k
由基于希尔伯特变换得到；hilbert(
·
)表示希尔伯特变换运算，具体运算公式如下：
[0027][0028]
由上式可知，希尔伯特变换运算的本质是将原始信号x(t)与1/πt之间进行卷积；
[0029]
步骤s2.2：计算设备的全寿命周期振动数据的解析信号e的平方包络幅值矩阵se∈rn×m如下：
[0030][0031]
其中，se
j,k
表示第j次采样中的第k个样本点所对应的平方包络幅值，由e
j,k
取幅值平方得到；
[0032]
步骤s2.3：基于快速傅里叶变换，计算设备的全寿命周期振动数据的平方包络谱幅值矩阵ses∈rn×m如下：
[0033][0034]
其中，ses
j,k
表示第j次采样中的第k个样本点所对应的平方包络谱幅值，由se
j,k
开展快速傅里叶变换得到。
[0035]
优选地，所述步骤s4中构建的权重w如下：
[0036][0037]
s.t.w'm
′
1＝r,mw≥0
[0038]
，dw≥0,hi＝sesw
[0039]
j＝1,2,
…
,n-1,εj≥0,
[0040]
其中，εj＝max(hi
j-hi
j 1
,0)是一个松弛变量，表示健康指数差分序列的单调性；cj是εj的权重，通过算术级数或者或几何级数来设定。
[0041]
优选地，所述cj的算术级数计算如下：
[0042][0043]cj
的几何级数计算如下：
[0044]cj
＝c1·qj-1
，j＝1,
…
,n-1
[0045]
其中，c1事先设定为一个常数来确定cj序列；因此，目标项用来表示健康指数序列非单调的数量；如果期望得到一个单调递增的健康指数序列，则εj应包含的零元素应越多越好，即越稀疏越好；r是平方包络谱幅值的权重和；λ是超参数；m∈rm×m是一个对角矩阵，其对角元素根据平方包络谱幅值的趋势确定；如果某一频率下的平方包络谱幅值具有递增的趋势，则m对应的对角元素为1，否则为-1；d∈r
(n-1)
×m表示平方包络谱幅值差分矩阵，具体表达式如下：
[0046][0047]
第一个约束条件w'm
′
1＝r表示平方包络谱幅值的权重和，具体表达式如下：
[0048][0049]
第二个约束条件mw≥0限制平方包络谱幅值的权重符号，具体表达式如下：
[0050][0051]
第三个约束条件dw≥0限制平方包络谱幅值差分与平方包络谱幅值的权重符号的一致，具体表达式如下：
[0052][0053]
本发明还提供一种基于稀疏退化建模的设备性能退化评估系统，所述系统包括如下模块：
[0054]
模块m1：以固定的时间间隔采集振动信号的方式，收集设备全寿命周期的时域振
动数据矩阵x∈rn×m，表示如下：
[0055][0056]
其中，n表示设备全寿命周期过程中采样次数；m表示每次采样的样本数；x
j，
.∈r1×m表示第j次从设备采集到的振动样本序列，j＝1,2
…
,n；x
j,k
表示第j次采样中的第k个样本点，k＝1,2,
…
,m；
[0057]
模块m2：基于希尔伯特变换和快速傅里叶变换将n次振动数据序列转化为平方包络谱幅值ses∈rn×m，表示如下：
[0058][0059]
其中，ses
j,k
表示第j次采样中的第k个样本点所对应的平方包络谱幅值，由x
j,k
由基于希尔伯特变换和快速傅里叶变换得到；
[0060]
模块m3：基于平方包络谱幅值融合的全寿命周期健康指数hi∈rn×1定义如下：
[0061][0062]
其中，hij(j＝1,2
…
,n)为第j次采样时间点的健康指数；w∈rm×1＝[w1,w2,
…
,wm]
′
为平方包络谱幅值的权重；wk(k＝1,2,
…
,m)为每次采样中的第k个平方包络谱幅值的权重；
[0063]
模块m4：基于全寿命周期振动数据的平方包络谱幅值矩阵ses，构建一种稀疏退化凸优化模型自动确定平方包络谱幅值的权重w；
[0064]
模块m5：根据模块m4构建的稀疏退化凸优化模型求解平方包络谱幅值的权重w；
[0065]
模块m6：根据模块m3的全寿命周期健康指数定义计算hi，对设备全寿命周期进行早期故障检测和单调退化评估。
[0066]
优选地，所述模块m1中的设备全寿命周期实验过程中共采集n次，每次采集m个振动数据序列，设备的全寿命周期振动数据矩阵x∈rn×m表示如下：
[0067][0068]
其中，n表示设备全寿命周期过程中采样次数；m表示每次采样的样本数；x
j，
.∈r1×m表示第j次从设备采集到的振动样本序列，j＝1,2
…
,n；x
j,k
表示第j次采样中的第k个样
本点，k＝1,2,
…
,m。
[0069]
优选地，所述模块m2包括如下模块：
[0070]
模块m2.1：首先利用希尔伯特变换运算，将n次的振动数据序列转化为解析信号矩阵e∈rn×m表示如下：
[0071][0072]
其中，|
·
|表示复数取模符号；x∈rn×m表示设备全寿命周期的时域振动数据矩阵；x
j，
.∈r1×m(j＝1,2
…
,n)表示第j次从设备采集到的振动样本序列；e
j,k
表示第j次采样中的第k个样本点所对应的解析信号，由x
j,k
由基于希尔伯特变换得到；hilbert(
·
)表示希尔伯特变换运算，具体运算公式如下：
[0073][0074]
由上式可知，希尔伯特变换运算的本质是将原始信号x(t)与1/πt之间进行卷积；
[0075]
模块m2.2：计算设备的全寿命周期振动数据的解析信号e的平方包络幅值矩阵se∈rn×m如下：
[0076][0077]
其中，se
j,k
表示第j次采样中的第k个样本点所对应的平方包络幅值，由e
j,k
取幅值平方得到；
[0078]
模块m2.3：基于快速傅里叶变换，计算设备的全寿命周期振动数据的平方包络谱幅值矩阵ses∈rn×m如下：
[0079][0080]
其中，ses
j,k
表示第j次采样中的第k个样本点所对应的平方包络谱幅值，由se
j,k
开展快速傅里叶变换得到。
[0081]
优选地，所述模块m4中构建的权重w如下：
[0082][0083]
s.t.w'm
′
1＝r,mw≥0
[0084]
，dw≥0,hi＝sesw
[0085]
j＝1,2,
…
,n-1,εj≥0,
[0086]
其中，εj＝max(hi
j-hi
j 1
,0)是一个松弛变量，表示健康指数差分序列的单调性；cj是εj的权重，通过算术级数或者或几何级数来设定。
[0087]
优选地，所述cj的算术级数计算如下：
[0088][0089]cj
的几何级数计算如下：
[0090]cj
＝c1·qj-1
，j＝1,
…
,n-1
[0091]
其中，c1事先设定为一个常数来确定cj序列；因此，目标项用来表示健康指数序列非单调的数量；如果期望得到一个单调递增的健康指数序列，则εj应包含的零元素应越多越好，即越稀疏越好；r是平方包络谱幅值的权重和；λ是超参数；m∈rm×m是一个对角矩阵，其对角元素根据平方包络谱幅值的趋势确定；如果某一频率下的平方包络谱幅值具有递增的趋势，则m对应的对角元素为1，否则为-1；d∈r
(n-1)
×m表示平方包络谱幅值差分矩阵，具体表达式如下：
[0092][0093]
第一个约束条件w'm
′
1＝r表示平方包络谱幅值的权重和，具体表达式如下：
[0094][0095]
第二个约束条件mw≥0限制平方包络谱幅值的权重符号，具体表达式如下：
[0096][0097]
第三个约束条件dw≥0限制平方包络谱幅值差分与平方包络谱幅值的权重符号的一致，具体表达式如下：
[0098][0099]
与现有技术相比，本发明具有如下的有益效果：
[0100]
1、本发明通过采用平方包络谱幅值融合的方法，来构造全寿命健康指数，解决设
备全寿命周期早期故障检测和单调退化评估的问题；
[0101]
2、本发明通过一种稀疏退化凸优化模型自动确定平方包络谱幅值的权重，解决了手动提取特征和筛选特征的问题；
[0102]
3、本发明所构建的健康指数具有明确数学定义和物理可解释性；
[0103]
4、本发明所优化得到的平方包络谱幅值的权重，能有效提取出稀疏的设备故障特征频率用于设备全寿命周期性能退化评估。
附图说明
[0104]
通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：
[0105]
图1为本发明的流程图；
[0106]
图2为本发明全寿命周期健康指数及其平方包络谱幅值的权重示例1的示意图；
[0107]
图3为本发明全寿命周期健康指数及其平方包络谱幅值的权重示例2的示意图。
具体实施方式
[0108]
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变化和改进。这些都属于本发明的保护范围。
[0109]
本发明利用故障特征信号在平方包络谱上的稀疏特性,提出一种基于稀疏退化建模和平方包络谱融合的设备性能退化评估方法。该方法通过构造稀疏退化评估凸优化模型，自动优化和确定平方包络谱的权重来构造全寿命周期健康指数。基于所构造全寿命周期健康指数以及优化的包络谱权重，实现设备早期故障检测和单调退化趋势评估。
[0110]
一种基于稀疏退化建模的设备性能退化评估方法包括以下步骤：
[0111]
步骤s1：假设以固定的时间间隔采集一次振动信号的方式，收集到一台设备全寿命周期的时域振动数据。设备的全寿命周期振动数据矩阵x∈rn×m表示如下：
[0112][0113]
其中，n表示设备全寿命周期过程中采样次数；m表示每次采样的样本数；x
j，
.∈r1×m表示第j次从设备采集到的振动样本序列，j＝1,2
…
,n；x
j,k
表示第j次采样中的第k个样本点，k＝1,2,
…
,m。
[0114]
步骤s2：基于希尔伯特变换和快速傅里叶变换将n次的振动数据序列转化为平方包络谱幅值ses∈rn×m，表示如下：
[0115][0116]
其中，ses
j,k
表示第j次采样中的第k个样本点所对应的平方包络谱幅值，由x
j,k
由基于希尔伯特变换和快速傅里叶变换得到。
[0117]
具体流程如下：
[0118]
步骤s2.1：首先利用希尔伯特变换运算，将n次的振动数据序列转化为解析信号矩阵e∈rn×m表示如下：
[0119][0120]
其中，|
·
|表示复数取模符号；x∈rn×m表示设备全寿命周期的时域振动数据矩阵；x
j，
.∈r1×m(j＝1,2
…
,n)表示第j次从设备采集到的振动样本序列；e
j,k
表示第j次采样中的第k个样本点所对应的解析信号，由x
j,k
由基于希尔伯特变换得到；hilbert(
·
)表示希尔伯特变换运算，具体运算公式如下：
[0121][0122]
由上式可知，希尔伯特变换运算的本质是将原始信号x(t)与1/πt之间进行卷积。
[0123]
步骤s2.2：计算设备的全寿命周期振动数据的解析信号e的平方包络幅值矩阵se∈rn×m如下：
[0124][0125]
其中，se
j,k
表示第j次采样中的第k个样本点所对应的平方包络幅值，由e
j,k
取幅值平方得到；
[0126]
步骤s2.3：基于快速傅里叶变换，计算设备的全寿命周期振动数据的平方包络谱幅值矩阵ses如下：
[0127][0128]
其中，ses
j,k
表示第j次采样中的第k个样本点所对应的平方包络谱幅值，由x
j,k
由
基于希尔伯特变换和快速傅里叶变换得到。由于每一次从设备采集的振动数据序列的平方包络谱幅值是对称的，所以一般取其前一半的平方包络谱幅值。
[0129]
步骤s3：基于平方包络谱幅值融合的全寿命周期健康指数hi∈rn×1定义如下：
[0130][0131]
其中，hij(j＝1,2
…
,n)表示第j次采样时间点的健康指数；w∈rm×1＝[w1,w2,
…
,wm]
′
表示平方包络谱幅值的权重；wk(k＝1,2,
…
,m)表示每次采样中的第k个平方包络谱幅值的权重。
[0132]
步骤s4：基于全寿命周期振动数据的平方包络谱幅值矩阵ses，构建了一种稀疏退化凸优化模型自动确定平方包络谱幅值的权重w，构建如下：
[0133][0134]
s.t.w'm
′
1＝r,mw≥0
[0135]
，dw≥0,hi＝sesw
[0136]
j＝1,2,
…
,n-1,εj≥0,
[0137]
其中，εj＝max(hi
j-hi
j 1
,0)是一个松弛变量，表示健康指数差分序列的单调性；cj是εj的权重，可通过算术级数或者或几何级数来设定；
[0138]cj
的算术级数可以计算如下：
[0139][0140]cj
的几何级数可以计算如下：
[0141]cj
＝c1·qj-1
，j＝1,
…
,n-1
[0142]
其中，c1事先设定为一个常数来确定cj序列；因此，目标项可以用来表示健康指数序列非单调的数量；如果期望得到一个单调递增的健康指数序列，则εj应包含的零元素应越多越好，即越稀疏越好；因此，||ε||1用以增强其稀疏性来提高健康指数序列的单调性；||w||1用以增强平方包络谱幅值权重的稀疏性；r是平方包络谱幅值的权重和，通常可以设定为1或者其它整数；λ是超参数，用来调整健康指数单调性和平方包络谱幅值权重的稀疏性；m∈rm×m是一个对角矩阵，其对角元素可以根据平方包络谱幅值的趋势来确定；如果某一频率下的平方包络谱幅值具有递增的趋势，则m对应的对角元素为1，否则为-1；d∈r
(n-1)
×m表示平方包络谱幅值差分矩阵，具体表达式如下：
[0143][0144]
第一个约束条件w'm
′
1＝r表示平方包络谱幅值的权重和，具体表达式如下：
[0145][0146]
第二个约束条件mw≥0限制平方包络谱幅值的权重符号，具体表达式如下：
[0147][0148]
第三个约束条件dw≥0限制平方包络谱幅值差分与平方包络谱幅值的权重符号的一致，用于进一步提升健康指数的单调退化趋势，具体表达式如下：
[0149][0150]
步骤s5：根据步骤s4构建的稀疏退化凸优化模型求解平方包络谱幅值的权重w。
[0151]
步骤s6：根据步骤s3的全寿命周期健康指数定义计算hi，用于设备全寿命周期的早期故障检测和单调退化评估。
[0152]
本发明还提供一种基于稀疏退化建模的设备性能退化评估系统，所述系统包括如下模块：
[0153]
模块m1：以固定的时间间隔采集振动信号的方式，收集设备全寿命周期的时域振动数据矩阵x∈rn×m，表示如下：
[0154][0155]
其中，n表示设备全寿命周期过程中采样次数；m表示每次采样的样本数；x
j，
.∈r1×m表示第j次从设备采集到的振动样本序列，j＝1,2
…
,n；x
j,k
表示第j次采样中的第k个样本点，k＝1,2,
…
,m。
[0156]
模块m1中的设备全寿命周期实验过程中共采集n次，每次采集m个振动数据序列，设备的全寿命周期振动数据矩阵x∈rn×m表示如下：
[0157][0158]
其中，n表示设备全寿命周期过程中采样次数；m表示每次采样的样本数；x
j，
.∈r1×m表示第j次从设备采集到的振动样本序列，j＝1,2
…
,n；x
j,k
表示第j次采样中的第k个样
本点，k＝1,2,
…
,m。
[0159]
模块m2：基于希尔伯特变换和快速傅里叶变换将n次振动数据序列转化为平方包络谱幅值ses∈rn×m，表示如下：
[0160][0161]
其中，ses
j,k
表示第j次采样中的第k个样本点所对应的平方包络谱幅值，由x
j,k
由基于希尔伯特变换和快速傅里叶变换得到。
[0162]
模块m2.1：首先利用希尔伯特变换运算，将n次的振动数据序列转化为解析信号矩阵e∈rn×m表示如下：
[0163][0164]
其中，|
·
|表示复数取模符号；x∈rn×m表示设备全寿命周期的时域振动数据矩阵；x
j，
.∈r1×m(j＝1,2
…
,n)表示第j次从设备采集到的振动样本序列；e
j,k
表示第j次采样中的第k个样本点所对应的解析信号，由x
j,k
由基于希尔伯特变换得到；hilbert(
·
)表示希尔伯特变换运算，具体运算公式如下：
[0165][0166]
由上式可知，希尔伯特变换运算的本质是将原始信号x(t)与1/πt之间进行卷积；
[0167]
模块m2.2：计算设备的全寿命周期振动数据的解析信号e的平方包络幅值矩阵se∈rn×m如下：
[0168][0169]
其中，se
j,k
表示第j次采样中的第k个样本点所对应的平方包络幅值，由e
j,k
取幅值平方得到；
[0170]
模块m2.3：基于快速傅里叶变换，计算设备的全寿命周期振动数据的平方包络谱幅值矩阵ses∈rn×m如下：
[0171][0172]
其中，ses
j,k
表示第j次采样中的第k个样本点所对应的平方包络谱幅值，由se
j,k
开展快速傅里叶变换得到。
[0173]
模块m3：基于平方包络谱幅值融合的全寿命周期健康指数hi∈rn×1定义如下：
[0174][0175]
其中，hij(j＝1,2
…
,n)为第j次采样时间点的健康指数；w∈rm×1＝[w1,w2,
…
,wm]
′
为平方包络谱幅值的权重；wk(k＝1,2,
…
,m)为每次采样中的第k个平方包络谱幅值的权重。
[0176]
模块m4：基于全寿命周期振动数据的平方包络谱幅值矩阵ses，构建一种稀疏退化凸优化模型自动确定平方包络谱幅值的权重w。
[0177]
模块m4中构建的权重w如下：
[0178][0179]
s.t.w'm
′
1＝r,mw≥0
[0180]
，dw≥0,hi＝sesw
[0181]
j＝1,2,
…
,n-1,εj≥0,
[0182]
其中，εj＝max(hi
j-hi
j 1
,0)是一个松弛变量，表示健康指数差分序列的单调性；cj是εj的权重，通过算术级数或者或几何级数来设定。
[0183]cj
的算术级数计算如下：
[0184][0185]cj
的几何级数计算如下：
[0186]cj
＝c1·qj-1
，j＝1,
…
,n-1
[0187]
其中，c1事先设定为一个常数来确定cj序列；因此，目标项用来表示健康指数序列非单调的数量；如果期望得到一个单调递增的健康指数序列，则εj应包含的零元素应越多越好，即越稀疏越好；r是平方包络谱幅值的权重和；λ是超参数；m∈rm×m是一个对角矩阵，其对角元素根据平方包络谱幅值的趋势确定；如果某一频率下的平方包络谱幅值具有递增的趋势，则m对应的对角元素为1，否则为-1；d∈r
(n-1)
×m表示平方包络谱幅值差分矩阵，具体表达式如下：
[0188][0189]
第一个约束条件w'm
′
1＝r表示平方包络谱幅值的权重和，具体表达式如下：
[0190][0191]
第二个约束条件mw≥0限制平方包络谱幅值的权重符号，具体表达式如下：
[0192][0193]
第三个约束条件dw≥0限制平方包络谱幅值差分与平方包络谱幅值的权重符号的一致，具体表达式如下：
[0194][0195]
模块m5：根据模块m4构建的稀疏退化凸优化模型求解平方包络谱幅值的权重w。
[0196]
模块m6：根据模块m3的全寿命周期健康指数定义计算hi，对设备全寿命周期进行早期故障检测和单调退化评估。
[0197]
本领域技术人员知道，除了以纯计算机可读程序代码方式实现本发明提供的系统及其各个装置、模块、单元以外，完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统及其各个装置、模块、单元以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同功能。所以，本发明提供的系统及其各项装置、模块、单元可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置、模块、单元也可以视为硬件部件内的结构；也可以将用于实现各种功能的装置、模块、单元视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
[0198]
以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变化或修改，这并不影响本发明的实质内容。在不冲突的情况下，本技术的实施例和实施例中的特征可以任意相互组合。