基于改进SVR复杂设备剩余使用寿命预测

摘要

本发明公开了基于改进SVR复杂设备剩余使用寿命预测，属于复杂设备寿命预测领域，预测结果中预测值与真实值很接近，波动趋势基本相同，这说明了半监督SVR保留了长期预测与利用历史数据的功能，且在性能方面表现更加出色，而相对于基于全局与局部预测的半监督SVR来说，半监督SVR在样本矩阵的组建方法上基本相同，但只需要进行一次全局预测的训练过程，大大地缩短了时间复杂度。实验验证了半监督SVR预测器同时具有几个优点：长期预测的能力，历史数据的充分利用，比较高的预测精度与少量的时间花费。其中在减少时间花费上的表现最为突出。

说明书

技术领域

本发明涉及复杂设备寿命预测领域，尤其涉及基于改进SVR复杂设备剩余使用寿命预测。

背景技术

随着工程设计、生产、组装、维修等技术的发展进步，工程中一些关键设备的可靠性得到了很大的提高，相应设备的寿命也大大延长。通常这类高可靠性、长寿命的设备的产量不高，并且价格昂贵。因此，应用传统的基于寿命数据的预测方法往往代价巨大。随着监测技术的迅速发展，状态监测数据(又叫退化数据)在可靠性评估中的应用越来越广泛，并且状态监测数据包含着更多的设备退化信息。因此，采用退化数据的寿命预测方法比采用寿命数据的方法更加受到专家学者的青睐。

设备的性能退化往往是由于其内部状态和外在工作环境共同作用的结果，并且往往表现出随机动态性。其中内部状态与使用材料有关，而外部工作环境则包括温度、适度、电磁环境、气压、振动等可能影响设备性能的因素。剩余寿命预测的结果是对设备进行检修、替换、部件订购等管理活动制定决策的基础，因此准确预测设备剩余寿命是最大化设备可靠性、最优化设备管理活动和最小化设备运行费用的关键。而要得到准确的剩余寿命预测结果，则应尽可能地将所有可能影响设备性能退化的因素考虑到所建立的退化模型中。

发明内容

本发明为解决上述问题，而提出的基于改进SVR复杂设备剩余使用寿命预测。

为了实现上述目的，本发明采用了如下技术方案：

基于改进SVR复杂设备剩余使用寿命预测，包括以下步骤；

S1、根据复杂设备可靠性、维修性和保障性函数建立复杂设备可靠性模型，其概率度量为“复杂设备完好率”；

S2、通过S1中建立的复杂设备可靠性模型，将时间序列数据变换成l×n维的历史数据样本矩阵X

S3、利用核函数K，将X

S4、计算

S5、设定参数C,ξ,ε,λ，并带入半监督SVR模型的优化函数中，求解二次规划问题，得到

S6、将

S7、将数据带入预测器f，得到预测值Y，Y为一个行向量，即最终结果是一个长期预测结果。

优选地，所述S1中建立的复杂设备可靠性模型时，首先建立复杂设备完成准备能进行运行的概率为：

P

式中：R(t)为前次运行过程中无故障概率，t为运行过程持续时间； P(t

其中从该系统执行前次运行后处于待修状态的瞬间开始计算，至下次运行命令到达前，该系统将充分恢复其规定状态的概率P(t

式中：m(t

将上述两式合并可得：

对于复杂设备的复杂系统，其故障可近似认为服从指数分布，而 R(t)为前次运行过程中无故障概率，t为运行过程持续时间，所以： R(t)＝e

对于指数分布，由于：

因此：

为计算复杂设备的完好率，假设前次运行发生故障且修理后执行第二次任务，再次发生故障于第二次任务结束以后，此时：

当δ

即完好率与可靠性参数MTBF、维修性参数MTTR及运行时间t 关系的表达式为：

优选地，所述S1中机械可靠性模型的边界条件是：该模型边界条件是：在先前运行任务中没有发生现在需要修复的故障，系统是随时可以再次立即使用的；即使是发生了这类故障，只要其维修时间不超过再次下达运行命令投入使用前的时间，系统将能够随时迎接下一任务，机械完好率是对机械运行状态的直接量度，假设保障维修系统其他因素一定，机械使用维修保障费用与其机械完好率的关系的重要度最大。

优选地，所述S4中的对称半正定矩阵式如下：

式中：Ω是一个(n×m)×(n×m)的对称半正定矩阵；K

优选地，所述S5中的半监督SVR模型的优化函数为：

式中，

与现有技术相比，本发明提供了基于改进SVR复杂设备剩余使用寿命预测，具备以下有益效果：

1.本发明预测结果中预测值与真实值很接近，波动趋势基本相同，这说明了半监督SVR保留了长期预测与利用历史数据的功能，且在性能方面表现更加出色，而相对于基于全局与局部预测的半监督SVR来说，半监督SVR在样本矩阵的组建方法上基本相同，但只需要进行一次全局预测的训练过程，大大地缩短了时间复杂度。实验验证了半监督SVR预测器同时具有几个优点：长期预测的能力，历史数据的充分利用，比较高的预测精度与少量的时间花费。其中在减少时间花费上的表现最为突出。

附图说明

图1为本发明提出的基于改进SVR复杂设备剩余使用寿命预测的资料实验结果图；

图2为本发明提出的基于改进SVR复杂设备剩余使用寿命预测的资料实验结果图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。

在本发明的描述中，需要理解的是，术语“上”、“下”、“前”、“后”、“左”、“右”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。

实施例1：

基于改进SVR复杂设备剩余使用寿命预测，包括以下步骤；

S1、根据复杂设备可靠性、维修性和保障性函数建立复杂设备可靠性模型，其概率度量为“复杂设备完好率”；

S2、通过S1中建立的复杂设备可靠性模型，将时间序列数据变换成l×n维的历史数据样本矩阵X

S3、利用核函数K，将X

S4、计算

S5、设定参数C,ξ,ε,λ，并带入半监督SVR模型的优化函数中，求解二次规划问题，得到

S6、将

S7、将数据带入预测器f，得到预测值Y，Y为一个行向量，即最终结果是一个长期预测结果。

进一步，优选地，S1中建立的复杂设备可靠性模型时，首先建立复杂设备完成准备能进行运行的概率为：

P

式中：R(t)为前次运行过程中无故障概率，t为运行过程持续时间； P(t

其中从该系统执行前次运行后处于待修状态的瞬间开始计算，至下次运行命令到达前，该系统将充分恢复其规定状态的概率P(t

式中：m(t

将上述两式合并可得：

对于复杂设备的复杂系统，其故障可近似认为服从指数分布，而 R(t)为前次运行过程中无故障概率，t为运行过程持续时间，所以： R(t)＝e

对于指数分布，由于：

因此：

为计算复杂设备的完好率，假设前次运行发生故障且修理后执行第二次任务，再次发生故障于第二次任务结束以后，此时：

当δ

即完好率与可靠性参数MTBF、维修性参数MTTR及运行时间t 关系的表达式为：

进一步，优选地，S1中机械可靠性模型的边界条件是：该模型边界条件是：在先前运行任务中没有发生现在需要修复的故障，系统是随时可以再次立即使用的；即使是发生了这类故障，只要其维修时间不超过再次下达运行命令投入使用前的时间，系统将能够随时迎接下一任务，机械完好率是对机械运行状态的直接量度，假设保障维修系统其他因素一定，机械使用维修保障费用与其机械完好率的关系的重要度最大。

进一步，优选地，S4中的对称半正定矩阵式如下：

式中：Ω是一个(n×m)×(n×m)的对称半正定矩阵；K

进一步，优选地，S5中的半监督SVR模型的优化函数为：

式中，

实施例2：

实验选取某一复杂设备的VI数据。训练资料采用了2019年8 月份前1000条数据做为当前数据，采用2018年1月份第1001-1400 共400条数据作为历史数据，将它们作为算法训练预测器过程的输入，得到的输出是2019年8月份第1001-1400条数据的预测值，然后与2019年8月份第1001-1400的400条实际数据进行比较分析；参数选择C＝4,ξ＝ξ

表1

实施例3：

实验选取某一复杂设备的VI数据。训练资料采用了2008年6 月份前1000条数据做为当前数据，采用2007年6月份第1001-1400 共400条数据作为历史数据，将它们作为算法训练预测器过程的输入，得到的输出是2008年6月份第1001-1400条数据的预测值，然后与2008年6月份第1001-1400的400条实际数据进行比较分析；参数选择C＝4,ξ＝ξ

表2

以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，根据本发明的技术方案及其发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。