法律状态公告日
法律状态信息
法律状态
2019-04-05
未缴年费专利权终止 IPC(主分类):G06F19/00 授权公告日:20151104 终止日期:20180422 申请日:20130422
专利权的终止
2015-11-04
授权
授权
2013-08-21
实质审查的生效 IPC(主分类):G06F19/00 申请日:20130422
实质审查的生效
2013-07-24
公开
公开
技术领域
本发明涉及一种指数分布的串联系统可靠度置信下限估计方法,它可以精确估计单元寿命服从指数分布的串联系统可靠度置信下限,适用于可靠性评估、可靠性验证等相关领域。
背景技术
电子产品的寿命一般服从指数分布,关于寿命服从指数分布的电子产品某一单元而言,工程上常表现为随机截尾样本,而对这类不完全样本估计可靠度函数较为困难,因此本发明针对随机截尾样本的可靠度函数估计具有非常重要的现实意义。
电子产品通常包含多个单元,目前评估含多个单元组成的电子产品的系统可靠性的方法有LM法、MML法,二者均属于近似方法,计算方便,在工程上得到广泛应用,但其评估结果都偏保守,评估结果的精度难以控制。
系统可靠性评估的精确方法得到的评估结果具有较好的统计性质,但这类方法计算困难,工程实用性不强,因此,探索计算方便的系统可靠性评估精确方法具有十分重要的意义,而目前关于这方面的研究较为薄弱,为此本发明提供了一种指数分布串联系统可靠度置信下限估计方法。
发明内容
本发明给出了寿命服从指数分布的随机截尾样本的系统可靠度置信下限估计方法,研究较为成熟的定时截尾样本与定数截尾样本是随机截尾样本的特殊形式。首先通过分位数填充算法将随机截尾数据补充成虚拟完全样本,然后基于信仰推断给出系统内每个单元失效率λi的信仰分布,最后根据指数分布的特点,从失效率角度出发,通过蒙特卡罗法得到在给定置信水平下的系统失效率上限,进而得到工程上特别关注的系统可靠度置信下限。
(1)发明目的
在实际研制生产过程中,寿命试验的时间过长会导致费用的增加,有时还会使试验过期而失去作用,因此缩短试验时间而获得的随机截尾样本在工程实际中较为常见,由于随机截尾样本属于不完全样本,故相关可靠度指标难以得到,其中可靠度置信下限是产品设计定型的重要指标之一,但考虑到计算方便,目前在工程上一般采用近似方法得到产品的可靠度置信下限, 基于此本发明提供了一种单元服从指数分布的随机截尾样本的串联系统可靠度置信下限估计方法,它是一种计算简便的精确方法。
(2)技术方案:
本发明一种指数分布的串联系统可靠度置信下限估计方法,其具体实施所涉及到的理论知识与前提条件如下:
①伽马分布及相关性质
形状参数为r,尺度参数为λ的Gamma分布Ga(r,λ)是一种连续概率函数,其概率密度函数为:
Gamma分布相关性质:
1.r=1时的伽马分布就是指数分布,即
Ga(1,λ)=Exp(λ)
2.的伽马分布是自由度为n的χ2分布,即
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利用极大似然估计得到
对于参数为λ指数分布F(t)=1-e-λt,(t≥0)而言,按下面公式可以得到上述指数法分布的参数λ的极大似然估计
其中,步骤一所述的失效时间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次,每次计算得到的该串联系统失效率分别记为式中符号说明如下:N为计算系统失效率的总循环次数;j是第j次计算系统失效率,故 j=1,2,...N;为第j次计算得到的系统失效率的值;
步骤六:对从小到大排列,分别记为
式中符号说明如下:是N个系统失效率的值从小到大排列在各位次上的值,例如代表N个系统失效率的值从小到大排在第2位次上的值;
步骤七:给定1-α置信水平下,该串联系统的失效率上界为
式中符号说明如下:[N(1-α)]代表向下取整,为该串联系统的失效率上界;
步骤八:得到串联系统在1-α置信水平下的可靠度置信下限为
式中符号说明如下:为该串联系统在1-α置信水平下的可靠度置信下限,t为时间,其它与步骤七一致,为该串联系统的失效率上界。
在步骤二里使用的分位数填充算法将删失样本补充成虚拟完全样本,再进行参数估计,其算法的收敛性及相合性,现有公知技术能够证明,使得系统可靠度函数估计的准确性得以提高。
在步骤三中,关于信仰推断结果的精确性在公知技术中都能证明在给定置信水平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)
:第j次循环得到的系统失效率;(见图3)
:目前观测到的随机截尾样本下单元失效率;(见图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寿命试验数据
解:
步骤一:根据已观测到的失效时间t1,...,tr利用极大似然估计得到
步骤二:根据式(2)将随机截尾数据补充成虚拟完全样本
单元A的初始分布函数估计(t≥0)由题意可得:
τ1=34,τ2=113,τ3=169,m1=3,m2=2,m3=2,
待补充数据为:τ11,τ12,τ13,τ21,τ22,τ31,τ32,
其中,
同理可以得到单元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
步骤四:分别从A/B/C三个分布λa,λb,λc中随机抽取一个样本计算系统失效率
步骤五:重复步骤四N=10000次,每次计算得到的系统失效率分别记为步骤六:对从小到大排列,分别记为
步骤七:根据式(6),取1-α=95%置信水平下,系统的失效率上界为:
步骤八:故根据式(7)得到该电子产品在1-α=95%置信水平下的可靠度函数置信下限为:
机译: 实例基推理产生的重合的可靠度估计系统和方法,故障机分析产生的错误日志的系统和方法
机译: 无线通信系统中基于可靠度概念的CFO估计方法
机译: 一种估计位置的方法,一种广播近邻基站的位置信息的方法以及一种协商估计位置能力的方法,可确定位置估计所需的参数