欢迎来到微分方程的世界！
在本章中，我们将跨越简单代数的范畴，开始探讨事物是如何变化的。微分方程 (Differential Equation, DE) 就是包含变化率（导数）的方程。为什么这很重要？因为现实世界中的几乎一切事物——从咖啡如何变凉，到蹦极者的摆动——都可以用这些方程来描述。读完这些笔记后，你将能够解开这些难题，并预测系统随时间演变的规律。如果起初觉得有点棘手，别担心；我们会将其拆解成简单、易于掌握的步骤！
1. 一阶微分方程：积分因子法
你已经知道如何解像 \(\frac{dy}{dx} = 2x\) 这样的简单方程。但如果方程长成这样呢：\(\frac{dy}{dx} + P(x)y = Q(x)\)？这时我们使用一个巧妙的技巧，称为积分因子 (Integrating Factor)。
检查“标准形式”
在做任何事之前，你的方程必须完全符合这个样子：
\(\frac{dy}{dx} + P(x)y = Q(x)\)
如果 \(\frac{dy}{dx}\) 前面有数字或 \(x\)，你必须先将整个方程除以该项！这是学生最常犯的错误。
循序渐进：积分因子法
- 找出积分因子，记作 \(I(x)\)。公式为：\(I(x) = e^{\int P(x) dx}\)。
- 将标准形式方程中的每一项都乘以这个 \(I(x)\)。
- 方程的左边现在会神奇地变成 \((I(x) \cdot y)\) 的导数。你可以将方程重写为：\(\frac{d}{dx}(I(x) \cdot y) = I(x) \cdot Q(x)\)。
- 将等式两边同时对 \(x\) 积分。
- 整理得到 \(y\)。别忘了你的积分常数 (\(+C\))！
小知识速览：
如果你的解包含 \(+C\)，它被称为通解 (General Solution)（它代表一族曲线）。如果你得到坐标（初始条件）来求出 \(C\) 的特定值，它则被称为特解 (Particular Solution)。
重点总结：积分因子法实际上是“还原”了乘积法则的微分过程。记住，计算 \(I(x)\) 之前一定要先整理成标准形式！
2. 二阶齐次微分方程
现在我们要探讨包含二阶导数的方程：\(a\frac{d^2y}{dx^2} + b\frac{dy}{dx} + cy = 0\)。因为等式右边为零，我们称之为齐次 (Homogeneous) 方程。
辅助方程 (Auxiliary Equation)
为了求解，我们“猜测”其解的形式为 \(y = e^{mx}\)。这引导我们得出辅助方程：
\(am^2 + bm + c = 0\)
你只需像解普通二次方程那样求解即可！你得到的根（解）类型决定了解的形式：
- 情况 1：两个相异实根 (\(m_1\) 和 \(m_2\))
解：\(y = Ae^{m_1x} + Be^{m_2x}\) - 情况 2：重实根 (\(m\))
解：\(y = (A + Bx)e^{mx}\) - 情况 3：复数根 (\(m = \alpha \pm \beta i\))
解：\(y = e^{\alpha x}(A\cos(\beta x) + B\sin(\beta x))\)
你知道吗？情况 3 发生的情境就是物体产生振荡时，例如吉他弦的振动！
重点总结：解二阶齐次微分方程就像解一个二次方程，然后套用正确的答案“模板”一样简单。
3. 二阶非齐次微分方程
如果方程右边不等于零呢？\(a\frac{d^2y}{dx^2} + b\frac{dy}{dx} + cy = f(x)\)。
完整解由两部分组成：\(y = \text{互补函数 (CF)} + \text{特解 (PI)}\)。
如何求解：
- 求 CF：将方程设为零，按上一节的方法求解。
- 求 PI：根据 \(f(x)\) 的形式进行“猜测”。
常见的 PI 猜测：
- 如果 \(f(x) = \text{多项式 (例如 } x^2)\)，试猜 \(y = px^2 + qx + r\)。
- 如果 \(f(x) = e^{kx}\)，试猜 \(y = \lambda e^{kx}\)。
- 如果 \(f(x) = \sin(kx)\) 或 \(\cos(kx)\)，试猜 \(y = p\cos(kx) + q\sin(kx)\)。
避免常见错误：如果你的 PI “猜测”已经存在于 CF 中，那是行不通的！你必须将猜测乘以 \(x\)（甚至 \(x^2\)）来使其与 CF 区分开来。
重点总结：你可以将 CF 视为系统的“自然”行为，而 PI 则视为对外力的“反应”。
4. 微分方程建模：简谐运动与阻尼
微分方程是运动学 (Kinematics) 的语言。记住这些术语：
- 位移 (Displacement) = \(x\)
- 速度 (Velocity) = \(\dot{x}\) 或 \(\frac{dx}{dt}\)
- 加速度 (Acceleration) = \(\ddot{x}\) 或 \(\frac{d^2x}{dt^2}\)
简谐运动 (SHM)
标准的 SHM 方程为 \(\ddot{x} = -\omega^2x\)。
其解永远是：\(x = A\cos(\omega t + \phi)\) 或 \(x = p\cos(\omega t) + q\sin(\omega t)\)。
阻尼振荡 (Damped Oscillations)
在现实世界中，摩擦力会减缓运动。我们将其建模为：\(a\ddot{x} + b\dot{x} + cx = 0\)。
辅助方程的根决定了它减速的方式：
- 过阻尼 (Over-damping)：两个实根。系统缓慢回到平衡状态而不产生振荡。（例如：门闭门器）。
- 临界阻尼 (Critical damping)：重根。回到平衡状态最快且不产生振荡的方式。（例如：汽车悬挂系统）。
- 欠阻尼 (Under-damping)：复数根。系统会来回摆动（振荡），但摆幅会越来越小。
重点总结：你对二次方程运算的结果，能直接预测一个物理系统会是弹跳、缓慢回归还是完美停留在平衡点。
5. 联立一阶微分方程
有时两个变量相互依赖。例如，森林中的兔子数量 (\(x\)) 和狐狸数量 (\(y\))。这些被称为耦合方程 (coupled equations)。
如何求解：
- 从两个方程开始，例如 \(\frac{dx}{dt} = ax + by\) 和 \(\frac{dy}{dt} = cx + dy\)。
- 将第一个方程再次微分，得到 \(\frac{d^2x}{dt^2}\) 项。
- 将第二个方程代入这个新方程中以“消去”\(y\)。
- 你最终会得到一个关于 \(x\) 的单一二阶微分方程，然后就可以使用我们之前学过的方法来求解！
鼓励的话：这算是本章的“最终 Boss”，但它其实只是把你已经学过的两个东西结合起来：代入法和二阶微分方程！
重点总结：联立微分方程通过将其转化为一个二阶微分方程来求解。先解出其中一个变量，再利用它求出第二个变量。
总结检查清单
- 你能识别一阶线性微分方程并找到积分因子吗？
- 你记得辅助方程的三种情况吗？
- 你能根据方程右侧选择合适的特解 (PI) 吗？
- 你能将二阶微分方程的数学与过阻尼、欠阻尼、临界阻尼联系起来吗？
- 你能通过消去变量来解联立微分方程吗？