CN103198223B - 一种电子产品实时可靠性的预测方法 - Google Patents
一种电子产品实时可靠性的预测方法 Download PDFInfo
- Publication number
- CN103198223B CN103198223B CN201310126195.XA CN201310126195A CN103198223B CN 103198223 B CN103198223 B CN 103198223B CN 201310126195 A CN201310126195 A CN 201310126195A CN 103198223 B CN103198223 B CN 103198223B
- Authority
- CN
- China
- Prior art keywords
- sigma
- time
- electronic product
- reliability
- average
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Complex Calculations (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种电子产品实时可靠性的预测方法基于贝叶斯方法和伪失效寿命,在更为恰当地估计先验分布中的未知参数和选择实时可靠性公式的基础上计算出当前电子产品实时可靠性。首先运用曲线拟合推出n个伪失效寿命,再选择正态分布来表示n个伪失效寿命数据的分布,接着借助时间序列样本生成法,估计出分布中的未知时变参数均值μj0及方差σj0,得到先验密度函数,然后根据现场数据xj可以更新时变参数并得到时变参数的后验密度函数的均值μcj和方差σcj 2,最后利用本发明设计的实时可靠性公式计算当前电子产品的实时可靠性。通过实验验证,本发明电子产品实时可靠性的预测方法对电子产品实时可靠性的预测精度高,能准确地对电子产品的实时可靠性进行预测。
Description
技术领域
本发明属于可靠性分析技术领域,更为具体地讲,涉及一种电子产品实时可靠性的预测方法。
背景技术
传统的可靠性分析方法假设可靠度是基于时间的概率分布,并通过对产品进行大量的实验并采样其寿命数据获取该分布。获得的分布往往反映了在相同测试环境下同一类产品的平均特征,常用于系统的设计阶段。但是对于运行系统中的部件,它的实时可靠性有着自身的个性特征,即同类产品特定个体之间的可靠性存在着差异,与同一类产品的可靠性也不尽相同。若用传统可靠性方法获取的失效分布进行实时定量可靠性分析,难免出现较大的误差,也会造成维护成本偏高和维护效率低下。因此,获取运行产品实时可靠性便成为安全性与可靠性工程的强烈需求。
贝叶斯方法提供的是一种计算假设概率的方法。这种方法是基于假设的先验概率、给定假设条件下观察到不同数据的概率以及观察到的数据本身而得出来的。贝叶斯方法可以最大程度上利用现场信息,融合现场数据。由于最后结果是用现场数据去更新先验分布后得到的,因此在推导先验概率时,需要采用恰当的方法去估计先验分布中的未知参数,以保证最后结果的准确性。由于曲线拟合总是会引入更多的误差,因此在实时可靠性计算公式上,应该直接选择反映总体实时可靠性的公式,而不是先计算单独时间点的实时可靠性,在通过曲线拟合得到总体的实时可靠性。
发明内容
本发明的目的在于克服现有技术的不足,提供一种电子产品实时可靠性的预测方法,以提高预测精度,准确地对电子产品的实时可靠性进行预测。
为实现以上目的,本发明电子产品实时可靠性的预测方法,其特征在于,包括以下步骤:
(1)、选取数量为n的同类电子产品,n大于等于10-15,对每个电子产品从最初运行开始,间隔一定时间ts提取测试点的信号幅值作为电子产品性能退化的历史数据,得到n组由m个信号幅值组成的历史数据;
(2)、分别对每一组历史数据进行曲线拟合,根据各组拟合结果来选择最理想的拟合模型,根据拟合模型来外推电子产品信号幅值在到达失效阈值d时的伪失效寿命xi,其中i=1,2,…,n,为电子产品序号;
(3)、计算电子产品的寿命T的方差σj 2:
方差:
电子产品的寿命T的均值μj为正态分布,均值μj的均值μj0及方差σj0 2通过以下步骤得到;
3.1)、根据伪失效寿命xi,得到时间序列样本yk:
其中,k=1,2,…,n
3.2)、根据时间序列样本yk计算出均值μj0及方差σj0 2;
(4)、计算用现场数据更新历史数据得到的包含后验信息的伪失效寿命T的均值μcj和方差
4.1)、对于当前电子产品也从最初运行开始,间隔一定时间ts提取测试点的信号幅值,得到当前电子产品性能退化的由m个信号幅值组成的现场数据;根据拟合模型来外推当前电子产品信号幅值在到达失效阈值d时的伪失效寿命xj;
4.2)、计算均值μcj和方差σcj 2:
μcj=μja
σcj 2=σja 2+σj 2
其中,
(5)、计算当前电子产品实时可靠性
将均值μcj和方差σcj 2代入可靠性计算公式,得到实时可靠度R(t):
其中,Φ(·)为标准正态分布函数。
本发明的目的是这样实现的:
本发明电子产品实时可靠性的预测方法基于贝叶斯方法和伪失效寿命,在更为恰当地估计先验分布中的未知参数和选择实时可靠性公式的基础上计算出当前电子产品实时可靠性。首先运用曲线拟合推出n个伪失效寿命,再选择正态分布来表示n个伪失效寿命数据的分布,接着借助时间序列样本生成法,估计出分布中的未知时变参数均值μj0及方差σj0,得到先验密度函数,然后根据现场数据xj可以更新时变参数并得到时变参数的后验密度函数的均值μcj和方差σcj 2,最后利用本发明设计的实时可靠性公式计算当前电子产品的实时可靠性。通过实验验证,本发明电子产品实时可靠性的预测方法对电子产品实时可靠性的预测精度高,能准确地对电子产品的实时可靠性进行预测。
附图说明
图1是本发明电子产品实时可靠性的预测方法一种具体实施方式流程图;
图2是历史数据拟合曲线图;
图3是拟合曲线与历史数据的残差曲线图;
图4是电子产品实时可靠性曲线图。
具体实施方式
下面结合附图对本发明的具体实施方式进行描述,以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是,在以下的描述中,当已知功能和设计的详细描述也许会淡化本发明的主要内容时,这些描述在这里将被忽略。
图1是本发明电子产品实时可靠性的预测方法一种具体实施方式流程图。
在本实施例中,如图2所示,本发明电子产品实时可靠性的预测方法包括以下步骤:
步骤ST1:选取数量为n的同类电子产品,对每个电子产品从最初运行开始,间隔一定时间ts提取测试点的信号幅值作为电子产品性能退化的历史数据,得到n组由m个信号幅值组成的历史数据;
步骤ST2:分别对每一组历史数据进行曲线拟合,根据各组拟合结果来选择最理想的拟合模型,根据拟合模型来外推电子产品信号幅值在到达失效阈值d时的伪失效寿命xi,其中i=1,2,…,n,为电子产品序号;
通常的退化轨迹有三种:线性即x=a+bt),凸形即log(x)=a+bt和凹形即log(x)=a+blog(t)。在本实施例中,利用Matlab的曲线拟合工具箱cftool来历史数据进行曲线拟合,根据拟合曲线和历史数据之间的残差的平方和最小的原则来选择最理想的拟合模型。
在本实施例中,每一组历史数据拟合后得到拟合曲线与历史数据之间的一系列残差当作一个离散随机过程,借助于自回归AR模型来估计预测残差,用估计出的残差结果来修正外推出的伪失效寿命xi。
步骤ST3:计算电子产品的寿命T的方差σj 2:
方差:
电子产品的寿命T的均值μj为正态分布,其均值μj0及方差σj0 2通过以下步骤得到;
3.1)、根据伪失效寿命xi,得到时间序列样本yk:
其中,k=1,2,…,n
3.2)、根据时间序列样本yk计算出均值μj0及方差σj0 2;
在通常情况下,产品的失效分布都为正态分布,即寿命T~N(μ,σ2)其概率密度函数为:
其中,μ、σ2分别为寿命T的均值和方差。
先验密度函数的确定
对于稳定的电子产品或者系统,在估计性能参数时,本发明电子产品的寿命T的均值μj未知,而方差σj 2已知的情况下进行。
在本发明中,σj 2可以用下面的公式估计:
均值μj服从正态分布,即μj~N(μj0,σj0 2),其中的未知参数即均值μj0和方差σj0 2的估计需要涉及到利用数量有限的历史数据来生成若干组时间序列样本,通过这些时间序列样本来计算未知参数即均值μj0和方差σj0 2。
通常生成时间序列样本的办法的两种:一种是分组样本容量小,分组数量大,即将历史数据尽量均匀分成较多的组,而每一个样本里的数据较少;另一种则是分组样本容量大,分组数量少,即将历史数据分成较少的组,而每一个样本里的数据较多。用各组样本的均值作为统计对象来计算均值μj0和方差σj0 2。
但是这两种方法都存在着一定的问题。对于第一种方法,由于样本内数据数量很少,存在着即使获得大量数据后均值μj0依然离散和趋近于0这样一个明显错误结果的风险。对于第二种方法,虽然样本内数据数量较多,不存在上面的风险,但是由于样本的数量很少,导致计算统计意义上的均值μj0和方差σj0 2时有时会精度不高,无法反映真实的情况。
针对以上两种方法存在的问题,本发明提出一种时序序列样本分割法。将第一个伪失效寿命x1作为第一个时间序列样本y1,前两个伪失效寿命x1、x2的均值作为第二个时间序列样本y2,前三个数据伪失效寿命x1、x2、x3的均值作为第三个时间序列样本y2,…,前n个数据伪失效寿命的均值作为第n个时间序列样本yn,这样就生成了n个时间序列样本,每个时间序列样本的数据数量分别为1,2,3,…,n。这样生成时间序列样本的方法既可以保证时间序列样本的数量,又可以在一定程度上保证时间序列样本内数据的数量,这样就不会出现上面两种方法所存在的问题和风险了,实现了更为优化的样本生成方法。
步骤ST4:计算用现场数据更新历史数据得到的包含后验信息的伪失效寿命T的均值μcj和方差
4.1)、对于当前电子产品也从最初运行开始,间隔一定时间ts提取测试点的信号幅值,得到当前电子产品性能退化的由m个信号幅值组成的现场数据;根据拟合模型来外推当前电子产品信号幅值在到达失效阈值d时的伪失效寿命xj;
4.2)、计算均值μcj和方差σcj 2:
μcj=μja
σcj 2=σja 2+σj 2
其中,
则寿命T的后验密度函数为f(T|xj),根据贝叶斯公式,寿命T的均值的先验均值μj的后验密度函数为:
由于正态分布具有共轭特性,即后验分布和先验分布具有相同的分布类型可以得到:f(μj|xj)=f(μj;μja,σja),其中:
然后可以得到f(T|xj),有
由上式可以看出,f(T|xj)是一个正态分布,均值μcj=μja,方差σcj 2=σja 2+σj 2。μcj和σcj 2分别用现场数据更新了历史数据后得到的包含了后验信息的伪失效寿命T的均值和方差。
步骤ST5:计算当前电子产品实时可靠性
将均值μcj和方差σcj 2代入可靠性计算公式,得到实时可靠度R(t):
其中,Φ(·)为标准正态分布函数。
这样就可以针对一类电子产品中的样品(即选取的数量为n的同类电子产品)计算出实时的总体可靠度了。这样也避免先计算各个时刻的实时可靠度,再通过曲线拟合来得出总体可靠度的方法存在的一个问题,就是会在拟合阶段引入额外的误差,降低最后结果的精度。
实例验证
表1是验证本发明预测方法的一种具体电子产品(GaAs激光器)性能(工作电流)退化数据表。一共提取了15个样品,做了15组退化数据;对于每一组退化数据,每间隔ts=250小时读取一次工作电流的幅值,一共读取了17次。其中将第13组退化数据选为现场数据,另外14组数据选为历史数据,即n=14。
表1
将每一组历史数据都进行曲线拟合,通过拟合检验,发现线性拟合最能符合退化数据的变化趋势,如图2所示。
拟合曲线与历史数据之间的残差,以第12组样本在0-3000小时之间的数据为例,如图3所示。将残差数据带入AR模型,得到第3250小时处的残差预测值,以此类推再一次得到第3500小时、3750小时、4000小时处的残差预测值,依次修正伪失效寿命。
n=14组历史数据得到电子产品信号幅值在到达失效阈值d=10.0时的伪失效寿命xi,其中i=1,2,…,n即14个伪失效寿命值,然后采样本发明中时间序列样本生成法,生成用来估计的时间序列样本yk。从而得到根据时间序列样本yk计算出均值μj0及方差σj0,然后可以将现场数据融合入历史数据中,得到更新后的参数均值μcj和方差σcj 2。最后将μcj和方差σcj 2带入本发明中提出的可靠性公式中就可以得到最终的当前电子产品的可靠性,如图4所示。
在图4中,从左向右的四条曲依次是融合了前11个现场数据即m=11、前13个现场数据即m=13、前15个现场数据即m=15和融合17个即全部现场数据即m=17的实时可靠性结果曲线。从图4可以清楚地看到,随着融合现场数据量即信号幅值点m的增加,可靠性曲线逐渐向融合全部现场数据时获得的产品可靠性曲线靠近,这也说明利用的现场信息越多,就越能得到反映当前电子产品的实时可靠性模型。
尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
Claims (2)
1.一种电子产品实时可靠性的预测方法,其特征在于,包括以下步骤:
(1)、选取数量为n的同类电子产品,n大于等于10,对每个电子产品从最初运行开始,间隔一定时间ts提取测试点的信号幅值作为电子产品性能退化的历史数据,得到n组由m个信号幅值组成的历史数据;
(2)、分别对每一组历史数据进行曲线拟合,根据各组拟合结果来选择最理想的拟合模型,根据拟合模型来外推电子产品在到达失效阈值d时的伪失效寿命xi,其中i为电子产品序号,i=1,2,…,n;
(3)、计算电子产品的寿命T的方差σj 2:
方差:
电子产品的寿命T的均值μj为正态分布,均值μj的均值μj0及方差σj0 2通过以下步骤得到;
3.1)、根据伪失效寿命xi,得到时间序列样本yk:
其中,k=1,2,…,n
3.2)、根据时间序列样本yk计算出均值μj0及方差σj0 2;
(4)、计算用现场数据更新历史数据得到的包含后验信息的伪失效寿命T的均值μcj和方差
4.1)、对于当前电子产品也从最初运行开始,间隔一定时间ts提取测试点的信号幅值,得到当前电子产品性能退化的由m个信号幅值组成的现场数据;根据拟合模型来外推当前电子产品在到达失效阈值d时的伪失效寿命xj;
4.2)、计算均值μcj和方差σcj 2:
μcj=μja
σcj 2=σja 2+σj 2
其中,
(5)、计算当前电子产品实时可靠性
将均值μcj和方差σcj 2代入可靠性计算公式,得到实时可靠度R(t):
其中,Φ(·)为标准正态分布函数。
2.根据权利要求1所述的实时可靠性的预测方法,其特征在于,在步骤(2)中,所述的选择最理想的拟合模型为根据拟合曲线和历史数据之间的残差的平方和最小的原则来选择最理想的拟合模型;
每一组历史数据拟合后得到拟合曲线与历史数据之间的一系列残差当作一个离散随机过程,借助于自回归AR模型来估计预测残差,用估计出的残差结果来修正外推出的伪失效寿命xi。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310126195.XA CN103198223B (zh) | 2013-04-12 | 2013-04-12 | 一种电子产品实时可靠性的预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310126195.XA CN103198223B (zh) | 2013-04-12 | 2013-04-12 | 一种电子产品实时可靠性的预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103198223A CN103198223A (zh) | 2013-07-10 |
CN103198223B true CN103198223B (zh) | 2015-12-09 |
Family
ID=48720777
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310126195.XA Active CN103198223B (zh) | 2013-04-12 | 2013-04-12 | 一种电子产品实时可靠性的预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103198223B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113406557A (zh) * | 2021-06-23 | 2021-09-17 | 深圳市计量质量检测研究院(国家高新技术计量站、国家数字电子产品质量监督检验中心) | 一种传播式充电桩远程检定方法 |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105354402A (zh) * | 2014-08-18 | 2016-02-24 | 鲍珂 | 一种车辆变速箱耗损型失效的可靠性评估方法 |
CN104765945B (zh) * | 2015-02-09 | 2016-08-24 | 广东电网有限责任公司电力科学研究院 | 一种基于数量置信度的设备缺陷率检测方法 |
CN104978459B (zh) * | 2015-07-08 | 2017-11-21 | 电子科技大学 | 一种igbt剩余寿命预测方法 |
CN106248381B (zh) * | 2016-10-11 | 2019-04-09 | 西安交通大学 | 一种基于多特征和相空间的滚动轴承寿命动态预测方法 |
CN106932708A (zh) * | 2017-02-10 | 2017-07-07 | 电子科技大学 | 电子封装焊点疲劳寿命分析方法 |
CN108038317B (zh) * | 2017-12-18 | 2021-10-29 | 广东科鉴检测工程技术有限公司 | 精密仪器性能参数保持期预测方法和系统 |
CN110457654A (zh) * | 2019-08-08 | 2019-11-15 | 哈尔滨理工大学 | 一种基于现场数据的机载设备可靠性预测方法 |
CN113743670B (zh) * | 2021-09-08 | 2023-05-09 | 电子科技大学 | 一种基于gru模型的电路故障实时预测方法及验证电路 |
CN117056692B (zh) * | 2023-10-09 | 2024-01-23 | 山东芯赛思电子科技有限公司 | 一种SiC基电机驱动器件的老化预测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101620045A (zh) * | 2009-07-31 | 2010-01-06 | 北京航空航天大学 | 基于时间序列的步进应力加速退化试验可靠性评估方法 |
US20110018575A1 (en) * | 2009-07-23 | 2011-01-27 | International Business Machines Corporation | Method and system for assessing reliability of integrated circuit |
CN102033182A (zh) * | 2010-12-10 | 2011-04-27 | 北京航空航天大学 | 一种固体钽电解电容器寿命预测方法 |
CN102629300A (zh) * | 2012-03-15 | 2012-08-08 | 北京航空航天大学 | 一种基于灰色预测模型的步进应力加速退化数据评估方法 |
-
2013
- 2013-04-12 CN CN201310126195.XA patent/CN103198223B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110018575A1 (en) * | 2009-07-23 | 2011-01-27 | International Business Machines Corporation | Method and system for assessing reliability of integrated circuit |
CN101620045A (zh) * | 2009-07-31 | 2010-01-06 | 北京航空航天大学 | 基于时间序列的步进应力加速退化试验可靠性评估方法 |
CN102033182A (zh) * | 2010-12-10 | 2011-04-27 | 北京航空航天大学 | 一种固体钽电解电容器寿命预测方法 |
CN102629300A (zh) * | 2012-03-15 | 2012-08-08 | 北京航空航天大学 | 一种基于灰色预测模型的步进应力加速退化数据评估方法 |
Non-Patent Citations (5)
Title |
---|
Comparison of Statistical Models for the Lumen Lifetime Distribution of High Power White LEDs;Jiajie Fan, K.c.Yung,Michael Pecht;《Prognostics & System Health Management Conference》;20121231;全文 * |
Real-Time Performance Reliability Prediction;Huitian Lu, William J. Kolarik, and Susan S. Lu;《IEEE TRANSACTIONS ON RELIABILITY》;20011231;第50卷(第4期);353-357 * |
一种基于伪失效寿命的LED可靠性快速评价方法;郭伟玲,樊星,崔德胜,吴国庆,俞鑫;《发光学报》;20130228;第34卷(第2期);213-217 * |
一种基于伪失效阂值的寿命分析方法;陈亮,胡昌华;《解放军理工大学学报》;20110430;第12卷(第2期);156-160 * |
基于贝叶斯方法的退化失效型产品实时可靠性评估;赵熠,董豆豆,周经伦,冯静;《国防科技大学学报》;20071231;第29卷(第6期);118-120 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113406557A (zh) * | 2021-06-23 | 2021-09-17 | 深圳市计量质量检测研究院(国家高新技术计量站、国家数字电子产品质量监督检验中心) | 一种传播式充电桩远程检定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103198223A (zh) | 2013-07-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103198223B (zh) | 一种电子产品实时可靠性的预测方法 | |
Young | Advances in real–time flood forecasting | |
Fitzgerald et al. | Monte Carlo transition dynamics and variance reduction | |
Garratt et al. | Real-time representations of the output gap | |
CN103310115B (zh) | 一种多目标跟踪的杂波估计方法 | |
CN103955892B (zh) | 一种目标跟踪方法及扩展截断无迹卡尔曼滤波方法、装置 | |
US20130297239A1 (en) | Estimating Loss Rates of Links in Smart Grids | |
CN102082560A (zh) | 一种基于集合卡尔曼滤波的粒子滤波方法 | |
CN105164675A (zh) | 增量式地更新统计 | |
CN109800915A (zh) | 一种基于缺失数据集的交通流量预测方法 | |
CN102508774A (zh) | 基于新环境因子函数的软件可靠性增长模型的建模方法 | |
CN105718642A (zh) | 一种基于门限自回归模型的参考时间尺度产生方法 | |
Warin | Variations on branching methods for non linear PDEs | |
Li et al. | Comparison on PPCA, KPPCA and MPPCA based missing data imputing for traffic flow | |
Shang et al. | On parameter estimation for the generalized gamma distribution based on left‐truncated and right‐censored data | |
CN102681936B (zh) | 一种针对金融系统测试结果的验证方法及装置 | |
Li et al. | A self‐exciting marked point process model for drought analysis | |
Garratt et al. | Real time representation of the UK output gap in the presence of model uncertainty | |
KR20230061338A (ko) | 오류 감소 기술 | |
Richards | A fractal forecasting model for financial time series | |
CN112396535A (zh) | 智能电网的管理方法、装置、设备及存储介质 | |
Vatamidou et al. | Corrected phase-type approximations of heavy-tailed risk models using perturbation analysis | |
Kunihama et al. | Efficient estimation and particle filter for max‐stable processes | |
Howe et al. | A total estimator of the Hadamard function used for GPS operations | |
CN113762685A (zh) | 异常用户的检测方法、装置及电子设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |