微分方程 – 二阶 (FP2.11)
你好,未来的数学家!欢迎来到二阶微分方程(DEs)的世界。别担心,尽管这个名字听起来很深奥,但它们仅仅是用来描述系统如何基于加速度发生变化的数学工具。这在模拟物理现象(如振动、电路,甚至经济模型)时非常有用。
在本章中,我们将学习如何在已知包含 \(y\)、\(\frac{dy}{dx}\) 和 \(\frac{d^2y}{dx^2}\) 的方程的情况下,求出 \(y\) 的完整精确公式。
快速了解:为什么叫“二阶”?
微分方程的“阶”由其中出现的最高导数决定。如果方程中包含 \(\frac{d^2y}{dx^2}\)(即加速度),它就是二阶的。这意味着我们期望通解中有两个任意常数(通常为 \(A\) 和 \(B\)),这就像在力学中,你需要初始位移和初始速度来完整描述一个抛体运动的轨迹一样!
1. 线性二阶微分方程的结构
我们主要关注常系数线性微分方程。这意味着导数项和 \(y\) 本身都仅以一次幂出现,且它们的系数(\(a, b, c\))均为常数(根据考纲要求,这些通常为整数)。
两大主要类型
根据等式右侧的函数 \(f(x)\),二阶线性微分方程分为两类:
-
1. 齐次方程 (Homogeneous Equations)
这是 \(f(x)=0\) 的“基本情况”。$$ a\frac{d^2y}{dx^2} + b\frac{dy}{dx} + cy = 0 $$
类比:这描述了系统的自然振荡,就像真空中的弹簧振动(没有外力作用)。 -
2. 非齐次方程 (Non-Homogeneous Equations)
这是当 \(f(x)\) 为非零函数的情况。$$ a\frac{d^2y}{dx^2} + b\frac{dy}{dx} + cy = f(x) $$
类比:这描述了在有外力(如引擎运转或持续推力)作用下的系统。
核心要点: 解非齐次方程需要两部分:齐次方程的解(称为互补函数,Complementary Function)和对应外力项的额外解(称为特解,Particular Integral)。
2. 求解齐次方程:互补函数 (CF)
齐次方程 \( a\frac{d^2y}{dx^2} + b\frac{dy}{dx} + cy = 0 \) 的解被称为互补函数(\(y_{CF}\))。
求解它们的秘诀在于辅助方程 (AE)。我们假设解的形式为 \(y = e^{\lambda x}\),代入后,微分方程就简化为关于 \(\lambda\) 的二次方程:
辅助方程 (AE)
对于 \( a\frac{d^2y}{dx^2} + b\frac{dy}{dx} + cy = 0 \),其辅助方程为:
$$ a\lambda^2 + b\lambda + c = 0 $$
CF 的形式完全取决于该二次方程根的性质。
情况 1:两个不同的实根 (\(\lambda_1 \neq \lambda_2\))
当判别式 \((b^2 - 4ac)\) 大于零时出现。
通解(CF)是两个指数解的线性组合: $$ y_{CF} = Ae^{\lambda_1 x} + Be^{\lambda_2 x} $$
类比:这通常代表“过阻尼”运动,系统平稳回到平衡位置而不产生振荡。
情况 2:一个重实根 (\(\lambda_1 = \lambda_2 = \lambda\))
当判别式 \((b^2 - 4ac)\) 等于零时出现。
通解(CF)必须包含一个额外的 \(x\) 因子,以确保两项是线性无关的: $$ y_{CF} = (A + Bx)e^{\lambda x} $$
常见错误提醒:千万别漏掉第二项的 \(x\)!如果你写成 \(Ae^{\lambda x} + Be^{\lambda x}\),它会合并为 \((A+B)e^{\lambda x}\),这样就只有一个任意常数,而不是两个。
情况 3:共轭复根 (\(\lambda = \alpha \pm i\beta\))
当判别式 \((b^2 - 4ac)\) 小于零时出现。由于 \(a, b, c\) 是实整数,根必然成对出现。
虽然原始的指数形式为 \(y = A'e^{(\alpha + i\beta)x} + B'e^{(\alpha - i\beta)x}\),但这通常使用欧拉公式 (\(e^{i\theta} = \cos\theta + i\sin\theta\)) 转化为更简洁的三角函数形式: $$ y_{CF} = e^{\alpha x} (A \cos(\beta x) + B \sin(\beta x)) $$
类比:这代表振荡运动。\(\alpha\) 决定了振荡是增强(\(\alpha > 0\))、衰减(\(\alpha < 0\),即阻尼),还是保持恒定(\(\alpha = 0\),即简谐运动)。
✅ 快速回顾:CF 公式
设 \(\lambda_1, \lambda_2\) 为 \( a\lambda^2 + b\lambda + c = 0 \) 的根。
- 不同实根: \( y_{CF} = Ae^{\lambda_1 x} + Be^{\lambda_2 x} \)
- 重实根 (\(\lambda\)): \( y_{CF} = (A + Bx)e^{\lambda x} \)
- 复根 (\(\alpha \pm i\beta\)): \( y_{CF} = e^{\alpha x} (A \cos(\beta x) + B \sin(\beta x)) \)
3. 求解非齐次方程:特解 (PI)
在求解 \( a\frac{d^2y}{dx^2} + b\frac{dy}{dx} + cy = f(x) \) 时,通解 (\(y_{GS}\)) 是互补函数 (CF) 与特解 (PI) 之和:
$$ y_{GS} = y_{CF} + y_{PI} $$
CF 处理的是自然行为(如第 2 部分所述)。PI 是代入原方程后能满足等式的任意函数。由于本考纲中 \(f(x)\) 形式有限,寻找 \(y_{PI}\) 涉及到合理的猜测。
寻找 \(y_{PI}\) 的分步指南
PI 必须与 \(f(x)\) 具有相同的通用形式。
第一步:根据 \(f(x)\) 猜测 \(y_{PI}\) 的形式。
| \( f(x) \) (外力项) | \( y_{PI} \) 的猜测 | 最高次数限制 |
|---|---|---|
| 多项式 (如 \(x^3\)) | 同次数的一般多项式 (如 \(Ax^3 + Bx^2 + Cx + D\)) | 4 次 |
| 指数函数 (如 \(e^{kx}\)) | \( Pe^{kx} \) | 不适用 |
| 三角函数 (如 \(\sin(kx)\) 或 \(\cos(kx)\)) | \( P \cos(kx) + Q \sin(kx) \) (必须同时包含 sin 和 cos!) | 不适用 |
第二步:求导并代入。
对你的猜测 (\(y_{PI}\)) 求两次导,得到 \(\frac{dy_{PI}}{dx}\) 和 \(\frac{d^2y_{PI}}{dx^2}\)。将这三项代入原始的非齐次微分方程。
第三步:待定系数法。
比较等式两边 \(e^{kx}\)、\(\cos(kx)\)、\(\sin(kx)\) 或 \(x\) 的幂次的系数。这将得到一组联立方程,解出 PI 猜测中的未知常数(\(A, B, P, Q\) 等)。
关键重叠规则(修正规则)
这是最容易出错的地方!
如果你选择的 \(y_{PI}\) 猜测(来自第一步)已经是互补函数 (\(y_{CF}\)) 中的一项,那么标准的猜测将失效(代入齐次部分时会直接变为零)。
如果出现重叠,你必须将标准猜测乘以最小可能的 \(x\) 的幂次(通常是 \(x\) 或 \(x^2\)),直到它与 CF 项线性无关。
例子:如果 \(f(x) = 5e^{3x}\),标准猜测是 \(y_{PI} = Pe^{3x}\)。但如果 CF 为 \(y_{CF} = Ae^{3x} + Be^{-x}\),则 PI 猜测与 \(Ae^{3x}\) 重叠了。
此时修正后的 PI 猜测应为:\( y_{PI} = Px e^{3x} \)。
你知道吗? 这种乘以 \(x\) 的操作是必要的,因为强制输入的频率(来自 \(f(x)\))与自然频率(来自 CF)相匹配,从而引发了称为共振的现象。系统的响应随时间线性增长,这就是 \(x\) 项的由来!
核心要点: 一定要先求出 \(y_{CF}\),然后将其与猜测的 \(y_{PI}\) 形式进行对比。如果两者匹配,将 \(y_{PI}\) 乘以 \(x\)。
4. 求特解(使用边界条件)
一旦你得到了通解 \(y_{GS} = y_{CF} + y_{PI}\),它仍然包含来自 CF 的两个任意常数 \(A\) 和 \(B\)。为了找到唯一的特解,你需要两个额外的信息,即边界条件或初始条件。
这些条件通常指定了在特定点 \(x\)(通常在 \(x=0\) 时)的 \(y\) 或 \(\frac{dy}{dx}\) 的值。
求特解的分步步骤
- 求 \(y_{GS}\): 计算 \(y_{CF}\) 和 \(y_{PI}\) 并合并:\(y = y_{CF} + y_{PI}\)。
- 对 \(y_{GS}\) 求导: 计算 \(\frac{dy}{dx}\)。(如果条件涉及导数,你将需要用到它)。
- 应用条件 1: 将第一个边界条件(例如 \(y(0)=5\))代入 \(y_{GS}\)。这会建立一个关于 \(A\) 和 \(B\) 的方程。
- 应用条件 2: 将第二个边界条件(例如 \(y'(0)=1\))代入 \(\frac{dy}{dx}\)。这会建立第二个关于 \(A\) 和 \(B\) 的方程。
- 求解: 解这个二元一次方程组,求出 \(A\) 和 \(B\)。
- 最终答案: 将 \(A\) 和 \(B\) 的值代回 \(y_{GS}\),写出最终的特解。
给同学的提示:求导时一定要反复检查计算!这里的错误不可避免地会导致 A 和 B 的值出错。请务必细心,特别是在处理复杂的 CF 或 PI 项而需要使用乘积法则时。
📝 最终核心总结
1. 识别 DE 类型:齐次 (\(=0\)) 或非齐次 (\(=f(x)\))。
2. 通解结构为:\( y_{GS} = y_{CF} + y_{PI} \)。
3. CF (互补函数): 使用辅助方程 \( a\lambda^2 + b\lambda + c = 0 \) 求得。形式取决于根是不同实根、重实根还是复根。
4. PI (特解): 通过猜测与 \(f(x)\) 相似的形式求得。关键点:检查是否与 CF 重叠,如有必要,乘以 \(x\)。
5. 特解 (Particular Solution): 使用给定的两个边界条件来确定常数 \(A\) 和 \(B\) 的具体数值。