性能退化产品可靠度置信区间估计数值仿真方法技术

一种性能退化产品可靠度置信区间估计数值仿真方法，其步骤如下。1）确定随机过程类型并估计特征参数；2）确定各时间点可靠度的估计值；3）产生一定数目的伪样本；4）计算各伪样本各时间点可靠度的估计值；5）得到可靠度置信区间。本发明专利技术为性能退化产品可靠度提供了一种符合适应性强、操作性强的置信区间估计方法，为产品的寿命设计、分析及验证提供依据和保障。

Numerical simulation method for reliability confidence interval estimation of performance degraded products

【技术实现步骤摘要】
性能退化产品可靠度置信区间估计数值仿真方法一、
本专利技术涉及性能退化产品在退化试验下的可靠度置信区间估计方法，具体是基于随机过程的性能退化产品可靠度置信区间估计数值仿真方法，属于可靠性建模技术与寿命预测分析领域，用于对产品的可靠度进行区间估计，为产品的寿命设计、分析及验证提供依据和保障。二、
技术介绍
目前对产品退化失效的研究中，大多采用随机系数或混合随机系数的回归模型描述产品的退化过程，退化量的随机性是利用模型中的随机系数表现的。但是，即便是在额定工作环境下工作，产品也可能会受到许多随机因素的影响，而随机系数或混合随机系数的回归模型中的随机系数，一般只能表征产品个体样品间的特性差异的随机性，它很难描述环境对产品状态变化的随机影响，因此，利用回归模型描述产品的退化规律时，将存在一定的局限性。然而，产品状态随时间的变化本质上是一个随机过程，因此使用随机过程进行退化建模将更符合工程实践。由于随机过程具有良好的计算和分析性质，我们很容易得到基于随机过程的产品可靠度的点估计，不过，仅仅得到点估计还远不能保证产品可靠度的计算精度，为了提升产品可靠度估计的可信性，对可靠度进行置信区间估计是必要的。由于不存在简单方法去估计性能退化产品可靠度的标准差，这时可利用数值仿真方法对基于随机过程的产品可靠度进行置信区间估计。本专利技术利用数值仿真方法对基于随机过程的产品可靠度进行区间估计，为提升产品可靠度估计的可信性提供保障。三、
技术实现思路
（一）目的本专利技术的目的是针对性能退化产品的特点和可靠度置信区间估计的需求以及现有技术的不足，提供了一种适应性强、操作性强的性能退化产品可靠度置信区间估计数值仿真方法，一方面降低了产品可靠度置信区间的估计难度，另一方面，为提升产品可靠度估计的可信性提供保障。（二）技术方案一种性能退化产品可靠度置信区间估计数值仿真方法，该方法是以性能退化产品的性能退化数据为基础，利用所属随机过程特征参数，生成一系列的Bootstrap伪样本，每个样本是单位时间退化量的有放回抽样。通过对伪样本的一系列计算，获得产品可靠度的置信区间。本专利技术一种性能退化产品可靠度置信区间估计数值仿真方法，其步骤如下。步骤1，确定随机过程类型并估计特征参数。通过性能退化数据，依据性能退化数据所属随机过程类型，确定随机过程特征参数，并计算得到特征参数估计值。特征参数估计值可采用矩估计和极大似然估计进行确定。步骤2，确定各时间点可靠度的估计值。通过步骤１得到的特征参数估计值，根据退化数据所属随机过程，计算得到各时间点可靠度的点估计值。步骤3，产生一定数目的伪样本。利用步骤１得到的特征参数估计值，根据退化数据所属随机过程，生成一定数目的伪样本。步骤4，计算各伪样本各时间点可靠度的估计值。利用伪样本计算得到对应退化轨道，其中每条退化轨道就像原始样本那样，计算各时间点可靠度的估计值。步骤5，得到可靠度置信区间。根据既定置信度，对各时间点作出可靠度的Bootstrap区间估计。其中，所述的“性能退化产品”，是指本专利技术所述的方法适用于下列产品，该性能退化产品包括：（1）长贮退化失效型产品，是指在产品全寿命周期过程中，其绝大多数时间处于贮存状态的产品，该产品随着服役时间的增加，将产生劣化、老化现象；（2）长时间连续工作退化失效型产品，是指在产品全寿命周期过程中，其绝大多数时间处于工作状态的产品，该产品随着服役时间的增加，将产生劣化、老化现象；（3）非连续工作退化失效型产品，是指在产品全寿命周期过程中，其一部分时间处于贮存状态、一部分时间处于工作状态、且贮存与工作状态不断交替的产品，该产品随着服役时间的增加，将产生劣化、老化现象。其中，在步骤1中所述的“性能退化数据”，是指本专利技术所述的方法是针对该类数据使用的，该类数据是通过对表征性能退化产品功能的某些量进行连续测量获得的。本专利技术的性能退化数据可以是一元的，也可以是多元的，若性能退化数据是多元的，可以采用仿真的方法进行可靠度的点估计求解。本专利技术需至少获得一个样本的性能退化数据，若同时获得了多个样本的性能退化数据，则需先根据各产品监测时刻点，采用插值方法对性能数据进行插值，将各样本的测量时刻对齐；然后得到各测量时刻的样本均值；再将样本均值随时间变化的序列看作是单样本性能变化数据，这样就将多样本数据转化成了单样本性能退化数据。其中，在步骤1中所述的“随机过程类型”，是指本专利技术所述的方法是在下列随机过程类型基础上进行的，该随机过程类型包括：维纳过程、复合泊松过程和伽马过程中的一种或数种。其中，在步骤1中所述的特征参数估计方法：是指矩估计法、极大似然估计法中的一种或数种。这些方法的具体情况如下。（1）矩估计法：假设参数定义阶矩为.而阶样本矩为.,。解出上述方程组，即可求得参数的矩估计。（2）极大似然估计法：。似然函数定义为.极大似然估计是使得其中，在步骤2中所述的“根据退化数据所属随机过程，计算得到各时间点可靠度的点估计值”，是指针对不同随机过程，其可靠度估计方法各不一样。针对常见随机过程类型具体情况如下。（1）维纳过程：由于其寿命的分布服从逆高斯分布，可以得到其产品可靠度的点估计为其中，为失效阈值，，为维纳过程扩散系数的估计值。（2）伽马过程：其产品可靠度的点估计为其中，为失效阈值，为不完全Gamma函数，,分布形状参数的估计值，分布尺度参数的估计值。（3）复合泊松过程：其产品可靠度的点估计为其中，为失效阈值，参数的估计值，为漂移系数的估计值，为扩散系数的估计值。其中，在步骤3中所述的产生伪样本的方法，是指产生一个数B,B是Bootstrap抽取的样本个数，具体如下。（a）利用，b=1,…,B，在每个Bootstrap样本下得到T个单位时间退化量，T与原始样本相同，于是得到B条退化轨道；（b）利用B条退化轨道，得出，b=1,…,B。其中，在步骤4中所述的计算各伪样本各时间点可靠度的估计值的方法，是指每个伪样本就像原始样本那样，用来作出的估计值。其中，在步骤5中所述的计算可靠度置信区间的方法，是指对每时刻t作出的Bootstrap置信区间，具体如下。（a）把在步骤４中获得的B个的估计值，b=1,…,B，按从小到大排列成，其中，t为定值；（b）分布函数的近似的置信区间为其中，（三）本专利技术的优点：本专利技术在考虑充分性、必要性、符合工程习惯、可论证性、可设计性、可验证性原则的基础上，根据性能退化产品的特点提出一种科学、规范化可靠度置信区间估计数值仿真方法，一方面填补了以往针对性能退化产品可靠度置信区间计算量大、算法复杂的缺陷；另一方面为性能退化产品的寿命设计、分析及验证提供依据和保障。四、附图说明本文档来自技高网...

【技术保护点】
1.一种性能退化产品可靠度置信区间估计数值仿真方法，该方法是以性能退化产品的性能退化数据为基础，利用所属随机过程特征参数，生成一系列的Bootstrap伪样本，每个样本是单位时间退化量的有放回抽样，通过对伪样本的一系列计算，获得产品可靠度的置信区间，方法的具体步骤如下：/n步骤1，确定随机过程类型并估计特征参数，通过性能退化数据，依据性能退化数据所属随机过程类型，确定随机过程特征参数，并计算得到特征参数估计值，特征参数估计值可采用矩估计和极大似然估计进行确定；/n步骤2，确定各时间点可靠度的估计值，通过步骤1得到的特征参数估计值，根据退化数据所属随机过程，计算得到各时间点可靠度的点估计值；/n步骤3，产生一定数目的伪样本，利用步骤1得到的特征参数估计值，根据退化数据所属随机过程，生成一定数目的伪样本；/n步骤4，计算各伪样本各时间点可靠度的估计值，利用伪样本计算得到对应退化轨道，其中每条退化轨道按照原始样本的处理方法，计算各时间点可靠度的估计值；/n步骤5，得到可靠度置信区间，根据既定置信度，对各时间点作出可靠度的Bootstrap区间估计。/n

【技术特征摘要】
1.一种性能退化产品可靠度置信区间估计数值仿真方法，该方法是以性能退化产品的性能退化数据为基础，利用所属随机过程特征参数，生成一系列的Bootstrap伪样本，每个样本是单位时间退化量的有放回抽样，通过对伪样本的一系列计算，获得产品可靠度的置信区间，方法的具体步骤如下：
步骤1，确定随机过程类型并估计特征参数，通过性能退化数据，依据性能退化数据所属随机过程类型，确定随机过程特征参数，并计算得到特征参数估计值，特征参数估计值可采用矩估计和极大似然估计进行确定；
步骤2，确定各时间点可靠度的估计值，通过步骤1得到的特征参数估计值，根据退化数据所属随机过程，计算得到各时间点可靠度的点估计值；
步骤3，产生一定数目的伪样本，利用步骤1得到的特征参数估计值，根据退化数据所属随机过程，生成一定数目的伪样本；
步骤4，计算各伪样本各时间点可靠度的估计值，利用伪样本计算得到对应退化轨道，其中每条退化轨道按照原始样本的处理方法，计算各时间点可靠度的估计值；
步骤5，得到可靠度置信区间，根据既定置信度，对各时间点作出可靠度的Bootstrap区间估计。


2.根据权利要求1所述的一种性能退化产品可靠度置信区间估计数值仿真方法，其特征在于，所述的“性能退化产品”包括：（1）长贮退化失效型产品，指在产品全寿命周期过程中，其绝大多数时间处于贮存状态的产品，该产品随着服役时间的增加，将产生劣化、老化现象；（2）长时间连续工作退化失效型产品，指在产品全寿命周期过程中，其绝大多数时间处于工作状态的产品，该产品随着服役时间的增加，将产生劣化、老化现象；（3）非连续工作退化失效型产品，指在产品全寿命周期过程中，其一部分时间处于贮存状态、一部分时间处于工作状态、且贮存与工作状态不断交替的产品，该产品随着服役时间的增加，将产生劣化、老化现象。


3.根据权利要求1所述的一种性能退化产品可靠度置信区间估计数值仿真方法，其特征在于，在步骤1中所述的“性能退化数据”是指：本发明所述的方法主要用于该类数据，该类数据是通过对表征性能退化产品功能的某些量进行连续测量获得的，本发明的性能退化数据可以是一元的，也可以是多元的，若性能退化数据是多元的，可以采用仿真的方法进行可靠度的点估计求解，本发明需至少获得一个样本的性能退化数据，若同时获得了多个样本的性能退化数据，则需先根据各产品监测时刻点，采用插值方法对性能数据进行插值，将各样本的测量时刻对齐；然后得到各测量时刻的样本均值；再将样本均值随时间变化的序列看作是单样本性能变化数据，这样就将多样本数据转化成了单样本性能退化数据。


4.根据权利要求1所述的一种性能退化产品可靠度置信区间估计数值仿真方法，其特征在于，在步骤...

【专利技术属性】
技术研发人员：冯静，孙权，
申请(专利权)人：湖南银杏可靠性技术研究所有限公司，
类型：发明
国别省市：湖南;43