CN105260304B - 一种基于qbgsa‑rvr的软件可靠性预测方法 - Google Patents
一种基于qbgsa‑rvr的软件可靠性预测方法 Download PDFInfo
- Publication number
- CN105260304B CN105260304B CN201510675965.5A CN201510675965A CN105260304B CN 105260304 B CN105260304 B CN 105260304B CN 201510675965 A CN201510675965 A CN 201510675965A CN 105260304 B CN105260304 B CN 105260304B
- Authority
- CN
- China
- Prior art keywords
- software
- kernel function
- qbgsa
- rvr
- particle
- 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.)
- Active
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明适用于软件可靠性预测的技术领域,公开了一种基于QBGSA‑RVR的软件可靠性预测方法,其基于量子衍生二进制引力搜索算法优化相关向量机,即QBGSA‑RVR,包括如下步骤:观测并记录软件累计失效时间序列;软件失效序列数据归一化、选择核函数用于相关向量机、选择m值;参数优化;更新每个粒子位置,直至满足终止条件;得到最优核函数参数值,使用相关向量机进行预测;最后数据回放,得到下一时间段软件失效时间预测值。本发明提供的一种基于QBGSA‑RVR的软件可靠性预测方法,步骤简单,QBGSA‑RVR较一些广泛使用的神经网络及核函数方法预测效果更优,在小样本情况下使用QBGSA‑RVR方法进行软件可靠性预测也可以保持很好的预测效果,大大提高了软件的可靠性预测效果。
Description
【技术领域】
本发明涉及软件可靠性预测的技术领域,特别涉及一种基于QBGSA-RVR的软件可靠性预测方法。
【背景技术】
随着互联网技术的飞速发展,应用软件***规模越做越大越复杂,其可靠性越来越难保证。应用本身对***运行的可靠性要求越来越高,在一些关键的应用领域,如航空、航天等,其可靠性要求尤为重要,在银行等服务性行业,其软件***的可靠性也直接关系到自身的声誉和生存发展竞争能力。特别是软件可靠性比硬件可靠性更难保证,会严重影响整个***的可靠性。在许多项目开发过程中,对可靠性没有提出明确的要求,开发商也不在可靠性方面花更多的精力,往往只注重速度、结果的正确性和用户界面的友好性等,而忽略了可靠性。在投入使用后才发现大量可靠性问题,增加了维护困难和工作量,严重时只有束之高阁,无法投入实际使用。
软件可靠性预测指的是使用测试或运行过程中收集的失效数据对软件未来失效情况做出预测,它在软件质量保障以及软件成本控制中起着决定性作用,然而由于软件失效过程的复杂性及非线性,软件失效的精确预测异常困难。近年来,支持向量回归、相关向量回归等核函数技术在软件可靠性预测中得到了较好的应用,然而核函数软件可靠性预测中的参数优化还有许多需要解决的问题。
量子衍生二进制引力搜索算法(Quantum-inspired Binary GravitationalSearch Algorithm),简称QBGSA,是一种基于量子计算原理的新型概率优化方法,它以量子计算的一些概念和理论为基础,用量子位编码来表示粒子,用量子门更新来完成进化搜索,利用旋转变换来更新粒子的位置,基于量子衍生二进制引力搜索算法优化的相关向量回归,简称QBGSA-RVR,使用量子衍生二进制引力搜索算法优化相关向量回归估计中核函数参数,具有较好的预测效果,因此提出一种基于QBGSA-RVR的软件可靠性预测方法。
【发明内容】
本发明的目的在于克服上述现有技术的不足,提供一种基于QBGSA-RVR的软件可靠性预测方法,其旨在解决现有技术中核函数软件可靠性预测中的参数优化效果较差、软件可靠性预测效果较低的技术问题。
为实现上述目的,本发明提出了一种基于QBGSA-RVR的软件可靠性预测方法,其基于量子衍生二进制引力搜索算法优化相关向量机,即QBGSA-RVR,包括如下步骤:
步骤一、观测并记录软件失效时间序列数据集t1,t2,…,tn,失效时间t与在它之前发生的m次失效时间之间存在映射关系,选择合适的核函数对软件失效数据进行归一化处理,归一化映射公式为将软件失效序列数据转换到(0,1)区间,得到n个软件失效序列数据,方便相关向量机进行学习、预测;
步骤二、若核函数宽度值r(i)∈(Maxr,Minr),随机取r(0)∈(Maxr,Minr),使用式子将r(0)映射得到初始化混沌变量x(0)∈(0,1),应用Logistic映射生成N个混沌变量x(i+1),通过变换r(i+1)=Minr+x(i+1)·(Maxr-Minr),i=0,1,2,…,N-1,得到N个的核函数宽度值,使用M维量子位对它们进行二进制编码,得到规模为N的初始化种群;
步骤三、使用已观测到的软件失效序列数据t1,t2,…,tn进行相关向量机向量回归学习,计算不同核函数宽度下的适应度,适应度函数为其中表示归一化后的软件失效时间预测值,t′i为归一化后的软件失效时间实测值;
步骤四、计算计算best(t)和worst(t),使用式子以及适应度函数计算每个粒子i的惯性质量mi,其中Tmax为最大循环次数,Mi(t)为粒子i的引力质量;
步骤五、根据公式计算每个粒子i在循环t时维数d的旋转因子Δθid(t),式中其值在(θmin,θmax)之间单调增长,用以控制旋转角度的大小,pkd表示粒子k在第d维上的位置,pid表示粒子i在第d维上的位置;
步骤六、进行旋转变换通过式子更新每个粒子位置上的值,其中rand(0,1)表示(0,1)之间的随机数;
步骤七、判断是否满足终止条件,若不满足,则回转至步骤三,若满足,则得到最优核函数参数值;
步骤八、使用最优核函数参数情况下的相关向量机对下一时间段的归一化软件失效时间进行预测,预测完成后使用映射将数据回放即可得到真实预测值。
作为优选,所述步骤一中合适的核函数包括高斯核函数、线性核函数、多项式核函数、柯西核函数、拉普拉斯核函数、对称三角核函数、双曲正割核函数、平方正弦基核函数。
作为优选,所述高斯核函数包括具有非线性特性的高斯核函数r表示核函数宽度。
作为优选,所述步骤一中m的取值为5~15之间。
作为优选,所述步骤四中best(t)的计算公式为worst(t)的计算公式为
作为优选,所述步骤五的旋转因子Δθid(t)计算公式中,其中fi(t)为适应度函数
作为优选,所述在的计算公式中,Mk(t)和Mi(t)分别为粒子k和粒子i的引力质量,Rik是粒子i和粒子k之间的Hamming距离。
作为优选,所述的公式中τ为常量。
本发明的有益效果:与现有技术相比,本发明提供的一种基于QBGSA-RVR的软件可靠性预测方法,步骤合理,使用量子衍生二进制引力搜索算法优化相关向量回归估计核函数参数,QBGSA-RVR较一些广泛使用的神经网络及核函数方法预测效果更优;其次,使用不同样本数量情况下,甚至在小样本情况下使用QBGSA-RVR方法进行软件可靠性预测也可以保持很好的预测效果,改善了核函数软件可靠性预测中参数优化效果较差的问题,且预测偏差在5%区间内,提高了软件可靠性的预测效果。
本发明的特征及优点将通过实施例结合附图进行详细说明。
【附图说明】
图1是本发明一种基于QBGSA-RVR的软件可靠性预测方法的步骤示意图。
【具体实施方式】
为使本发明的目的、技术方案和优点更加清楚明了,下面通过附图中及实施例,对本发明进行进一步详细说明。但是应该理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限制本发明的范围。此外,在以下说明中,省略了对公知结构和技术的描述,以避免不必要地混淆本发明的概念。
参阅图1,本发明实施例提供一种基于QBGSA-RVR的软件可靠性预测方法,其基于量子衍生二进制引力搜索算法优化相关向量机,即QBGSA-RVR,包括如下步骤:
步骤一、观测并记录软件失效时间序列数据集t1,t2,…,tn,失效时间t与在它之前发生的m次失效时间之间存在映射关系,选择合适的核函数对软件失效数据进行归一化处理,归一化映射公式为将软件失效序列数据转换到(0,1)区间,得到n个软件失效序列数据,方便相关向量机进行学习、预测;
步骤二、若核函数宽度值r(i)∈(Maxr,Minr),随机取r(0)∈(Maxr,Minr),使用式子将r(0)映射得到初始化混沌变量x(0)∈(0,1),应用Logistic映射生成N个混沌变量x(i+1),通过变换r(i+1)=Minr+x(i+1)·(Maxr-Minr),i=0,1,2,…,N-1,得到N个的核函数宽度值,使用M维量子位对它们进行二进制编码,得到规模为N的初始化种群;
步骤三、使用已观测到的软件失效序列数据t1,t2,…,tn进行相关向量机向量回归学习,计算不同核函数宽度下的适应度,适应度函数为其中表示归一化后的软件失效时间预测值,t′i为归一化后的软件失效时间实测值;
步骤四、计算计算best(t)和worst(t),best(t)的计算公式为worst(t)的计算公式为使用式子以及适应度函数计算每个粒子i的惯性质量mi,其中Tmax为最大循环次数,Mi(t)为粒子i的引力质量;
步骤五、根据公式计算每个粒子i在循环t时维数d的旋转因子Δθid(t),式中其值在(θmin,θmax)之间单调增长,用以控制旋转角度的大小,pkd表示粒子k在第d维上的位置,pid表示粒子i在第d维上的位置,Mk(t)和Mi(t)分别为粒子k和粒子i的引力质量,Rik是粒子i和粒子k之间的Hamming距离,fi(t)为适应度函数
步骤六、进行旋转变换通过式子更新每个粒子位置上的值,其中rand(0,1)表示(0,1)之间的随机数;
步骤七、判断是否满足终止条件,若不满足,则回转至步骤三,若满足,则得到最优核函数参数值;
步骤八、使用最优核函数参数情况下的相关向量机对下一时间段的归一化软件失效时间进行预测,预测完成后使用映射将数据回放即可得到真实预测值。
其中,步骤一中合适的核函数包括高斯核函数、线性核函数、多项式核函数、柯西核函数、拉普拉斯核函数、对称三角核函数、双曲正割核函数、平方正弦基核函数。
进一步地,为达到更好的预测效果,在本发明实施例中核函数采用具有非线性特性的高斯核函数r表示核函数宽度。
更进一步地,步骤一中m的取值为5~15之间,能取得相对较好的预测性能。在本发明实施例中m取8。
实验:
实验中用到的软件失效数据集如表1所示,表1中列出了所使用数据集的类型、记录的失效次数等。它们是从实际软件工程项目测试与运行过程中收集来的,经过分析和整理,精度高,适用于比较、验证和评价软件可靠性模型,包含累计执行时间数据和累计失效次数数据。本文实验过程中,使用软件累计执行时间数据进行学习预测,取所有数据集的前三分之二作为学习数据,对后面三分之一数据进行预测后与真实数据进行比较以评价模型的预测能力。
表1 实验分析所使用的数据集
数据集 | 代码行数 | 发生失效数 | 软件类型 | 用于学习的失效数据 | 学习组数 |
DATA1 | 21700 | 136 | 实时命令控制应用*** | 1~90 | 82 |
DATA2 | 10000 | 118 | 动态飞行控制*** | 1~81 | 73 |
DATA3 | 22500 | 180 | 动态飞行控制*** | 1~120 | 112 |
DATA4 | 38500 | 213 | 动态飞行控制*** | 1~144 | 136 |
使用广泛采用的平均相对预测误差对QBGSA-RVR的拟合及预测性能进行评价,计算式子分别如下所示:
式中,n表示失效数据集中发生失效的总数,分别表示上下取整运算,ti表示失效时间实际值。在AEFitting式中,表示失效时间拟合值,AEFitting的值越小,说明模型的拟合能力越好;在AEPredicting中,表示失效时间预测值,AEPredicting的值越小,说明模型的预测能力越强。
从而,用于QBGSA-RVR学习的顺序失效数据为{t1,t2,…,tp},其中如表1所示,四个数据集DATA1至DATA4上p的值分别为90、81、120和144,用于QBGSA-RVR进行学习的样本数为P-8,则DATA1至DATA4上的进行学习的样本数分别为82、73、112以及136,如表2所示。
表2 用于QBGSA-RVR的软件失效数据
QBGSA-RVR中初始值设置:种群规模N=30,最大循环次数Tmax=200,α=(0.5,0.5,…,0.5),核函数宽度初始范围r∈(0,10),θmin=0.001π,θmax=0.050π,τ=8%。
实验结果数据中,DATA1至DATA4上使用QBGSA-RVR的失效时间预测值与实际值的对比,绝大部分预测值落在实际值的5%区间内,这显示出QBGSA-RVR具有极优的软件可靠性预测能力,较现有的软件预测具有更好的预测效果,更接近实际值。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换或改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种基于QBGSA-RVR的软件可靠性预测方法,其基于量子衍生二进制引力搜索算法优化相关向量机,即QBGSA-RVR,其特征在于:包括如下步骤:
步骤一、观测并记录软件失效时间序列数据集t1,t2,···,tn,失效时间t与在它之前发生的m次失效时间之间存在映射关系,选择合适的核函数对软件失效数据进行归一化处理,归一化映射公式为将软件失效序列数据转换到(0,1)区间,得到n个软件失效序列数据,方便相关向量机进行学习、预测;
步骤二、若核函数宽度值r(i)∈(Maxr,Minr),随机取r(0)∈(Maxr,Minr),使用式子将r(0)映射得到初始化混沌变量x(0)∈(0,1),应用Logistic映射生成N个混沌变量x(i+1),通过变换r(i+1)=Minr+x(i+1)·(Maxr-Minr),i=0,1,2,···,N-1,得到N个的核函数宽度值,使用M维量子位对它们进行二进制编码,得到规模为N的初始化种群;
步骤三、使用已观测到的软件失效序列数据t1,t2,···,tn进行相关向量机向量回归学习,计算不同核函数宽度下的适应度,适应度函数为其中表示归一化后的软件失效时间预测值,ti'为归一化后的软件失效时间实测值;
步骤四、计算计算best(t)和worst(t),best(t)的计算公式为worst(t)的计算公式为使用式子以及适应度函数计算每个粒子i的惯性质量mi,其中Tmax为最大循环次数,Mi(t)为粒子i的引力质量;
步骤五、根据公式计算每个粒子i在循环t时维数d的旋转因子Δθid(t),式中其值在(θmin,θmax)之间单调增长,用以控制旋转角度的大小,其中τ为常量,fi(t)为适应度函数Mk(t)和Mi(t)分别为粒子k和粒子i的引力质量,Rik是粒子i和粒子k之间的Hamming距离,pkd表示粒子k在第d维上的位置,pid表示粒子i在第d维上的位置;
步骤六、进行旋转变换通过式子更新每个粒子位置上的值,其中rand(0,1)表示(0,1)之间的随机数;
步骤七、判断是否满足终止条件,若不满足,则回转至步骤三,若满足,则得到最优核函数参数值;
步骤八、使用最优核函数参数情况下的相关向量机对下一时间段的归一化软件失效时间进行预测,预测完成后使用映射将数据回放即可得到真实预测值。
2.如权利要求1所述的一种基于QBGSA-RVR的软件可靠性预测方法,其特征在于:所述步骤一中合适的核函数包括高斯核函数、线性核函数、多项式核函数、柯西核函数、拉普拉斯核函数、对称三角核函数、双曲正割核函数、平方正弦基核函数。
3.如权利要求2所述的一种基于QBGSA-RVR的软件可靠性预测方法,其特征在于:所述高斯核函数包括具有非线性特性的高斯核函数r表示核函数宽度。
4.如权利要求1所述的一种基于QBGSA-RVR的软件可靠性预测方法,其特征在于:所述步骤一中m的取值为5~15之间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510675965.5A CN105260304B (zh) | 2015-10-19 | 2015-10-19 | 一种基于qbgsa‑rvr的软件可靠性预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510675965.5A CN105260304B (zh) | 2015-10-19 | 2015-10-19 | 一种基于qbgsa‑rvr的软件可靠性预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105260304A CN105260304A (zh) | 2016-01-20 |
CN105260304B true CN105260304B (zh) | 2018-03-23 |
Family
ID=55100004
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510675965.5A Active CN105260304B (zh) | 2015-10-19 | 2015-10-19 | 一种基于qbgsa‑rvr的软件可靠性预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105260304B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106021097B (zh) * | 2016-05-09 | 2019-02-26 | 北京航空航天大学 | 基于测试特征的软件可靠性指标区间估计方法 |
CN106354640A (zh) * | 2016-08-29 | 2017-01-25 | 江苏科技大学 | 基于改进pso的软件可靠性模型参数估计方法 |
CN112000706B (zh) * | 2020-04-30 | 2021-05-18 | 中南民族大学 | 一种基于rvm-l模型的互联网舆情预警机制分析方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6330491B1 (en) * | 1999-07-21 | 2001-12-11 | Nicholas Lion | Integrated system and method of vending prescription medications using a network of remotely distributed, automated dispensing units |
US6879944B1 (en) * | 2000-03-07 | 2005-04-12 | Microsoft Corporation | Variational relevance vector machine |
CN101520652A (zh) * | 2009-03-03 | 2009-09-02 | 华中科技大学 | 一种数控装备服役可靠性的评估方法 |
CN101546290A (zh) * | 2009-04-30 | 2009-09-30 | 上海交通大学 | 提高面向对象型软件中类阶层质量预测准确度的方法 |
CN102004661A (zh) * | 2010-06-09 | 2011-04-06 | 电子科技大学 | 一种通用的数据驱动的软件和***可靠性模型及参数优化方法 |
CN102411339A (zh) * | 2011-11-30 | 2012-04-11 | 华中科技大学 | 一种数控装备性能可靠性评估方法 |
CN103093094A (zh) * | 2013-01-14 | 2013-05-08 | 湖州师范学院 | 基于核偏最小二乘回归算法的软件失效时间预测方法 |
CN103093095A (zh) * | 2013-01-14 | 2013-05-08 | 湖州师范学院 | 基于核主成分回归算法的软件失效时间预测方法 |
CN103106139A (zh) * | 2013-01-14 | 2013-05-15 | 湖州师范学院 | 基于相关向量回归估计的软件失效时间预测方法 |
-
2015
- 2015-10-19 CN CN201510675965.5A patent/CN105260304B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6330491B1 (en) * | 1999-07-21 | 2001-12-11 | Nicholas Lion | Integrated system and method of vending prescription medications using a network of remotely distributed, automated dispensing units |
US6879944B1 (en) * | 2000-03-07 | 2005-04-12 | Microsoft Corporation | Variational relevance vector machine |
CN101520652A (zh) * | 2009-03-03 | 2009-09-02 | 华中科技大学 | 一种数控装备服役可靠性的评估方法 |
CN101546290A (zh) * | 2009-04-30 | 2009-09-30 | 上海交通大学 | 提高面向对象型软件中类阶层质量预测准确度的方法 |
CN102004661A (zh) * | 2010-06-09 | 2011-04-06 | 电子科技大学 | 一种通用的数据驱动的软件和***可靠性模型及参数优化方法 |
CN102411339A (zh) * | 2011-11-30 | 2012-04-11 | 华中科技大学 | 一种数控装备性能可靠性评估方法 |
CN103093094A (zh) * | 2013-01-14 | 2013-05-08 | 湖州师范学院 | 基于核偏最小二乘回归算法的软件失效时间预测方法 |
CN103093095A (zh) * | 2013-01-14 | 2013-05-08 | 湖州师范学院 | 基于核主成分回归算法的软件失效时间预测方法 |
CN103106139A (zh) * | 2013-01-14 | 2013-05-15 | 湖州师范学院 | 基于相关向量回归估计的软件失效时间预测方法 |
Non-Patent Citations (3)
Title |
---|
基于相关向量回归的仿真元建模方法;吴冰等;《计算机工程》;20100215;第36卷(第3期);第24-27页 * |
软件可靠性预测中不同核函数的预测能力评估;楼俊钢等;《计算机学报》;20130630;第36卷(第6期);第1303-1311页 * |
软件可靠性预测的相关向量机模型;楼俊钢等;《计算机研究与发展》;20130731;第50卷(第7期);第1542-1550页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105260304A (zh) | 2016-01-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Song et al. | Short-term wind speed forecasting with Markov-switching model | |
Wang et al. | Improved LSTM-based time-series anomaly detection in rail transit operation environments | |
Le Son et al. | Remaining useful life estimation based on stochastic deterioration models: A comparative study | |
Lu et al. | Software defect prediction using semi-supervised learning with dimension reduction | |
WO2020237729A1 (zh) | 一种基于模式转移的虚拟机混合备用动态可靠性评估方法 | |
Sunahase et al. | Pairwise hits: Quality estimation from pairwise comparisons in creator-evaluator crowdsourcing process | |
CN111523778A (zh) | 基于粒子群算法和梯度提升树的电网运行安全评估方法 | |
Su et al. | Gaussian process machine-learning method for structural reliability analysis | |
JP6783002B2 (ja) | 企業の債務不履行予測システム及びその動作方法 | |
CN105260304B (zh) | 一种基于qbgsa‑rvr的软件可靠性预测方法 | |
CN105260786A (zh) | 一种电力推进***仿真可信度评估模型综合优化方法 | |
CN103617259A (zh) | 一种基于有社会关系和项目内容的贝叶斯概率矩阵分解推荐方法 | |
Guo et al. | Time series prediction method based on LS-SVR with modified gaussian RBF | |
CN109472048A (zh) | 基于稀疏多项式混沌扩展评估智能电表结构可靠度的方法 | |
CN108205713A (zh) | 一种区域风电功率预测误差分布确定方法和装置 | |
CN114266384A (zh) | 一种基于融合模型的快速预测台风路径的方法及其装置 | |
Park et al. | Source term estimation using deep reinforcement learning with Gaussian mixture model feature extraction for mobile sensors | |
Yang et al. | Fast flow field prediction of three-dimensional hypersonic vehicles using an improved Gaussian process regression algorithm | |
Li et al. | Learning to infer voltage stability margin using transfer learning | |
CN107104747A (zh) | 无线时变信道中的多径分量的分簇方法 | |
Wu et al. | An MCMC based EM algorithm for mixtures of Gaussian processes | |
Zhang et al. | Fuzzy stochastic Petri nets and analysis of the reliability of multi‐state systems | |
Arab et al. | Enhancing weighted uniform simulation for structural reliability analysis | |
Chen et al. | Distributed hierarchical deep reinforcement learning for large-scale grid emergency control | |
CN113642785A (zh) | 基于先验信息的空间碎片轨道长期预报方法、***及设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |