[发明专利]一种动态环境下的可靠性预计方法及系统

权利要求书

1.一种动态环境下的可靠性预计方法，其特征在于，所述方法包括以下步骤：

步骤1，利用GJB299C-2006或者MIL-HDBK-217方法获取静态环境下的可靠性指标；

步骤2，开展马尔可夫过程的模拟，获取马尔可夫过程的一条轨道；

步骤3，开展马尔可夫过程参数的模拟，获取马尔可夫过程的转移率矩阵；

步骤4，开展基于马尔可夫过程的动态环境建模与仿真，获取一次产品失效时间；

步骤5，反复使用基于马尔可夫过程的动态环境建模与仿真算法若干次，根据所得结果获得产品的故障率即MTBF值。

2.根据权利要求1所述的动态环境下的可靠性预计方法，其特征在于，步骤2所述开展马尔可夫过程的模拟，获取马尔可夫过程的一条轨道，具体包括：

对于马氏过程X＝{X(t),t≥0}，当给定X(t_n)＝i_n时，停留在i_n状态的时间{θ_n,n≥1}条件独立，且θ_n服从参数为的指数分布，而状态转移概率

初始分布为{P_i(0),i∈S}即系统在0时刻处于状态i，转移率矩阵Q＝{q_ij}的马氏链X＝{X(t),t≥0}轨道的模拟可分为以下三部分：首先进行初始分布P_i(0)的模拟；然后模拟停留时间{θ_n,n≥1}；最后对状态转移进行模拟；

具体模拟步骤如下：

步骤一：取两串相互独立的同为(0,1)上的均匀分布随机变量序列{U_n,n≥0}及{V_n,n≥1}；

步骤二：若U₀满足：

则取X(0)＝i₀作为初始状态；

步骤三：取出U₁，令作为X(t)停留在X(0)＝i₀的时间；

步骤四：取出V₁，若i₁∈S使V₁满足：

则取X(t₁)＝X(θ₁)＝i₁；

步骤五：取出U₂，令作为X(t)停留在X(t₁)＝i₁状态的停留时间，再取t₂＝θ₁+θ₂，作为第二次跳跃时刻；

步骤六：取出V₂，若i₂∈S使V₂满足：

则取X(t₂)＝i₂；

......继续以上步骤

步骤七：取出U_n，令作为X(t)停留在X(t_n-1)＝i_n-1状态的停留时间，再取

步骤八：取出V_n，若i_n∈S使V_n满足

则取X(t_n)＝i_n；

如此继续重复如上步骤，即可模拟X的一条轨道。