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

权利要求书

1.一种估计对数正态型单元可靠性参数的方法，其特征在于，针对寿命服从对数正态分布的单元开展可靠性参数估计，包括以下步骤：

步骤一：确定候选的寿命分布参数，根据已有对数正态型单元可靠性数据的分布规律，初步确定待估对数正态分布单元的寿命对数均值参数下限值μ_min和对数均值参数上限值μ_max及待估对数正态分布的对数标准差参数下限值σ_min和对数标准差参数上限值σ_max；

在确定的对数正态分布单元的寿命对数均值参数区间内，等间隔生成n1个候选对数正态分布的对数均值参数，候选参数中，相邻对数均值参数的步长相等为d1，根据对数均值参数的步长依次计算出n1个对数正态分布的对数均值参数μ1_j1，其中1≤j1≤n1；

在确定的对数正态分布单元的寿命对数标准差参数区间内，等间隔生成n2个候选对数正态分布的对数标准差参数，其中，相邻对数标准差参数的步长相等为d2，根据对数标准差参数的步长依次计算出n2个对数正态分布的对数标准差参数σ1_j2，其中1≤j2≤n2；

取n＝n1×n2，对μ1_j1和σ1_j2进行遍历组合获得n组候选的分布参数(μ_j,σ_j)，其中1≤j≤n；

步骤二：遍历寿命对数均值参数μ_j和对数标准差参数σ_j的组合计算似然值，对于每个候选参数组合，针对一组包含m个对数正态型单元检测信息的数据组，根据第i个检测信息所包含的在T_i时刻的单元状态信息F_i确定其对应的计算系数W_i，依据m个检测数据不断迭代更新似然值L_j，迭代之初，设定每个候选参数组合对应的似然值初始值为0，在每个候选参数组合所对应的迭代完毕后的似然值中，寻找最大值即L_M，则最大似然值L_M对应的μ_M为对数均值参数的估计值，σ_M为对数标准差参数的估计值。

2.根据权利要求1所述的对数正态型单元可靠性参数的方法，其特征在于：所述步骤一中生成候选分布参数的具体计算过程如下：

(1)确定对数正态分布的对数均值参数μ1_j1和步长d1的具体计算过程如下：

其中，μ_max表示对数正态分布的对数均值参数上限，μ_min表示对数正态分布的对数均值参数下限，n1为正整数，且n1≥2；

(2)确定对数正态分布的对数标准差参数σ1_j2和步长d2的具体计算过程如下：

其中，σ_max表示对数正态分布的对数标准差参数上限，σ_min表示对数正态分布的对数标准差参数下限，n2为正整数，且n2≥2；

(3)μ1_j1和σ1_j2遍历组合的计算方式为：

令j＝1；第一层循环中遍历j1＝1:n1，第二层循环中遍历j2＝1:n2，令μ_j＝μ1_j1，σ_j＝σ1_j2，j＝j+1；其中，μ_max≥μ1_j1≥μ_min，σ_max≥σ1_j2≥σ_min。

3.根据权利要求1所述的对数正态型单元可靠性参数的方法，其特征在于：所述步骤二中计算W_i和似然值L_j的计算公式如下：

其中，W_i为计算系数，log(*)为自然对数函数，σ_j为对数正态分布的对数标准差参数，μ_j为对数正态分布的对数均值参数。

4.根据权利要求1或3所述的一种估计对数正态型单元可靠性参数的方法，其特征在于：所述步骤二中似然值L_j的遍历计算过程如下：

(1)令j＝1；

(2)令i＝1，L_j＝0；

(3)计算系数

其中，μ_j为对数正态分布的均值参数，σ_j为对数正态分布的对数标准差参数，μ_j为对数正态分布的对数均值参数，L_j为似然值，T_i为第i个产品的检测时刻；

(4)更新i＝i+1，若i≤m则转(3)，否则转(5)；

(5)令j＝j+1，若j≤n则转(2)，否则(6)；

(6)在L_j中找到最大值，其中j的取值范围为1≤j≤n，记其为L_M，则μ_M为对数均值参数的估计值，σ_M为对数标准差参数的估计值。