如何求解二阶常系数齐次微分方程的通解？

Question

如何求解形如 $ y'' + py' + qy = 0 $ 的二阶常系数齐次微分方程的通解？当特征方程 $ r^2 + pr + q = 0 $ 的判别式小于零时，出现共轭复根，此时通解应如何构造？是否只需将复根代入指数函数形式即可？常见误区是忽略实部与虚部的分离，导致解的形式错误。此外，如何从复数解还原为实值通解？是否必须引入正弦与余弦函数？这些问题在振动系统与电路分析中尤为关键。

Answer 1

一、二阶常系数齐次微分方程的通解求解方法

对于形如 $ y'' + py' + qy = 0 $ 的二阶常系数齐次线性微分方程，其求解的核心在于特征方程的分析。我们首先引入特征方程：

$$
r^2 + pr + q = 0
$$

该方程的根决定了微分方程通解的形式。根据判别式 $ \Delta = p^2 - 4q $ 的符号，可分为三种情况：

1. Δ > 0：两个不等实根 $ r_1, r_2 $，通解为 $ y(x) = C_1 e^{r_1 x} + C_2 e^{r_2 x} $
2. Δ = 0：重根 $ r $，通解为 $ y(x) = (C_1 + C_2 x)e^{rx} $
3. Δ < 0：共轭复根 $ r = \alpha \pm i\beta $，此时进入本文重点讨论的情形。

二、判别式小于零时的复根情形分析

当 $ \Delta < 0 $ 时，特征方程的根为一对共轭复数：

$$
r_{1,2} = \frac{-p}{2} \pm i \frac{\sqrt{|\Delta|}}{2} = \alpha \pm i\beta
$$

其中 $ \alpha = -p/2 $ 是实部，$ \beta = \sqrt{|\Delta|}/2 $ 是虚部。若直接将复根代入指数形式，得到两个复值解：

• $ y_1(x) = e^{(\alpha + i\beta)x} $
• $ y_2(x) = e^{(\alpha - i\beta)x} $

但这并不意味着通解就是 $ y(x) = C_1 e^{(\alpha + i\beta)x} + C_2 e^{(\alpha - i\beta)x} $（其中 $ C_1, C_2 \in \mathbb{C} $），因为通常我们要求的是实值函数解，尤其是在物理建模如振动系统或RLC电路中。

三、从复数解还原为实值通解的方法

利用欧拉公式：

$$
e^{i\theta} = \cos\theta + i\sin\theta
$$

可将复指数展开为三角函数形式：

$$
e^{(\alpha + i\beta)x} = e^{\alpha x}(\cos(\beta x) + i\sin(\beta x))
$$

$$
e^{(\alpha - i\beta)x} = e^{\alpha x}(\cos(\beta x) - i\sin(\beta x))
$$

将这两个解进行线性组合，并取实部与虚部分离，可得两个线性无关的实值解：

• $ y_1(x) = e^{\alpha x}\cos(\beta x) $
• $ y_2(x) = e^{\alpha x}\sin(\beta x) $

因此，实值通解为：

$$
y(x) = e^{\alpha x}(C_1 \cos(\beta x) + C_2 \sin(\beta x)), \quad C_1, C_2 \in \mathbb{R}
$$

Δ 值	根类型	通解形式
Δ > 0	两不同实根	$ C_1 e^{r_1 x} + C_2 e^{r_2 x} $
Δ = 0	重根	$ (C_1 + C_2 x)e^{r x} $
Δ < 0	共轭复根	$ e^{\alpha x}(C_1 \cos\beta x + C_2 \sin\beta x) $

四、常见误区与技术细节解析

在实际应用中，开发者或工程师常犯以下错误：

1. 误认为复根直接代入即可得到实解，忽略了解空间必须是实向量空间的要求。
2. 未分离实部与虚部，导致后续数值仿真或控制系统设计中出现非物理结果。
3. 在编程实现时（如Python或MATLAB），使用复数解而未提取实部，造成数据可视化异常。

示例代码（Python）展示如何构造实值解：

import numpy as np
import matplotlib.pyplot as plt

# 参数设置：α = -0.5, β = 2
alpha = -0.5
beta = 2.0
x = np.linspace(0, 10, 500)
C1, C2 = 1.0, 0.5

# 构造实值通解
y = np.exp(alpha * x) * (C1 * np.cos(beta * x) + C2 * np.sin(beta * x))

plt.plot(x, y)
plt.title("Damped Oscillation: Real-valued Solution")
plt.xlabel("x")
plt.ylabel("y(x)")
plt.grid(True)
plt.show()

五、工程应用场景与延伸思考

在IT相关领域，尤其是嵌入式系统、信号处理和控制系统中，此类微分方程广泛存在于：

• 机械振动模型（如悬架系统）
• RLC电路中的电流响应
• 滤波器设计中的阶跃响应分析
• 机器人运动控制中的阻尼振荡行为建模

以下是典型RLC串联电路的微分方程模型：

$$
L\frac{d^2q}{dt^2} + R\frac{dq}{dt} + \frac{1}{C}q = 0
$$

对应特征方程为：

$$
Lr^2 + Rr + \frac{1}{C} = 0
$$

当 $ R^2 < 4L/C $ 时，出现欠阻尼状态，即复根情形，电流呈现衰减振荡。

graph TD A[微分方程 y''+py'+qy=0] --> B{判别式 Δ = p²-4q} B -->|Δ > 0| C[两实根 → 指数解] B -->|Δ = 0| D[重根 → 多项式×指数] B -->|Δ < 0| E[复根 → 三角×指数] E --> F[使用欧拉公式展开] F --> G[提取实部与虚部] G --> H[构造实值通解] H --> I[应用于振动/电路系统]