常用软件可靠性模型推导

常用软件可靠性模型推导

本章针对软件可靠性IEEE P1633标准给出的模型参数的极大似然估计和最小二乘估计的详细推导，给出了求解公式。

随机过程类的软件可靠性数学模型主要包括马尔科夫过程模型(Markov Process Model)和非齐次泊松过程模型(NHPP).一般假定错误出现率在软件无改动的区间内是常数，并且随着错误数据的减少而下降，这样的模型数据马尔科夫过程模型，如Jelinski_Moranda 模型.另外，排错过程中的累积错误数作为时间的函数N(t)在一定的条件下可以近似为一个非齐次泊松过程，这一类的数学模型属于非齐次泊松过程模型。如Goel_Okumoto 模型，Schneiwind 模型等.

另外本章还讨论了一个非随机过程类模型Littlewood －Verrall 模型，L_V 模型应用贝叶斯方法研究软件可靠性。

对于大的样本，极大似然法是非常有效的估计方法，但只针对小样本或者中等大小的样本，用最小二乘法比较好。下面将针对各个模型给出具体的参数估计推到过程。 1

指数模型

1.1 指数模型简介与假设 1.2 指数模型推导 ,t R e λλ-=其中为常数

则有

1

λλ

m(t)=t,MTBF=

1.3 指数模型参数估计 (1) 数据要求：

测试时间： i t ,00=t ；累计失效数：i n (i t 时刻对应失效数)。 (2) 参数点估计：

测试终止时刻测试时间为f t ，累积失效数为f n ，则参数估计值为：

f f

n t λ=

2 Jelinski_Moranda(J_M)模型

2.1 J_M 模型简介与假设

由Jelinski －Moranda 开发的可靠性模型是最早建立且现在仍然使用着的模型之一，该模型现在正用在麦克唐奈道格拉斯海军工程中。它是最具代表性的早期软件可靠性马尔可夫过程的数学模型。随后的许多工作，都是在它的基础上，

对其中与软件开发实际不相适合的地方进行改进而提出的。因此，在这个意义上来说，JM 模型又是对后面的工作有着广泛影响的模型之一。

模型假设为：

(1)缺陷检测率（即单位时间内被发现的软件缺陷数，简称缺陷率）与程序当前的残留缺陷数成正比；

(2)所有缺陷导致失效的可能性是相同的，每个缺陷的等级相同 (3)缺陷被查出时，失效是相互独立的； (4)缺陷率在相邻失效时间间隔内保持不变

(5)软件（测试）运行方式与预测的（实际）运行剖面相同

(6)每一个缺陷一旦经检测发现即被瞬时剔除，且不引入新的缺陷； 2.2 J_M 模型推导

假设(2)确保不同的失效强度具有相同的分布特性，假设(3)可以简化推导过程，假设(5)确保在某一特定的环境下使用数据采集进行模型评价的正确性在软件的第i-1次失效与第i 次失效的时间间隔内，失效强度函数相同:

))1((/10--=i N t φ (1-1)

式中，0N -----比率常数, 0N ----观测开始之后软件中的缺陷总数，t 为平均失效间隔我们先推导时间t 所经历的故障数M(t)。由于J-M 模型为二项型式模型所以

M(t)＝m 的概率服从二项式分布，即

()[]()[]()[]m

u a m

a m

u t F t F C m t M P --==00

1(1-2)

其中()t F a 是累计分布函数，0u 是初始时刻，软件潜在的总失效数。它的均值函数为

()()t F u np t u a 0==，

(1-3)

对()t u 微分可得故障强度

()()t f u t a 0=λ

(1-4)

由于

()()[]()

[]()[]

11

11

11|1|1|0-+--+-+---=-==i j

a

i j i j a i j i u i j t t

F t t F C i t M i t M P 其中

()1|-i j a t t F ＝

()()()

111----i a i a j a t F t F t F

所以

[]()()[]

∑=----====≤0

1||111u i

j i j i i i i i t M j t M P t T t T P

()[]1

10

|11+----=i u i i a t t F

(1-5)

又

()()()⎥⎦⎤⎢⎣⎡--=-=⎰t

a a dx x Z t R t F 0exp 11

(1-6)

所以

()()⎥⎦⎤⎢⎣⎡--=⎰t

t a e a e dx x Z t t F exp 1| (1-7)

()()()()⎥⎦

⎤⎢⎣⎡+--=⎥⎦

⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡-==>⎰⎰--+---i

i i

i t t a i u t t a i i i i dx x Z i u dx x Z t T t T P 1011ex p ex p |01

11

(1-8)

又

()

[]()()⎥⎦⎤⎢⎣⎡+--==>=⎰+-----'

11

1exp ||011

1'i

i i t t t a i i i i i i

dx x Z i u t T t T P t t R (1-9)

所以

(

)

()

()

()()

'10'

1'1'1

'

1||1

|i i a i i i i i i i

t t Z i u dt t t dR t t R t t Z ++-=-=---- (1-10)

由J-M 模型的假设知，J-M 模型是指数类模型，即它的单位错误故障率是常数，用Φ表示，()()Φ+-=∴-1|01'i u t t Z i i J-M 模型的概率密度函数

()()[]()[]()

1

''

'

1

'

'

1

'1'||1||-----=-==i i

i i

i i a i i a t t

R t t

R t t F t t f (1-11)

又

()[][]()()[]

'00'

'01'1exp 11(exp |i i i i a t i u i u t i u t t f Φ+--Φ+-=Φ+---=-

()[][]

'

'

01(exp i i t i u t f Φ+---=

()()[]

'

001ex p )1(i t i u i N Φ+--Φ--=

由失效率函数为：0((1))N i λφ=-- 则累积失效均值函数

0()((1))m t N i t

φ=--

2.3 J_M 模型参数估计-最大似然法 (3) 数据要求：

软件发生失效时刻： i t ,00=t ；软件失效间隔序列：i x =i t -1-i t