NHPP软件可靠性模型的推广

  • 格式:pdf
  • 大小:190.84 KB
  • 文档页数:3

第22卷第1期湖 北 工 业 大 学 学 报2007年02月Vol.22No.1 JournalofHubeiUniversityofTechnologyFeb.2007

[收稿日期]2006-10-30[作者简介]钱丹青(1984-),女,江苏常州人,华中科技大学硕士研究生,研究方向:应用概率统计.

[文章编号]1003-4684(2007)0120004203NHPP软件可靠性模型的推广

钱丹青,刘次华(华中科技大学数学系,湖北武汉430074)

[摘 要]在一般的意义上对一类非齐次泊松过程模型作了推广,用多项式函数来拟合纠错经验函数,在一定程度上简化了模型的选择.给出了上述推广的可行性的证明,并对模型参数的性质进行了讨论.

[关键词]软件可靠性;非齐次泊松过程;G-ONHPP模型;极大似然估计[中图分类号]O213.2[文献标识码]:A

软件可靠性模型可以用随机过程的理论来描述,人们通常可以将软件可靠性或与软件可靠性直接有关的量,如:可靠度或故障率等,表示成时间以及软件产品的特性,或者开发过程的函数.软件可靠性模型通常给出了软件可靠性对上述各变量的一种依赖关系.非齐次泊松过程(NHPP)模型就是用非齐次泊松过程来描述软件的失效过程.NHPP类软件可靠性增长模型已成为可靠性工程实践中非常重要的工具[1].软件的可靠度增长[2]是指消耗的测试时间与累积错误数m(t)的关系.此对应关系可以用一条增长曲线来描述,在实际观测中根据曲线形状,可分为两种类别,指数型增长和S字母型增长,由此得出NHPP模型分别称为指数型可靠度和S字母型可靠度增长模型.两类模型的区别在于指数模型认为软件故障检出率仅决定于系统中残存的错误数,而S字母模型还认为检出率是时间的函数.G-ONHPP模型[3]的m(t)=a(1-e-bt),就是典型的指数模型.而威布尔模型(m(t)=a(1-e-btc))以及延迟S字母型模型(m(t)=a(1-(1+bt)e-bt))就是S字母型增长模型.在实际应用中,错误检出率与时间的关系通常是未知的,笔者给出了在这种情况下建立可靠度模型的一个一般方法.1 NHPP软件可靠性模型的推广用经验函数f(t)来表示错误检出率与时间的关系.则上面介绍的模型即是根据不同经验函数得出的.笔者选用多项式函数来拟合经验函数,可以在对经验函数没有充分的认识的情况下对这一类问题建立可靠性模型,从而简化了NHPP可靠性模型的选择.

1.1 模型假设1)在任何不相交的时间区间中检测到的错误数

相互独立.

2)t时刻累积错误数N(t)是服从参数为m(t)

的非齐次泊松过程,t时刻错误检出率与剩余错误期望值成正比,且是f(t)的函数.

3)累积错误数的期望函数m(t)满足边界条件:m(0)=0,limt→∞m(t)=a.其中,a是最终可能被检测出的错误总数的期望值.

4)纠错不会引入新的错误.由假设可知下列微分方程成立:

m′(t)=b[a-m(t)]f(t) (b>0). 笔者采用多项式f

n(t)

=anxn+an-1xn-1+…+

a0来近似可积函数f(t),

通过极大似然法估计多项

式的系数,从而得到一个近似的失效模型.

1.2 经验函数f(t)的多项式逼近引理1[4] 若f∈L(E)(即f在可测集E上Lebesgue可积),则对任给的ε>0,存在R

n上具有

紧支集的连续函数g(x)(即函数的非零定义域的子集是紧集),使得

∫E|f(x)-g(x)|

dx)

定理1 存在一个一致收敛于可积函数f(t)的多项式序列{fn(t)},使得:mn(t)一致收敛于m(t).其中m(t),mn(t)分别是下列微分方程的解:m′(t)=b(a-m(t))f(t),(2)m′n(t)=b(a-mn(t))fn(t).(3)其中z(t)是式(1)中在积分意义下逼近f(t)的连续函数.证明 由引理1,可知存在连续函数z(t)在积分意义下逼近可积函数f(t),故对Πε>0,∫dc|f(t)-z(t)|dt0,当n>N时,对Πt有|fn(t)-z(t)|N时∫dc|f(t)-fn(t)|dt<∫dc|f(t)-z(t)|dt+∫dc|z(t)-fn(t)|dt0,ϖN,当n>N时,|m(t)-mn(t)|错误累积函数所得的参数估计值为^

θ

n,多项式逼近

所得的参数估计值为^

θ

nm,其中m为多项式次数,n

为样本容量.

定理2 存在M,当m>M时,在大样本条件下,

^θnm

a.sθ (m,n→∞)

.

证明 两似然函数相减得

|lnLn-lnL|≤∑nk=1(yk-yk-1)・

|exp(-b∫tk0f(s)ds)-exp(-b∫tk-10f(s)ds)exp(-b∫tk0fm(s)ds)-exp(-b∫tk-10fm(s)ds)|+

a|exp(-b∫tn0f(s)ds)-exp(-b∫tn0fm(s)ds)|.因为exp(-b∫tk0f(s)ds)

exp(-b∫tk0fm(s)ds)

=

exp(-b∫tk0(f(s)-fm(s))ds).

由f

m(s)→f(s)(m→∞)可知

exp(-b∫tk0f(s)ds)

exp(-b∫tk0fm(s)ds)

→1 (m→∞)]

exp(-b∫tk0f(s)ds)-

exp(-b∫tk0fm(s)ds)→0 (m→∞)

对Πk成立.所以

lnexp(-b∫tk0f(s)ds)-exp(-b∫tk-10f(s)ds)exp(-b∫tk0fm(s)ds)-exp(-b∫k-10fm(s)ds)=ln

[1+exp(-b∫tk0f(s)ds)-exp(-b∫tk0fm(s)ds)exp(-b∫tk0fm(s)ds)-exp(-b∫tk-10fm(s)ds)+

exp(-b∫tk-10fm(s)ds)-exp(-b∫tk-10f(s)ds

)

exp(-b∫tk0fm(s)ds)-exp(-b∫tk-10fm(s)ds

)

]

→0 (m→∞)

,

exp(-b∫tn0f(s)ds)-

exp(-b∫tn0fm(s)ds)→0 (m→∞)

.

故lnLm-lnL~o(1/n),由似然函数的一致收敛性显然可以得出^

θ

mn→^θn(m→∞).

由于失效率函数

5 第22卷第1期 钱丹青等 NHPP软件可靠性模型的推广λ(t,θ)=m′(t,θ)可识别,故极大似然估计有相合性,即^θna.sθ(m→∞).所以^θnma.sθ(m→∞).上述证明可知用多项式逼近经验函数在大样本情形下的参数估计也具有相合性,此外,在大多数情况下MLE还具有渐近正态性等优良性质,是大样本情况下统计中常用的一种参数估计方式.2 模拟设g(t)=ct2+dt+e.由极大似然估计的不变性,m(t)=a(1-exp(-bt3-ct2-dt)).lnLn=∑ni=1(ni-ni-1)[lna+ln(exp(-bt3i-1-ct2i-1-dti-1)-exp(-bt3i-ct2i-dti))]-a(1-exp(-bt3n-ct2n-dtn))-∑ni=1ln[(ni-ni-1)!].用数值算法解出参数,得到的拟合曲线(图1).图1 拟合曲线图其回归曲线的相关指数R2=1-∑ni=1(ni-^ni)2∑ni=1(ni-󰁦ni)2=0.9964[5]而用威布尔模型

m(t)=a(1-exp(-βtα))[6]

拟和的相关指数R

2

=0.9961,可知笔者所选用的

模型具有更好的精度.

3 结束语软件可靠性模型的数量多,应用时具有一定的灵活性,但是却使得模型选择成为了一个重大的难题,此外模型的复杂性也给参数估计提出了不小的挑战.笔者所给的模型,利用多项式在更一般的意义下推广了NHPP模型,从而简化了模型选择和参数估计.

[ 参 考 文 献 ][1] HoangPham,LarsNordmann,XuemeiZhang.AGeneralImperfect2Software2DebuggingModelwithS2ShapedFault2DetectionRate[J].IEEETransonReli2ability,1999,48(2):169-175.[2] 赵晓华.计算机软件可靠性与质量管理[M].北京:中国经济出版社,1992:111-164.

[3] GoelAL,OkumotoK.Time2dependentError2detectionRateModelforSoftwareReliabilityandotherPer2formanceMeasures[J].IEEETrans.Rel,1979,R228(3):206-211.

[4] 周民强.实变函数论[M].北京:北京大学出版社,2001:181-188.[5] 刘宏伟,杨孝宗,曲 峰,等.非齐次泊松过程类软件可靠性增长模型[J].同济大学学报(自然科学版),2004,

32(8):1037-1038.[6] 黄锡滋.软件的可靠性与安全性[M].北京:科学出版社,1993:120-123.

TheGeneralizationofNHPPSoftwareReliabilityModelingQIANDan2qing,LIUCi2hua(Dep.ofMathematics,HuazhongUniv.ofSci.&Tech.,Wuhan430074,China)

Abstract:NHPP(Non2homogeneousPoissonProcess)softwarereliabilitymodelisanimportantsoftwarereliabilitymodel.ThispapergivesageneralizationofanNHPPmodel,usesapolynomialtosimulateanexperientialfunction,sopredigeststheselectionofdifferentmodelstoacertainextent.Thefeasibilityofthegeneralizationhasbeenproved.Thecharactersofitsparametersarediscussedattheend.Keywords:softwarereliability;non2homogeneousPoissonprocess;G2ONHPPmodel;maximumlikeli2hoodestimate[责任编辑:张培炼]

6湖 北 工 业 大 学 学 报2007年第1期