微分方程 (DEs) 简介
你好！欢迎来到进阶纯数学 2 (Further Pure Mathematics 2) 的最后一个章节：微分方程 (Differential Equations)。如果这一章初看起来有点复杂，请别担心；微分方程可以说是高等数学中最强大的工具，它将纯微积分直接与现实世界的物理、工程学和人口模型联系在了一起。
在 A-Level 数学 (9709) 中，你已经学习了如何求解简单的可分离变量微分方程。在进阶数学中，我们将为你提供处理两种主要复杂微分方程的技巧：一阶线性方程和二阶线性方程，并利用积分因子 (IF)、通解函数 (CF) 和特解 (PI) 等专业工具来解决它们。
让我们开始吧！
1. 一阶线性微分方程
一阶微分方程仅涉及 \(\frac{dy}{dx}\)、\(y\) 以及关于 \(x\) 的函数。求解线性形式的关键技巧是使用积分因子 (Integrating Factor, IF)。
1.1 标准形式与积分因子
一阶线性微分方程必须整理为标准形式：
$$ \frac{dy}{dx} + P(x)y = Q(x) $$
其中 \(P(x)\) 和 \(Q(x)\) 是关于 \(x\) 的函数（或常数）。
关键术语：积分因子 (\(\mu\))
积分因子是一个函数 \(\mu\)，我们将整个方程乘以它，使得方程左侧 (LHS) 恰好成为一个乘积的导数：\(\frac{d}{dx} (\mu y)\)。
$$ \mu = e^{\int P(x) dx} $$
类比：你可以把积分因子看作“万能胶水”。它能确保当你将其与微分方程左侧相乘时，所有项都能整齐地粘合在一起，转化为一个易于积分的导数形式。
1.2 使用积分因子的逐步求解方法
- 标准化：确保方程形式为 \(\frac{dy}{dx} + P(x)y = Q(x)\)。
- 确定 \(P(x)\)：找出乘以 \(y\) 的函数。
- 计算积分因子：求出 \(\mu = e^{\int P(x) dx}\)。（注意：此处无需添加任意常数 \(C\)。）
- 相乘：将标准化后的微分方程两边同时乘以积分因子 \(\mu\)。
- 简化左侧：左侧必须简化为 \(\frac{d}{dx} (\mu y)\)。如果不成立，请检查你的积分因子计算是否有误！
- 积分：对等式两边关于 \(x\) 进行积分： $$ \mu y = \int Q(x) \mu dx + C $$
- 求解 \(y\)：分离出 \(y\)，从而得到通解 (General Solution)。
快速回顾：积分因子法
| 形式 | 积分因子 (\(\mu\)) | 积分后的形式 | | :--- | :--- | :--- | | \(\frac{dy}{dx} + P(x)y = Q(x)\) | \(e^{\int P(x) dx}\) | \(\mu y = \int Q(x) \mu dx + C\) |
2. 二阶线性微分方程（常系数）
二阶线性微分方程的形式为： $$ a\frac{d^2y}{dx^2} + b\frac{dy}{dx} + cy = f(x) $$
通解 (GS) 总是由两部分之和组成：
$$ y_{GS} = y_c + y_p $$
$$ y_{GS} = \text{补余函数 (CF)} + \text{特解 (PI)} $$
2.1 第一部分：补余函数 (\(y_c\))
补余函数 (CF) 是齐次方程（即 \(f(x)=0\)）的解： $$ a\frac{d^2y}{dx^2} + b\frac{dy}{dx} + cy = 0 $$
我们通过建立辅助方程 (Auxiliary Equation, AE) 来求解，这是一个二次方程，通过将导数替换为变量（通常为 \(m\)）得到：
$$ am^2 + bm + c = 0 $$
辅助方程的根的性质决定了补余函数的形式。
情况 1：两个不同的实根 (\(m_1 \neq m_2\))
如果判别式 (\(b^2 - 4ac\)) 大于零，根为实数且不相等。
$$ y_c = Ae^{m_1 x} + Be^{m_2 x} $$
例如：若根为 \(m=3\) 和 \(m=-1\)，则 \(y_c = Ae^{3x} + Be^{-x}\)。
情况 2：重实根 (\(m_1 = m_2 = m\))
如果判别式等于零，则存在一个重根 \(m\)。
$$ y_c = (A + Bx)e^{mx} $$
记忆小贴士：当出现重根时，在第二项中引入因子 \(x\)，以确保两个解在数学上是线性无关的。
情况 3：共轭复根 (\(m = \alpha \pm i\beta\))
如果判别式小于零，根为复数。
$$ y_c = e^{\alpha x}(A \cos \beta x + B \sin \beta x) $$
其中 \(\alpha\) 是实部，\(\beta\) 是虚部（使用 \(\beta\) 的正值）。
例如：若根为 \(m = 2 \pm 3i\)，则 \(\alpha=2, \beta=3\)，故 \(y_c = e^{2x}(A \cos 3x + B \sin 3x)\)。
关键点 (CF)：补余函数包含两个任意常数 (\(A\) 和 \(B\))，因为它源于二阶方程。其形式完全取决于辅助方程的根。
2.2 第二部分：特解 (\(y_p\))
特解 (PI) 是满足非齐次方程 \(a\frac{d^2y}{dx^2} + b\frac{dy}{dx} + cy = f(x)\) 的任意解。它不含任意常数。
我们通过根据 \(f(x)\) 的形式提出一个试探函数 (Trial Function) 来寻找特解。
| 若 \(f(x)\) 为... | 尝试 \(y_p\) (试探函数) | 示例 |
|---|---|---|
| \(n\) 次多项式 | \(n\) 次通项多项式 | 若 \(f(x) = x^2 - 1\)，尝试 \(y_p = Px^2 + Qx + R\) |
| 指数函数 \(ke^{bx}\) | \(Pe^{bx}\) | 若 \(f(x) = 5e^{3x}\)，尝试 \(y_p = Pe^{3x}\) |
| 三角函数 \(a \cos px + b \sin px\) | \(P \cos px + Q \sin px\) | 若 \(f(x) = \sin 2x\)，尝试 \(y_p = P \cos 2x + Q \sin 2x\) |
求解常数 \(P, Q, R\) 的步骤：
- 提出试探函数 \(y_p\)。
- 对 \(y_p\) 进行一阶导 (\(\frac{dy_p}{dx}\)) 和二阶导 (\(\frac{d^2y_p}{dx^2}\)) 计算。
- 将其代回原微分方程。
- 通过比较等式两边的系数或项来解出这些常数。
2.3 关键情况：重合（共振）
这是最常见的陷阱！
避免常见的错误：如果你的试探函数 \(y_p\) 已经包含在补余函数 \(y_c\) 中（即它与 \(y_c\) 中的某一项重复），代入后将无法解出常数。
解决方法：如果发生重复，将试探函数 \(y_p\) 乘以 \(x\)。如果新的试探函数仍然与 \(y_c\) 中的项重复，则乘以 \(x^2\)。
示例场景：
- CF: \(y_c = Ae^{2x} + Be^{-x}\)。
- \(f(x)\) 为 \(4e^{2x}\)。
- 最初尝试的 PI: \(P e^{2x}\)。这与 \(Ae^{2x}\) 重复。
- 新的试探 PI: \(y_p = Px e^{2x}\)。
教学大纲关联：大纲明确提到了当题目给出如 \(kx \cos 2x\) 这种特殊形式作为 PI 时，需要求出系数 \(k\)。这种情况正好发生在标准的 PI (\(P \cos 2x + Q \sin 2x\)) 因与 CF 中的项重复而失效时（此时辅助方程的根为 \(\pm 2i\)）。
关键点 (PI)：根据 \(f(x)\) 选择试探函数。一定要检查它是否与 CF 中的任何项重复。如果重复，请乘以 \(x\)。
3. 使用代换法简化微分方程
有时，微分方程在给定的形式下既不是线性的，也不是可分离变量的。教学大纲要求使用特定的代换法，将复杂的方程转化为已知如何求解的形式（例如一阶线性微分方程或二阶常系数线性微分方程）。
3.1 代换法以化为常系数方程 (\(x = e^t\))
这种代换用于包含 \(x^2\frac{d^2y}{dx^2}\) 和 \(x\frac{dy}{dx}\) 等项的方程。这通常被称为欧拉-柯西方程 (Euler-Cauchy equation)。
令 \(x = e^t\)，则 \(t = \ln x\)。利用链式法则进行导数转换：
一阶导数： $$ \frac{dy}{dx} = \frac{dy}{dt} \cdot \frac{dt}{dx} = \frac{dy}{dt} \cdot \frac{1}{x} $$ $$ x\frac{dy}{dx} = \frac{dy}{dt} $$
二阶导数： $$ \frac{d^2y}{dx^2} = \frac{d}{dx} \left( \frac{1}{x}\frac{dy}{dt} \right) $$ （在右侧使用乘法法则） $$ x^2\frac{d^2y}{dx^2} = \frac{d^2y}{dt^2} - \frac{dy}{dt} $$
将这些表达式代回，即可将原方程变为关于 \(t\) 的常系数微分方程，你可以使用第 2 节中的 CF/PI 方法求解。
3.2 代换法以化为可分离变量形式 (\(y = ux\))
此代换通常用于形式为 \(\frac{dy}{dx} = F\left(\frac{y}{x}\right)\) 的一阶齐次微分方程。
教学大纲示例形式：将 \(\frac{dy}{dx} = \frac{x+y}{x-y}\) 化为可分离变量形式。
令 \(y = ux\)。我们需要利用乘法法则找出 \(u, x\) 和 \(\frac{du}{dx}\) 表示的 \(\frac{dy}{dx}\)：
$$ \frac{dy}{dx} = u(1) + x\frac{du}{dx} $$ $$ \frac{dy}{dx} = u + x\frac{du}{dx} $$
将 \(y=ux\) 及新的导数代入原方程：
$$ u + x\frac{du}{dx} = \frac{x + ux}{x - ux} = \frac{x(1 + u)}{x(1 - u)} = \frac{1 + u}{1 - u} $$
现在，重排方程以分离变量 \(u\) 和 \(x\)：
$$ x\frac{du}{dx} = \frac{1 + u}{1 - u} - u = \frac{1 + u - u(1 - u)}{1 - u} = \frac{1 + u^2}{1 - u} $$
$$ \frac{1 - u}{1 + u^2} du = \frac{1}{x} dx $$
这个新方程现在是可分离变量的，可以通过直接积分求解，最后将 \(u\) 换回 \(\frac{y}{x}\) 即可。
关键点 (代换法)：代换法使我们能够通过将非标准微分方程转换为已知的线性标准形式，从而解决难题。
4. 初始条件与解释
你求出的通解 (GS) 包含任意常数（一阶方程为 \(C\)，二阶为 \(A\) 和 \(B\)）。这代表了一个无限的解族。
4.1 寻找特解 (Particular Solution)
要找到唯一的解——即特解，必须使用题目给出的初始条件 (ICs)。
对于一阶微分方程，你需要一个条件（例如 \(y(0)=5\)）。
对于二阶微分方程，你需要两个条件（例如 \(y(0)=2\) 和 \(y'(0)=1\)）。
应用步骤：
- 求出通解 \(y\)。
- 如有必要（对于二阶方程），对 \(y\) 求导以得到 \(\frac{dy}{dx}\) 或 \(y'\)。
- 将第一个初始条件代入 \(y\)，求出一个常数（或常数间的关系）。
- 如果需要，将第二个初始条件代入 \(y'\)，求出剩余的常数。
- 将这些值代回通解，得到特解。
4.2 环境背景下的解释
微分方程常用于建模动态系统（运动、人口、热传导）。教学大纲要求你根据所建模的问题对解进行解释。
- 补余函数 (CF)： 代表系统的自然响应或瞬态行为。如果根为负，这些项会随时间衰减（例如，逐渐停止的震荡）。
- 特解 (PI)： 代表系统的强制响应或稳态行为。这是系统长期的表现，通常由外部输入 \(f(x)\) 决定。
- 特解 (Particular Solution)： 一旦确定了常数，这就描述了系统从给定初始状态开始的特定行为。
你知道吗？
如果你正在求解一个弹簧振子系统的微分方程，辅助方程中的复根意味着系统会发生震荡（类似正弦波），而实根则意味着它在不震荡的情况下回到平衡状态（即“过阻尼”）。数学能让你精确预见物理系统的表现！
关键点 (ICs 和解释)：初始条件固定了常数，从而定义了一条唯一的路径。始终将数学成分（CF, PI）与问题中数量的物理意义（如时间、位移、电荷）联系起来。