2017-09-10 63 views
1

我正在尋找一種方式來產生柯爾莫哥洛夫 - 查普曼方程Mathcad來解決馬爾可夫鏈的問題。 問題是要找到系統的狀態之一的概率。系統有N個組件。 我有2^N個節點(狀態),和2 * N參數的曲線圖:N a的,這是概率第N個成分會分解和N b的,這是破組件將成爲健康的概率再次。生成柯爾莫哥洛夫 - 查普曼方程馬爾可夫過程

N可以是接近10,這意味着,將有1024個方程式至少,所以我在尋找一種方式來產生這些方程。

的Mathcad是不是必需的,任何其他的數學引擎是OK(甚至Python的,如果再能與水木清華可以解決像scipy

+0

首先寫一個方程,然後推廣。 「尋找生成方式」是什麼意思?您是否在編寫方程式或編程時遇到困難?前者不是一個好的SO問題;後者可能是如果你發佈你寫的代碼。 – duffymo

+0

您使用馬可夫鏈的概率矩陣形式嗎?如果是這樣,你只需要矩陣乘法。 [矩陣形式(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我不認爲你需要動力迭代計算的冪(或矩陣的限制,或它們的特徵向量)2×2矩陣。所有這些都可以通過分析來完成。 – Anton