微分方程 – 二阶 (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}\) 的值。

求特解的分步步骤

  1. 求 \(y_{GS}\): 计算 \(y_{CF}\) 和 \(y_{PI}\) 并合并:\(y = y_{CF} + y_{PI}\)。
  2. 对 \(y_{GS}\) 求导: 计算 \(\frac{dy}{dx}\)。(如果条件涉及导数,你将需要用到它)。
  3. 应用条件 1: 将第一个边界条件(例如 \(y(0)=5\))代入 \(y_{GS}\)。这会建立一个关于 \(A\) 和 \(B\) 的方程。
  4. 应用条件 2: 将第二个边界条件(例如 \(y'(0)=1\))代入 \(\frac{dy}{dx}\)。这会建立第二个关于 \(A\) 和 \(B\) 的方程。
  5. 求解: 解这个二元一次方程组,求出 \(A\) 和 \(B\)。
  6. 最终答案: 将 \(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\) 的具体数值。