@@ -172,13 +172,13 @@ plt.show()
172
172
173
173
## 边际分布
174
174
175
- 现在让我们来看看某个固定时间点 $T$ 时 $X_T$ 的边际分布 $\psi_T$。
175
+ 现在让我们来看看某一固定时间点 $T$ 时 $X_T$ 的边际分布 $\psi_T$。
176
176
177
- 我们将通过在给定初始条件 $X_0$ 的情况下生成多个 $X_T$ 的样本来实现这一点 。
177
+ 我们将通过在给定初始条件 $X_0$ 的情况下,生成多个 $X_T$ 的样本来实现 。
178
178
179
179
通过这些 $X_T$ 的样本,我们可以构建其分布 $\psi_T$ 的图像。
180
180
181
- 这里是一个可视化示例,其中 $T=50$。
181
+ 下面是 $T=50$的情况下,一个可视化示例 。
182
182
183
183
``` {code-cell} ipython3
184
184
T = 50
@@ -219,7 +219,7 @@ axes[1].hist(sample,
219
219
plt.show()
220
220
```
221
221
222
- 通过绘制更多样本 ,我们可以得到一个更清晰的图像
222
+ 通过抽取更多样本 ,我们可以得到一个更清晰的图像
223
223
224
224
``` {code-cell} ipython3
225
225
T = 50
@@ -241,18 +241,18 @@ ax.hist(sample,
241
241
plt.show()
242
242
```
243
243
244
- 请注意分布呈双峰
244
+ 注意到分布呈双峰
245
245
246
- * 大多数公司已经补货两次,但少数公司只补货一次 (见上图路径)。
247
- * 第二类公司的库存较低 。
246
+ * 大多数公司已经补了两次货,但也有少部分公司只补货一次 (见上图路径)。
247
+ * 第二种公司的库存较少 。
248
248
249
- 我们也可以使用 [ 核密度估计] ( https://en.wikipedia.org/wiki/Kernel_density_estimation ) 来近似这个分布。
249
+ 我们还可以使用 [ 核密度估计] ( https://baike.baidu.com/item/核密度估计/10349033 ) 来近似这个分布。
250
250
251
251
核密度估计可以被理解为平滑的直方图。
252
252
253
- 当被估计的分布可能是平滑的时候 ,核密度估计比直方图更可取。
253
+ 当被估计的分布很可能是平滑的时候 ,核密度估计比直方图更可取。
254
254
255
- 我们将使用[ scikit-learn] ( https://scikit-learn.org/stable/ ) 中的核密度估计器
255
+ 我们将使用[ scikit-learn] ( https://scikit-learn.org/stable/ ) 中的核密度估计量
256
256
257
257
``` {code-cell} ipython3
258
258
from sklearn.neighbors import KernelDensity
@@ -282,30 +282,30 @@ plt.show()
282
282
283
283
这个模型是渐近平稳的,具有唯一的平稳分布。
284
284
285
- (有关平稳性的背景讨论 ,请参见{doc}`我们关于AR(1)过程的讲座 <intro:ar1_processes>`——基本概念是相同的。)
285
+ (作为背景知识,有关平稳性的讨论 ,请参见{doc}`我们关于AR(1)过程的讲座 <intro:ar1_processes>`——基本概念是相同的。)
286
286
287
- 特别是,边际分布序列$\{\psi_t\}$正在收敛到一个唯一的极限分布,该分布不依赖于初始条件 。
287
+ 特别是,边际分布序列$\{\psi_t\}$正在收敛到一个唯一的极限分布,且该分布不依赖于初始条件 。
288
288
289
- 虽然我们在这里不会证明这一点,但我们可以通过模拟来研究它 。
289
+ 虽然我们不会在此证明这一点,但我们可以通过模拟来研究这一性质 。
290
290
291
- 你的任务是根据上述讨论 ,在时间点$t = 10, 50, 250, 500, 750$生成并绘制序列$\{\psi_t\}$。
291
+ 你的任务是,根据上述讨论 ,在时间点$t = 10, 50, 250, 500, 750$生成并绘制序列$\{\psi_t\}$。
292
292
293
- (核密度估计器可能是呈现每个分布的最佳方式 。)
293
+ (核密度估计量可能是呈现每个分布最佳的方式 。)
294
294
295
- 你应该能看到收敛性,体现在连续分布之间的差异越来越小 。
295
+ 你应该能看到收敛性,体现在两个连续分布之间的差异越来越小 。
296
296
297
- 尝试不同的初始条件来验证,从长远来看,分布在不同初始条件下是不变的 。
297
+ 尝试不同的初始条件来验证,长期来看,不同初始条件下分布是不变的 。
298
298
```
299
299
300
300
``` {solution-start} id_ex1
301
301
:class: dropdown
302
302
```
303
303
304
- 以下是一个可能的解决方案 :
304
+ 以下是一种可能的解法 :
305
305
306
- 这些计算涉及大量的CPU周期,所以我们试图高效地编写代码 。
306
+ 因为其中的计算涉及大量的CPU周期,所以我们试图更高效地编写代码 。
307
307
308
- 这意味着编写一个专门的函数 ,而不是使用上面的类。
308
+ 也就是编写一个专门的函数 ,而不是使用上面的类。
309
309
310
310
``` {code-cell} ipython3
311
311
s, S, mu, sigma = firm.s, firm.S, firm.mu, firm.sigma
@@ -354,11 +354,11 @@ ax.legend()
354
354
plt.show()
355
355
```
356
356
357
- 注意到在 $t=500$ 或 $t=750$ 时密度几乎不再变化。
357
+ 注意到,在 $t=500$ 或 $t=750$ 时密度几乎不再变化。
358
358
359
- 我们已经得到了平稳密度的合理近似 。
359
+ 我们得到了平稳密度的合理近似 。
360
360
361
- 你可以通过测试几个不同的初始条件来确信初始条件并不重要 。
361
+ 你可以通过测试多个不同的初始条件,来确定初始条件确实不重要 。
362
362
363
363
例如,尝试用所有公司从 $X_0 = 20$ 或 $X_0 = 80$ 开始重新运行上面的代码。
364
364
@@ -368,7 +368,7 @@ plt.show()
368
368
``` {exercise}
369
369
:label: id_ex2
370
370
371
- 使用模拟计算从 $X_0 = 70$ 开始的公司在前50个周期内需要订货两次或更多次的概率 。
371
+ 使用模拟的方式,计算从 $X_0 = 70$ 开始的公司在前50个周期内需要订货两次或更多的概率 。
372
372
373
373
你需要一个较大的样本量来获得准确的结果。
374
374
```
@@ -377,9 +377,9 @@ plt.show()
377
377
:class: dropdown
378
378
```
379
379
380
- 这是一个解决方案 。
380
+ 这里是一种解法 。
381
381
382
- 同样 ,由于计算量相对较大,我们编写了一个专门的函数而不是使用上面的类。
382
+ 同样地 ,由于计算量相对较大,我们编写了一个专门的函数而不是使用上面的类。
383
383
384
384
我们还将使用跨公司的并行化处理。
385
385
@@ -420,7 +420,7 @@ print(f"至少发生两次缺货的频率 = {freq}")
420
420
421
421
根据你的系统配置,运行速度的差异可能会很大。
422
422
423
- (在我们的台式机上 ,速度提升了5倍。)
423
+ (在台式机上 ,速度提升了5倍。)
424
424
425
425
``` {solution-end}
426
426
```
0 commit comments