マルコフ連鎖の母期待値の時間発展の数値計算(markovexpect01) - 課題p052

情報

• 出題:2019-05-14
• 実行/提出期限:2019-06-04 23:55:00
• 提出
  • markovexpect01.c プログラム
  • markovexpect01.xlsx 結果のデータとグラフ

学習目標

• マルコフ連鎖で母期待値の時間変化を計算するプログラムを書ける.

課題

状況の説明

課題markov01と同じ状況, 同じ転置推移確率行列\(M\)を考える.

課題のタスク

• 課題markov01で作成したプログラム, または解答のプログラムを改造して, 時刻ごとに母期待値 \(\mathrm{E}[X(t)]\),母分散\(\mathrm{V}[X(t)]\), \(X(t)>1\)の母比率\(p\)を計算する, 仕様の通りのプログラムを作成しよう.
• T=30で実行し, 横軸t, 縦軸母期待値or母分散or母比率でグラフを描こう.

プログラムの入力

• 最終時刻T(0以上の整数)
• t=0での初期分布ベクトルの成分p[0] p[1] ... p[m-1](m個の小数)

プログラムの入力例

3
1.0 0.0 0.0 0.0  

プログラムの出力

• 1行目に#T=に続いて 最終時刻\(T\)
• 2行目にラベル t,p0,p1,p2,p3,EX,VX,P
• 2行目にコンマで区切って, 時刻t=0の\(t, p_0,p_1,p_2,\ldots,p_{m-1},\mathrm{E}[X(T)]\),母分散\(\mathrm{V}[X(T)]\),母比率\(p\)
• 3行目にコンマで区切って, 時刻t=1の\(t, p_0,p_1,p_2,\ldots,p_{m-1},\mathrm{E}[X(T)]\),母分散\(\mathrm{V}[X(T)]\),母比率\(p\)
• …
• T+2行目にコンマで区切って, 時刻t=Tの\(t, p_0,p_1,p_2,\ldots,p_{m-1},\mathrm{E}[X(T)]\),母分散\(\mathrm{V}[X(T)]\),母比率\(p\)

プログラムの出力例

#T=3
0,1.0,0.0,0.0,0.0,0.0,1,1,0
1,1.0,0.0,0.0,0.0,0.0,1,1,0
2,1.0,0.0,0.0,0.0,0.0,1,1,0
3,1.0,0.0,0.0,0.0,0.0,1,1,0

アドバイス

公式\(\mathrm{V}[X(t)]=\mathrm{E}[X(t)^2]-\mathrm{E}[X(t)]^2\)を使うのがいいでしょう.

標本期待値じゃなく母期待値なのでご注意. \(m=4\)って標本サイズじゃなくて母分布の量だよね. 標本抽出してるわけじゃないからもちろんそうなんですけど.

print_distを書き替えるのがいちばん簡単かも. 関数プロトタイプ宣言も同時に変更しないと.