11月14日作业

1、

实验代码：

import numpy as np
matrix = np.matrix([[0.9,0.075,0.025],[0.15,0.8,0.05],[0.25,0.25,0.5]], dtype=float)
vector1 = np.matrix([[0.3,0.4,0.3]], dtype=float)
for i in range(100):
    vector1 = vector1*matrix
    print （"Current round:" , i+1）
    print （vector1）

实验结果：

 

 2、

实验代码：

matrix = np.matrix([[0.9,0.075,0.025],[0.15,0.8,0.05],[0.25,0.25,0.5]], dtype=float)
vector1 = np.matrix([[0.7,0.1,0.2]], dtype=float)
for i in range(100):
    vector1 = vector1*matrix
    print （"Current round:" , i+1）
    print （vector1）

实验结果：

 3、实验总结：

通过课堂的学习和编写运行实验代码，对马尔科夫链有了初步的理解认识，学习到马尔科夫链在矩阵上的应用的优势。