CN109325289B - 一种估计电子件可靠性参数的方法 - Google Patents
一种估计电子件可靠性参数的方法 Download PDFInfo
- Publication number
- CN109325289B CN109325289B CN201811083876.1A CN201811083876A CN109325289B CN 109325289 B CN109325289 B CN 109325289B CN 201811083876 A CN201811083876 A CN 201811083876A CN 109325289 B CN109325289 B CN 109325289B
- Authority
- CN
- China
- Prior art keywords
- parameter
- value
- parameters
- mean
- equal
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Landscapes
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种估计电子件可靠性参数的方法,包括以下步骤:步骤一:确定候选的寿命分布参数,根据已有电子件可靠性数据分布规律,初步确定待估计电子件寿命均值参数的下限值μmin和上限值μmax,在确定的参数区间内,等间隔生成n个候选分布参数,其中,候选分布参数中,相邻参数之间的步长相等为d1,根据待测电子件均值参数的步长依次计算出n个候选均值参数μj,其中1≤j≤n;步骤二:遍历寿命均值参数μj计算似然值Lj,寻找最大似然值记为LM,则最大似然值LM对应的寿命分布参数μM即为均值参数的估计值。本发明利用“少量的可靠性试验数据+在产品研制、生产、使用等阶段产生的大量数据”,估算产品寿命的分布规律。
Description
技术领域
本发明属于可靠性测试技术领域,特别涉及一种估计电子件可靠性参数的方法。
背景技术
可靠性是描述产品质量的核心属性,通常用寿命的分布规律(分布类型和参数)来定量描述可靠性。理论上,针对产品开展大量的可靠性试验,可以获得足够数量的产品寿命数据,然后可采用成熟的数理统计方法来估计出产品寿命的分布类型和参数。但在实际工作中,针对产品开展大量可靠性试验,往往意味着高昂的经济成本和漫长的试验耗时,因此更常见的做法是利用“少量的可靠性试验数据+在产品研制、生产、使用等阶段产生的大量数据”,估算产品寿命的分布规律。在产品的可靠性试验中,一般配备有专门的在线检测设备,用于实时监测产品的完好性状态,及时记录产品的故障时刻,因此可以获得产品寿命的数值。但在产品研制、生产、使用等这些非可靠性试验场景下,不一定配备有专门的在线检测设备,只能定期或不定期地对产品进行完好性检查,因而不能准确获知产品的故障时刻,也就无法获得寿命的数值信息。
发明内容
为了克服背景技术中存在的缺陷,本发明提供一种估计电子件可靠性参数的方法。
为了实现上述目的,本发明采用的技术方案为:一种估计电子件可靠性参数的方法,包括以下步骤:
步骤一:确定候选的寿命分布参数,根据已有电子件可靠性数据分布规律,初步确定待估计电子件寿命均值参数的下限值μmin和上限值μmax,在确定的参数区间内,等间隔生成n个候选分布参数,候选分布参数中,相邻参数之间的步长相等为d1,根据待测电子件均值参数的步长依次计算出n个候选均值参数μj,其中1≤j≤n;
步骤二:遍历寿命均值参数μj计算似然值Lj,对于每个寿命均值参数μj,针对一组包含m个电子元件检测信息的数据组,根据第i个检测信息所包含的在Ti时刻的状态信息Fi确定其对应的计算系数Wi,依据m个检测数据不断迭代更新似然值Lj,迭代之初,设定每个候选参数对应的似然值初始值为0,在每个候选参数所对应的迭代完毕后的似然值中,寻找最大值即LM,则最大似然值LM对应的寿命分布参数μM即为均值参数的估计值。
在上述方案中,所述步骤一中均值参数μj和步长d1的具体计算过程如下:
其中,μmax表示指数分布的均值参数上限,μmin表示指数分布的均值参数下限,n为正整数,且n≥2。
在上述方案中,所述步骤二中Wi和似然值Lj的计算公式如下:
其中,log(*)为自然对数函数,μj为指数分布的均值参数,Lj为似然值,Ti为第i个产品的检测时刻。
在上述方案中,所述步骤二中似然值Lj的遍历计算过程如下:
(1)令j=1;
(2)令i=1,Lj=0;
(4)更新i=i+1,若i≤m则转(3),否则转(5);
(5)令j=j+1,若j≤n则转(2),否则(6);
(6)在Lj(1≤j≤n)中找到最大值,记其为LM,则μM为电子件寿命均值参数的估计值。
与现有技术相比,本发明的有益效果是:利用“少量的可靠性试验数据+在产品研制、生产、使用等阶段产生的大量数据”,估算产品寿命的分布规律,避免针对产品开展大量的可靠性试验造成的人力、物力及财力的消耗。
具体实施方式
以下结合某电子件的案例,对本发明作进一步的详细描述。
本发明一种估计电子件可靠性参数的方法,包括以下步骤:
步骤一:确定候选的寿命分布参数,根据已有电子件可靠性数据分布规律,初步确定待估计电子件寿命均值参数的下限值μmin和上限值μmax,在确定的参数区间内,等间隔生成n个候选分布参数,其中,候选分布参数中,相邻参数之间的步长相等为d1,根据待测电子件均值参数的步长依次计算出n个候选均值参数μj,其中1≤j≤n;
其中,均值参数μj和步长d1的具体计算过程如下:
其中,μmax表示指数分布的均值参数上限,μmin表示指数分布的均值参数下限,n为正整数,且n≥2;
步骤二:遍历寿命均值参数μj计算似然值Lj,对于每个寿命均值参数μj,针对一组包含m个电子元件检测信息的数据组,根据第i个检测信息所包含的在Ti时刻的状态信息Fi确定其对应的计算系数Wi,依据m个检测数据不断迭代更新似然值Lj,迭代之初,设定每个候选参数对应的似然值初始值为0,在每个候选参数所对应的迭代完毕后的似然值中,寻找最大值记为LM,则最大似然值LM对应的寿命分布参数μM即为均值参数的估计值。
其中,似然值Lj的遍历计算过程如下:
(1)令j=1;
(2)令i=1,Lj=0;
(4)更新i=i+1,若i≤m则转(3),否则转(5);
(5)令j=j+1,若j≤n则转(2),否则(6);
(6)在Lj(1≤j≤n)中找到最大值,记其为LM,则μM为电子件寿命均值参数的估计值。
实施例,某电子件的[F T]型可靠性数据如下表,从工程经验可知该电子件的寿命服从指数分布,试估计其均值参数。
从以往经验得知,该电子件的均值参数在100~2600范围内,以500为步长,生成6个候选的分布参数μj,1≤j≤6,计算结果如下表:
从表中可以看出,在Lj(1≤j≤6)中的最大值为L4,则μ4=1600为该电子件寿命均值参数的估计值。
为进一步验证本发明方法的可行性,建立了如下仿真模型。
假定某电子件的寿命服从指数分布Exp(μ)。
(1)产生k1个随机数simTi(1≤i≤k1),simTi服从指数分布Exp(μ),用于模拟电子件的寿命值。令Fi=0,Ti=simTi得到k1组[Fi Ti],1≤i≤k1。
(2)产生k2个随机数simTi(k1+1≤i≤k1+k2),simTi服从指数分布Exp(μ),用于模拟电子件的寿命值。
(3)产生k2个均匀随机数simTci(k1+1≤i≤k1+k2),用于模拟检查时刻。
利用以上仿真得到的k1+k2组寿命数据[Fi Ti]后,再应用本文方法可得到分布参数的估计值。以μ=1600,k1=5,k2=15为例,大量仿真后用本文方法得到的该电子件寿命分布均值参数统计结果的均值为1581.8、根方差为417.8。如果使用k1+k2组寿命数据simTi的话,采用理论方法计算得到的均值参数统计结果的均值为1564.1、根方差为376.9,二者的差异在工程允许范围以内。
以上仅为本发明的实施方式,并非因此限制本发明的专利范围,凡是利用本发明说明书内容所作的等效结构或等效流程变换,或直接或间接运用在其他的相关技术领域,均同理包括在本发明的专利保护范围内。
Claims (3)
1.一种估计电子件可靠性参数的方法,其特征在于,针对寿命服从指数分布的电子部件开展可靠性参数估计,包括以下步骤:
步骤一:确定候选的寿命分布参数,根据已有电子件可靠性数据分布规律,初步确定待估计电子件寿命均值参数的下限值μmin和上限值μmax,在确定的参数区间内,等间隔生成n个候选分布参数,其中,n个候选分布参数中,相邻参数之间的步长相等为d1,根据待估电子件均值参数的步长依次计算出n个候选均值参数μj,其中1≤j≤n;
步骤二:遍历寿命均值参数μj计算似然值Lj,对于每个寿命均值参数μj,针对一组包含m个电子元件检测信息的数据组,根据第i个检测信息所包含的在Ti时刻的状态信息Fi确定其对应的计算系数Wi,依据m个检测数据不断迭代更新似然值Lj,迭代之初,设定每个候选参数对应的似然值初始值为0,在每个候选参数所对应的迭代完毕后的似然值中,寻找最大值即LM,则最大似然值LM对应的寿命分布参数μM即为均值参数的估计值;
所述步骤二中Wi和似然值Lj的计算公式如下:
其中Wi为计算系数,log(*)为自然对数函数,μj为均值参数,Ti为第i个产品的检测时刻。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811083876.1A CN109325289B (zh) | 2018-09-17 | 2018-09-17 | 一种估计电子件可靠性参数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811083876.1A CN109325289B (zh) | 2018-09-17 | 2018-09-17 | 一种估计电子件可靠性参数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109325289A CN109325289A (zh) | 2019-02-12 |
CN109325289B true CN109325289B (zh) | 2023-03-10 |
Family
ID=65264817
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811083876.1A Active CN109325289B (zh) | 2018-09-17 | 2018-09-17 | 一种估计电子件可靠性参数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109325289B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110309487B (zh) * | 2019-05-16 | 2023-05-16 | 中国人民解放军海军工程大学 | 基于备件保障数据的指数型单元寿命分布参数估计方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102323509A (zh) * | 2011-10-10 | 2012-01-18 | 上海电力学院 | 有机电致发光显示器加速寿命试验参数范围的确定方法 |
CN103218534A (zh) * | 2013-04-22 | 2013-07-24 | 北京航空航天大学 | 一种右截尾型寿命数据分布选择方法 |
CN104182377A (zh) * | 2014-09-02 | 2014-12-03 | 北京航空航天大学 | 一种基于β似然函数的参数估计方法 |
CN105718722A (zh) * | 2016-01-18 | 2016-06-29 | 中国人民解放军国防科学技术大学 | 基于定时截尾寿命试验数据的产品可靠度估计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8275575B2 (en) * | 2009-03-20 | 2012-09-25 | International Business Machines Corporation | Determining component failure rates using accelerated life data |
-
2018
- 2018-09-17 CN CN201811083876.1A patent/CN109325289B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102323509A (zh) * | 2011-10-10 | 2012-01-18 | 上海电力学院 | 有机电致发光显示器加速寿命试验参数范围的确定方法 |
CN103218534A (zh) * | 2013-04-22 | 2013-07-24 | 北京航空航天大学 | 一种右截尾型寿命数据分布选择方法 |
CN104182377A (zh) * | 2014-09-02 | 2014-12-03 | 北京航空航天大学 | 一种基于β似然函数的参数估计方法 |
CN105718722A (zh) * | 2016-01-18 | 2016-06-29 | 中国人民解放军国防科学技术大学 | 基于定时截尾寿命试验数据的产品可靠度估计方法 |
Non-Patent Citations (1)
Title |
---|
数据驱动的可靠性评估与寿命预测研究进展:基于协变量的方法;喻勇等;《自动化学报》;20180228(第2期);第216-227页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109325289A (zh) | 2019-02-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111008502B (zh) | 一种数字孪生驱动的复杂装备故障预测方法 | |
CN109740687B (zh) | 一种基于dlae的发酵过程故障监测方法 | |
US11003738B2 (en) | Dynamically non-gaussian anomaly identification method for structural monitoring data | |
CN109271319B (zh) | 一种基于面板数据分析的软件故障的预测方法 | |
CN110147323A (zh) | 一种基于生成对抗网络的变更智能检查方法及装置 | |
CN110297852B (zh) | 一种基于pca-粗糙集的船舶涂装缺陷知识获取方法 | |
CN108830492B (zh) | 一种基于大数据的确定抽检商家的方法 | |
CN116611747B (zh) | 一种基于云计算的资源信息智能管理系统 | |
CN107704704A (zh) | 一种基于贝叶斯理论的继电器类单机贮存可靠性评估方法 | |
CN109388861B (zh) | 一种指数型单元寿命分布参数估计方法 | |
CN109325289B (zh) | 一种估计电子件可靠性参数的方法 | |
CN110717273B (zh) | 一种工艺过程仿真边界条件构建方法 | |
CN109284478B (zh) | 一种估计对数正态型单元可靠性参数的方法 | |
CN109325287B (zh) | 一种估计机械件可靠性参数的方法 | |
WO2019000891A1 (zh) | Pcb加投率计算模型构建方法和装置 | |
CN109323876B (zh) | 一种估计伽玛型单元可靠性参数的方法 | |
CN109344472B (zh) | 一种估计机电件可靠性参数的方法 | |
CN109145502B (zh) | 一种威布尔型单元寿命分布参数估计方法 | |
CN109446481B (zh) | 一种对数正态型单元寿命分布参数估计方法 | |
Bala et al. | Use of the multiple imputation strategy to deal with missing data in the ISBSG repository | |
CN109376334B (zh) | 一种正态型单元的寿命分布参数估计方法 | |
CN109857581A (zh) | 一种优化选择软件可靠性增长模型的方法 | |
CN109388860B (zh) | 一种伽玛型单元寿命分布参数估计方法 | |
CN110309533B (zh) | 基于备件保障数据的正态型单元寿命分布参数估计方法 | |
CN113821419A (zh) | 一种基于svr和高斯函数的云服务器老化预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |