CN103218533A - 一种指数分布的串联系统可靠度置信下限估计方法 - Google Patents
一种指数分布的串联系统可靠度置信下限估计方法 Download PDFInfo
- Publication number
- CN103218533A CN103218533A CN201310139717XA CN201310139717A CN103218533A CN 103218533 A CN103218533 A CN 103218533A CN 201310139717X A CN201310139717X A CN 201310139717XA CN 201310139717 A CN201310139717 A CN 201310139717A CN 103218533 A CN103218533 A CN 103218533A
- Authority
- CN
- China
- Prior art keywords
- sample
- lambda
- distribution
- crash rate
- formula
- 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.)
- Granted
Links
Images
Landscapes
- Complex Calculations (AREA)
Abstract
本发明给出了一种指数分布的串联系统可靠度置信下限估计方法,它包括八个步骤。它首先基于信仰推断给出系统内每个单元失效率λi的信仰分布,再根据指数分布的特点,从失效率角度出发得到了系统的组成结构关系,最后通过蒙特卡罗法得到在给定置信水平下的系统失效率上限,进而得到工程上特别关注的系统可靠度置信下限。本发明提出一种求解单元寿命服从指数分布的串联系统可靠度置信下限的精确方法,计算简便,对于随机截尾数据给出了数据填充算法,扩展了使用范围。
Description
技术领域
本发明涉及一种指数分布的串联系统可靠度置信下限估计方法,它可以精确估计单元寿命服从指数分布的串联系统可靠度置信下限,适用于可靠性评估、可靠性验证等相关领域。
背景技术
电子产品的寿命一般服从指数分布,关于寿命服从指数分布的电子产品某一单元而言,工程上常表现为随机截尾样本,而对这类不完全样本估计可靠度函数较为困难,因此本发明针对随机截尾样本的可靠度函数估计具有非常重要的现实意义。
电子产品通常包含多个单元,目前评估含多个单元组成的电子产品的系统可靠性的方法有LM法、MML法,二者均属于近似方法,计算方便,在工程上得到广泛应用,但其评估结果都偏保守,评估结果的精度难以控制。
系统可靠性评估的精确方法得到的评估结果具有较好的统计性质,但这类方法计算困难,工程实用性不强,因此,探索计算方便的系统可靠性评估精确方法具有十分重要的意义,而目前关于这方面的研究较为薄弱,为此本发明提供了一种指数分布串联系统可靠度置信下限估计方法。
发明内容
本发明给出了寿命服从指数分布的随机截尾样本的系统可靠度置信下限估计方法,研究较为成熟的定时截尾样本与定数截尾样本是随机截尾样本的特殊形式。首先通过分位数填充算法将随机截尾数据补充成虚拟完全样本,然后基于信仰推断给出系统内每个单元失效率λi的信仰分布,最后根据指数分布的特点,从失效率角度出发,通过蒙特卡罗法得到在给定置信水平下的系统失效率上限,进而得到工程上特别关注的系统可靠度置信下限。
(1)发明目的
在实际研制生产过程中,寿命试验的时间过长会导致费用的增加,有时还会使试验过期而失去作用,因此缩短试验时间而获得的随机截尾样本在工程实际中较为常见,由于随机截尾样本属于不完全样本,故相关可靠度指标难以得到,其中可靠度置信下限是产品设计定型的重要指标之一,但考虑到计算方便,目前在工程上一般采用近似方法得到产品的可靠度置信下限, 基于此本发明提供了一种单元服从指数分布的随机截尾样本的串联系统可靠度置信下限估计方法,它是一种计算简便的精确方法。
(2)技术方案:
本发明一种指数分布的串联系统可靠度置信下限估计方法,其具体实施所涉及到的理论知识与前提条件如下:
①伽马分布及相关性质
形状参数为r,尺度参数为λ的Gamma分布Ga(r,λ)是一种连续概率函数,其概率密度函数为:
Gamma分布相关性质:
1.r=1时的伽马分布就是指数分布,即
Ga(1,λ)=Exp(λ)
3.伽马分布的可加性。设随机变量X~Ga(r1,λ),Y~Ga(r2,λ),,且X与Y独立,则
Z=X+Y~Ga(r1+r2,λ)。
4.伽马分布的线性性质。设Xi~Ga(αi,aiλ)(αi>0,ai>0),i=1,2,…,n,X1,X2,…,Xn相互独立,则
5.伽马分布与卡方分布的关系。设X1,X2,…,Xn是指数分布Exp(λ)的随机样本,则
②卡方分布及相关性质
若k个随机变量Z1、……、Zk是相互独立,符合标准正态分布的随机变量(数学期望为0、方差为1),则随机变量X的平方和
被称为服从自由度为k的卡方分布,记作
X~X2(k).
卡方分布的概率密度函数为:
其中x≥0,当x≤0时fk(x)=0。这里Γ代表Gamma函数。
(3)实施方案:
假设某复杂系统由m个单元串联组成串联系统,如图1所示:其寿命分别服从失效率为λ1,...λm的指数分布,则该串联系统失效率λ*为
式中m代表串联系统中的单元个数,λ1,...λm表示第1,……m个单元的失效率,λ*为该串联系统的失效率;
从指数分布F(t)=1-e-λt,(t≥0)的总体中随机抽取n个样本进行寿命试验,t为时间,在试验中观察其工作状态,若某一时刻开始某一样本不能继续工作,则该时刻记为该样本的失效时间。试验结束后,有r个样本失效,记录其失效时间依次为t1,...ti...,tr,剩下的n-r个样本
未失效,分别于τj时间有mj个撤离试验,且∑mj=n-r,其中,τj记为第j次撤离样本的时间,称为删失时间。
本发明一种指数分布的串联系统可靠度置信下限估计方法,该方法具体步骤如下:
步骤一:根据已观测到的失效时间t1,...,tr及删失时间τk利用极大似然估计得到
其中,步骤一所述的失效时间t1,...,tr是指从指数分布F(t)=1-e-λt,(t≥0)的总体中随机抽取n个样本进行寿命试验,在试验中观察其工作状态,若某一时刻开始某一样本不能继续工作,则该时刻记为该样本的失效时间,试验结束后,有r个样本失效,记录其失效时间依次为t1,...,tr,剩下的n-r个样本未失效,分别于τj时间有mj个撤离试验,且∑mj=n-r,其中,τj记为第j次撤离样本的时间,称为删失时间,k为τj的个数,mj是在τj时刻撤离试验的个数。是指数分布F(t)=1-e-λt,(t≥0)的参数λ的极大似然估计。
式中符号说明如下:由得到F(t)的反函数j为在τi时刻 撤离的mi个样本中的第j个,其他符号与步骤一一致,即k为τj的个数,mj是在τj时刻撤离试验的个数,是指数分布F(t)=1-e-λt,(t≥0)的参数λ的极大似然估计;
步骤三:基于虚拟完全样本通过信仰推断得到每个单元的失效率λi的分布Fi。其中λi为第i个单元的失效率,Fi为第i个单元的失效率的分布,具体算法如下:
故
进而根据信仰推断得到:
式中符号说明如下:寿命服从指数分布F(t)=1-e-λt,(t≥0)的单元的失效率分布,该n个试验单元的总试验时间,χ2(2n)是自由度为2n的卡方分布;
步骤四:对于每一个单元而言,利用步骤三的式子(5)可以得到各单元失效率的分布,从该失效率分布中随机抽取一个样本λi,根据式子(1)计算该串联系统失效率
式中符号说明如下:λ*为该串联系统失效率,λi是从第i个单元的失效率分布函数中随机抽取的一个失效率样本,m为该串联系统的单元总个数;
步骤七:给定1-α置信水平下,该串联系统的失效率上界为
步骤八:得到串联系统在1-α置信水平下的可靠度置信下限为
在步骤二里使用的分位数填充算法将删失样本补充成虚拟完全样本,再进行参数估计,其算法的收敛性及相合性,现有公知技术能够证明,使得系统可靠度函数估计的准确性得以提高。
在步骤三中,关于信仰推断结果的精确性在公知技术中都能证明在给定置信水平1-α下,由信仰推断得到的单元可靠度函数,具有频率意义下的实际置信水平1-α。这与经典方法得到的精确置信水平的置信界一致。因此基于信仰推断得到的失效率分布函数是精确的。
(3)本发明的优点:
i本发明针对工程实际中常见又较难处理的随机截尾寿命试验样本,提出一种解决单元寿命服从指数分布的串联系统可靠度置信下限估计的方法,是基于传统意义下定时截尾样本与定数截尾样本进行可靠性评估的一种泛化推广。
ii本发明提出的基于信仰推断与蒙特卡罗方法的系统可靠度置信下限估计是一种精确方法,根据指数分布的特点,从失效率角度出发,得到串联系统可靠度置信下限估计,简化了传统精确方法的计算困难。
附图说明
图1本发明串联系统结构图
图2随机截尾试验单元A/B/C及系统可靠度函数示意图
图3随机截尾试验系统可靠度置信下限估计流程图
图中符号说明如下:
m:串联系统的单元个数(见图1)
x1,x2,...xm:串联系统中第1,2…m个单元(见图1)
i:单元编号;(见图3)
j:循环次数编号;(见3)
N:给定需循环次数;(见图3)
Fi:单元i的失效率的信仰分布;(见图3)
λi:单元i的一个失效率样本;(见图3)
具体实施方式
见图1-图3,本发明针对工程应用较常见的随机截尾样本,提出一种单元服从指数分布的串联系统可靠度置信下限估计方法,具体实施方式如下:
已知电子产品由3个电子单元串联组成,记为单元A、B、C,分别对与3个单元同批生产出来的单元进行随机截尾试验,单元A随机抽取10个,单元B抽12个,单元C抽15个,对于单元A的寿命试验,在发生第一部失效的时候,撤下3部未失效的单元A,在第二部和第三部失效时,分别撤下2部未失效的单元,试验过程依次记录的失效及删失数据见表1-3,设三种单元寿命都服从指数分布,试计算该系统在95%的置信水平下的可靠度。
表1单元A寿命试验数据
表2单元B寿命试验数据
表3单元C寿命试验数据
解:
步骤二:根据式(2)将随机截尾数据补充成虚拟完全样本
τ1=34,τ2=113,τ3=169,m1=3,m2=2,m3=2,
待补充数据为:τ11,τ12,τ13,τ21,τ22,τ31,τ32,
其中, (j=1,...,m1)其他类似得到:
同理可以得到单元B/C的虚拟完全样本
单元B:
τ11≈1972.1,τ12≈3927.4,
τ21≈2978.3,
τ31≈4069.3,
τ41≈2983.1,τ42≈5118.4,
τ51≈3944.3
单元C:
τ11≈6346,
τ21≈4027.7,τ22≈6474,τ23≈10655.9
τ41≈7358,
τ51≈7520,
τ61≈6112.3,τ62≈10294.2
步骤三:基于以上虚拟完全样本通过信仰推断得到每个单元的失效率λi的分布Fi
步骤七:根据式(6),取1-α=95%置信水平下,系统的失效率上界为:
步骤八:故根据式(7)得到该电子产品在1-α=95%置信水平下的可靠度函数置信下限为:
Claims (1)
1.一种指数分布的串联系统可靠度置信下限估计方法,其特征在于:该方法具体步骤如下:
其中,所述的失效时间t1,...,tr是指从指数分布F(t)=1-e-λt,(t≥0)的总体中随机抽取n个样本进行寿命试验,在试验中观察其工作状态,若某一时刻开始某一样本不能继续工作,则该时刻记为该样本的失效时间,试验结束后,有r个样本失效,记录其失效时间依次为t1,...,tr,剩下的n-r个样本未失效,分别于τj时间有mj个撤离试验,且∑mj=n-r;其中,τj记为第j次撤离样本的时间,称为删失时间;k为τj的个数,mj是在τj时刻撤离试验的个数,是指数分布F(t)=1-e-λt,(t≥0)的参数λ的极大似然估计;
步骤二:通过分位数填充算法将随机截尾数据按下式补充成虚拟完全样本其中t1,...,tr即在步骤一中表示的失效时间,τij表示在τi时刻撤离的第j个样本经补充后的虚拟失效时间,按下式得到:
式中符号说明如下:由得到F(t)的反函数j为在τi时刻撤离的mi个样本中的第j个,其他符号与步骤一一致,即k为τj的个数,mj是在τj时刻撤离试验的个数,是指数分布F(t)=1-e-λt,(t≥0)的参数λ的极大似然估计;
步骤三:基于虚拟完全样本通过信仰推断得到每个单元的失效率λi的分布Fi;其中λi为第i个单元的失效率,Fi为第i个单元的失效率的分布,具体算法如下:
故
进而根据信仰推断得到:
式中符号说明如下:寿命服从指数分布F(t)=1-e-λt,(t≥0)的单元的失效率分布,该n个试验单元的总试验时间,χ2(2n)是自由度为2n的卡方分布;
步骤四:对于每一个单元而言,利用步骤三的式子(5)得到各单元失效率的分布,从该失效率分布中随机抽取一个样本λi,根据式子(1)计算该串联系统失效率
式中符号说明如下:λ*为该串联系统失效率,λi是从第i个单元的失效率分布函数中随机抽取的一个失效率样本,m为该串联系统的单元总个数;
式中符号说明如下:N为计算系统失效率的总循环次数;j是第j次计算系统失效率,故j=1,2,...N;为第j次计算得到的系统失效率的值;
步骤七:给定1-α置信水平下,该串联系统的失效率上界为
步骤八:得到串联系统在1-α置信水平下的可靠度置信下限为
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310139717.XA CN103218533B (zh) | 2013-04-22 | 2013-04-22 | 一种指数分布的串联系统可靠度置信下限估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310139717.XA CN103218533B (zh) | 2013-04-22 | 2013-04-22 | 一种指数分布的串联系统可靠度置信下限估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103218533A true CN103218533A (zh) | 2013-07-24 |
CN103218533B CN103218533B (zh) | 2015-11-04 |
Family
ID=48816312
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310139717.XA Expired - Fee Related CN103218533B (zh) | 2013-04-22 | 2013-04-22 | 一种指数分布的串联系统可靠度置信下限估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103218533B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104462757A (zh) * | 2014-11-03 | 2015-03-25 | 南京航空航天大学 | 基于监测数据的Weibull分布可靠性序贯验证试验方法 |
CN104899423A (zh) * | 2015-05-06 | 2015-09-09 | 同济大学 | 一种动车组子系统关键部件运用可靠性评估方法 |
CN105069214A (zh) * | 2015-07-31 | 2015-11-18 | 北京航空航天大学 | 一种基于非线性相关分析的工艺可靠性评估方法 |
CN105373687A (zh) * | 2014-08-18 | 2016-03-02 | 鲍珂 | 一种基于多阶段试验的重型车辆动力系统可靠性评估方法 |
CN105718722A (zh) * | 2016-01-18 | 2016-06-29 | 中国人民解放军国防科学技术大学 | 基于定时截尾寿命试验数据的产品可靠度估计方法 |
CN106202938A (zh) * | 2016-07-14 | 2016-12-07 | 熊猫电子集团有限公司 | 一种通信装备可靠性评估方法 |
CN107132834A (zh) * | 2017-04-12 | 2017-09-05 | 北京航空航天大学 | 基于统计量排序的和联系统可靠度精确置信下限估计方法 |
CN108334667A (zh) * | 2018-01-09 | 2018-07-27 | 北京航空航天大学 | 一种基于环境折合的加速可靠性试验方案设计方法 |
CN111007418A (zh) * | 2019-12-30 | 2020-04-14 | 电子科技大学 | 一种基于扩展性指数分布的锂电池剩余寿命预测方法 |
CN112001077A (zh) * | 2020-08-21 | 2020-11-27 | 中国石油大学(华东) | 一种基于非完全数据的石化安全关键设备寿命的评估方法 |
CN112528505A (zh) * | 2020-12-14 | 2021-03-19 | 西南交通大学 | 一种指数分布型产品可靠性评估方法 |
CN113064112A (zh) * | 2021-03-23 | 2021-07-02 | 广东电网有限责任公司计量中心 | 一种智能电能表使用寿命评估方法及系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102243695A (zh) * | 2011-08-05 | 2011-11-16 | 北京航空航天大学 | 零膨胀泊松分布参数的一种稳健估计方法 |
-
2013
- 2013-04-22 CN CN201310139717.XA patent/CN103218533B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102243695A (zh) * | 2011-08-05 | 2011-11-16 | 北京航空航天大学 | 零膨胀泊松分布参数的一种稳健估计方法 |
Non-Patent Citations (4)
Title |
---|
DAVID W. COIT ET AL.: "System-Reliability Confidence-Intervals for Complex-Systems With Estimated Component-Reliability", 《IEEE TRANSACTIONS ON RELIABILITY》, vol. 46, no. 4, 31 December 1997 (1997-12-31) * |
TENG-FA TSAO ET AL.: "Composite Reliability Evaluation Model for Different Types of Distribution Systems", 《IEEE TRANSACTIONS ON POWER SYSTEMS》, vol. 18, no. 2, 31 May 2003 (2003-05-31), XP011096395, DOI: doi:10.1109/TPWRS.2002.807095 * |
杨军 等: "复杂系统平均剩余寿命综合评估方法", 《航空学报》, vol. 28, no. 6, 15 November 2007 (2007-11-15) * |
胡思贵: "指数分布区间型删失数据的可靠度最优置信下限", 《贵州大学学报(自然科学版)》, vol. 24, no. 6, 30 November 2007 (2007-11-30) * |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105373687A (zh) * | 2014-08-18 | 2016-03-02 | 鲍珂 | 一种基于多阶段试验的重型车辆动力系统可靠性评估方法 |
CN104462757B (zh) * | 2014-11-03 | 2017-10-31 | 南京航空航天大学 | 基于监测数据的Weibull分布可靠性序贯验证试验方法 |
CN104462757A (zh) * | 2014-11-03 | 2015-03-25 | 南京航空航天大学 | 基于监测数据的Weibull分布可靠性序贯验证试验方法 |
CN104899423A (zh) * | 2015-05-06 | 2015-09-09 | 同济大学 | 一种动车组子系统关键部件运用可靠性评估方法 |
CN104899423B (zh) * | 2015-05-06 | 2017-12-15 | 同济大学 | 一种动车组子系统关键部件运用可靠性评估方法 |
CN105069214A (zh) * | 2015-07-31 | 2015-11-18 | 北京航空航天大学 | 一种基于非线性相关分析的工艺可靠性评估方法 |
CN105069214B (zh) * | 2015-07-31 | 2018-02-16 | 北京航空航天大学 | 一种基于非线性相关分析的工艺可靠性评估方法 |
CN105718722A (zh) * | 2016-01-18 | 2016-06-29 | 中国人民解放军国防科学技术大学 | 基于定时截尾寿命试验数据的产品可靠度估计方法 |
CN106202938A (zh) * | 2016-07-14 | 2016-12-07 | 熊猫电子集团有限公司 | 一种通信装备可靠性评估方法 |
CN107132834A (zh) * | 2017-04-12 | 2017-09-05 | 北京航空航天大学 | 基于统计量排序的和联系统可靠度精确置信下限估计方法 |
CN107132834B (zh) * | 2017-04-12 | 2019-03-08 | 北京航空航天大学 | 基于统计量排序的和联系统可靠度精确置信下限估计方法 |
CN108334667A (zh) * | 2018-01-09 | 2018-07-27 | 北京航空航天大学 | 一种基于环境折合的加速可靠性试验方案设计方法 |
CN108334667B (zh) * | 2018-01-09 | 2019-12-20 | 北京航空航天大学 | 一种基于环境折合的加速可靠性试验方案设计方法 |
CN111007418A (zh) * | 2019-12-30 | 2020-04-14 | 电子科技大学 | 一种基于扩展性指数分布的锂电池剩余寿命预测方法 |
CN111007418B (zh) * | 2019-12-30 | 2020-12-29 | 电子科技大学 | 一种的锂电池剩余寿命预测方法 |
CN112001077A (zh) * | 2020-08-21 | 2020-11-27 | 中国石油大学(华东) | 一种基于非完全数据的石化安全关键设备寿命的评估方法 |
CN112528505A (zh) * | 2020-12-14 | 2021-03-19 | 西南交通大学 | 一种指数分布型产品可靠性评估方法 |
CN112528505B (zh) * | 2020-12-14 | 2022-03-25 | 西南交通大学 | 一种指数分布型产品可靠性评估方法 |
CN113064112A (zh) * | 2021-03-23 | 2021-07-02 | 广东电网有限责任公司计量中心 | 一种智能电能表使用寿命评估方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN103218533B (zh) | 2015-11-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103218533A (zh) | 一种指数分布的串联系统可靠度置信下限估计方法 | |
Canty et al. | Package ‘boot’ | |
CN103942457B (zh) | 基于关联向量机回归的水质参数时间序列预测方法 | |
Zheng et al. | Shifted Gamma-Generalized Pareto Distribution model to map the safety continuum and estimate crashes | |
CN107301118B (zh) | 一种基于日志的故障指标自动标注方法与系统 | |
CN102360332B (zh) | 一种软件可靠性加速测试与评估方法及其计算机辅助工具 | |
Tshimanga | Hydrological uncertainty analysis and scenario-based streamflow modelling for the Congo River Basin | |
CN105843733B (zh) | 一种大数据平台的性能检测方法及装置 | |
CN102945310A (zh) | 基于自治计算的流行病传播网络建模与推断 | |
CN106442271A (zh) | 岩心渗透率模拟方法及装置 | |
CN103679540A (zh) | 一种基于储层构型要素界面的储层模型的粗化方法及系统 | |
CN103793599A (zh) | 一种基于隐马尔科夫模型的出行异常检测方法 | |
CN106324672A (zh) | 一种岩石物理建模方法及系统 | |
CN109284477A (zh) | 一种水文序列的丰枯组合概率计算方法及装置 | |
CN105093351B (zh) | 识别储层微裂缝的方法 | |
CN108205713A (zh) | 一种区域风电功率预测误差分布确定方法和装置 | |
CN104712328A (zh) | 快速评价复杂油藏中单个流动单元动用状况的方法 | |
CN108734218A (zh) | 一种多传感器系统的信息融合方法和装置 | |
CN106569981A (zh) | 一种适用于大规模数据集的统计参量确定方法及系统 | |
CN106125713A (zh) | 一种区间删失情况下可靠性增长的评估与预测方法 | |
CN105046089A (zh) | 一种预测强降雨及洪涝灾害的方法 | |
CN106777928A (zh) | 面向正态分布数据样本的Bayes验前分布构造方法 | |
CN100420209C (zh) | 自动进行方案对比的可信网络仿真系统 | |
Verdon et al. | Long-term drought risk assessment in the Lachlan River Valley–A paleoclimate perspective | |
Mohaghegh | Surrogate reservoir model |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20151104 Termination date: 20180422 |
|
CF01 | Termination of patent right due to non-payment of annual fee |