[发明专利]基于统计量排序的和联系统可靠度精确置信下限估计方法

权利要求书

1.一种基于统计量排序的和联系统可靠度精确置信下限估计方法，其特征在于：具体的实施步骤如下：

步骤一：确定样本空间

若一和联系统由N个支任务单元与1个中心任务单元C构成，其中，支任务单元i(i＝1,2,…,N)完成任务的比重为β_i(i＝1,2,…,N)，

假设支任务单元i(i＝1,2,…,N)的寿命服从指数分布，其失效分布函数为F_i(t)＝1-exp(-t/θ_i)，中心任务单元C的寿命也服从指数分布，其失效分布函数为F_C(t)＝1-exp(-t/θ_C)；

现针对和联系统的支任务单元开展随机截尾试验，支任务单元i(i＝1,2,…,N)的试验结果为用(t_i,y_i)(i＝1,2,…,N)简化表示，其中，t_ij表示第i(i＝1,2,…,N)个支任务单元第j(j＝1,2,…,M_i)次随机截尾试验的截尾时间，y_ij表示第i(i＝1,2,…,N)个支任务单元第j(j＝1,2,…,M_i)次随机截尾试验的结果，当y_ij＝0时，表示到截尾时间t_ij时，故障已经发生；当y_ij＝1时，表示到截尾时间t_ij时，故障没有发生；同时，针对系统的中心任务单元有随机截尾试验的结果用(t_C,y_C)简化表示，其中，t_Ck表示中心任务单元第k(k＝1,2,…,M_C)次随机截尾试验的截尾时间，y_Ck表示中心任务单元第k(k＝1,2,…,M_C)次随机截尾试验的结果，当y_Ck＝0时，表示到截尾时间t_Ck时，故障已经发生；当y_Ck＝1时，表示到截尾时间t_Ck时，故障没有发生；

那么，(t_i,y_i)(i＝1,2,…,N)与(t_C,y_C)构成了和联系统所有单元随机截尾试验的全部试验结果，令(T_i,Y_i)＝(t_i,y_i)(i＝1,2,…,N)，(T_N+1,Y_N+1)＝(t_C,y_C)，用(T,Y)简化表示系统所有组成单元随机截尾试验的全部试验结果，称之为样本空间，该空间中包含种试验情况；

步骤二：统计量排序

在试验样本空间的基础上，对样本空间中的所有试验结果按照“好坏”进行排序，把试验结果“好”的情况排在试验结果“坏”的情况之前，选取系统可靠度点估计这一统计量作为衡量结果“好坏”的标准；那么，基于系统可靠度点估计得到试验样本空间排序的主要计算过程分为单元产品寿命分布参数估计、系统可靠度计算和统计量排序三部分：

1)单元产品寿命分布参数估计：

对于第i(i＝1,2,…,N)个支任务单元的某一确定试验结果其产品寿命分布参数估计方法如下：

(1)若试验结果中包含失效样本，即存在则采用极大似然估计方法给出单元级产品寿命分布参数估计，构造

使得L_i取得最大值情况的参数值即为指数分布失效率参数θ_i的极大似然估计值；式中：t_ij表示第i(i＝1,2,…,N)个支任务单元第j(j＝1,2,…,M_i)次随机截尾试验的截尾时间，y_ij表示第i(i＝1,2,…,N)个支任务单元第j(j＝1,2,…,M_i)次随机截尾试验的结果，当y_ij＝0时，表示到截尾时间t_ij时，故障已经发生；当y_ij＝1时，表示到截尾时间t_ij时，故障没有发生；

(2)若试验样本均为截尾样本，即则使用在置信水平0.5下的最优置信下限估计，来代替指数分布失效率参数的点估计，即

式中：t_ij表示第i(i＝1,2,…,N)个支任务单元第j(j＝1,2,…,M_i)次随机截尾试验的截尾时间；

同理，对系统所有支任务单元及中心单元寿命分布的参数进行估计，从而得到与

2)系统可靠度计算

根据和联系统可靠度的定义知，系统可靠度

式中：R_c(t)为中心单元的可靠度点估计，R_i(t)为第i个支任务单元的可靠度点估计，和为上一步骤中得到的单元产品寿命分布参数估计值；

因此，由上述公式能计算样本空间中所有试验结果下系统可靠性的点估计R_S(t)；

3)统计量排序

将系统可靠度R_S(t)这一统计量按照从大到小的顺序进行排列，同时将其对应的试验情况进行排序，将系统可靠度最高的一组试验结果排在第一位，记作(T,Y)₍₁₎，其中，下标“(·)”中的数字表示排序后的序号，从而，对样本空间中的全部样本进行排序，得到对应的序号ω；

步骤三：计算系统可靠度的精确置信下限

现在有一组关于和联系统各组成单元的随机截尾试验，其试验结果为为根据步骤二的结果找到在样本空间中的排序，那么在该组试验条件下，置信水平1-α下和联系统可靠度的精确置信下限的解算过程如下：

其中，

式中：t_ij表示第i(i＝1,2,…,N)个支任务单元第j(j＝1,2,…,M_i)次随机截尾试验的截尾时间，t_Cj表示中心单元第j(j＝1,2,…,M_C)次随机截尾试验的截尾时间，y_ij表示第i(i＝1,2,…,N)个支任务单元第j(j＝1,2,…,M_i)次随机截尾试验的结果，y_Cj表示中心单元第j(j＝1,2,…,M_C)次随机截尾试验的结果，结果为0表示到截尾时间时故障已经发生；结果为1表示到截尾时间时故障没有发生；

采用内点法求解这一非线性不等式约束下的优化问题，得到置信水平1-α下和联系统可靠度的精确置信下限为R_L(t)。