@@ -106,9 +106,9 @@ import random
106106```
107107
108108(lss_pgs)=
109- ### 基本要素
109+ ### 基本参数
110110
111- 模型的基本要素是
111+ 模型的基本参数是
112112
1131131 . 矩阵 $A, C, G$
1141141 . 冲击分布,我们将其特定为 $N(0,I)$
@@ -121,7 +121,7 @@ import random
121121
122122稍后我们将看到如何计算这些分布及其矩。
123123
124- #### 鞅差异冲击
124+ #### 鞅差分冲击
125125
126126``` {index} single: Linear State Space Models; Martingale Difference Shocks
127127```
139139\mathbb{E} [w_{t+1} | x_t, x_{t-1}, \ldots ] = 0
140140$$
141141
142- 这个条件比 $\{ w_t\} $ 是独立同分布且 $w_ {t+1} \sim N(0,I)$ 的条件要弱 。
142+ 这是一个比 $\{ w_t\} $ 是独立同分布且 $w_ {t+1} \sim N(0,I)$ 还更弱的条件 。
143143
144144### 示例
145145
146- 通过适当选择基本要素 ,各种动态系统都可以用线性状态空间模型来表示。
146+ 通过适当选择基本参数 ,各种动态系统都可以用线性状态空间模型来表示。
147147
148148以下示例有助于突出这一点。
149149
150- 这些示例也说明了 "找到状态是一门艺术"这一智慧格言。
150+ 这些示例也阐明了 "找到状态是一门艺术"这一智慧格言。
151151
152152(lss_sode)=
153153#### 二阶差分方程
@@ -225,7 +225,7 @@ G = [0, 1, 0]
225225plot_lss(A, C, G)
226226```
227227
228- 稍后将要求您重新创建这个图形 。
228+ 后续将请你重现该图示 。
229229
230230#### 单变量自回归过程
231231
@@ -265,7 +265,7 @@ C = \begin{bmatrix}
265265 \end{bmatrix}
266266$$
267267
268- 矩阵$A$具有向量$\begin{bmatrix}\phi_1 & \phi_2 & \phi_3 & \phi_4 \end{bmatrix}$的 * 伴随矩阵* 形式 。
268+ 矩阵$A$具有向量$\begin{bmatrix}\phi_1 & \phi_2 & \phi_3 & \phi_4 \end{bmatrix}$* 伴随矩阵* 的形式 。
269269
270270下图显示了当
271271
274274y_{-3} = 1
275275$$
276276
277- 时该过程的动态变化
277+ 该过程的动态变化
278278
279279(lss_uap_fig)=
280280``` {code-cell} ipython3
@@ -367,7 +367,7 @@ A = \begin{bmatrix}
367367 \end{bmatrix}
368368$$
369369
370- 容易验证$A^4 = I$,这意味着$x_t$是严格周期的,周期为4 :[ ^ foot1 ]
370+ 容易验证$A^4 = I$,这意味着$x_t$是一个周期严格为 4 的序列 :[ ^ foot1 ]
371371
372372$$
373373x_{t+4} = x_t
505505
506506其中 $x_ {1t}$ 是 $x_t$ 的第一个元素。
507507
508- 右边的第一项是鞅差的累积和,因此是一个[ 鞅] ( https://en.wikipedia.org/wiki/Martingale_%28probability_theory%29 ) 。
508+ 右边的第一项是鞅差的累积和,因此是一个[ 鞅] ( https://baike.baidu.com/item/%E9%9E%85/56830519 ) 。
509509
510510第二项是时间的平移线性函数。
511511
@@ -581,7 +581,7 @@ $y_t$ 的方差-协方差矩阵可以很容易地证明为
581581``` {index} single: Linear State Space Models; Distributions
582582```
583583
584- 一般来说,知道随机向量的均值和方差-协方差矩阵并不如知道完整的分布那么好 。
584+ 一般来说,知道一个随机向量的均值和方差-协方差矩阵,并不完全等同于知道其完整分布 。
585585
586586然而,在某些情况下,仅仅这些矩就能告诉我们所需要知道的一切。
587587
@@ -594,7 +594,7 @@ $y_t$ 的方差-协方差矩阵可以很容易地证明为
5945941 . 我们对基本量的高斯分布假设
5955951 . 正态性在线性运算下得以保持
596596
597- 事实上,[ 众所周知] ( https://en.wikipedia.org/wiki/Multivariate_normal_distribution#Affine_transformation )
597+ 事实上,[ 众所周知] ( https://baike.baidu.com/item/%E7%9F%A9%E9%98%B5%E6%AD%A3%E6%80%81%E5%88%86%E5%B8%83/22800089?fromModule=lemma_inlink )
598598
599599``` {math}
600600:label: lss_glig
@@ -628,7 +628,7 @@ x_t \sim N(\mu_t, \Sigma_t)
628628y_t \sim N(G \mu_t, G \Sigma_t G')
629629```
630630
631- ### 集合解释
631+ ### 系综解释
632632
633633我们应该如何解释{eq}` lss_mgs_x ` --{eq}` lss_mgs_y ` 定义的分布?
634634
@@ -744,23 +744,23 @@ plt.show()
744744
745745通过观察图形并尝试不同的参数,你可以理解总体分布是如何依赖于{ref}` 上面列出的<lss_pgs> ` 模型基本要素的,这种依赖关系是通过分布的参数体现的。
746746
747- #### 集合均值
747+ #### 系综均值
748748
749749在前面的图中,我们通过以下方式近似了$y_T$的总体分布:
750750
7517511 . 生成$I$条样本路径(即时间序列),其中$I$是一个很大的数
7527521 . 记录每个观测值$y^i_T$
7537531 . 对这个样本制作直方图
754754
755- 正如直方图近似总体分布一样,* 集合 * 或* 横截面* 平均值
755+ 正如直方图近似总体分布一样,* 系综 * 或* 横截面* 平均值
756756
757757$$
758758\bar y_T := \frac{1}{I} \sum_{i=1}^I y_T^i
759759$$
760760
761- 近似期望值$\mathbb{E} [ y_T] = G \mu_T$(这是由大数定律所暗示的 )。
761+ 近似期望值$\mathbb{E} [ y_T] = G \mu_T$(由大数定律所示 )。
762762
763- 这里是一个模拟,比较了在时间点$t=0,\ldots,50$处的集合平均值和总体均值 。
763+ 这里是一个模拟,比较了在时间点$t=0,\ldots,50$处的系综平均值和总体均值 。
764764
765765参数与前面的图表相同,样本量相对较小($I=20$)。
766766
@@ -799,7 +799,7 @@ ax.legend(ncol=2)
799799plt.show()
800800```
801801
802- $x_t$ 的集成均值为
802+ $x_t$ 的系综均值为
803803
804804$$
805805\bar x_T := \frac{1}{I} \sum_{i=1}^I x_T^i \to \mu_T
831831
832832要计算 $x_0, x_1, \ldots, x_T$ 的联合分布,回想
833833
834- 联合密度和条件密度通过以下规则相关联
835-
834+ 联合密度与条件密度存在如下关系
836835$$
837836p(x, y) = p(y \, | \, x) p(x)
838837\qquad \text{(联合密度 }=\text{ 条件密度 }\times\text{ 边际密度)}
846845p(x_0, x_1, \ldots, x_T) = p(x_0) \prod_{t=0}^{T-1} p(x_{t+1} \,|\, x_t)
847846$$
848847
849- 边际密度 $p(x_0)$ 就是原始的 $N(\mu_0, \Sigma_0)$。
848+ 边缘密度 $p(x_0)$ 就是原始的 $N(\mu_0, \Sigma_0)$。
850849
851850根据{eq}` st_space_rep ` ,条件密度为
852851
864863\Sigma_{t+j, t} := \mathbb{E} [ (x_{t+j} - \mu_{t+j})(x_t - \mu_t)' ]
865864```
866865
867- 基本计算表明
866+ 通过基本计算可知
868867
869868``` {math}
870869:label: eqnautocov
10331032
10341033由于状态向量中的第一个分量是常数,我们永远不会有 $\mu_t \to 0$。
10351034
1036- 我们如何找到满足常数状态分量的平稳解 ?
1035+ 我们如何求解具有常数状态分量的平稳解 ?
10371036
10381037#### 具有常数状态分量的过程
10391038
@@ -1116,11 +1115,11 @@ $\Sigma_t$ 和 $\Sigma_{t+j,t}$ 的平稳值满足
11161115
11171116假设我们正在处理一个协方差平稳过程。
11181117
1119- 在这种情况下,我们知道当样本量 $I$ 趋向无穷时,整体均值将收敛到 $\mu_ {\infty}$。
1118+ 在这种情况下,我们知道当样本量 $I$ 趋向无穷时,系综均值将收敛到 $\mu_ {\infty}$。
11201119
11211120#### 时间平均
11221121
1123- 理论上,跨模拟的整体平均很有趣,但在现实生活中,我们通常只观察到 * 单个* 实现 $\{ x_t, y_t\} _ {t=0}^T$。
1122+ 虽然模拟的系综平均在理论上很有意义,但在现实中,我们通常只能观测到 * 单个* 实现 $\{ x_t, y_t\} _ {t=0}^T$。
11241123
11251124因此现在让我们取一个单独的实现并形成时间序列平均值
11261125
11511150
11521151在某些情况下,观测方程 $y_t = Gx_t$ 会被修改以包含一个误差项。
11531152
1154- 这个误差项通常表示真实状态只能被不完美地观测到这一概念 。
1153+ 通常这个误差项反映了真实状态只能被不完全观测的特性 。
11551154
11561155为了在观测中引入误差项,我们引入
11571156
11721171
11731172序列$\{ v_t\} $被假定与$\{ w_t\} $相互独立。
11741173
1175- 过程$\{ x_t\} $不会被观测方程中的噪声所修改 ,其矩、分布和稳定性特征保持不变 。
1174+ 过程$\{ x_t\} $不受观测方程噪声的影响 ,其矩、分布及稳定性性质均保持不变 。
11761175
11771176$y_t$的无条件矩从{eq}` lss_umy ` 和{eq}` lss_uvy ` 现在变为
11781177
@@ -1240,7 +1239,7 @@ x_{t+j} = A^j x_t + A^{j-1} C w_{t+1} + A^{j-2} C w_{t+2} +
12401239\cdots + A^0 C w_{t+j}
12411240$$
12421241
1243- 根据IID特性,当前和过去的状态值不能提供关于未来冲击值的任何信息 。
1242+ 根据IID性质,当前及过去的状态值无法提供关于冲击未来值的任何信息 。
12441243
12451244因此$\mathbb{E}_ t[ w_ {t+k}] = \mathbb{E}[ w_ {t+k}] = 0$。
12461245
@@ -1303,15 +1302,15 @@ $V_j$收敛的一个充分条件是$A$的特征值的模都严格小于1。
13031302(lm_fgs)=
13041303## 代码
13051304
1306- 我们之前的模拟和计算都基于 [ QuantEcon.py] ( http://quantecon.org/quantecon-py ) 包中的[ lss.py] ( https://github.com/QuantEcon/QuantEcon.py/blob/master/quantecon/lss.py ) 文件。
1305+ 我们之前的模拟和计算都是基于 [ QuantEcon.py] ( http://quantecon.org/quantecon-py ) 包中的[ lss.py] ( https://github.com/QuantEcon/QuantEcon.py/blob/master/quantecon/lss.py ) 文件。
13071306
13081307该代码实现了一个用于处理线性状态空间模型的类(包括模拟、计算矩等功能)。
13091308
13101309你可能不太熟悉的一个Python结构是在` moment_sequence() ` 方法中使用生成器函数。
13111310
13121311如果你忘记了生成器函数是如何工作的,请回去[ 阅读相关文档] ( https://python-programming.quantecon.org/python_advanced_features.html#generators ) 。
13131312
1314- 使用示例在练习的解答中给出 。
1313+ 具体示例详见习题解答部分 。
13151314
13161315## 练习
13171316
0 commit comments