软件可靠性模型算法分析与评价

尹晶杰：软件可靠性模型的算法分析与评价

软件可靠性模型算法分析与评价

尹晶杰

摘要：本文首先对三个经典的软件可靠性模型(J-M、G-O、S-W)进行参数计算的数值算法设计，在此基础上通过可靠性数学关系得到失效间隔时间的密度函数、分布函数、可靠性函数以及失效率函数。其次，分别采用未确知模型、J-M模型、G-O模型、S-W模型针对具体实例的失效时间进行预测评估，并对各模型的可靠性评估曲线进行描绘。第三，利用Delphi开发软件设计并开发完成了一个简易的软件可靠性评估工具。该工具嵌入了包括未确知模型在内的四个软件可靠性模型（J-M模型、G-O模型、S-W模型），能够输出模型评估结果和评估曲线，并具有计算各模型评价准则值（KS值、PL值、模型噪声）和绘制用于模型评价的PLR图、u?结构图、关键词：软件可靠性；软件可靠性模型

Abstract: Here originally in the paper, model in three software reliability model (J-M, G-O , S-W ) at first ask the parameter algorithm to be designed , draw the parameter of each model, receive invalid density function of spacing interval , distribute function , reliability function and software failure rate function through dependability mathematics relation on this basis.Secondly, on the basis of the above function, including software reliability based on unascertained theory model , J-M model , G-O model , S-W model predict the assessment to the failure time of the concrete instance separately, design through Delphi one simple interface describe to every reliability assessment curve of model.Moreover , utilize Delphi to designed and development a simple software reliability estimation tool. This tool inlayed three traditional software reliability models: J-M model, G-O model, S-W model and the new model put forward in this paper. It is not only can exports the estimation results but also can provides assessment curve , including calculating every model appraise criterion value (KS value , PL value , model noise ) and PLR chart that is used to model comparison, u-plot and the y- plot. Keywords: Software reliability Software reliability model

y?结构图的功能。

1. 基本概念

1.1 软件可靠性的定义

关于软件可靠性的确切含义，学术界有过长期的争论，经过长期的争论和研究，1983年美国IEEE计算机学会对 “软件可靠性”一次正式做出如下定义：

(1) 在规定条件下，在规定的时间内，软件不引起系统失效的概率，该概率是系统输入和系统使用的函数，也是软件中存在的错误的函数；系统输入将确定是否会遇到已存在的错误（如果错误存在的话）；

(2) 在规定的时间周期内，在所述条件下程序执行所要求的功能的能力。

1.2 软件可靠性参数

下面对几个主要的软件可靠性参数进行介绍： (1) 可靠度

软件可靠度R是指软件在规定的条件下、规定的时间段内完成预定的功能的概率。或者说是软件在规定时间内无失效发生的概率。

用随机变数?表示从软件运行开始到系统失效所经历的时间，用F?(t)表示?的分布函数，用t表示任意给定的时刻，用R?(t)表示软件在t时刻的可靠度，则数学公式如下：

R?(t)?Pr???t??1?F?(t) (1-1)

(2) 失效率

139

河北工程学院优秀毕业设计（论文）

失效率是指软件在t时刻没有发生失效的条件下，在(t??t)区间内，当?t很小时，单位时间内发生失效的概率。

用随机变数?表示从软件运行开始到系统失效所经历的时间，用t表示任意给定的时刻，用?(t)表示失效率，则失效率数学公式为：

?(t)?limPr?t??t???t/??t??t

?t?0 (1-2)

(3) 平均失效间隔时间MTBFMTBF是指两次相邻失效时间间隔的均值。假设当两次相邻失效时间间隔为?，?具有累计概

率密度函数F(t)?P(??t)，即可靠度函数

R(t)?1?F(t)?P(??t)

?则 MTTF??R(t)dt (1-4)

02. 系统设计

2.1设计思想

整个软件可靠性评估工具的设计采取模块化思想，一个模型设计为一个模块，都具有数据输入和评估结果（包括模型中的参数估计值、软件可靠性参数、模型的可靠性评估曲线等）输出。

系统的数据流图如图2-1所示：

用 户 用户命令 用户命令 可靠性评估 可靠性模型 评估结果 结果显示 可靠性模型库 失效数据 模型选择 可靠性模型 图2-1 系统的数据流图

2.2 数值算法设计

2.2.1 Goel-Okumoto模型参数计算

Goel-Okumoto模型参数的最大似然估计：

在应用G-O模型进行实际观测、分析软件故障数据时，首先要推定模型中的均值函数m(t)中的未定参数a和b。下面我们用最大似然估计法估计a和b的值。

按假设（4），若t时刻累积故障数为y，则得到N(t)的概率密度为：

140

尹晶杰：软件可靠性模型的算法分析与评价

Pr{N(t)?y}?m(t)yy!e?m(t)?(a(1?e?bt))yy!exp{?a(1?e?bt)} (2-1)

从式(2-1)可得出参数的似然函数为：

L(a,b)??i?1m(a(e?bti?1?e?bti))ni?ni?1(ni?ni?1)!?exp{?a(1?e?btm)} (2-2)

通过极大似然参数估计解的特性的分析实现了该计算方法的计算机编程，具体步骤如下： 步骤1：令给定的误差允许值为D??ti?1NiNtN。如果0?D?12，则xl?(1?2D)/2，故

xr?1D，转步骤2；如果D?12，则参数估计无解，转步骤5。

步骤2：计算xm?(xr?xl)/2，如果xr?xl??y，则转步骤4。 步骤3：f?(1?Dxm)e则xr?xm，转步骤2。

步骤4：计算b?xm?(D?1)xm?1，如果f??y，则xl?xm转步骤2；如果f???y，

?btNxmtN和a?N/(1?e)。

步骤5：停止计算。

2.2.2 Jelinski-Moranda模型参数计算

Jelinski-Moranda模型参数的最大似然估计：

在模型公式中有两个未知参数N0和?，如果在测试过程中测得的失效间隔时间为x1,x2,?,xn，则参数N0和?的似然函数为:

nL(N0,?)???(N0?i?1)exp{??(N0?i?1)xi} (2-5)

i?1通过极大似然参数估计解的特性的分析实现了该计算方法的计算机编程，具体步骤如下： 设ex为针对自变量x而给定的误差精度控制值，ey为针对自变量y而给定的误差精度控制值。计算

nnf(N)??N?i?1i?11?nN?p , p??(i?1)(ti?1i?ti?1)

tn步骤1：如果P?n?12，则left?n?1,right?n，转步骤2；如果P?n?12，则终止计算。

步骤2：如果f(right)?ey，则left?right和right?right?1，重复步骤2；若

?ex?f(right)?ey，则root?right，转步骤5；若f(right)??ex，转步骤3。

right?left步骤3：如果right?left?ex，则root?，转步骤5；如果right?left?ex，则

2right?left，转步骤4。 root?2步骤4：如果f(root)?ey，则left?root，转步骤3；如果?ey?f(root)?ey，则转步骤5；如果f(root)??ey，则right?root，转步骤3。

141