@@ -12,7 +12,9 @@ kernelspec:
12
12
---
13
13
14
14
# 就业的湖泊模型
15
+
15
16
## 大纲
17
+
16
18
除了 Anaconda 中包含的内容外,本讲座还需要以下库:
17
19
18
20
``` {code-cell} ipython3
@@ -26,6 +28,7 @@ plt.rcParams['font.family'] = ['Source Han Serif SC']
26
28
```
27
29
28
30
## 湖泊模型
31
+
29
32
这个模型有时被称为** 湖泊模型** ,因为存在两个工人群体:
30
33
1 . 当前就业的工人。
31
34
2 . 当前失业但正在寻找工作的工人。
@@ -45,8 +48,11 @@ plt.rcParams['font.family'] = ['Source Han Serif SC']
45
48
```
46
49
47
50
## 动态变化
51
+
48
52
设 $e_t$ 和 $u_t$ 分别表示在时间 $t$ 时就业和失业工人的数量。
53
+
49
54
工人总人口为 $n_t = e_t + u_t$。
55
+
50
56
因此,失业和就业工人的数量按以下方式变化:
51
57
52
58
``` {math}
@@ -59,7 +65,7 @@ plt.rcParams['font.family'] = ['Source Han Serif SC']
59
65
\end{aligned}
60
66
```
61
67
62
- 我们可以将{eq}` lake_model ` 安排为矩阵形式的线性方程组 $x_ {t+1} = Ax_t$,其中
68
+ 我们可以将{eq}` lake_model ` 整理为矩阵形式的线性方程组 $x_ {t+1} = Ax_t$,其中
63
69
64
70
$$
65
71
x_{t+1} =
87
93
88
94
因此,该系统的长期结果可能取决于初始条件 $x_0$ 和矩阵 $A$。
89
95
90
- 我们关心 $u_t$ 和 $e_t$ 随时间如何演变 。
96
+ 我们关心 $u_t$ 和 $e_t$ 如何随时间如何演变 。
91
97
92
- 我们应该预期什么样的长期失业率和就业率 ?
98
+ 失业率和就业率从长期来看会如何变化 ?
93
99
94
100
长期结果是否取决于初始值 $(u_0, e_o)$?
95
101
96
102
### 可视化长期结果
103
+
97
104
让我们首先绘制失业率 $u_t$、就业率 $e_t$ 和劳动力 $n_t$ 的时间序列图。
98
105
99
106
``` {code-cell} ipython3
@@ -105,7 +112,7 @@ class LakeModel:
105
112
------------
106
113
λ:标量
107
114
当前失业工人的找到工作的比率
108
- α:标量
115
+ α:标量
109
116
当前就业工人的解雇率
110
117
b:标量
111
118
进入劳动力市场的比率
@@ -143,7 +150,7 @@ class LakeModel:
143
150
x:迭代器
144
151
包含就业率和失业率的序列
145
152
"""
146
- x0 = np.atleast_1d(x0) # 以防万一, 重新转换为数组
153
+ x0 = np.atleast_1d(x0) # 重新转换为数组
147
154
x_ts= np.zeros((2, T))
148
155
x_ts[:, 0] = x0
149
156
for t in range(1, T):
@@ -180,16 +187,15 @@ for ax in axes:
180
187
plt.tight_layout()
181
188
plt.show()
182
189
```
183
- 将劳动力 $n_t$ 以恒定速率增长的观察结果并不令人惊讶 。
190
+ 不出所料,我们观察到劳动力 $n_t$以恒定的速率增长 。
184
191
185
- 这与失业和就业池只有一个流入源 (新进入者池)的事实相吻合 。
192
+ 这与只有一个流入源 (新进入者池)流向失业和就业池的事实相吻合 。
186
193
187
- 劳动力市场系统的流入和流出
188
- 在长期内由劳动力市场的固定退出率和进入率决定。
194
+ 从长期来看,劳动力市场系统的流入和流出是由劳动力市场的恒定退出率和进入率决定的。
189
195
190
196
具体来说,令 $\mathbb{1}=[ 1, 1] ^\top$ 为一个全1向量。
191
197
192
- 观察到
198
+ 我们可以得出
193
199
194
200
$$
195
201
\begin{aligned}
207
213
208
214
### 佩龙-弗罗贝尼乌斯定理的应用
209
215
210
- 直观上 ,如果我们将失业池和就业池视为一个封闭系统,其增长应与劳动力相似 。
216
+ 从直觉上讲 ,如果我们将失业池和就业池视为一个封闭系统,其增长应该与劳动力相似 。
211
217
212
- 接下来我们询问 $e_t$ 和 $u_t$ 的长期增长率
213
- 是否也由 $1+b-d$ 主导,就像劳动力一样。
218
+ 接下来我们要问, $e_t$ 和 $u_t$ 的长期增长率
219
+ 是否也像劳动力一样主要由 $1+b-d$ 决定?
214
220
215
- 如果我们求助于 {ref}` 佩龙-弗罗贝尼乌斯定理<perron-frobe> ` ,答案将更加清晰。
221
+ 如果我们应用 {ref}` 佩龙-弗罗贝尼乌斯定理<perron-frobe> ` ,答案将更加清晰。
216
222
217
223
佩龙-弗罗贝尼乌斯定理的重要性源于以下事实:
224
+
218
225
首先,在现实世界中,我们遇到的大多数矩阵都是非负矩阵。
219
226
220
- 其次,许多重要模型都是简单的线性迭代模型,
221
- 从初始条件 $x_0$ 开始,然后按规则 $x_ {t+1} = Ax_t$ 或简写为 $x_t = A^tx_0$ 递归演变。
227
+ 其次,许多重要模型只是线性迭代模型,它们从初始条件$x_0$开始,然后通过规则$x_ {t+1} = Ax_t$或简写为$x_t = A^tx_0$递归演化。
222
228
223
- 这个定理有助于表征主导特征值 $r(A)$,它
224
- 决定了这个迭代过程的行为。
229
+ 这个定理有助于描述主特征值$r(A)$,它决定了这个迭代过程的行为。
225
230
226
231
#### 主导特征向量
227
- 现在我们通过展示佩龙-弗罗贝尼乌斯定理如何帮助我们分析湖泊模型来说明它的力量。
232
+
233
+ 现在我们学习如何运用佩龙-弗罗贝尼乌斯定理来帮助我们分析湖泊模型。
228
234
229
235
由于 $A$ 是非负且不可约的矩阵,佩龙-弗罗贝尼乌斯定理意味着:
230
236
241
247
A \phi = r(A) \phi, \quad \psi A = r(A) \psi
242
248
$$
243
249
244
- - 如果进一步 $A$ 是正的,那么当 $<\psi, \phi> = \psi \phi=1$ 时,我们有
250
+ - 如果进一步假设 $A$ 是正的,那么当 $<\psi, \phi> = \psi \phi=1$ 时,我们有
245
251
246
252
$$
247
253
r(A)^{-t} A^t \to \phi \psi
248
254
$$
249
255
250
- 最后一个陈述意味着长期来看,$A^t$ 的量级与 $r(A)^t$ 的量级相同,其中 $r(A)$ 在本讲座中可被视为主导特征值 。
256
+ 最后一个陈述意味着长期来看,$A^t$ 的量级与 $r(A)^t$ 的量级相同,其中 $r(A)$ 在本讲座中可被视为主特征值 。
251
257
252
258
因此,长期来看,$x_t = A^t x_0$ 的量级也由 $r(A)^t$ 主导。
253
259
259
265
\min_j \text{colsum}_j (A) \leq r(A) \leq \max_j \text{colsum}_j (A)
260
266
```
261
267
262
- 注意,对于 $j=1,2$,$\text{colsum}_ j(A) = 1 + b - d$,根据 {eq}` PF_bounds ` ,我们因此可以得出主特征值为 $r(A) = 1 + b - d$。
268
+ 注意,对于 $j=1,2$,有 $\text{colsum}_ j(A) = 1 + b - d$,且根据 {eq}` PF_bounds ` ,我们可以得出主特征值为 $r(A) = 1 + b - d$。
263
269
264
270
令 $g = b - d$ 表示总劳动力的整体增长率,因此 $r(A) = 1 + g$。
265
271
@@ -304,11 +310,11 @@ def plot_time_paths(lm, x0=None, T=1000, ax=None):
304
310
305
311
if ax is None:
306
312
fig, ax = plt.subplots(figsize=(10, 8))
307
- # Plot line D
313
+ # 绘制线$D$
308
314
s = 10
309
- ax.plot([0, s * ū], [0, s * ē], "k--", lw=1, label='set $D$')
315
+ ax.plot([0, s * ū], [0, s * ē], "k--", lw=1, label='设置 $D$')
310
316
311
- # Set the axes through the origin
317
+ # 设置通过原点的坐标轴
312
318
for spine in ["left", "bottom"]:
313
319
ax.spines[spine].set_position("zero")
314
320
for spine in ["right", "top"]:
@@ -364,6 +370,7 @@ $D := \{ x \in \mathbb{R}^2 : x = \alpha \bar{x} \; \text{对某些} \; \alpha >
364
370
这表明所有这样的序列在长期内都具有很强的相似性,由主特征向量 $\bar{x}$ 决定。
365
371
366
372
#### 负增长率
373
+
367
374
在上面说明的例子中,我们考虑的参数使得劳动力的总体增长率 $g>0$。
368
375
369
376
现在假设我们面临 $g<0$ 的情况,即劳动力呈负增长。
@@ -389,7 +396,7 @@ plot_time_paths(lm, x0=x0)
389
396
390
397
由于 $A$ 的列和为 $r(A)=1$,左特征向量是 $\mathbb{1}^\top=[ 1, 1] $。
391
398
392
- 佩龙-弗罗贝尼乌斯理论意味着
399
+ 佩龙-弗罗贝尼乌斯理论告诉我们
393
400
394
401
$$
395
402
r(A)^{-t} A^{t} \approx \bar{x} \mathbb{1}^\top = \begin{bmatrix} \bar{u} & \bar{u} \\ \bar{e} & \bar{e} \end{bmatrix}.
414
421
415
422
后者意味着 $\bar{u}$ 和 $\bar{e}$ 分别是长期失业率和就业率。
416
423
417
- 具体来说,我们有失业率和就业率: 当 $t \to \infty$ 时,$x_t / n_t = A^t n_0 / n_t \to \bar{x}$。
424
+ 具体来说,当 $t \to \infty$ 时,我们有失业率和就业率: $x_t / n_t = A^t n_0 / n_t \to \bar{x}$。
418
425
419
426
为了说明这些比率的动态,令 $\hat{A} := A / (1+g)$ 为 $r_t := x_t/ n_t$ 的转移矩阵。
427
+
420
428
比率的动态遵循
421
429
422
430
$$
423
431
r_{t+1} = \frac{x_{t+1}}{n_{t+1}} = \frac{x_{t+1}}{(1+g) n_{t}} = \frac{A x_t}{(1+g)n_t} = \hat{A} \frac{x_t}{n_t}
424
432
=\hat{A} r_t.
425
433
$$
426
434
427
- 注意到 $\hat{A}$ 的列和都为 1,因此 $r(\hat{A})=1$。
435
+ 注意 $\hat{A}$ 的列和都为 1,因此 $r(\hat{A})=1$。
428
436
429
- 可以验证 $\bar{x}$ 也是 $\hat{A}$ 对应于 $r(\hat{A})$ 的右特征向量,即 $\bar{x} = \hat{A} \bar{x}$。
437
+ 我们可以验证 $\bar{x}$ 也是 $\hat{A}$ 对应于 $r(\hat{A})$ 的右特征向量,即 $\bar{x} = \hat{A} \bar{x}$。
430
438
431
439
此外,对于任何 $r_0 = x_0 / n_0$,当 $t \to \infty$ 时,$\hat{A}^t r_0 \to \bar{x}$,因为上述讨论意味着
432
440
433
441
$$
434
442
r_t = \hat{A}^t r_0 = (1+g)^{-t} A^t r_0 = r(A)^{-t} A^t r_0 \to \begin{bmatrix} \bar{u} & \bar{u} \\ \bar{e} & \bar{e} \end{bmatrix} r_0 = \begin{bmatrix} \bar{u} \\ \bar{e} \end{bmatrix}.
435
443
$$
436
444
437
- 这在下面有所说明。
445
+ 我们用下面的代码来说明
438
446
439
447
``` {code-cell} ipython3
440
448
lm = LakeModel()
@@ -471,7 +479,7 @@ plt.tight_layout()
471
479
plt.show()
472
480
```
473
481
474
- 为了更直观地理解收敛性 ,我们在下面不使用佩龙-弗罗贝尼乌斯定理进一步解释收敛过程。
482
+ 为了能更直观地理解收敛性 ,我们在下面不使用佩龙-弗罗贝尼乌斯定理进一步解释收敛过程。
475
483
476
484
假设 $\hat{A} = P D P^{-1}$ 是可对角化的,其中 $P = [ v_1, v_2] $ 由 $\hat{A}$ 的特征向量 $v_1$ 和 $v_2$ 组成,
477
485
分别对应于特征值 $\gamma_1$ 和 $\gamma_2$,
0 commit comments