CN110309487B - 基于备件保障数据的指数型单元寿命分布参数估计方法 - Google Patents

Abstract

本发明涉及一种基于备件保障数据的指数型单元寿命分布参数估计方法，该方法首先根据保障任务的相关数据生成n组候选分布参数，然后初始化似然值；再根据k次保障任务的相关数据计算保障任务概率，更新似然值；最后在更新后的似然值中找到最大似然值，则最大似然值对应的指数分布的均值参数即为估计结果。本发明提出的寿命分布参数估计方法，能在因缺少专门的可靠性试验导致缺少寿命值数据时，有效利用实际工作产生的备件保障数据，较为准确地估计寿命分布参数。

Description

基于备件保障数据的指数型单元寿命分布参数估计方法

技术领域

本发明涉及产品可靠性预测技术领域，具体涉及一种基于备件保障数据的指数型单元寿命分布参数估计方法。

背景技术

产品由各种单元组成。通常以产品的寿命分布函数来定量描述产品的可靠性。在本文中，指数型单元是指寿命服从指数分布的单元。一般来说，电子零部件等单元的寿命服从指数分布，如：印制电路板插件、电子部件、电阻、电容、集成电路等。寿命X服从指数分布记作X～Exp(μ)，μ为指数分布的均值参数，X的密度函数为

寿命分布函数由函数类别和分布参数决定，本质上是一种描述统计规律的函数。因此，只有当寿命样本值达到一定数量规模时，该函数才能较为准确的描述产品寿命的分布规律。在实际工作中，工作环境往往对产品可靠性有较大影响，因此需要关注工作环境下产品的寿命分布规律。在实际工作中，产品工况数据的记录质量与专门的可靠性试验相比往往较差，常常由于记录不及时或记录不全，导致故障发生时间、产品累积工作时间等直接反映产品寿命值的相关数据丢失，只留有保障任务开始前的备件数量、任务成功与否、任务计划时间等数据。这使得基于寿命值数据的分布参数常规估计方法难以从这种备件保障数据中分析产品的可靠性。

发明内容

本发明针对现有技术中存在的技术问题，提供一种基于备件保障数据的指数型单元寿命分布参数估计方法，该方法约定：共有k次保障任务的相关数据，第i次保障任务时该单元的计划工作时间记为Tw_i(以下简称为任务时间)，该单元在装备内的装机数量记为Nzj_i，任务开始前的备件数量记为Nbj_i，保障任务结果记为F_i，1≤i≤k。在保障任务期间，Nzj_i个该型单元之间视为串联关系，以备件更换故障件的方式来排除故障，若发生的故障次数不大于备件数量，则本次保障任务成功(F_i记为1)，否则视为失败(F_i记为0)。即：以形如[Tw_i Nzj_i Nbj_i F_i]的备件保障数据来描述一次保障任务。在积累了k组备件保障数据后，估计寿命分布参数。

本发明解决上述技术问题的技术方案如下：

一种基于备件保障数据的指数型单元寿命分布参数估计方法，包括以下步骤：

步骤1，根据保障任务的相关数据生成n组候选分布参数μ_j,1≤j≤n，其中，μ_j表示指数分布的均值参数，n为正整数；

步骤2，初始化似然值L_j,1≤j≤n，令L_j＝0；

步骤3，根据k次保障任务的相关数据计算保障任务概率P1，更新似然值L_j；

步骤4，在更新后的似然值L_j(1≤j≤n)中找到最大似然值，记为L_M，则似然值L_M对应的μ_M即为指数分布的均值参数的估计结果。

进一步的，步骤1中所述的候选分布参数的生成方法包括：

参数μ_j＝μ_min+(j-1)d1,1≤j≤n，其中，

μ_max表示指数分布的均值参数上限，μ_min表示指数分布的均值参数下限，n为正整数，且n≥2。

进一步的，所述步骤3包括：

步骤3.1，令候选参数序号j＝1；

步骤3.2，令任务序号i＝1，1≤i≤k，k为保障任务最大次数；

步骤3.3，对于第i次保障任务，利用下式计算保障任务概率P1_i；

其中，Tw_i为第i次保障任务时该单元的计划工作时间，Nzj_i为该单元在装备内的装机数量，Nbj_i为任务开始前的备件数量；F_i表示保障任务结果，保障任务成功则F_i＝1，保障任务失败则F_i＝0；

步骤3.4，更新i＝i+1，若i≤k则跳转至步骤3.3，否则跳转至步骤3.5；

步骤3.5，令

L_j反映了估计值μ_j对真值μ的似然程度；

步骤3.6，令j＝j+1，若j≤n则跳转至步骤3.2，否则跳转至步骤4。

本发明的有益效果是：本发明提出的寿命分布参数估计方法，能在因缺少专门的可靠性试验导致缺少寿命值数据时，有效利用实际工作产生的备件保障数据，较为准确地估计寿命分布参数。

具体实施方式

以下结合实施例对本发明的原理和特征进行描述，所举实例只用于解释本发明，并非用于限定本发明的范围。

产品由各种单元组成。在本文中，指数型单元是指寿命服从指数分布的单元。一般来说，电子零部件等单元的寿命服从指数分布，如：印制电路板插件、电子部件、电阻、电容、集成电路等。寿命X服从指数分布记作X～Exp(μ)，μ为指数分布的均值参数，X的密度函数为

约定：共有k次保障任务的相关数据，第i次保障任务时该单元的计划工作时间记为Tw_i(以下简称为任务时间)，该单元在装备内的装机数量记为Nzj_i，任务开始前的备件数量记为Nbj_i，保障任务结果记为F_i，1≤i≤k。在保障任务期间，Nzj_i个该型单元之间视为串联关系，以备件更换故障件的方式来排除故障，若发生的故障次数不大于备件数量，则本次保障任务成功(F_i记为1)，否则视为失败(F_i记为0)。即：以形如[Tw_i Nzj_i Nbj_i F_i]的备件保障数据来描述一次保障任务。在积累了k组备件保障数据后，具体估计寿命分布参数的方法如下：

1、确定候选的寿命分布参数

生成n个候选的分布参数μ_j,1≤j≤n，其中，μ_j表示指数分布的均值参数，n为正整数。生成候选的分布参数的具体实现方法为：

参数μ_j＝μ_min+(j-1)d1,1≤j≤n，其中，

μ_max表示指数分布的均值参数上限，μ_min表示指数分布的均值参数下限，n为正整数，且n≥2。n决定了对均值参数的估计精度。

2、初始化似然值

初始化似然值L_j,1≤j≤n，令L_j＝0。

3、遍历计算似然值

3.1令候选参数序号j＝1

3.2令任务序号i＝1

3.3对于第i次保障任务，计算保障任务概率P1_i：

3.4更新i＝i+1，若i≤k则转3.3，否则转3.5。

3.5令

L_j反映了估计值μ_j对真值μ的似然程度。

3.6令j＝j+1，若j≤n则转3.2，否则转4。

4、输出结果

在L_j(1≤j≤n)中找到最大值，记其为L_M，则μ_M为电子件寿命均值参数的估计值。

实施例

某电子件的[Tw_i Nzj_i Nbj_i F_i]型备件保障数据如表1所示，从工程经验可知该电子件的寿命服从指数分布，试估计其均值参数。

表1某电子件的备件保障数据

计算过程如下：

1、确定候选的寿命分布参数

从以往经验得知，该电子件的均值参数在100～1900范围内，以200为步长，生成10个候选的分布参数μ_j,1≤j≤10。

2、初始化似然值

初始化似然值L_j,1≤j≤10，令L_j＝0。

3、遍历计算似然值

3.1令j＝1

3.2令i＝1

3.3对于第i次保障任务，计算保障任务概率P1_i：

3.4更新i＝i+1，若i≤15则转3.3，否则转3.5。

表2是j＝1～10时，15次保障任务概率结果。

表2保障任务概率P1_i

3.5令

L_j反映了估计值μ_j对真值μ的似然程度。候选分布参数μ_j及对应的似然值L_j计算结果如表3。

表3似然值

3.6令j＝j+1，若j≤10则转3.2，否则转4。

4、输出结果

从表3可知，μ₆＝1100对应的似然值最大，因此将其作为该单元的寿命分布参数估计值。

可建立以下仿真模型模拟备件保障过程：

1)产生Nzj个随机数simT_i,1≤i≤Nzj用于模拟装备中Nzj个单元的寿命，simT_i服从该类单元的寿命分布规律；

2)初始化当前备件数量N1＝Nbj，该类单元的累积工作时间simTw＝0；

3)把simT_i,1≤i≤Nzj按照从小到大重新排序；

4)更新simTw＝simTw+simT₁，simTw为故障发生时刻，然后遍历更新simT_i＝simT_i-simT₁,1≤i≤Nzj；

5)若N1＞0，则按照该类单元的寿命分布规律产生随机数t1，并令simT₁＝t1用于模拟换件维修，更新N1＝N1-1后，转3)；

若N1≤0，则转6)；

6)若simTw≥Tw，则本次保障任务成功，令F＝1；否则，本次保障任务失败，令F＝0。

运行上述仿真模型，可以得到大量形如[Tw_i Nzj_i Nbj_i F_i]的备件保障数据，对本发明的方法进行仿真验证。以μ＝1000,k＝15，装机数1～5为例进行大量仿真。如果基于备件保障数据，用本文方法得到的该单元的寿命分布均值参数统计结果的均值为1017.2、根方差为105.2。如果基于寿命数据simT_i，采用理论方法计算得到的均值参数统计结果的均值为1003.1、根方差为53.4。二者的差异在工程允许范围以内。

以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (2)

1.一种基于备件保障数据的指数型单元寿命分布参数估计方法，其特征在于，包括以下步骤：

步骤1，根据保障任务的相关数据生成n组候选分布参数μ_j,1≤j≤n，其中，μ_j表示指数分布的均值参数，n为正整数；

步骤2，初始化似然值L_j,1≤j≤n，令L_j＝0；

步骤3，根据k次保障任务的相关数据计算保障任务概率P1，更新似然值L_j；

步骤4，在更新后的似然值L_j中找到最大似然值，记为L_M，则似然值L_M对应的μ_M即为指数分布的均值参数的估计结果；

所述步骤3包括：

步骤3.1，令候选参数序号j＝1；

步骤3.2，令任务序号i＝1，1≤i≤k，k为保障任务最大次数；

步骤3.3，对于第i次保障任务，利用下式计算保障任务概率P1_i；

其中，Tw_i为第i次保障任务时该单元的计划工作时间，Nzj_i为该单元在装备内的装机数量，Nbj_i为任务开始前的备件数量；F_i表示保障任务结果，保障任务成功则F_i＝1，保障任务失败则F_i＝0；

步骤3.4，更新i＝i+1，若i≤k则跳转至步骤3.3，否则跳转至步骤3.5；

步骤3.5，令

L_j反映了估计值μ_j对真值μ的似然程度；

步骤3.6，令j＝j+1，若j≤n则跳转至步骤3.2，否则跳转至步骤4。

2.根据权利要求1所述的方法，其特征在于，步骤1中所述的候选分布参数的生成方法包括：

参数μ_j＝μ_min+(j-1)d1,1≤j≤n，其中，

μ_max表示指数分布的均值参数上限，μ_min表示指数分布的均值参数下限，n为正整数，且n≥2。