KNOWECON WIKI

ARTICLE

差分方程

差分方程 (Difference Equation) 差分方程 (Difference Equation) 是描述一个序列中某一项与其前若干项之间关系的方程。与微分方程刻画连续时间动态不同，差分方程是分析离散时间动态系统的核心数学工具，广泛应用于宏观经济学、金融学、计量经济学和数理经济学等领域。任何涉及离散时间跨期决策、动态优化或时间序列建模的问题，几乎都离

浏览 0 更新 2025-11-08

差分方程 (Difference Equation)

差分方程 (Difference Equation) 是描述一个序列中某一项与其前若干项之间关系的方程。与微分方程刻画连续时间动态不同，差分方程是分析离散时间动态系统的核心数学工具，广泛应用于宏观经济学、金融学、计量经济学和数理经济学等领域。任何涉及离散时间跨期决策、动态优化或时间序列建模的问题，几乎都离不开差分方程。

差分方程的本质在于，它不直接给出序列的通项 yt y_t ，而是将不同时点的值关联起来。例如，一阶差分方程的形式为 yt=f(yt−1,t) y_t = f(y_{t-1}, t) ，而微分方程的对应形式为 y˙(t)=f(y(t),t) \dot{y}(t) = f(y(t), t) 。前者刻画跳跃式的逐期演化，后者刻画瞬时变化。

基本概念与分类

差分算子

为表述方便，首先定义差分算子 (Difference Operator)。对于序列 {yt}t=0∞ \{y_t\}_{t=0}^{\infty} ：

一阶前向差分 (Forward Difference)：

Δyt=yt+1−yt\Delta y_t = y_{t+1} - y_t

一阶后向差分 (Backward Difference)：

∇yt=yt−yt−1\nabla y_t = y_t - y_{t-1}

高阶差分可递归定义，如二阶前向差分：

Δ2yt=Δ(Δyt)=yt+2−2yt+1+yt\Delta^2 y_t = \Delta(\Delta y_t) = y_{t+2} - 2y_{t+1} + y_t

在经济模型中，后向差分更为常见——大多数经济变量在时刻 t t 的值取决于 t−1 t-1 期的值。但通过指标平移，前向与后向形式等价。

差分方程的分类

一个 k k 阶差分方程的一般形式为：

F(t,yt,yt+1,…,yt+k)=0F(t, y_t, y_{t+1}, \ldots, y_{t+k}) = 0

其中 k k 为方程的阶数 (Order)，即方程中涉及的时期跨度的最大值。按性质可分为：

  1. 线性 vs 非线性：若 F F 关于 yt y_t 及其各阶滞后（或超前）是线性的，则为线性差分方程；否则为非线性的。例如 yt=ayt−1+b y_t = a y_{t-1} + b 是线性一阶差分方程，而 yt=ryt−1(1−yt−1) y_t = r y_{t-1}(1 - y_{t-1}) （Logistic方程）是非线性的。
  1. 齐次 vs 非齐次：对于线性差分方程，若不含与 yt y_t 无关的外生项（即常数为零），则为齐次方程；否则为非齐次方程。例如，yt−2yt−1=0 y_t - 2y_{t-1} = 0 是齐次的，而 yt−2yt−1=3 y_t - 2y_{t-1} = 3 是非齐次的。
  1. 自治 vs 非自治：若方程中不显含时间变量 t t ，则为自治方程；否则为非自治方程。自治系统的动力学行为不随时间平移而改变，分析更为简便。

线性常系数差分方程

这是经济学中应用最广泛的一类差分方程。一个 k k 阶线性常系数差分方程的形式为：

yt=a1yt−1+a2yt−2+⋯+akyt−k+bty_t = a_1 y_{t-1} + a_2 y_{t-2} + \cdots + a_k y_{t-k} + b_t

其中 a1,…,ak a_1, \ldots, a_k 为常数系数，bt b_t 为驱动项（外生项或输入项）。

一阶线性差分方程

一阶情形是最基础的：

yt=ayt−1+by_t = a y_{t-1} + b

迭代求解是理解一阶方程最直观的方法。从初始条件 y0 y_0 出发：

y1=ay0+by2=ay1+b=a2y0+(a+1)b⋮yt=aty0+b∑i=0t−1ai\begin{aligned} y_1 &= a y_0 + b \\ y_2 &= a y_1 + b = a^2 y_0 + (a+1)b \\ &\vdots \\ y_t &= a^t y_0 + b \sum_{i=0}^{t-1} a^i \end{aligned}

若 a≠1 a \neq 1 ，利用等比数列求和公式：

yt=aty0+b⋅1−at1−ay_t = a^t y_0 + b \cdot \frac{1 - a^t}{1 - a}

若 a=1 a = 1 ，则 yt=y0+bt y_t = y_0 + bt （等差数列）。

通解结构：一般解 = 齐次解 + 特解。齐次方程 yt=ayt−1 y_t = a y_{t-1} 的解为 yth=Cat y_t^h = C a^t （C C 为待定常数）。特解（稳态）通过设 yt=yt−1=yˉ y_t = y_{t-1} = \bar{y} 求得：若 a≠1 a \neq 1 ，yˉ=b/(1−a) \bar{y} = b/(1-a) 。

稳定性：当 ∣a∣<1 |a| < 1 时，齐次解随 t→∞ t \to \infty 而衰减至零，系统从任意初始点收敛到稳态 yˉ \bar{y} ，此时方程是渐近稳定的。当 ∣a∣>1 |a| > 1 时，系统发散。当 a=1 a = 1 时，系统无稳态，yt y_t 线性增长。另外，若 a<0 a < 0 ，系统表现为振荡收敛（∣a∣<1 |a|<1 ）或振荡发散（∣a∣>1 |a|>1 ）。

二阶线性差分方程

二阶齐次方程的标准形式为：

yt+a1yt−1+a2yt−2=0y_t + a_1 y_{t-1} + a_2 y_{t-2} = 0

特征方程法是最标准的解法。假设解具有指数形式 yt=λt y_t = \lambda^t ，代入得特征方程：

λ2+a1λ+a2=0\lambda^2 + a_1 \lambda + a_2 = 0

令其根为 λ1,λ2 \lambda_1, \lambda_2 （可以是实根或共轭复根）：

情形一：两个不等的实根 (λ1≠λ2 \lambda_1 \neq \lambda_2 )

通解为：

yt=C1λ1t+C2λ2ty_t = C_1 \lambda_1^t + C_2 \lambda_2^t

情形二：重实根 (λ1=λ2=λ \lambda_1 = \lambda_2 = \lambda )

通解为：

yt=(C1+C2t)λty_t = (C_1 + C_2 t) \lambda^t

情形三：共轭复根 (λ1,2=re±iθ \lambda_{1,2} = r e^{\pm i\theta} )

通解为：

yt=rt(C1cos⁡θt+C2sin⁡θt)y_t = r^t (C_1 \cos \theta t + C_2 \sin \theta t)

其中 r=∣λ∣=a2 r = |\lambda| = \sqrt{a_2} 为模长，θ \theta 为辐角。

稳定性条件：当且仅当所有特征根的模长（绝对值）都严格小于 1 时，系统是渐近稳定的——即无论从何种初始条件出发，yt y_t 最终趋于零（对自治系统而言，趋于唯一的稳态）。对一阶方程，条件为 ∣a∣<1 |a| < 1 ；对二阶方程，条件归结为 ∣a2∣<1 |a_2| < 1 且 1+a1+a2>0 1 + a_1 + a_2 > 0 且 1−a1+a2>0 1 - a_1 + a_2 > 0 。复根时出现阻尼振荡（r<1 r<1 ），这是经济学中周期现象的重要来源。

对于非齐次方程 yt+a1yt−1+a2yt−2=bt y_t + a_1 y_{t-1} + a_2 y_{t-2} = b_t ，通解为齐次通解加上一个特解。当 bt b_t 为常数时，特解为 yˉ=b/(1+a1+a2) \bar{y} = b/(1 + a_1 + a_2) （分母非零时）。

线性差分方程组

许多经济模型涉及多个相互依存的变量，需要用差分方程组来描述。一阶线性差分方程组的标准形式为：

yt=Ayt−1+b\mathbf{y}_t = A \mathbf{y}_{t-1} + \mathbf{b}

其中 yt \mathbf{y}_t 为 n×1 n \times 1 向量，A A 为 n×n n \times n 系数矩阵，b \mathbf{b} 为常数向量。

求解思路与标量情形类似：稳态满足 yˉ=Ayˉ+b \bar{\mathbf{y}} = A \bar{\mathbf{y}} + \mathbf{b} ，即 yˉ=(I−A)−1b \bar{\mathbf{y}} = (I - A)^{-1}\mathbf{b} （当 I−A I-A 可逆时）。齐次部分 yth=Aty0h \mathbf{y}_t^h = A^t \mathbf{y}_0^h 的动态行为由 A A 的特征值决定：若所有特征值的模长都小于 1，系统渐近稳定。

经济学应用

蛛网模型

蛛网模型是差分方程在微观经济学中的经典应用。假设农产品的供给取决于上一期价格（生产时滞），需求取决于当期价格：

Qts=−γ+δPt−1,Qtd=α−βPtQ_t^s = -\gamma + \delta P_{t-1}, \quad Q_t^d = \alpha - \beta P_t

市场出清条件 Qts=Qtd Q_t^s = Q_t^d 导出：

Pt=α+γβ−δβPt−1P_t = \frac{\alpha + \gamma}{\beta} - \frac{\delta}{\beta} P_{t-1}

这是一阶线性差分方程。均衡价格 Pˉ=(α+γ)/(β+δ) \bar{P} = (\alpha + \gamma)/(\beta + \delta) 。若 δ/β<1 \delta/\beta < 1 （供给弹性小于需求弹性绝对值），价格收敛到均衡（蛛网向内收敛）；若 δ/β>1 \delta/\beta > 1 ，价格发散；若相等则产生持续振荡。

萨缪尔森乘数-加速数模型

萨缪尔森将乘数效应与加速原理结合，解释经济周期：

Yt=Ct+It+G0Ct=cYt−1(0<c<1)It=v(Ct−Ct−1)(v>0)\begin{aligned} Y_t &= C_t + I_t + G_0 \\ C_t &= c Y_{t-1} \quad (0 < c < 1) \\ I_t &= v (C_t - C_{t-1}) \quad (v > 0) \end{aligned}

消去 Ct C_t 和 It I_t 得到关于 Yt Y_t 的二阶差分方程：

Yt=c(1+v)Yt−1−cvYt−2+G0Y_t = c(1+v) Y_{t-1} - c v Y_{t-2} + G_0

根据参数 c c 和 v v 的取值，模型可以产生单调收敛、阻尼振荡、爆炸振荡等多种动态行为，成功刻画了经济波动的内生机理。

离散时间索洛增长模型

索洛增长模型在离散时间下的资本积累方程为：

Kt+1=sF(Kt,Lt)+(1−δ)KtK_{t+1} = s F(K_t, L_t) + (1-\delta) K_t

人均形式（kt=Kt/Lt k_t = K_t/L_t ，设人口增长率 n n ）：

kt+1=sf(kt)+(1−δ)kt1+nk_{t+1} = \frac{s f(k_t) + (1-\delta) k_t}{1+n}

这是一阶非线性差分方程。对于柯布-道格拉斯生产函数 f(k)=kα f(k) = k^\alpha ，稳态满足 skα=(n+δ)k s k^\alpha = (n+\delta)k ，即 kˉ=(s/(n+δ))1/(1−α) \bar{k} = (s/(n+\delta))^{1/(1-\alpha)} 。在稳态附近线性化可以分析收敛速度：kt k_t 以速率 λ=(1−α)(n+δ)/(1+n) \lambda = (1-\alpha)(n+\delta)/(1+n) 向稳态收敛。

ARMA模型与时间序列

在计量经济学中，ARMA模型本质上就是随机差分方程。AR(1)过程：

yt=ϕyt−1+εt,εt∼WN(0,σ2)y_t = \phi y_{t-1} + \varepsilon_t, \quad \varepsilon_t \sim \text{WN}(0, \sigma^2)

当 ∣ϕ∣<1 |\phi| < 1 时，过程是协方差平稳的（这是差分方程稳定性条件在随机环境中的对应）。更一般的 ARMA(p,q) 过程涉及 p p 阶自回归差分结构，其平稳性条件由对应特征方程的所有根落在单位圆外（等价于差分方程特征根模长小于 1）来保证。

与微分方程的联系

差分方程与微分方程之间存在深刻的对应：微分方程 y˙=f(y) \dot{y} = f(y) 的欧拉离散化 yt+1−yth=f(yt) \frac{y_{t+1} - y_t}{h} = f(y_t) 恰为一阶差分方程。反过来，当时间间隔趋近于零时，适当的差分方程趋近于对应的微分方程。这种联系意味着，许多连续时间经济模型可以数值求解为差分方程；同理，差分方程的稳定性理论（特征根分析）也与微分方程的稳定性理论（特征值实部符号分析）形成镜像：微分方程要求特征值实部为负以保证稳定，差分方程要求特征根模长小于 1。