对香港天气预报的马尔科夫链性质分析及python模拟

查阅香港天文台的数据后发现，香港的雨量预报包括9天的预报情况，下雨的机率大致分为3个等级：高=0.7/中=0.5/低=0.3。

因此，我把 3 状态建立为如下图所示的马尔可夫链进行模拟。

然后我在 python 中对这个马尔可夫链进行了 100 次模拟。我发现不管start_matrix是什么，也就是初始状态的矩阵，得到的矩阵都会收敛到一个固定的矩阵。这完全符合马尔可夫链的无记忆特性。所以我认为天气预报确实是由一个马尔可夫链系统组成的。

import matplotlib.pyplot as plt
import numpy as np
transfer_matrix = np.array([[0.5,0.25,0.25],[