• 検索結果がありません。

隠れマルコフモデルによる不完全デバッグ環境でのソフトウェア信頼性評価法に関する考察

N/A
N/A
Protected

Academic year: 2021

シェア "隠れマルコフモデルによる不完全デバッグ環境でのソフトウェア信頼性評価法に関する考察"

Copied!
2
0
0

読み込み中.... (全文を見る)

全文

(1)

1−D−2

2001年度日本オペレーションズ・リサーチ学会

秋季研究発表会

隠れマルコフモデルによる不完全デバッグ環境での

ソフトウェア信頼性評価法に関する考察

01108985 鳥取大学 *木村光宏 KIMURAMitsuhiro

O1702425 鳥取大学 山田茂 mMADAShigeru

1 はじめに

ソフトウェアの定量的な信頼性評価手法は開発されたソフ

トウェア成果物の信頼性や可用性を計測し予測するために必要 であり,ソフトウェア工学における重要な問題の1つとなって いる.そのようなソフトウェア信頼性評価問題に取り組むため に,確率・統計論に基づいて構築された数多くのソフトウェア

信頼性評価モデルがこの30年間に多くの研究者らによって開

発されており.これらのモデルの中には,実用に供されている

モデルもある叶 しかしながら一般に,確率モデルはそのモデ

ルの成立のためにいくつかの仮定を必要とするが,その仮定の いくつかは実際の現象を記述するには厳しすぎるという批判が ある,

ソフトウェア信頼性評価モデルの研究分野において,完全

デバッグの仮定はモデルを簡単化するのに非常に役立つことが よく知られているが,多くのソフトウェア品質/信頼性研究者 らと専門家らは,ソフトウェアテスト環境に非現実的な完全デ バッグ環境を仮定することに批判的である.そのような背景の 下,この仮定を緩めたいくつかの不完全デバッグソフトウェア 信頼性評価モデルが開発された.しかしながら,これらのモデ

ルを実際のテスト工程などで適用する際に必要となる,不完全

デバッグ率を与える未知パラメータを推定する実際的な方法は 提案されていないのが現状である.

本研究では,不完全デバッグ環境を考慮した簡単なソフト

ウェア信頼性評価モデルをとりあげ,このモデルに含まれる未 知パラメーータの実用的な推定方法に焦点を当てる.具体的に

は,OkumotoaIldGoel【2】により提案された不完全デバッグモ

デルとテスト工程において実測された発見フォールト数データ から,彼らの研究では実用的に推定することができなかった不 完全デバッグ率の未知パラメータを推定する方法について考察 する.

2 モデルの記述

2.1 モデルの仮定

本研究が対象とするソフトウェアのテスト及びデバッグ環境に 対して次のように仮定する. ● ソフトウェアは要求仕様に基づきあらかじめ規定されたテ ストケースを使ってテストされる. ● ひとつのソフトウェア故障はひとつのソフトウェアフォー ルトによって引き起こされる. ● デバッグ過程において,ソフトウニ1ニアフォーールトは正しく 修正されるか,あるいは正しく修正されるとともに新しい フォールトがプログラムの別の部分(当該テストケース〝) 実行によってテストされるプログラムパス以外の部分)に 作り込まれる. 一般に,ソフトウェアのテスト工程では,テストケースの実 行によりソフトウェア内にフォールトが潜在していることが判 明すると,そのソフトウェアは詳しく調べられ,最終的に当該 テストケ一子が首尾まく処理されるまでデバッグされる.その 結果として,あるテストケースがソフトウェアシステムに対し て正しく処理されたとき,そのテストケ山スが処理される際に 通過するプログラムパスは少なくともそのテストケースに関し てはフォールトを含んではいない.言い換えれば,ソフトウェ アデバッグ作業者があるテストケースで発見されたフォールト

の修正に失敗するということは,デバッグ作業が当該テスト

ケースが正しく処理されるまで行われることを考えると,新し いフォールトがプログラムのある部分に作りこまれることにな り,この作り込まれたフォールトは,以後に適用されるテスト ケースによって発見されるか,あるいはテスト工程が終了する までもはや発見されないことになる.本研究ではこのようなデ バッグ環境を仮定する.

2.2 従来のモデル

前節の仮定に基づき,OkumotoandGoel【2]は不完全デバッグ

環境を考慮したソフトウェア信頼性評価モデルを提案した.こ のモデルはマルコフ過程としてソフトウェア内の潜在フォール ト数の振舞いを表したものであり,その状態遷移図を図1に表

す.この図におけるパラメータpは完全デバッグ率を,qは不

完全デバッグ率をそれぞれ表す(p+q=1).また肌まテスト 開始時の潜在フォールト数を表す.

8−グー&8

区= 状態遷移図

二のモデルの状態Jにおける滞在時間分布関数,すなわち,

ソフトウェア内に潜在する真のフォールト数がノ個であるとき,

次のソフトウェア故障が発生するまで〝)時間分布を, ち(り =1−eXI)ト入ノり(ノ=0,1,・・,〃), (1) により表す・ここで,入ノはソフトウエア故障発生に対するハ ザードレート(11孔Zar〔lrat(、)を表す・特に入.ノ=入・Jとすれば, こ〝)分布関数は最初にJelillSkiall(川Iorall(laトi】によって提案さ れたもび)を表す. OkllmOtOall〔1Goel[2]は最尤法とベイズ型推定によって統計 的に推定する方法を議論したが,この方法にはデータとしてそ れぞれのデバッグが完全であったか否かに関するデータが得ら −62− © 日本オペレーションズ・リサーチ学会. 無断複写・複製・転載を禁ず.

(2)

れることを前提としている.しかしながらそのようなデータを

得ることは実際には非常に難しい.したがって,テスト工程に

おけるソフトウェアフォールト発見時間間隔の様な,観測する ことが可能なデー タからモデルに含まれる未知パラメータを推 定することが必要になる.

2.3 本研究のモデル

前節に挙げた従来のモデルの未知パラメータの推定に関する問

題点を解決するため,本研究では隠れマルコフモデル(hidden_

Markovmodel)を用いた不完全デバッグモデルを考える.隠

れマルコフモデルはマルコフモデルの一種であるが,その状態

空間は隠れた状態,すなわち外部から観測できない状態を含ん

でいるところに特徴がある.このモデル化技法の有利な点は隠

れた状態に含まれる未知パラメータを推定することができるこ とにある.このような隠れマルコフモデルは主に音声認識の研

究分野で用いられてきた【4,5】

以下に,隠れマルコフモデルとしてモデル化する際に必要と なる緒量を定義する.

〃:ソフトウェア開発のテスト段階の開始時におけるソフト

ウェアシステム内の潜在フォールト数,すなわち初期状態

を表す. p‥完全デ/1ッグ率(0<p≦1) ヴ‥不完全デバッグ率(ヴ=1−p) 入j‥状態jにおけるソフトウェア故障発生に対するハザード レート(11aZardrate) jニ ソフトウェア内の潜在フォールト数を表す状態変数 巧(り:状態Jにおける滞在時間分布 汀j=状態Jの初期確率分布(∑J打ノ=1) αfJ‥状態才から状態ノへの推移確率

3 未知パラメータの推定法

本節では,前節で定義された隠れマルコフモデルに含まれる

未知パラメータ入J,ヴおよび〃を推定する方法について述べ

る.まず,分析するために得たフォールト発見時間間隔データ

を(豆,エ?)の形式に加工する.ここで,=ま五番目のフォールト

発見を表し,∬fは(五−1)番目からメ番目のフォールト発見時間

間隔を表す.

ニのデータに対して隠れマルコフモデルにおけるBaum_

Ⅵ硯ch再推定手続き[4,5]を用し、ることにより,不完全デバッグ

率り(=トp)をはじめとする未知パラメーータ入ノ,およびⅣを 推定することができる.

ここで,rは観測されたデータの個数を表す.また,

1(豆=0,j=0) ヴ(盲=j)

p(哀−j=1)’

0(otherwise)

(5) αり = dろ(た) dた J叫・) =入JeXp卜入Jた]

(J=Ⅳ+,Ⅳ+−1,‥.,1,0),(6)

打J恒(ごl)(f=1;メ=〃十,Ⅳ+−1,…;0) Ⅳ+

(∑αi(f−1)αiJ)軸木)(f=2,3,…,r;,(7)

i=O

j=Ⅳ+,〃+−1,‥.,1,0)

αJ(り =

頼=r;五=Ⅳ+,Ⅳ+−1,.‥,1,0)

〃+ ∑αiJ頃勘+1)卿+1) j=0

(f=r−1,r−2,...,1;

虞=Ⅳ十,Ⅳ+−1,…,1,0)

βi(t)=

,(8)

Ⅳ ̄十

榊)=αi(梱(f)/∑αi(r),

I=0 〃+

E誹)= αi(小町軸木+1)卿+1)/∑叫(r),

i=0

となる.ここに,Ⅳ+は繰返し推定のために必要となる潜在

フォールト数の初期値を表す.また,計算に先立って未知パラ メータヴ,入∫,および打=(汀0,れ,‥・,打Ⅳ+)に適当な初期値

を与える.これらの再評価手順を繰返した後,収束した値とし

て最終的な推定値を得ることができる.また,式(9)に与えら れた7i(りを用いて,最尤な状態推移の順序を推定することが

でき,吏番目のデバッグ作業が行われた状態の推定値giは

gi= argmaX[「′J(用(1≦豆≦r), (0≦J≦Ⅳ+) (11)

により得られる.さらに,ルにより表されるパラメータⅣの

推定値は式(11)からル=glによって与えられる.

謝辞 本研究の一部は文部科学省科学研究費補助金奨励研究(A)

(課題番号13780364)の下で行われた.

参考文献

【1]J・D・Musa,A・Iannino,andK.Okumot.0:SoPwart]Reliabili一

旬・■〟eαβ†▲僧n即もど,Predよc如†l,Ap〆血如†l,九l(:GraⅥ′−Hill,1eⅥ′ 1brk,1987. 【2]K・OkulnOtOandA・L・Goel‥“Classicalan〔1BaycsiaJliIlkrencc forthesoftwal,eilrlperfectdfll)ugglnglnOdel,”Tec山一icalllel)Ort To・78T2,Depart・mentOfIE&OR,Syracusel:niversit.y,1978. 【3]山田茂‥ ソフトウニ】ニア信頼性モデ′レm・基礎と〟甜ト,【】科技連, 東京,1994. 刷今井聖‥音声認識,共立出版,東京,1995.

f5]L・RL RahiherandB.H.Juang‥“Arlintroduct.iontohidden

MarkovModels,”IEEEASSPMagazine,1bl・3,To・1,pp・

4−16,1986.

3.1Iiaum−Ⅵね1ch再推定公式

以下にBauIn−11klぐ11再推定公式を示す. r】 r1

¢ = ∑如(り/∑「′0(り,

/=1 ′=1 (2) T 71

も = ∑棚]卿]/∑り伸帥れ

J=1 ′=1 (3) 勺 =「)(1) (4) −63− © 日本オペレーションズ・リサーチ学会. 無断複写・複製・転載を禁ず.

参照

関連したドキュメント

暑熱環境を的確に評価することは、発熱のある屋内の作業環境はいう

(ページ 3)3 ページ目をご覧ください。これまでの委員会における河川環境への影響予測、評

[r]

(評議員) 東邦協会 東京大学 石川県 評論家 国粋主義の立場を主張する『日

優越的地位の濫用は︑契約の不完備性に関する問題であり︑契約の不完備性が情報の不完全性によると考えれば︑

本稿で取り上げる関西社会経済研究所の自治 体評価では、 以上のような観点を踏まえて評価 を試みている。 関西社会経済研究所は、 年

実効性 評価 方法. ○全社員を対象としたアンケート において,下記設問に関する回答

今までの少年院に関する筆者の記述はその信瀝性が一気に低下するかもしれ