常微分方程（三）：高阶线性微分方程

一阶 ODE 只需记住一个数；二阶 ODE 则要记住两个。 正是这多出来的一点自由度，让同一个方程既能描述拨动的吉他弦、汽车的悬挂系统、FM 收音机里的 LC 谐振电路，也能刻画大风中摇晃的高楼。在所有这些情形中，系统行为总归为三种模式：持续振荡、略带过冲地回归平衡，或缓慢爬回原位——而决定具体走向的，始终是那个代数工具：特征方程。

本章将一次性构建完整的理论工具箱。我们会先推导核心方法、证明结构定理，随后不断在不同物理场景中重逢同一幅图景。

本章学习目标#

• 为什么牛顿定律遇上回复力就自然导出二阶方程
• 三大结构定理：叠加原理、Wronskian 判据、通解结构 $y = y_h + y_p$
• 常系数齐次方程的求解：特征多项式法（实根、重根、复根）
• 阻尼比 $\zeta$ 与欠阻尼 / 临界阻尼 / 过阻尼的三分法
• 待定系数法：适用于“简单”外力函数（如指数、三角、多项式）
• 参数变易法：处理任意连续外力的通用方法
• 共振现象及其振幅曲线解读
• 同一方程的电路实现（RLC 电路），以及如何将任意 $n$ 阶 ODE 转化为一阶系统

前置知识#

• 第二章：一阶方法 —— 变量分离法、积分因子法
• 线性代数基础：掌握 $2\times 2$ 行列式与复数运算即可

为什么二阶是“自然单位”#

其中 自然频率 $\omega_0 = \sqrt{k/m}$ ，阻尼比 $\zeta = b/(2\sqrt{mk})$ 。本章几乎所有例子都是该方程的特例。

二阶究竟带来了什么？ 一阶 ODE $\dot x = F(x,t)$ 在时刻 $t$ 的状态仅由单个数值 $x(t)$ 决定；而二阶 ODE 需要两个初始条件 $x(0)$ 和 $\dot x(0)$ ——你可以独立指定初始位置和速度，因此系统能够存储并交换两种形式的“能量”（如动能与势能、电压与电流等）。振荡正是这种能量交换最直观的表现。

解空间的结构#

标准形式#

当 $g \equiv 0$ 时称为齐次方程，否则为非齐次方程。

三条组织全局的定理#

(T1) 叠加原理： 若 $y_1, y_2$ 都是齐次方程的解，则对任意常数 $c_1, c_2$ ，线性组合 $c_1 y_1 + c_2 y_2$ 也是解。（这源于导数的线性性——动手验证一下便知。）

这实际上是存在唯一性定理的另一种表述：给定 $n$ 个初始条件 $y(x_0), y'(x_0), \dots, y^{(n-1)}(x_0)$ ，就能唯一确定这个 $n$ 维向量空间中的一个元素。

其中 $y_p$ 是任意一个特解。原因在于：若 $y$ 和 $\tilde y$ 都是非齐次方程的解，则它们的差 $y - \tilde y$ 满足齐次方程 $L[y - \tilde y] = 0$ 。

Wronskian：用行列式检验线性无关性#

对于 $n=2$ 的情形：$W(y_1, y_2) = y_1 y_2' - y_2 y_1'$ 。

判别准则： 只要在区间内存在某一点 $x_0$ 使得 $W(x_0) \neq 0$ ，这些解就线性无关；根据 Abel 恒等式，此时 Wronskian 在整个区间上处处非零。

Wronskian 恒为 $-1$ ，说明二者处处线性无关，构成方程 $y'' + y = 0$ 的一组基。

恒等于零。因为二者互为标量倍数，张成的空间维数仅为 1，而非 2。

常系数情形：特征方程#

核心技巧#

每个根 $r$ 都对应一个基础解。从根到基础解的映射构成了本节的全部内容。

三种典型情形#

示例： $y'' - 5y' + 6y = 0$ 的特征方程为 $(r-2)(r-3) = 0$ ，故通解为 $y = c_1 e^{2x} + c_2 e^{3x}$ 。

其原理在于：若 $P(r)$ 有二重根，则微分算子可分解为 $(D - r)^2 \cdot Q(D)$ ，而直接计算可得 $(D - r)^2[x e^{rx}] = 0$ 。

示例： $y'' - 4y' + 4y = 0$ 对应 $(r-2)^2 = 0$ ，故通解为 $y = (c_1 + c_2 x)e^{2x}$ 。

为何成立？ 复根对原本贡献 $C_1 e^{(\alpha+i\beta)x} + C_2 e^{(\alpha-i\beta)x}$ ，利用欧拉公式 $e^{i\beta x} = \cos\beta x + i\sin\beta x$ 重新组合，即可得到上述实形式。直观来看，$\alpha$ 控制指数包络的衰减或增长，$\beta$ 决定振荡频率。

示例： $y'' + 2y' + 5y = 0$ 的根为 $r = -1 \pm 2i$ ，故通解为 $y = e^{-x}(c_1\cos 2x + c_2\sin 2x)$ ——一个被指数衰减包络“困住”的正弦波。

在复平面上“读”特征根#

特征根在复平面上的位置，直接揭示了解的定性行为：

• 左半平面（$\Re(r) < 0$ ）：解随时间衰减，系统稳定。
• 右半平面（$\Re(r) > 0$ ）：解发散增长，系统不稳定。
• 虚轴上（$\Re(r) = 0$ ）：纯振荡，处于临界稳定状态。
• 偏离实轴（$\Im(r) \neq 0$ ）：解包含频率为 $|\Im(r)|$ 的振荡成分。

阻尼振动：三分法的图解#

其特征方程 $r^2 + 2\zeta\omega_0 r + \omega_0^2 = 0$ 的判别式为 $4\omega_0^2(\zeta^2 - 1)$ 。$\zeta^2 - 1$ 的符号决定了系统落入哪一类：

工程小注： 汽车悬挂通常调至 $\zeta \approx 0.6$ –$0.7$ ，略低于临界阻尼：车辆能快速平稳，且微小过冲不易察觉。若严格取 $\zeta = 1$ ，末段响应会显得“迟钝”；而 $\zeta = 0.1$ 则会导致过减速带后持续弹跳。相比之下，门吸器则故意设计为过阻尼（$\zeta > 1$ ），以避免关门时发出撞击声。

非齐次方程：待定系数法#

求解流程#

对于常系数非齐次方程 $L[y] = f(x)$ ，步骤如下：

1. 求解齐次方程 $L[y] = 0$ ，得到通解 $y_h$ 。
2. 根据 $f(x)$ 的形式，猜测一个特解 $y_p$ 的结构。
3. 将 $y_p$ 代入原方程，通过比较系数确定未知常数。
4. 合并结果：$y = y_h + y_p$ 。

猜测形式对照表#

共振修正规则： 如果所猜的 $y_p$ 形式本身已是齐次方程的解，则需将其乘以 $x$ （若为双重共振，则乘 $x^2$ ）。否则代入后会得到 $0 = f(x)$ 的矛盾。

完整示例#

求解 $y'' + y' + y = e^{-x/2}\cos x$ 。

解得 $A = \tfrac{2}{37}$ , $B = \tfrac{24}{37}$ 。

共振：当试探解落入齐次解空间#

此时振幅随时间线性增长——能量在每个周期被持续注入却无处耗散。若加入哪怕微小的阻尼，这种无界增长就会转变为有限峰值，峰值频率为 $\omega_r = \omega_0\sqrt{1 - 2\zeta^2}$ ，峰值高度与 $1/\zeta$ 成正比。

参数变易法：当“猜”不再可行#

为何需要它#

待定系数法仅适用于指数函数、多项式、正弦余弦及其乘积构成的外力。面对 $f(x) = \sec x, \tan x, \ln x, e^{x^2}$ 等更复杂的函数，我们需要一种对任意连续外力都有效的通用方法——这就是参数变易法。

公式（二阶情形）#

示例：$y'' + y = \sec x$ #

何时选用哪种方法？ 待定系数法在适用时更快捷：只需猜测形式并求解小型线性方程组。参数变易法虽对任意连续 $f$ 都有效，但需计算两个积分——这些积分本身可能难以解析求解。对于指数、三角或多项式外力，优先使用待定系数法；其余情况则选用参数变易法。

RLC 电路：同一方程的电气化身#

其中 $q(t)$ 是电容上的电荷，$\dot q = i(t)$ 为回路电流。将其与力学标准形式逐项对比：

相应地，标准参数变为：谐振角频率 $\omega_0 = 1/\sqrt{LC}$ ，阻尼比 $\zeta = \tfrac{R}{2}\sqrt{C/L}$ ，以及品质因数 $Q = 1/(2\zeta)$ （$Q$ 越高，共振峰越尖锐）。欠阻尼、临界阻尼、过阻尼的三分法在此完全适用。

化为一阶系统#

该系数矩阵 $\bigl(\begin{smallmatrix} 0 & 1 \\ -b & -a \end{smallmatrix}\bigr)$ 的特征值，恰好是特征多项式 $r^2 + ar + b$ 的根。换言之，本章关于特征根的所有讨论，本质上都是二维线性代数的另一种表述。这一点将在第六章中被充分运用。

总结#

方法选择指南#

根 → 基础解 速查表#

核心思想#

• 二阶 ODE 是描述“惯性 + 回复力”系统的自然语言。
• 特征根在复平面 $\mathbb{C}$ 中的位置，已蕴含了解的全部定性信息。
• Wronskian 是一个行列式工具，用于检验候选解是否张成完整的 $n$ 维解空间。
• 共振发生在试探解与齐次解重合之时；乘以 $x$ 不仅修复了代数矛盾，也从物理上解释了振幅的无界增长。
• 力学振动系统与 RLC 电路遵循完全相同的数学方程，因此物理直觉可无缝迁移。

练习题#

热身。

1. 求 $y'' - 3y' + 2y = 0$ 的通解。
2. 求 $y'' - 4y' + 4y = 0$ 的通解（重根情形）。
3. 解初值问题 $y'' + y = 0$ ，$y(0) = 1$ ，$y'(0) = 0$ 。
4. 解 $y'' + 4y = \sin 2x$ ，并指出所用的共振修正项。

练习。

5. 使用参数变易法，求 $y'' + y = \tan x$ 在区间 $(-\pi/2, \pi/2)$ 上的一个特解。
6. 某 RLC 电路参数为：$R = 100\,\Omega$ ，$L = 0.5\,\text{H}$ ，$C = 10\,\mu\text{F}$ 。计算 $\omega_0$ 、$\zeta$ 和 $Q$ ，并判断其阻尼状态。
7. 证明方程 $\ddot x + 2\zeta\omega_0\dot x + \omega_0^2 x = F_0\cos\omega t$ 的稳态振幅为 $A(\omega) = \frac{F_0}{\sqrt{(\omega_0^2 - \omega^2)^2 + (2\zeta\omega_0\omega)^2}},$ 并验证当 $\zeta < 1/\sqrt{2}$ 时，振幅峰值出现在 $\omega_r = \omega_0\sqrt{1 - 2\zeta^2}$ 。

编程。

8. 复现阻尼状态图：使用 scipy.integrate.solve_ivp 绘制 $\zeta = 0.1, 0.3, 0.7, 1.0, 2.0$ 、$\omega_0 = 2\pi$ 时的阶跃响应。
9. 模拟一个双摆（非线性耦合二阶系统），数值验证两条初始条件相差 $10^{-6}$ 的轨迹在几秒后显著分离。这是你对第九章 （混沌）的初次体验。

下一步#

下一章引入拉普拉斯变换——一个可以把微分方程整体翻译成代数方程的“魔法盒子”。变换之后，ODE 的初值条件、强迫项、阻尼项都变成代数操作，求解过程变成因式分解+查表。当你第一次用它解一个 RLC 电路时，会有“之前为什么不告诉我这个”的感觉。

参考文献#

• Boyce & DiPrima, Elementary Differential Equations and Boundary Value Problems, Wiley (2012)。按定理逐条推进的标准教材。
• Kreyszig, Advanced Engineering Mathematics, Wiley (2011)。工程应用尤为扎实，RLC 电路与机械振动讲解透彻。
• Nagle, Saff & Snider, Fundamentals of Differential Equations, Pearson (2017)。参数变易法部分叙述清晰。
• Strogatz, Nonlinear Dynamics and Chaos (1994)。第 4、5 章从几何角度解读线性稳定性，可直接衔接第七章 。
• MIT OpenCourseWare 18.03，Differential Equations —— Arthur Mattuck 主讲的视频课程。