产品可靠性综合评估方法与系统、设备、介质与流程

1.本发明涉及数据分析与评估技术领域，具体地，涉及一种产品可靠性综合评估方法与系统、设备、介质，尤其涉及一种融合多阶段多类型的产品可靠性综合评估方法与系统、设备、介质。


背景技术：

2.根据经典的可靠性理论，评估产品在实际使用环境下的可靠性水平，须依据产品在现场环境下采集到的数据进行统计推断，这意味着我们需要大量的现场数据。然而，很多产品并不具备大量现场数据，因此，充分利用产品寿命周期各阶段收集到的各类型可靠性信息并加以融合评估，是可靠性分析的一个重要的研究方向。
3.bayes理论是解决小样本统计推断的有效手段，可以融合多种数据进行可靠性评估。目前，国内外文献在利用bayes理论进行可靠度评估时，往往只是对产品的两种或少量几种数据进行融合，例如加速寿命试验数据与现场故障数据的融合、专家信息与现场成败型数据的融合等。然而，实际工程中存在着大量不同阶段不同类型的可靠性数据，这些数据从不同侧面不同角度反映了产品的可靠性水平。如何将现场信息进行融合也是亟需解决的问题。
4.经过检索，专利文献cn110874502a公开了一种基于多阶段试验数据折合的航天产品可靠性评估方法，包括如下步骤：进行航天产品的两阶段可靠性试验，获取航天产品第一阶段以及第二阶段的产品试验寿命数据；对第一阶段的产品试验寿命数据进行折合；确定第二阶段weibull寿命分布参数的验前分布，并得到联合分布数据下的bayes验后分布核函数；采用基于bayes的m
‑
h抽样方法进行可靠度的区间估计。该现有技术的不足之处在于需要有现场数据的情况下才能得到可靠性数据的分析，虽然也用到了bayes理论进行可靠度的评估，但是仅针对寿命数据进行分析，并不全面。
5.因此，亟需设计一套简便有效的处理方法在尽可能保留多阶段多类型信息所包含内容的前提下，将多阶段多类型可靠性数据进行融合，给出产品可靠性水平的精确估计。


技术实现要素：

6.针对现有技术中的缺陷，本发明的目的是提供一种产品可靠性综合评估方法与系统、设备、介质，是在缺乏足够的现场数据的情况下，充分利用产品在设计、制造、使用过程中各阶段产生的各类可靠性数据，在确保数据一致性的前提下，对各阶段各类型可靠性数据进行分类、预处理、融合和评估，以提高可靠性评估的精度。解决了在尽可能保留多阶段多类型信息所包含内容的前提下，将多阶段多类型可靠性数据进行融合，给出产品可靠性水平的精确估计。
7.根据本发明提供的一种产品可靠性综合评估方法，包括如下步骤：
8.步骤s1：将多阶段多类型数据进行分类与预处理；
9.步骤s2：将分类与预处理之后的多阶段多类型数据进行一致性检验；
10.步骤s3：一致性检验后对多阶段多类型数据进行折算和综合；
11.步骤s4：多阶段多类型数据进行折算和综合之后根据数据类型确定先验分布和超参数；
12.步骤s5：确定先验分布和超参数之后，根据bayes理论及相关模型，进行数据融合和可靠性评估。
13.优选地，步骤s1包括将多阶段多类型数据进行分类，数据类型如下：
14.产品现场类数据：产品实际使用或服役环境条件下的寿命型、成败型和退化型数据；
15.产品异环境类数据：产品在寿命周期各阶段开展的不同试验条件下不同于现场使用环境条件下的寿命型、成败型和退化型数据；
16.相似产品数据：又称为异总体数据，包括相似产品现场类数据和相似产品异环境类数据，相似产品现场类数据是指相似产品实际使用或服役环境条件下的寿命型、成败型和退化型数据，相似产品异环境类数据是指相似产品在寿命周期各阶段开展的不同试验条件下不同于现场使用环境条件下的寿命型、成败型和退化型数据；
17.部件级数据：产品各组成部件的现场类数据和异环境类数据，部件级现场类数据是指产品各组成部件实际使用或服役环境条件下的寿命型、成败型和退化型数据，部件级异环境类数据是指产品各组成部件在寿命周期各阶段开展的不同试验条件下不同于现场使用环境条件下的寿命型、成败型和退化型数据；
18.专家型数据：利用专家对产品可靠度的经验估计值来表示的数据。
19.优选地，步骤s1包括对多阶段多类型数据进行预处理：
20.针对寿命型数据的预处理：利用寿命分布类型检验和参数评估方法进行预处理，得到寿命类型及其参数估计；
21.针对异环境类数据的预处理：利用环境因子估算方法进行预处理，给出异环境类数据相对于现场数据的“环境因子”；
22.针对部件级数据的预处理：利用系统可靠性综合方法进行预处理，将部件级数据折合成产品数据；
23.针对退化型数据的预处理：利用退化建模方法进行预处理，将退化型数据转换成伪寿命型数据；
24.针对相似产品数据的预处理：利用相似因子估算方法进行预处理，给出相似产品数据的“相似因子”。
25.优选地，步骤s2包括：
26.步骤s2.1：对寿命分布类型进行一致性检验，对于寿命型数据，根据寿命类型，判断多阶段多类型寿命数据与产品现场类数据是否服从同一种寿命分布；
27.步骤s2.2：对失效机理进行一致性检验，判断异环境寿命数据是否满足失效机理不变条件，即环境因子不变原则。
28.优选地，步骤s3中对多阶段多类型数据的折算和综合包括：
29.对异环境数据的折算和综合：对于环境因子k
i0
≠1的异环境数据，根据数据类型将异环境数据折算成现场数据，再将折算后的现场数据进行综合；
30.对相似产品数据的折算和综合：对于相似因子f
i0
≠1的相似产品数据，根据数据类
型将相似产品数据折算成被评估产品数据，再将折算后的被评估产品数据进行综合。
31.优选地，步骤s4包括：
32.步骤s4.1：根据总体分布中未知参数的取值范围选取先验分布，基本方法如下：
33.对于[0,1]取值的参数，选用贝塔分布be(a,b)作为先验分布，
[0034]
对于(
‑
∞, ∞)取值的参数，选用正态分布n(μ,σ2)作为先验分布，
[0035]
对于(0, ∞)取值的参数，选用伽马分布ga(z,η)或倒伽马分布iga(z,η)作为先验分布；
[0036]
步骤s4.2：根据验前数据，确定验前分布的超参数，基本方法为利用验前信息，产生自助样本，得到各样本对应分布的参数的极大似然估计值样本，计算出参数的均值和方差，该均值和方差应等于对应参数的验前分布的均值和方差，联立方程组，得到寿命分布参数的验前分布的超参数。
[0037]
优选地，步骤s5包括：
[0038]
步骤s5.1：在已知总体分布类型和现场数据的情况下，确定各分布的似然函数；
[0039]
步骤s5.2：根据bayes理论，确定验后分布的核或分布密度函数；
[0040]
步骤s5.3：根据验后分布的核，判断验后分布的类型，得到参数的后验分布密度函数
[0041]
步骤s5.4：在已知参数的后验分布密度函数的情况下，根据总体分布，计算各可靠性评估指标的后验期望估计值，作为其贝叶斯估计值。
[0042]
根据本发明提供的一种产品可靠性综合评估系统，包括：
[0043]
模块m1：将多阶段多类型数据进行分类与预处理；
[0044]
模块m2：将分类与预处理之后的多阶段多类型数据进行一致性检验；
[0045]
模块m3：一致性检验后对多阶段多类型数据进行折算和综合；
[0046]
模块m4：多阶段多类型数据进行折算和综合之后根据数据类型确定先验分布和超参数；
[0047]
模块m5：确定先验分布和超参数之后，根据bayes理论及相关模型，进行数据融合和可靠性评估。
[0048]
根据本发明提供的一种存储有计算机程序的计算机可读存储介质，计算机程序被处理器执行时实现上述方法的步骤。
[0049]
根据本发明提供的一种产品可靠性综合评估设备，包括上述的产品可靠性综合评估系统或者上述的存储有计算机程序的计算机可读存储介质。
[0050]
与现有技术相比，本发明具有如下的有益效果：
[0051]
1、本发明通过从各阶段各类型数据的特点入手，对实际工程中存在的可靠性数据进行分类，对多阶段多类型数据进行分类和预处理，结合bayes及其相关模型，实现了多阶段多类型数据的产品可靠性综合评估。
[0052]
2、本发明能够充分利用产品在设计、制造、使用过程中的各类可靠性信息，特别是在缺乏足够的现场信息的情况下，可以提高可靠性评估的精度。
[0053]
3、本发明通过充分利用产品在设计、制造、使用过程中各阶段产生的各类可靠性数据，在确保数据一致性的前提下，对各阶段各类型可靠性数据进行分类、预处理、融合和评估，进一步提高了可靠性评估的精度。
附图说明
[0054]
通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：
[0055]
图1为本发明中的产品可靠性综合评估方法的基本流程图；
[0056]
图2为本发明中的多阶段多类型数据分类情况的示意图；
[0057]
图3为本发明中部件级异环境类成败型数据的预处理流程示意图；
[0058]
图4为本发明中部件级异环境类寿命型数据的预处理流程示意图；
[0059]
图5为本发明中部件级异环境类退化型数据的预处理流程示意图；
[0060]
图6为本发明中相似产品异环境类成败型数据的预处理流程示意图；
[0061]
图7为本发明中相似产品异环境类寿命型数据的预处理流程示意图；
[0062]
图8为本发明中相似产品异环境类退化型数据的预处理流程示意图；
[0063]
图9为本发明中产品现场类寿命型数据的预处理流程示意图。
具体实施方式
[0064]
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变化和改进。这些都属于本发明的保护范围。
[0065]
如图1
‑
图2所示，本发明提供了一种产品可靠性综合评估方法，从各阶段各类型数据的特点入手，对实际工程中存在的可靠性数据进行分类，给出了一种多阶段多类型数据的分类和预处理方法，结合bayes及其相关模型，提出了一套融合多阶段多类型数据的产品可靠性综合评估流程。包括如下步骤：
[0066]
步骤s1：将多阶段多类型数据进行分类与预处理，具体包括以下步骤：
[0067]
步骤s1.1：多阶段多类型数据进行分类，数据类型如下：
[0068]
产品现场类数据：产品实际使用或服役环境条件下的寿命型、成败型和退化型数据；
[0069]
产品异环境类数据：产品在寿命周期各阶段开展的不同试验条件下不同于现场使用环境条件下的寿命型、成败型和退化型数据；各阶段包括可靠性试验、寿命试验、加速寿命试验、功能试验、环境试验、定期试验等阶段。
[0070]
相似产品数据：又称为异总体数据，包括相似产品现场类数据和相似产品异环境类数据，相似产品现场类数据是指相似产品实际使用或服役环境条件下的寿命型、成败型和退化型数据，相似产品异环境类数据是指相似产品在寿命周期各阶段开展的不同试验条件下不同于现场使用环境条件下的寿命型、成败型和退化型数据；
[0071]
部件级数据：产品各组成部件的现场类数据和异环境类数据，部件级现场类数据是指产品各组成部件实际使用或服役环境条件下的寿命型、成败型和退化型数据，部件级异环境类数据是指产品各组成部件在寿命周期各阶段开展的不同试验条件下不同于现场使用环境条件下的寿命型、成败型和退化型数据。
[0072]
专家型数据：利用专家对产品可靠度的经验估计值来表示的数据。
[0073]
步骤s1.2：对多阶段多类型数据进行预处理，为了便于数据的表达和预处理，对上
述数据进一步归类，多阶段多类型数据可分为四类基本数据，即寿命型数据、成败型数据、退化型数据和专家型数据。四类基本数据类型的统一格式分别见表1至4：
[0074]
表1寿命型数据的统一格式
[0075]
样品编号寿命类型故障时间截尾时间相似因子环境因子
ꢀꢀꢀꢀꢀꢀ
[0076]
表2成败型数据的统一格式
[0077]
样品编号试验次数失败次数相似因子环境因子
ꢀꢀꢀꢀꢀ
[0078]
表3退化型数据的统一格式
[0079]
样品编号观测时刻性能退化量相似因子环境因子
ꢀꢀꢀꢀꢀ
[0080]
表4专家型数据的统一格式
[0081][0082]
表1中“寿命类型”由寿命分布类型检验预处理模块确定，“故障时间”和“截尾时间”用于记录不同截尾方式(包括定时、定数、定总、随机截尾)对应的样品故障时间和截尾时间。表2中的“试验次数”和“失败次数”用于记录样品成败型数据的总试验次数和失败次数。表3中的“观测时刻”和“性能退化量”用于记录样品退化型数据的观测时间点和性能退化量。表4中的可靠度的“点估计”和“置信下限估计”(给定置信度)为专家对可靠度给出的主观估计值。
[0083]
表1、表2和表3中的“相似因子”用于反映被评估产品与样品之间的相似程度，由相似产品数据预处理模块确定。当相似因子＝1时，表示样品技术状态与被评估产品完全相同；当0<相似因子<1时，表示样品的可靠性较被评估产品低；当相似因子>1时，表示样品的可靠性较被评估产品高。
[0084]
表1、表2和表3中的“环境因子”是异环境数据与现场数据的环境折算因子或加速系数。当环境因子＝1时，表示该数据为现场数据；当环境编号>0且≠1时，则表示该数据为异环境数据。
[0085]
进一步明确分类之后，分别针对寿命型数据、异环境类数据、部件级数据、退化型数据、相似产品数据进行预处理。
[0086]
针对寿命型数据的预处理：利用寿命分布类型检验和参数评估方法进行预处理，得到寿命类型及其参数估计；寿命类型包括指数寿命型、威布尔寿命型、正态寿命型、对数正态型等。
[0087]
针对异环境类数据的预处理：利用环境因子估算方法进行预处理，给出异环境类数据相对于现场数据的“环境因子”；异环境类数据包括产品异环境类数据、部件级异环境类数据、相似产品异环境类数据，环境因子估算方法是对不同分布类型的环境因子进行估算。
[0088]
针对部件级数据的预处理：利用系统可靠性综合方法进行预处理，将部件级数据折合成产品数据；部件级数据包括部件级现场类数据和异环境数据，系统可靠性综合方法包括mml、sr、l
‑
m、cmsr、bayes法等。
[0089]
针对退化型数据的预处理：利用退化建模方法进行预处理，将退化型数据转换成伪寿命型数据；退化型数据包括产品现场类退化型数据、部件级现场类退化型数据、部件级异环境类退化型数据，退化建模方法包括基于退化量分布、基于wiener、gamma等随机过程进行退化量分布。
[0090]
针对相似产品数据的预处理：利用相似因子估算方法进行预处理，给出相似产品数据的“相似因子”；相似产品数据包括相似产品寿命型数据、相似产品成败型数据、相似产品异环境类寿命型数据和相似产品异环境类成败型数据，相似因子估算方法包括综合评判法、专家打分法等。
[0091]
各类数据的处理流程如图3至图9所示。
[0092]
步骤s2：将分类与预处理之后的多阶段多类型数据进行一致性检验；由于现场数据能够真实地反映产品的可靠性水平，因此在一致性检验时，以产品现场数据为一致性检验基准，将其他多阶段多类型数据与现场数据进行寿命分布类型和失效机理一致性检验。只有通过一致性检验的数据才能用于产品的可靠性综合评估。基本方法如下：
[0093]
步骤s2.1：对于寿命型数据，根据寿命类型，判断多阶段多类型寿命数据(包括异环境寿命数据、相似产品寿命数据、伪寿命数据)与产品现场数据是否服从同一种寿命分布。
[0094]
步骤s2.2：对失效机理进行一致性检验，对于异环境寿命数据(环境因子≠1)，还需要进行失效机理的一致性检验。根据寿命分布的类型及其失效机理不变条件，见表5，判断异环境寿命数据是否满足失效机理不变条件，即环境因子不变原则。只有满足失效机理不变条件，才能将异环境数据和现场数据进行融合。
[0095]
表5常见寿命分布的环境因子与失效机理不变的条件
[0096][0097]
步骤s3：一致性检验后对多阶段多类型数据进行折算和综合。具体包括如下步骤：
[0098]
步骤s3.1：对异环境数据的折算和综合：对于环境因子k
i0
≠1的异环境数据，根据数据类型将异环境数据折算成现场数据，再将折算后的现场数据进行综合。基本方法如下：
[0099]
①
指数寿命型数据
[0100]
保持故障数z不变，试验总时间τ乘以加速因子k
i0
或其置信下限即
[0101][0102]
将所有异环境折算后的数据进行综合，其办法是：
[0103]
(z
′0,τ
′0)＝(z1,τ1·
k
10
) (z2,τ2·
k
20
)
…
(z
q
,τ
q
·
k
q0
)
ꢀꢀ
(2)
[0104]
②
其他寿命型数据(包括威布尔、正态、对数正态寿命型)
[0105]
所有故障时间数据和截尾数据乘以加速因子k
i0
或其置信下限即
[0106][0107]
其中：t
i
为环境s
i
下的试验数据，包括每个样本的故障时间或截尾时间。
[0108]
对所有异环境折算后的数据进行综合，其办法是：
[0109]
t
′0＝(t
′
10
,t
′
20
,
…
,t
′
q0
)
ꢀꢀ
(4)
[0110]
③
成败型数据
[0111]
保持失败数f
i
不变，试验次数n
i
乘以加速因子k
i0
[0112][0113]
对所有异环境折算后的数据进行综合，其办法是：
[0114]
(f
′0,n
′0)＝(f
′1,n
′1) (f
′2,n
′2)
…
(f
′
q
,n
′
q
)
ꢀꢀ
(6)
[0115]
步骤s3.2：对相似产品数据的折算和综合：对于相似因子f
i0
≠1的相似产品数据，根据数据类型将相似产品数据折算成被评估产品数据，再将折算后的被评估产品数据进行综合。基本方法与异环境数据的折算和综合方法相似，区别在于将环境因子换成相似因子。
[0116]
步骤s3.3：根据数据类型，将折算和综合后的异环境数据和相似产品数据再进行综合，得到寿命型数据和成败型数据的验前信息。
[0117]
步骤s4：多阶段多类型数据进行折算和综合之后根据数据类型确定先验分布和超参数。具体包括如下步骤：
[0118]
步骤s4.1：根据总体分布中未知参数的取值范围选取先验分布π(θ)，基本方法如下：
[0119]
对于[0,1]取值的参数，选用贝塔分布b(a,b)作为先验分布，超参数为(a,b)；对于(
‑
∞, ∞)取值的参数，选用正态分布n(μ,σ2)作为先验分布，超参数为(μ,σ2)；对于(0, ∞)取值的参数，选用伽马分布γ(z,τ)或倒伽马分布iγ(z,τ)作为先验分布，超参数为(z,τ)。
[0120]
各类验前分布的密度函数、期望和方差如表6所示。
[0121]
表6常见验前分布的概率密度函数、核、期望和方差
[0122][0123]
步骤s4.2：根据验前数据，确定验前分布的超参数，基本方法为利用验前信息，产生自助样本，得到各样本对应分布的参数的极大似然估计值样本，计算出参数的均值和方差，该均值和方差应等于对应参数的验前分布的均值和方差，联立方程组，得到寿命分布参数的验前分布的超参数。
[0124]
具体地，
①
对指数寿命型数据选取先验分布和超参数：
[0125]
指数分布的参数失效率λ,λ∈(0, ∞)的共轭验前分布为伽马分布，即π(λ)～γ(z,τ)，超参数(z,τ)分别表示指数型验前试验的故障数和试验总时间，其值由验前信息确定。
[0126]
确定超参数(z,τ)的方法如下：
[0127]
利用先验信息和bootstrap法获得失效率的若干估计值，记为计算先验均值和先验方差并令它们等于伽马分布γ(z,τ)的期望和方差，联立方程组解之，得到超参数的估计值：
[0128][0129]
②
对正态或对数正态型数据选取先验分布和超参数：
[0130]
正态或对数正态分布的参数μ和σ2未知且相关时，其参数(μ,σ2)的联合共轭分布是正态—逆伽马分布，即v0,μ0,k0为超参数，由验前数据确定。
[0131]
确定超参数的方法如下：
[0132]
对于均值，先利用验前信息和bootstrap法得到均值μ的验前均值和验前方差分别令其等于μ的验前分布的期望和方差
[0133][0134]
同理，对于方差，利用验前信息和bootstrap法得到方差σ2的验前均值和验前方差分别令其等于σ2的验前分布的期望和方差，即
[0135][0136]
联立上面两个方程组，解得
[0137][0138]
③
对威布尔寿命型数据选取先验分布和超参数：
[0139]
威布尔分布的参数有(m,η)。根据m的取值范围，分两种情况讨论：
[0140]
情况一：若已知失效率是递减的，则0<m<1，选取贝塔分布为m的先验分布，即取m～β(a,b)，λ＝(1/η)
m
～γ(z,τ)，则在m和λ是相互独立的假设下，它们的先验分布为：
[0141]
π(λ,m)
∝
λ
z
‑1e
‑
τλ
m
a
‑1(1
‑
m)
b
‑1ꢀꢀ
(10)
[0142]
超参数(a,b)和(z,τ)的确定方法如下：
[0143]
对于参数m，利用验前信息和bootstrap法得到m的验前均值和方差分别令其等于m的验前分布m～β(a,b)的期望和方差，即
[0144][0145]
求上述方程组，得到超参数(a,b)：
[0146][0147]
同样，对于参数λ，利用验前信息和bootstrap法的λ的验前均值和方差分别令其等于λ的验前分布λ＝(1/η)
m
～γ(z,τ)的期望和方差，即
[0148][0149]
求上述方程组，得到超参数(z,τ)：
[0150][0151]
情况二：已知失效率是递增的，则m>1，则可选取如下先验分布，使m
′
＝m
‑
1服从伽马分布。由于m>1，则令m
′
＝m
‑
1的先验分布为伽马分布，即取m
′
＝m
‑
1～γ(a,b)，λ＝(1/η)
m
～γ(d,τ)，则在m和λ是相互独立的假设下，它们的先验分布为
[0152]
π(λ,m)
∝
λ
z
‑1e
‑
τλ
(m
‑
1)
a
‑1e
‑
b(m
‑
1)
ꢀꢀ
(14)
[0153]
超参数(a,b)和(z,τ)的确定方法与情况一相同。
[0154]
④
对成败型数据选取先验分布和超参数：
[0155]
二项分布的参数为成功率或可靠度r：共轭验前分布为贝塔分布，即π(r)～b(a0,b0)，超参数a0,b0分别表示验前试验的成功次数和失败次数，其值由验前信息确定。验前分布密度函数：
[0156][0157]
确定超参数(a0,b0)的方法如下：
[0158]
利用先验信息和bootstrap法获得成功概率(可靠度)的若干估计值，记为
计算先验均值和先验方差并令它们等于贝塔分布b(α,β)的期望和方差，利用如下方程组得到超参数(a0,b0)：
[0159][0160]
⑤
当验前信息为专家信息时，选取先验分布方法如下：
[0161]
当现场数据为成败型时，选取贝塔分布为先验分布；
[0162]
当现场数据为寿命型时，选取负对数伽马分布为先验分布。
[0163]
确定超参数的方法：
[0164]
已知专家对可靠度的点估计和置信下限r
0l
(置信度1
‑
β)，当r的先验分布为贝塔分布时，由下列方程组确定超参数(a,b)
[0165][0166]
当r的先验分布为负对数伽马分布时，令任务时间t0，τ为试验总时间，则等效任务数为η＝τ/t0。由下列方程组确定超参数(z,η)
[0167][0168]
步骤s5：确定先验分布π(θ)和超参数之后，根据bayes理论及相关模型，进行数据融合和可靠性评估。具体包括如下步骤：
[0169]
步骤s5.1：在已知总体分布类型和现场数据的情况下，确定各分布的似然函数p(x|θ)。
[0170]
各类总体分布的似然函数如表7所示。
[0171]
表7常见总体分布的似然函数
[0172][0173]
步骤s5.2：根据bayes理论，确定验后分布的核或分布密度函数。
[0174]
验后分布密度函数为：
[0175][0176]
其中，c＝1/∫
θ
p(x|θ)π(θ)dθ为一个与θ无关的正则常数，p(x|θ)π(θ)为验后分布的核，即
[0177]
π(θ|x)
∝
p(x|θ)π(θ)
[0178]
步骤s5.3：根据验后分布的核，判断验后分布的类型，得到参数的后验分布密度函数。
[0179]
具体地，
①
对于二项分布，试验数据：(n,s)，n为试验次数，s为成功次数。
[0180]
参数：可靠度r,r∈[0,1]，其共轭验前分布为贝塔分布，即π(r)～be(s0,f0)，超参数s0,f0分别表示验前试验的成功次数和失败次数，其值由验前信息确定。
[0181]
验前分布密度函数
[0182][0183]
似然函数为
[0184][0185]
验后分布的核
[0186]
π(r|n,s)
∝
p(n,s|r)π(r)
[0187]
r的验后分布：也为贝塔分布，即π(r|n,s)～be(s0 s,f0 f)，其pdf为
[0188][0189]
②
指数分布
[0190]
试验数据：(z,τ)，z＝r为失效数，τ＝t
r
为试验总时间
[0191]
参数：失效率λ,λ∈(0, ∞)，其共轭验前分布为伽马分布，即π(λ)～ga(z0,τ0)，超参数z0,τ0分别表示指数型验前试验的失效数和试验总时间，其值由验前信息确定。
[0192]
似然函数为：
[0193]
p(z,τ|λ)＝
[0194]
λ验后分布：也为伽马分布，即，其pdf为
[0195][0196]
③
正态分布和对数正态分布
[0197]
试验数据：
[0198]
正态x
(1)
≤x
(2)
≤
…
≤x
(n)
[0199]
对数正态lnx
(1)
≤lnx
(2)
≤
…
≤lnx
(n)
[0200]
若为寿命数据，则均应大于0。这里只考虑完全样本的情况。
[0201]
参数：存在以下三种情况：
[0202]
情况(一)：μ未知(σ2已知)
[0203]
μ的验前共轭分布为正态分布，即超参数μ0,分别表示验前数据的均值和方差，由验前数据估计得到。
[0204]
验前分布的核为：
[0205][0206]
似然函数为：
[0207][0208]
μ的共轭验后分布仍为正态分布，即其中：
[0209][0210][0211][0212]
情况(二)：σ2未知(μ已知)
[0213]
σ2的验前共轭分布为逆伽马分布，即π(σ2)～iga(a,b)，超参数a,b由验前数据估计
得到。
[0214]
验前分布密度函数：
[0215][0216]
似然函数：
[0217][0218]
σ2的共轭验后分布仍为逆伽马分布
[0219][0220][0221]
情况(三)：μ和σ2未知且相关
[0222]
参数(μ,σ2)的联合共轭分布是正态—逆伽马分布，即v0,μ0,k0为超参数，由验前数据确定。
[0223]
其中：即
[0224][0225][0226]
则联合验前分布
[0227][0228]
似然函数为
[0229][0230]
则联合后验分布也为正态－逆伽马分布
[0231][0232]
其中：
[0233][0234]
μ的后验分布为自由度为v
n
的学生t分布，即
[0235]
σ2的后验分布为逆伽马分布，即
[0236]
④
威布尔分布
[0237]
试验数据：
[0238]
t1≤t2≤
…
≤t
r
，即n个样本的前r个故障时间
[0239]
两参数威布尔分布的分布函数：
[0240][0241]
参数：分以下三种情况讨论：
[0242]
情况(一)：形状参数m已知，η未知
[0243]
试验数据可以表示成则这些数据为来自指数分布exp(λ)，λ＝(1/η)
m
的前r个次序统计量，则参数λ的验后分布可以按照指数分布给出。
[0244]
情况(二)：失效率是递减的，即0＜m＜1，选取贝塔分布为m的先验分布。
[0245]
取m～β(a,b)，λ＝(1/η)
m
～γ(d,τ)，则在m和λ是相互独立的假设下，先验分布为
[0246]
π(λ,m)
∝
λ
d
‑1e
‑
τλ
m
a
‑1(1
‑
m)
b
‑1[0247]
给定试验样本数据条件下，似然函数为
[0248][0249]
其中：表示以为失效时间的总试验时间。
[0250]
则m和λ的联合后验分布为
[0251][0252]
即
[0253][0254]
其中：c1为与m和λ无关的常数，其值为
[0255][0256]
由此可得，m的后验分布为
[0257][0258]
λ的后验分布为
[0259][0260]
情况(三)：失效率是递增的，即m＞1，则可选取如下先验分布，使m
′
＝m
‑
1服从伽马分布，即m
′
＝m
‑
1～γ(a,b)，选取λ＝(1/η)
m
～γ(d,τ)，则在m和λ是相互独立的假设下，先验分布为
[0261]
π(λ,m)
∝
λ
d
‑1e
‑
τλ
(m
‑
1)
a
‑1e
‑
b(m
‑
1)
[0262]
给定试验样本数据条件下，似然函数为
[0263][0264]
其中：表示以为失效时间的总试验时间。
[0265]
则m和λ的联合后验分布为
[0266][0267]
即
[0268][0269]
其中：c2为与m和λ无关的常数，其值为
[0270][0271]
由此可得，m的后验分布为
[0272][0273]
λ的后验分布为
[0274][0275]
步骤s5.4：在已知参数的后验分布密度函数的情况下，根据总体分布，计算各可靠性评估指标的后验期望估计值，作为其贝叶斯估计值。
[0276]
具体地，
①
对于二项分布，r的贝叶斯估计(后验期望)：
[0277][0278]
②
指数分布
[0279]
λ的贝叶斯估计(后验期望)：
[0280][0281]
平均寿命的贝叶斯估计(后验期望)：
[0282][0283]
任务时间t0时的可靠度的贝叶斯估计(后验期望)：
[0284][0285]
③
正态分布和对数正态分布
[0286]
情况(一)
[0287]
μ的贝叶斯估计：
[0288][0289]
任务时间t0时的可靠度的点估计：
[0290][0291]
情况(二)
[0292]
σ2的贝叶斯估计：
[0293][0294]
任务时间t0时的可靠度的贝叶斯估计(后验期望)：
[0295][0296]
其中：c1为与s2有关的归一化常数。
[0297]
情况(三)
[0298]
μ的期望后验估计：
[0299][0300]
σ2的期望后验估计：
[0301][0302]
r的bayes估计：
[0303][0304]
其中：c2为一个与s2有关的常数。
[0305]
④
威布尔分布
[0306]
情况(一)
[0307]
同指数分布。
[0308]
情况(二)
[0309]
m的贝叶斯估计：
[0310][0311]
λ的贝叶斯估计：
[0312][0313]
可靠度的贝叶斯估计：
[0314][0315]
情况(三)
[0316]
m的贝叶斯估计：
[0317][0318]
λ的贝叶斯估计：
[0319][0320]
可靠度的贝叶斯估计：
[0321][0322]
本发明还提供了一种产品可靠性综合评估系统，包括：
[0323]
模块m1：将多阶段多类型数据进行分类与预处理；
[0324]
模块m2：将分类与预处理之后的多阶段多类型数据进行一致性检验；
[0325]
模块m3：一致性检验后对多阶段多类型数据进行折算和综合；
[0326]
模块m4：多阶段多类型数据进行折算和综合之后根据数据类型确定先验分布和超参数；
[0327]
模块m5：确定先验分布和超参数之后，根据bayes理论及相关模型，进行数据融合和可靠性评估。
[0328]
本发明又提供了一种存储有计算机程序的计算机可读存储介质，计算机程序被处理器执行时实现上述方法的步骤。
[0329]
本发明继续提供的一种产品可靠性综合评估设备，包括上述的产品可靠性综合评估系统或者上述的存储有计算机程序的计算机可读存储介质。
[0330]
本领域技术人员知道，除了以纯计算机可读程序代码方式实现本发明提供的系统
及其各个装置、模块、单元以外，完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统及其各个装置、模块、单元以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同功能。所以，本发明提供的系统及其各项装置、模块、单元可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置、模块、单元也可以视为硬件部件内的结构；也可以将用于实现各种功能的装置、模块、单元视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
[0331]
以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变化或修改，这并不影响本发明的实质内容。在不冲突的情况下，本技术的实施例和实施例中的特征可以任意相互组合。