👋 欢迎来到二阶微分方程的世界！
嘿，未来的数学家！欢迎来到进阶纯数学中最强大的章节之一：二阶微分方程 (Second Order Differential Equations, 简称 S.O.D.E.s)。
别担心，虽然名字听起来很唬人，但简单来说，这些方程能帮我们模拟现实世界中涉及加速度的现象，比如弹簧的振动、电子电路的稳定过程，甚至是钟摆的运动。
本章的目标是寻找一个满足包含 \(y\)、\(\frac{dy}{dx}\) 和 \(\frac{d^2y}{dx^2}\) 方程的函数 \(y\)。学完这些笔记后，你将掌握如何一步步解出它们！
第一节：二阶线性方程的结构
1.1 标准形式与常系数
在 FP2 中，我们只处理具有常系数的线性二阶微分方程。这意味着导数项前的系数都是常数，而不是 \(x\) 的函数。
其一般形式如下：
\[ a\frac{d^2y}{dx^2} + b\frac{dy}{dx} + cy = f(x) \]
- \(a, b, c\) 是常数（系数）。
- \(f(x)\) 是驱动系统的函数（如果是齐次方程，则为 0）。
🔑 核心概念：通解
S.O.D.E. 的最终解始终由两部分组成：
\[ y = y_{CF} + y_{PI} \]
- \(y_{CF}\)：互补函数 (Complementary Function, CF)。描述系统在忽略驱动函数 \(f(x)\) 时的自然行为。
- \(y_{PI}\)：特解 (Particular Integral, PI)。描述由驱动函数 \(f(x)\) 引起的特定响应。
让我们先从求解 CF 开始。
第二节：求互补函数 (\(y_{CF}\))
要求出 CF，我们必须先求解齐次方程 (Homogeneous Equation)，即令等式右侧为零：
\[ a\frac{d^2y}{dx^2} + b\frac{dy}{dx} + cy = 0 \]
2.1 辅助方程 (Auxiliary Equation, AE)
我们假设解具有指数形式 \(y = e^{mx}\)。将其代入齐次方程后，我们得到一个简单的二次方程，称为辅助方程 (AE)：
\[ am^2 + bm + c = 0 \]
该方程根 \(m\) 的性质决定了 CF 的形式。根据判别式 \(b^2 - 4ac\)，共有三种情况。
2.2 情况 1：实数且不相等的根 (\(m_1 \neq m_2\))
如果 AE 有两个不同的实根 \(m_1\) 和 \(m_2\)。
CF 的解：
\[ y_{CF} = A e^{m_1 x} + B e^{m_2 x} \]
(A 和 B 是由初始条件决定的任意常数。)
例题：解 \(\frac{d^2y}{dx^2} + 5\frac{dy}{dx} + 6y = 0\)。
AE: \(m^2 + 5m + 6 = 0 \Rightarrow (m+2)(m+3) = 0\)。根为 \(m_1 = -2\) 和 \(m_2 = -3\)。
CF: \(y_{CF} = A e^{-2x} + B e^{-3x}\)。
2.3 情况 2：实数且相等的根 (\(m_1 = m_2 = m\))
如果 AE 只有一个实根（二次方程为完全平方式）。
CF 的解：
\[ y_{CF} = (A x + B) e^{m x} \]
💡 为什么要乘 \(x\)？ 如果只写成 \(Ae^{mx} + Be^{mx}\)，我们就能将常数合并成一项，这无法构成二阶方程的通解（我们需要两个任意常数）。将其中一项乘以 \(x\) 确保了我们拥有两个真正独立的解。
例题：解 \(\frac{d^2y}{dx^2} - 4\frac{dy}{dx} + 4y = 0\)。
AE: \(m^2 - 4m + 4 = 0 \Rightarrow (m-2)^2 = 0\)。根为 \(m = 2\)。
CF: \(y_{CF} = (A x + B) e^{2x}\)。
2.4 情况 3：复共轭根 (\(m = \alpha \pm i\beta\))
如果判别式小于零，根为复共轭形式，其中 \(\alpha\) 为实部，\(\beta\) 为虚部。
CF 的解：
\[ y_{CF} = e^{\alpha x} (A \cos(\beta x) + B \sin(\beta x)) \]
你知道吗？ 这种解结构源自欧拉公式 (Euler's Identity)，它连接了指数函数与三角函数。这种形式完美地模拟了振荡系统（如摆动的钟摆或弹簧）。
例题：解 \(\frac{d^2y}{dx^2} + 2\frac{dy}{dx} + 5y = 0\)。
AE: \(m^2 + 2m + 5 = 0\)。使用求根公式：
\(m = \frac{-2 \pm \sqrt{4 - 4(5)}}{2} = \frac{-2 \pm \sqrt{-16}}{2} = \frac{-2 \pm 4i}{2} = -1 \pm 2i\)。
此处 \(\alpha = -1\)，\(\beta = 2\)。
CF: \(y_{CF} = e^{-x} (A \cos(2x) + B \sin(2x))\)。
快速回顾：CF 形式
- 不同实根 \(m_1, m_2\)： \(A e^{m_1 x} + B e^{m_2 x}\)
- 相等实根 \(m\)： \((A x + B) e^{m x}\)
- 复根 \(\alpha \pm i\beta\)： \(e^{\alpha x} (A \cos(\beta x) + B \sin(\beta x))\)
第三节：求特解 (\(y_{PI}\))
特解 (PI) 用于处理非齐次方程右侧的 \(f(x)\)：
\[ a\frac{d^2y}{dx^2} + b\frac{dy}{dx} + cy = f(x) \]
该方法包括根据 \(f(x)\) 的形式“猜测”一个 PI 的函数形式，然后将其代回微分方程以求解待定系数。
3.1 猜测 \(y_{PI}\) 的形式
猜测形式必须与 \(f(x)\) 具有相同的结构，但需包含待定系数（如 \(P, Q, R\)）。
情况 A：多项式 \(f(x)\)
如果 \(f(x)\) 是 \(n\) 次多项式，则猜测一个 \(n\) 次的一般多项式。
| 如果 \(f(x)\) 是： | 猜测 \(y_{PI}\) 为： |
|---|---|
| \(5\) (常数) | \(P\) |
| \(2x - 1\) (1次) | \(Px + Q\) |
| \(x^2 + 3\) (2次) | \(P x^2 + Q x + R\) |
情况 B：指数函数 \(f(x)\)
如果 \(f(x)\) 是 \(k e^{kx}\)。
| 如果 \(f(x)\) 是： | 猜测 \(y_{PI}\) 为： |
|---|---|
| \(4 e^{5x}\) | \(P e^{5x}\) |
情况 C：三角函数 \(f(x)\)
如果 \(f(x)\) 涉及正弦或余弦，猜测项必须同时包含正弦和余弦项（因为对正弦求导会得到余弦，反之亦然）。
| 如果 \(f(x)\) 是： | 猜测 \(y_{PI}\) 为： |
|---|---|
| \(7 \sin(3x)\) | \(P \cos(3x) + Q \sin(3x)\) |
| \(\cos(x)\) | \(P \cos(x) + Q \sin(x)\) |
PI 求解步骤：
- 构建猜测形式 \(y_{PI}\)。
- 求出 \(\frac{dy_{PI}}{dx}\) 和 \(\frac{d^2y_{PI}}{dx^2}\)。
- 将 \(y_{PI}\)、\(\frac{dy_{PI}}{dx}\) 和 \(\frac{d^2y_{PI}}{dx^2}\) 代回原方程。
- 对比等式左侧与右侧 \(f(x)\) 的系数。
- 解出 \(P, Q, R, ...\) 的方程组。
如果起初觉得有些麻烦也不要担心，代入和系数对比是这一步最难的地方，需要细心的代数运算！
3.2 关键法则：共振条件（重合）
这是学生最常犯的错误！如果你猜测的 \(y_{PI}\) 已经是 \(y_{CF}\) 的一部分，意味着该项代入齐次方程的结果为零，代入它将无法平衡 \(f(x)\) 部分。
类比：想象推秋千。如果你以秋千的自然频率（CF）去推，秋千摆动会变得巨大，但驱动函数 \(f(x)\) 并不会像预想的那样工作。为了修正，你需要调整解的结构。
修正方法：如果猜测的 \(y_{PI}\) 项已经存在于 \(y_{CF}\) 中，必须将整个猜测乘以 \(x\)。
重合示例（情况 B - 指数）：
假设 AE 的根为 \(m=3\) 和 \(m=-1\)。\(y_{CF} = A e^{3x} + B e^{-x}\)。
如果 \(f(x) = 5e^{3x}\)，初始猜测应为 \(y_{PI} = P e^{3x}\)。
因为 \(e^{3x}\) 已经在 CF 中了，所以该猜测无效。
修正后的猜测： \(y_{PI} = P x e^{3x}\)。
如果乘以 \(x\) 后依然重合怎么办？
这种情况只发生在重根（情况 2）下。如果 AE 有重根 \(m\)，CF 为 \((Ax+B)e^{mx}\)。如果 \(f(x) = k e^{mx}\)，则 \(e^{mx}\) 和 \(x e^{mx}\) 都是齐次方程的解。在这种极罕见的情况下，必须乘以 \(x^2\)。
修正后的猜测（双重重合）： \(y_{PI} = P x^2 e^{mx}\)。
🔥 PI 的关键点
一定要先求 \(y_{CF}\)，再构建 \(y_{PI}\)。检查是否重合。若有重合，将猜测乘以 \(x\)。
第四节：通解与边界条件
4.1 形成通解
当你求出包含任意常数 \(A\) 和 \(B\) 的互补函数 \(y_{CF}\) 以及求出待定系数 \(P, Q, ...\) 的特解 \(y_{PI}\) 后，将它们相加即可得到通解：
\[ y = y_{CF} + y_{PI} \]
这个通解包含两个任意常数 \(A\) 和 \(B\)。
4.2 利用初始条件和边界条件
为了求出特定物理问题的唯一解，我们必须利用边界条件或初始条件来确定 \(A\) 和 \(B\) 的值。
由于我们处理的是二阶方程，始终需要两个条件。它们通常涉及 \(y\) 在某一点的值（如 \(y(0)=5\)）和/或导数 \(\frac{dy}{dx}\) 在某一点的值（如 \(y'(0)=1\)）。
求最终解的步骤：
- 写出通解 \(y = y_{CF} + y_{PI}\)。
- 对通解求导得到 \(\frac{dy}{dx}\)。
- 应用第一个条件（例如，将 \(x=0, y=5\) 代入 \(y\) 方程）得到关于 \(A\) 和 \(B\) 的方程。
- 应用第二个条件（例如，将 \(x=0, y'=1\) 代入 \(\frac{dy}{dx}\) 方程）得到第二个方程。
- 联立方程组求解 \(A\) 和 \(B\)。
- 将 \(A\) 和 \(B\) 代回通解，得到唯一解。
常见错误警告！
应用涉及 \(\frac{dy}{dx}\) 的条件时，学生经常忘记包含 \(y_{PI}\) 的导数部分。记住：在代入初始条件计算斜率之前，务必对整个通解求导！
第五节：总结与回顾
解二阶微分方程是一项流程化工作。只要按照步骤操作，你一定能成功！
S.O.D.E.s 求解流程图
1. 确定齐次部分： \(a\frac{d^2y}{dx^2} + b\frac{dy}{dx} + cy = 0\)
2. 解辅助方程： 从 \(am^2 + bm + c = 0\) 求出根 \(m\)。
3. 确定 \(y_{CF}\)： 根据根的性质（不同实根、重根、复根）使用正确的公式。
4. 确定 \(y_{PI}\)： 观察 \(f(x)\) 并写出合适的猜测形式（多项式、三角函数、指数函数）。
5. 检查重合： 如果 \(y_{PI}\) 的猜测与 \(y_{CF}\) 中的项重合，将猜测乘以 \(x\)。
6. 求系数： 将 \(y_{PI}\) 及其导数代入原非齐次方程，解出系数 (P, Q, ...)。
7. 写出通解： \(y = y_{CF} + y_{PI}\)。
8. 应用条件： 利用关于 \(y\) 和 \(\frac{dy}{dx}\) 的两个条件解出常数 \(A\) 和 \(B\)。
你一定行！二阶微分方程将二次方程、微分和代数融合成一个精巧的过程。熟能生巧，加油！