融合产品成败型数据和故障时间数据的可靠性评估方法
【技术领域】
[0001 ] 本发明是一种利用贝叶斯理论，针对产品的可靠性数据中既存在其成功/失败数 据又存在其故障时间数据的情况，对其可靠性进行融合评估，以期达到利用现有的数据来 获取高精度可靠性评估结果的效果，属于可靠性评估技术领域。可用于复杂运行系统的可 靠度等指标的评估、考核。
【背景技术】
[0002] 在产品可靠性指标的评估过程中，产品故障时间数据和成功/失败数据都是其主 要的评估依据，但是这两种数据只能分别采用以进行评估，即仅使用故障时间数据或仅使 用成功/失败数据。一般情况下，当获得充足的故障时间数据或充足的成功/失败数据时， 可以依据相应的评估方法评估得到产品的可靠度。但当故障时间数据与成功/失败数据 都不充足时，评估工作便难以开展。在某些情况下，在产品的试验与运行过程中，在不同的 时段分别记录了其故障时间数据与成功/失败数据，虽然其中的单独一类数据难以开展评 估工作，但两类数据的总量却是充足而满足评估数据量要求的；如果将上述数据弃之不用， 无疑是一种损失，这就需要一种能够融合产品成败型数据和故障时间数据的可靠性评估方 法，来解决上述问题，获得满足精度要求的可靠性评估结果。
[0003] 在基于多源数据融合的可靠性评估方面，国内外已开展了一定的研宄，主要是同 类型数据间的融合和故障时间数据与退化数据间的融合，在故障时间数据与成功/失败数 据间的融合方面还鲜有研宄。因此提出融合产品成败型数据和故障时间数据的可靠性评估 方法是具有一定独创性的。
【发明内容】
[0004] 本发明的目的是为了解决上述问题，提出一种能够有效的融合产品成败型数据和 故障时间数据的可靠性评估方法，基于成功/失败数据和故障时间数据来获得具有较高精 度的可靠性评估结果。
[0005] 本发明的具体步骤为：
[0006] 步骤一、搜集成功/失败数据与故障时间数据；
[0007] 步骤二、构建融合模型；
[0008] 步骤三、构建贝叶斯评估模型；
[0009] 步骤四、评估可靠性指标；
[0010] 本发明的优点在于：
[0011] (1)本发明能够融合成败型数据与故障时间数据，评估得到产品的可靠性指标；
[0012] (2)本发明建立了成败型数据与故障时间数据间的联系，从全局的角度进行求解， 对不同类型数据间的关系进行了量化；
[0013] (3)本发明能够在产品数据稀缺的情况下，融合各方面信息，提高产品可靠度评估 精度。
【附图说明】
[0014] 图1是本发明的流程图；
【具体实施方式】
[0015] 下面将结合附图和实施例对本发明作进一步的详细说明。
[0016] 方法的流程图如图1所示，包括以下几个步骤：。
[0017] 步骤一、搜集成功/失败数据与故障时间数据
[0018] 对产品相关的成功/失败数据与故障时间数据进行搜集。
[0019] 若产品在n个样本中成功了 r次，依据成功发生的顺序将其表示为伯努利过程中 的数据序列Xi，其中i = 1…n，Xi等于0或1，若第i次失败则Xi= 0,若第i次成功则Xi =1，Xi= 1的次数为r。如某产品在5次试验中成功了 4次失败了 1次，如果是在第3次 试验时失败的话，那么成败/失败数据的伯努利过程表达为[1，1，0, 1，1]。
[0020] 同时搜集产品的故障时间数据I，j = 1…m，m为故障次数。
[0021] 步骤二、构建融合模型
[0022] (1).数据的概率模型及其关系
[0023] 对于成功/失败数据而言，每一次任务成功或失败的概率可以表示为：
[0024] P(Xi) =pxA! (1)
[0025] 其中p为任务成功的概率。
[0026] 对于故障时间数据而言，本专利假设产品故障时间服从指数分布，其可靠度可表 示为：
[0027] R(t) = e_At (2)
[0028] 其概率密度函数为：
[0029] f(t) = (3)
[0030] 若每次任务的平均时间为h，那么任务成功的概率可表示为：
[0031] p= ^eAt〇 (4)
[0032] (2).数据融合模型
[0033] 在上述内容的基础上，构建融合成功/失败数据与故障时间数据的模型。假设一 个模型的对数函数为W k= l〇gf(Zk| 0)，可利用伯努利分布将其似然函数表示为：
[0035] 那么故障时间的对数函数可表示为％=1〇琪1//〇，利用伯努利分布可将其似然 函数表示为：
[0037] 成功/失败数据服从伯努利分布，其似然函数可以表示为：
[0039] 定义cs为状态参数（当数据为成功/失败数据时，c s= 0。当数据为故障时间数 据时，cs= 1)，那么定义：
[0040] pBs =c\ -c'' +(l-c;)/; (8)
[0041] 因此实现了以伯努利分布为媒介，将成功/失败数据与故障时间数据融合在一个 模型之中，即融合模型：
[0043] 其中当数据为成功/失败数据时us=Xs，当数据为故障时间数据时u s= 1。
[0044] 步骤三、构建贝叶斯评估模型
[0045] 首先确定贝叶斯模型的总体分布为（9)，其次确定模型中参数A的先验分布，这 里选择伽玛分布作为X的先验分布，即：
[0046] 入~Gamma (a, b) (10)
[0047] 其中a，b为先验分布中的超参数。
[0048] 那么，可知未知参数的后验分布为：
[0050] 进而确定（11)为贝叶斯评估模型。
[0051] 步骤四、评估可靠性指标
[0052] 利用马尔可夫链蒙特卡罗（Markov Chain Monte Carlo，MCMC)方法对贝叶斯模型 (11)进行抽样模拟，获得未知变量的后验分布及未知变量的评估值。即获得参数A的评估 值i，从而得到产品在时刻t的可靠度评估值：
[0053] R(t) =e^h
[0054] 实施例：
[0055] 已知，某产品在试验与使用过程中的不同阶段，分别搜集得到成功/失败数据和 故障时间数据，现利用本专利提出的方法开展相应的融合评估，以获取产品的可靠度评估 值。
[0056] 步骤一、搜集成功/失败数据与故障时间数据
[0057] 搜集产品成功/失败数据共50组，其中成功45次，失败5次，将其表示为伯努利 过程中的数据序列为Xi，其中i = 1-50。
[0058] 已知每次任务的平均工作时间h为10小时。
[0059] 搜集得到产品的故障时间数据 20 组，T = [117,8,84,169,10, 2,82, 220,135,89， 52,134, 51，34,151，214,122,114,86,68](小时），表示为 T」，j = 1…20。
[0060] 步骤二、构建融合模型
[0061] 故障时间的对数函数可表示为uV/ =l〇g(Ae <)，依据公式（6)将其似然函数表示 为：
[0063] 依据公式（7)将成功/失败数据的似然函数表示为：
[0065] 定义cs为状态参数（当数据为成功/失败数据时，c s= 0。当数据为故障时间数 据时，cs= 1)，那么：
[0066] pBs =cs-eWfi +(l-cjp
[0067] 那么依据公式（9)获得融合模型：
[0069] 其中当数据为成功/失败数据时us=Xs，当数据为故障时间数据时u s= 1。
[0070] 步骤三、构建贝叶斯评估模型
[0071] 将（12)确定为总体分布，将分布参数A作为贝叶斯模型中的参数。同时选择 Gamma分布作为A的先验分布，将其具体形式设为：
[0072] 人~Gamma (1. 3, 1. 3)
[0073] 从而构建贝叶斯模型为：
[0075] 步骤四、评估可靠性指标
[0076] 利用马尔可夫链蒙特卡罗（Markov Chain Monte Carlo，MCMC)方法对贝叶斯模型 进行抽样模拟，获得A的后验分布及未知变量的评估值，如表1所示。
[0077] 表1参数评估值
[0078]
[0079] 从而得到产品工作20小时的可靠度为0. 819。该结果同产品的实际情况相符。
【主权项】
1.融合产品成败型数据和故障时间数据的可靠性评估方法，其特征在于，包括w下几 个步骤： 步骤一、捜集成功/失败数据与故障时间数据 对产品相关的成功/失败数据与故障时间数据进行捜集。 若产品在n个样本中成功了r次，依据成功发生的顺序将其表示为伯努利过程中的数 据序列Xi，其中i=1…n，Xi等于0或1，若第i次失败则Xi=0,若第i次成功则Xi=1， Xi=1的次数为r。 同时捜集产品的故障时间数据Tj.，j=1…m，m为故障次数。 步骤二、构建融合模型 (1) .数据的概率模型及其关系 对于成功/失败数据而言，每一次任务成功或失败的概率可W表示为：(1) 其中P为任务成功的概率。 对于故障时间数据而言，本专利假设产品故障时间服从指数分布，其可靠度可表示 为： R(t)=e_ 人t (2) 其概率密度函数为： f(t) =Ae_\t (3) 若每次任务的平均时间为t。，那么任务成功的概率可表示为：(4) (2) .数据融合模型 在上述内容的基础上，构建融合成功/失败数据与故障时间数据的模型。假设一个模 型的对数函数为Wk=logf(zj0)，可利用伯努利分布将其似然函数表示为：那么故障时间的对数函数可表示为=bg(Ac'勺，利用伯努利分布可将其似然函数 表示为：成功/失败数据服从伯努利分布，其似然函数可W表示为：(7) 定义C,为状态参数（当数据为成功/失败数据时，C,=0。当数据为故障时间数据 时，c,=1)，那么定义；(8) 因此实现了W伯努利分布为媒介，将成功/失败数据与故障时间数据融合在一个模型 (9) 之中，即融合模型： 其中当数据为成功/失败数据时u,=X,，当数据为故障时间数据时u,= 1。 步骤=、构建贝叶斯评估模型 首先确定贝叶斯模型的总体分布为巧），其次确定模型中参数A的先验分布，该里选 择伽玛分布作为A的先验分布，即； 入~Gamma(a,b) (10) 其中a,b为先验分布中的超参数。 那么，可知未知参数的后验分布为：进而确定（11)为贝叶斯评估模型。 步骤四、评估可靠性指标 利用马尔可夫链蒙特卡罗（Markov化ainMonteCarlo，MCMC)方法对贝叶斯模型（11) 进行抽样模拟，获得未知变量的后验分布及未知变量的评估值。即获得参数A的评估值i， 从而得到产品可靠度的评估值。2.根据权利要求1所述的融合产品成败型数据和故障时间数据的可靠性评估方法，其 特征在于，所述的步骤四具体为： 产品在时刻t的可靠度评估值为：
【专利摘要】本发明公开了一种融合产品成败型数据和故障时间数据的可靠性评估方法，包括以下几个步骤：步骤一、搜集成功/失败数据与故障时间数据；步骤二、构建融合模型；步骤三、构建贝叶斯评估模型；步骤四、评估可靠性指标；本发明能够融合成败型数据与故障时间数据，评估得到产品的可靠性指标；建立了成败型数据与故障时间数据间的联系，从全局的角度进行求解，对不同类型数据间的关系进行了量化；在产品数据稀缺的情况下，融合各方面信息，提高产品可靠度评估精度。
【IPC分类】G06F19/00
【公开号】CN104933323
【申请号】CN201510406033
【发明人】王立志, 王晓红, 范文慧, 宣金泉
【申请人】北京航空航天大学
【公开日】2015年9月23日
【申请日】2015年7月10日