[发明专利]一种估计伽玛型单元可靠性参数的方法

摘要

本发明公开了一种估计伽玛型单元可靠性参数的方法，包括以下步骤：步骤一：确定候选的寿命分布参数，根据已有伽玛型单元可靠性数据分布规律，确定待估伽玛型单元的寿命形状参数上限αmax和下限αmin及尺度参数上限bmax和下限bmin，根据上下限值分别确定形状参数的步长d1和尺度参数的步长d2，计算出α1j1和b1j2，然后对α1j1和b1j2进行遍历组合；步骤二：遍历寿命形状参数为αj和尺度参数为bj的组合计算似然值。寻找最大似然值记为LM，则该最大值对应的αM为伽玛型单元寿命形状参数的估计值，bM为寿命尺度参数的估计值。本发明利用“少量的可靠性试验数据+在产品研制、生产、使用等阶段产生的大量数据”，估算产品寿命的分布规律。

主权项

1.一种估计伽玛型单元可靠性参数的方法，其特征在于，针对寿命服从伽玛分布的单元开展可靠性参数估计，包括以下步骤：步骤一：确定候选的寿命分布参数，根据已有伽玛型单元可靠性数据分布规律，初步确定待估伽玛型单元寿命伽玛分布的形状参数上限αmax和形状参数下限αmin及待估伽玛型单元寿命伽玛分布的尺度参数上限bmax和尺度参数下限bmin；在确定的伽玛型单元伽玛形状参数区间内，等间隔生成n1个候选形状参数，候选参数中每个相邻参数之间的步长相等为d1，根据伽玛分布形状参数的步长依次计算出n1个伽玛分布的形状参数α1j1，其中1≤j1≤n1；在确定的伽玛型单元伽玛尺度参数区间内，等间隔生成n2个候选尺度参数，候选参数中每个相邻参数之间的步长相等为d2，根据伽玛分布尺度参数的步长依次计算出n2个伽玛分布的尺度参数b1j2，其中1≤j2≤n2；取n＝n1×n2，由α1j1和b1j2进行遍历组合获得n组候选的分布参数(αj,bj),1≤j≤n；步骤二：遍历寿命形状参数为αj和尺度参数为bj的组合计算似然值，对于每个候选参数组合，针对一组包含m个伽玛型单元检测信息的数据组，根据第i个检测信息所包含的在Ti时刻的单元状态信息Fi确定其对应的计算系数Wi，依据m个检测数据不断迭代更新似然值Lj，迭代之初，设定每个候选参数组合对应的似然值初始值为0，在每个候选参数组合所对应的迭代完毕后的似然值中，寻找最大值即LM，则该最大值对应的αM为伽玛型单元形状参数的估计值，bM为尺度参数的估计值。