📚 FP2 学习笔记:二阶微分方程
你好,未来的数学家!这一章我们要学习的二阶微分方程(Second Order Differential Equations, 简称 S.O.D.E.s)是进阶纯数(Further Pure Mathematics)中非常精彩的内容。别担心名字听起来很深奥——本质上,我们是在学习一套描述运动、振动和电路的数学语言。
掌握了 S.O.D.E.s,我们就能对那些“变化率既取决于当前值,又取决于变化率本身”的系统进行建模。想象一下弹簧、钟摆,甚至是谣言传播的速度!
本章学习目标:
- 求解具有常系数的二阶微分方程。
- 理解齐次方程的通解部分——互补函数 (Complementary Function, C.F.)。
- 求非齐次方程的特解——特解 (Particular Integral, P.I.)。
- 利用初始条件和边界条件确定唯一的特解。
1. 二阶微分方程简介
1.1. 基本形式与术语
二阶微分方程包含二阶导数,即 \(\frac{\text{d}^2y}{\text{d}x^2}\)。在 FP2 中,我们主要关注常系数线性二阶微分方程。
我们要研究的通式为:
$$ a \frac{\text{d}^2y}{\text{d}x^2} + b \frac{\text{d}y}{\text{d}x} + cy = f(x) $$
- \(a, b, c\):这些是固定的常数(数字),这让我们的计算轻松不少!
- \(\frac{\text{d}^2y}{\text{d}x^2}\):二阶导数(代表加速度或曲率)。
- \(f(x)\):这是驱动函数或强制项。
我们将这些方程分为两类:
- 齐次方程 (Homogeneous Equation):当 \(f(x) = 0\) 时。它描述的是系统自然、无外力干扰的行为。
- 非齐次方程 (Non-Homogeneous Equation):当 \(f(x) \neq 0\) 时。它描述的是系统在外部影响下的状态。
目标:求通解 (General Solution, G.S.)
最终解 \(y\) 总是由两部分组成的:
$$ \mathbf{y} = \mathbf{y_{CF}} + \mathbf{y_{PI}} $$
$$ \text{通解 (G.S.)} = \text{互补函数 (C.F.)} + \text{特解 (P.I.)} $$
类比:想象汽车的运动。互补函数 (C.F.) 是汽车自然运动的方式(例如,熄火后靠惯性滑行)。特解 (P.I.) 是由外力引起的运动(例如,一直踩着油门)。
2. 第一步:求解齐次方程(互补函数,C.F.)
要找到 C.F. (\(y_{CF}\)),我们必须解该方程:
$$ a \frac{\text{d}^2y}{\text{d}x^2} + b \frac{\text{d}y}{\text{d}x} + cy = 0 $$
2.1. 辅助方程 (Auxiliary Equation)
我们假设解的形式为 \(y = e^{mx}\)。将其代入齐次方程,可以直接得到一个关于 \(m\) 的简单二次方程,即辅助方程(或特征方程):
$$ \mathbf{am^2 + bm + c = 0} $$
解这个二次方程得到根 (\(m\)),这些根决定了 C.F. 的形式。由于是二次方程,根有三种可能性,对应三种不同的解法形式。
2.2. 情况 1:两个不同的实根 (\(m_1 \neq m_2\))
如果判别式 \((b^2 - 4ac)\) 大于零,你会得到两个不同的实数 \(m_1\) 和 \(m_2\)。
示例:\(m^2 + 5m + 6 = 0 \implies (m+2)(m+3) = 0\)。根为 \(m_1 = -2\), \(m_2 = -3\)。
C.F. 形式: $$ \mathbf{y_{CF} = Ae^{m_1x} + Be^{m_2x}} $$
注:\(A\) 和 \(B\) 是任意常数,我们稍后会通过初始条件求出它们。
2.3. 情况 2:一个重实根 (\(m_1 = m_2 = m\))
如果判别式 \((b^2 - 4ac)\) 等于零,你只会得到一个实根 \(m\)。这意味着辅助方程是一个完全平方式。
示例:\(m^2 - 4m + 4 = 0 \implies (m-2)^2 = 0\)。根为 \(m = 2\)。
C.F. 形式: $$ \mathbf{y_{CF} = (Ax + B)e^{mx}} $$
💡 记忆小贴士:当出现重根时,你需要添加额外的 \(Ax\) 项,以确保解是真正的“二阶”解(即包含两个独立的常数 \(A\) 和 \(B\))。
2.4. 情况 3:共轭复根 (\(\alpha \pm i\beta\))
如果判别式 \((b^2 - 4ac)\) 小于零,你将得到两个复根,它们总是以共轭对的形式出现:\(\alpha + i\beta\) 和 \(\alpha - i\beta\)。
示例:\(m^2 + 2m + 5 = 0\)。使用求根公式可得 \(m = -1 \pm 2i\)。此处 \(\alpha = -1\), \(\beta = 2\)。
C.F. 形式: $$ \mathbf{y_{CF} = e^{\alpha x}(A \cos(\beta x) + B \sin(\beta x))} $$
你知道吗?这种形式描述的是振动(如钟摆摆动),可以是阻尼振动 (\(\alpha < 0\)),也可以是发散振动 (\(\alpha > 0\))。
3. 第二步:求特解 (P.I.)
特解 (\(y_{PI}\)) 是满足非齐次方程的一个特定解:
$$ a \frac{\text{d}^2y}{\text{d}x^2} + b \frac{\text{d}y}{\text{d}x} + cy = f(x) $$
这里的关键是根据函数 \(f(x)\) 猜测 \(y_{PI}\) 的形式。
3.1. 猜测策略
我们选择一个与 \(f(x)\) 形式类似的试探解,并包含未知常数系数(例如 \(P, Q, R\))。
试探解表 (P.I. 形式)
设 \(f(x)\) 为等式右侧的项:
| 如果 \(f(x)\) 是... | 猜测 \(y_{PI}\) 为... | 示例 |
|---|---|---|
| \(n\) 次多项式(例如常数、一次、二次) | \(n\) 次通项多项式 | 如果 \(f(x) = 5x^2 - 3\),猜测 \(y_{PI} = Px^2 + Qx + R\)。 |
| 指数函数 (\(Ke^{kx}\)) | 指数函数的常数倍 | 如果 \(f(x) = 7e^{-x}\),猜测 \(y_{PI} = Pe^{-x}\)。 |
| 三角函数 (\(K \cos(\omega x)\) 或 \(K \sin(\omega x)\)) | 正弦和余弦的组合 | 如果 \(f(x) = 4 \sin(2x)\),猜测 \(y_{PI} = P \cos(2x) + Q \sin(2x)\)。 |
3.2. 求 P.I. 的步骤
- 猜测:根据 \(f(x)\) 选择正确的试探形式 \(y_{PI}\)。
- 求导:计算试探解的一阶导数 (\(y'_{PI}\)) 和二阶导数 (\(y''_{PI}\))。
- 代入:将 \(y_{PI}\), \(y'_{PI}\), 和 \(y''_{PI}\) 代入原二阶微分方程。
- 比较系数:令等式两侧对应项的系数相等(例如:先令 \(x^2\) 的系数相等,再令 \(x\) 的系数相等,最后令常数项相等),从而解出 \(P, Q, R, \dots\)。
如果计算出的方程很长也不要慌;慢慢来,整理好你的项!
3.3. 特殊情况:重叠(共振规则)
这是最常见的陷阱!当你的 P.I. 试探解中包含的某一项已经在互补函数 (C.F.) 中出现时,就会发生重叠。
如果在这种情况下不调整 P.I. 的猜测,代入步骤会导致 \(0 = f(x)\),这是不可能的!
重叠规则:
如果标准的 P.I. 猜测与 C.F. 中的项重叠,你必须将整个试探解乘以 \(\mathbf{x}\)。
重叠示例:
假设 \(y_{CF} = Ae^{2x} + Be^{-3x}\)。
如果 \(f(x) = 5e^{2x}\),最初的猜测是 \(y_{PI} = Pe^{2x}\)。
由于 \(e^{2x}\) 已经在 C.F. 中存在,该猜测失效。
正确的调整后猜测应为: \(y_{PI} = Px e^{2x}\)。
复杂重叠(重根情况):
如果 C.F. 含有重根,即 \(y_{CF} = (Ax + B)e^{2x}\),且 \(f(x) = 5e^{2x}\),此时乘以一次 \(x\) (\(Pxe^{2x}\)) 仍然会产生重叠(因为 \(Axe^{2x}\) 已经在 C.F. 中)。在这种特殊情况下,你必须乘以 \(\mathbf{x^2}\)。
正确调整后的猜测: \(y_{PI} = Px^2 e^{2x}\)。
4. 第三步 & 第四步:通解与特解
4.1. 构成通解 (G.S.)
一旦你成功求出 C.F.(包含常数 \(A\) 和 \(B\))以及 P.I.(确定了常数 \(P, Q, \dots\)),就可以将它们组合起来:
$$ y = y_{CF} + y_{PI} $$
这个 G.S. 包含两个任意常数 \(A\) 和 \(B\),它代表了一整个解族。
4.2. 求特解 (Particular Solution)
在应用问题中,我们通常需要一个特定的唯一解。我们会利用边界条件或初始条件来找出 \(A\) 和 \(B\) 的准确值。
由于方程是二阶的,你需要两个独立的条件:
- 初始条件:两个条件均在同一起始点给出(例如 \(x=0\))。
示例: \(y(0) = 5\)(起始位置)且 \(y'(0) = 1\)(起始速度)。 - 边界条件:在两个不同的点给出条件。
示例: \(y(0) = 5\) 且 \(y(1) = 10\)。
求 \(A\) 和 \(B\) 的步骤
- 写出 G.S.: \(y = y_{CF} + y_{PI}\)。
- 对 G.S. 求导: 算出 \(y' = y'_{CF} + y'_{PI}\)。(如果初始条件涉及 \(y'\),则需要这一步)。
- 代入条件: 将已知值(如 \(x=0, y=5\))代入 \(y\) 和 \(y'\) 的方程。
- 解联立方程: 现在你会得到一对关于 \(A\) 和 \(B\) 的线性联立方程。解出它们的具体值。
- 写出特解: 将 \(A\) 和 \(B\) 的值代回通解公式中。
提示:尽量利用 \(x=0\) 时的初始条件,因为 \(e^0 = 1\), \(\sin(0)=0\), 且 \(\cos(0)=1\),这会极大地简化联立方程!
1. 齐次 \(\rightarrow\) 辅助方程 \(\rightarrow\) C.F.(得到 A, B)。
2. 非齐次 \(\rightarrow\) 猜测 P.I.(若与 C.F. 重叠则调整) \(\rightarrow\) 代入并解出常数 P, Q, R。
3. G.S. = C.F. + P.I。
4. 代入初始/边界条件求出具体的 A 和 B。