@@ -363,7 +363,7 @@ mc.simulate_indices(ts_length=4)
363363
364364考虑任意状态 $y \in S$,我们想计算系统在下一时刻处于状态 $y$ 的概率。
365365
366- 根据[ 全概率公式] ( https://en.wikipedia.org/wiki/Law_of_total_probability ) ,我们可以将这个概率表示为:
366+ 根据[ 全概率公式] ( https://baike.baidu.com/item/%E5%85%A8%E6%A6%82%E7%8E%87%E5%85%AC%E5%BC%8F/9980676 ) ,我们可以将这个概率表示为:
367367
368368$$
369369\mathbb P \{X_{t+1} = y \}
@@ -1031,13 +1031,16 @@ ax.grid()
10311031ax.hlines(0, 0, N, lw=2, alpha=0.6) # 在零处画水平线
10321032
10331033for x0, col in ((0, 'blue'), (1, 'green')):
1034+
10341035 # 生成从x0开始的工人的时间序列
10351036 X = mc.simulate(N, init=x0)
1037+
10361038 # 计算每个n的失业时间比例
10371039 X_bar = (X == 0).cumsum() / (1 + np.arange(N, dtype=float))
1040+
10381041 # 绘图
10391042 ax.fill_between(range(N), np.zeros(N), X_bar - p, color=col, alpha=0.1)
1040- ax.plot(X_bar - p, color=col, label=f '$X_0 = \, {x0} $')
1043+ ax.plot(X_bar - p, color=col, label=fr '$X_0 = \, {x0} $')
10411044 ax.plot(X_bar - p, 'k-', alpha=0.6)
10421045
10431046ax.legend(loc='upper right')
@@ -1232,19 +1235,22 @@ Q = np.zeros((n, n), dtype=int)
12321235with open(infile) as f:
12331236 edges = f.readlines()
12341237for edge in edges:
1235- from_node, to_node = re.findall('\w', edge)
1238+ from_node, to_node = re.findall(r '\w', edge)
12361239 i, j = alphabet.index(from_node), alphabet.index(to_node)
12371240 Q[i, j] = 1
1241+
12381242# 创建相应的马尔可夫矩阵P
12391243P = np.empty((n, n))
12401244for i in range(n):
12411245 P[i, :] = Q[i, :] / Q[i, :].sum()
12421246mc = MarkovChain(P)
1243- # 计算稳态分布r
1247+
1248+ # 计算平稳分布r
12441249r = mc.stationary_distributions[0]
12451250ranked_pages = {alphabet[i] : r[i] for i in range(n)}
1246- # 打印解决方案,从最高排名到最低排名排序
1247- print('排名\n ***')
1251+
1252+ # 输出解,从最高排名到最低排名排序
1253+ print(r'排名\n ***')
12481254for name, rank in sorted(ranked_pages.items(), key=itemgetter(1), reverse=1):
12491255 print(f'{name}: {rank:.4}')
12501256```
@@ -1322,7 +1328,7 @@ $P$,如上所述。
13221328:class: dropdown
13231329
13241330可以在[QuantEcon.py](http://quantecon.org/quantecon-py)库中找到解决方案,
1325- 具体见 [这里](https://github.com/QuantEcon/QuantEcon.py/blob/master/quantecon/markov/approximation.py)。
1331+ 具体可见 [这里](https://github.com/QuantEcon/QuantEcon.py/blob/master/quantecon/markov/approximation.py)。
13261332
13271333```
13281334
0 commit comments