2017-09-10 33 views
1

マルコフ連鎖問題を解くためにMathCadのKolmogorov-Chapman方程式を生成する方法を探しています。 問題は、システムがいずれかの状態にある確率を見つけることです。システムにはN個のコンポーネントがあります。 私は2^N個のノード(状態)と2 * Nパラメータ:N aのグラフを持っています。これはN番目のコンポーネントが故障する可能性があり、N b ''が壊れたコンポーネントが正常になる確率です再び。マルコフ過程のKolmogorov-Chapman方程式を生成する

Nは10に近いことがあります。つまり、少なくとも1024式があることを意味します。これらの式を生成する方法を探しています。

(その後scipyようななめらかで解決できる場合もパイソン)Mathcadのは、要件、任意の他の数学エンジンがOKでない

+0

まず、1つの方程式を書いてそれを一般化します。 「生成する方法を探している」とはどういう意味ですか?方程式を書いたりプログラミングしたりするのに問題がありますか?前者は良いことではない。後者は、あなたが書いたコードを投稿した場合かもしれません。 – duffymo

+0

Makov連鎖の確率行列形式を使用していますか?もしそうなら、行列の乗算が必要です。 [行列形式](https://www.youtube.com/watch?v=W5P4kCpdhho) –

答えて

1

問題のいずれかの状態にあるシステムの確率を見つけることです。

あなたの確率過程はN 独立しマルコフ連鎖の組み合わせであるため、問題は2 2のシステムを解くことよりも、別途各鎖のための方程式を解くことによりアプローチする方がはるかに簡単です* Nのため方程式全プロセス。

個々の鎖のための遷移行列は次のとおりです。

enter image description here

ので、コルモゴロフ・チャップマン方程式は以下のとおりです。

enter image description here

enter image description here

そして、あなたがまだ必要な場合は、システム全体の方程式:

enter image description here

+0

力の反復が良いヒューリスティックである可能性がありますか? – wildplasser

+1

@wildplasser 2x2の行列をべき乗(またはそのような行列や固有ベクトルの限界)を計算するためには、パワーリテイリングは必要ないと思います。それはすべて分析的に行うことができます。 – Anton