状态转移矩阵 和 马尔科夫链的平稳分布

状态转移矩阵
https://blog.csdn.net/gdp12315_gu/article/details/49949669
马尔科夫链的平稳分布
http://www.docin.com/p-2033406117.html

第二个链接代码

import numpy as np
import matplotlib.pyplot as plt


s0=(0.52, 0.3, 0.18)
s0=np.array(s0)
p=np.array([[0.6, 0.3, 0.1],
           [0.35, 0.35, 0.3],
            [0.1, 0.2, 0.7]])
a=[]
b=[]
c=[]
for i in range(60):
    print(np.dot(s0, p))
    s0=np.dot(s0,p)
    a.append(s0[0])
    b.append(s0[1])
    c.append(s0[2])

plt.plot(a)
plt.plot(b)
plt.plot(c)
plt.show()