如何确定二阶微分方程待定特解的形式？

Question

在求解非齐次二阶线性微分方程时，如何根据自由项（如多项式、指数函数、三角函数及其组合）正确设定待定特解的形式是一个关键问题。常见困惑是：当自由项与齐次解重复时，为何要将特解形式乘以 $ x $ 或 $ x^2 $？例如，对于方程 $ y'' - 2y' + y = e^x $，其齐次解含 $ e^x $，此时为何特设解为 $ y_p = Ax^2 e^x $ 而非 $ Ae^x $？若不乘以足够幂次的 $ x $，会导致代入后无法匹配原方程右边项。掌握这一规则的逻辑依据与适用边界，是准确构造特解形式的技术难点。

Answer 1

非齐次二阶线性微分方程特解构造的深度解析

1. 问题背景与基本概念回顾

在求解形如：

\[ y'' + p(x)y' + q(x)y = f(x) \]

的非齐次二阶线性微分方程时，通解结构为：

\[ y = y_h + y_p \]

其中 \( y_h \) 是对应齐次方程的通解，\( y_p \) 是一个特解。关键挑战在于如何根据自由项 \( f(x) \) 正确设定 \( y_p \) 的形式。

常见的 \( f(x) \) 包括多项式、指数函数 \( e^{\alpha x} \)、三角函数 \( \sin(\beta x), \cos(\beta x) \)，以及它们的乘积或线性组合。

标准方法是“待定系数法”，其核心是：基于 \( f(x) \) 的类型预设一个含待定系数的形式，代入原方程求解系数。

然而，当 \( f(x) \) 的某部分恰好是齐次解的一部分时，直接代入会导致矛盾——即无法匹配右边项。

例如，考虑方程：

\[ y'' - 2y' + y = e^x \]

其特征方程为 \( r^2 - 2r + 1 = 0 \)，解得重根 \( r = 1 \)，故齐次解为：

\[ y_h = (C_1 + C_2 x)e^x \]

2. 冲突出现：为何不能直接设 \( y_p = Ae^x \)？

若尝试设特解为 \( y_p = Ae^x \)，则计算导数：

• \( y_p' = Ae^x \)
• \( y_p'' = Ae^x \)

代入原方程：

\[ Ae^x - 2Ae^x + Ae^x = 0 \neq e^x \]

左边恒为0，无法等于右边的 \( e^x \)。这意味着 \( Ae^x \) 实际上是齐次解的一部分，不能作为非齐次方程的有效特解。

因此，必须对原始猜测形式进行“修正”——引入额外因子 \( x^k \)，使得新形式不再属于齐次解空间。

3. 修正机制：乘以 \( x^k \) 的数学逻辑

设原猜测形式为 \( y_p^{(0)} = e^x \)，但该函数满足齐次方程，即：

\[ L[y_p^{(0)}] = 0 \]

其中 \( L \) 是微分算子 \( L = D^2 - 2D + 1 \)（\( D = d/dx \)）。

为了打破这种“零响应”，我们尝试更高阶的试探函数。令：

\[ y_p = A x^k e^x \]

目标是找到最小的 \( k \)，使得 \( L[y_p] \neq 0 \) 且能匹配 \( e^x \)。

k 值	形式	是否属于齐次解	L[y_p] 是否可非零
0	Ae^x	是	否（恒为0）
1	Axe^x	是（含于 y_h）	仍为0
2	Ax²e^x	否	可非零

可见，由于特征根 \( r=1 \) 是二重根，需乘以 \( x^2 \) 才能跳出齐次解空间。

4. 一般规则与分类表

设非齐次项为：

\[ f(x) = P_m(x) e^{\alpha x} \cos(\beta x) \quad \text{或} \quad P_m(x) e^{\alpha x} \sin(\beta x) \]

则特解形式为：

\[ y_p = x^k e^{\alpha x} \left[ Q_m(x) \cos(\beta x) + R_m(x) \sin(\beta x) \right] \]

其中 \( k \) 是“重数”——即 \( \alpha + i\beta \) 作为特征根的重数（若不是根，则 \( k=0 \)）。

1. 若 \( \alpha + i\beta \) 不是特征根 → \( k = 0 \)
2. 若是单根 → \( k = 1 \)
3. 若是重根 → \( k = 2 \)

例如，在本例中 \( \alpha = 1, \beta = 0 \)，\( r=1 \) 是二重根 ⇒ \( k=2 \)，故设 \( y_p = A x^2 e^x \)。

5. 验证：代入正确形式求解

设 \( y_p = A x^2 e^x \)，计算导数：

y_p'  = A(2x e^x + x^2 e^x) = A e^x (2x + x^2)
y_p'' = A[(2 + 2x + 2x + x^2) e^x] = A e^x (x^2 + 4x + 2)

代入原方程：

\[ y_p'' - 2y_p' + y_p = A e^x (x^2 + 4x + 2 - 4x - 2x^2 + x^2) = A e^x (2) \]

令其等于 \( e^x \) ⇒ \( 2A = 1 \) ⇒ \( A = 1/2 \)

因此，特解为 \( y_p = \frac{1}{2} x^2 e^x \)

6. 流程图：特解构造决策路径

graph TD
    A[开始: 给定非齐次方程] --> B[求齐次解 y_h]
    B --> C[分析自由项 f(x) 类型]
    C --> D{f(x) 形式?}
    D -->|P(x)| E[设多项式形式]
    D -->|e^{αx}| F[设 Ae^{αx}]
    D -->|sin/cos| G[设组合形式]
    D -->|混合| H[设复合形式]
    H --> I[检查 α+iβ 是否为特征根]
    I --> J{k = 重数}
    J --> K[乘以 x^k]
    K --> L[代入原方程求系数]
    L --> M[得到 y_p]

7. 常见错误与边界情况

• 误判重数：将单根当作重根，导致过度乘以 \( x \)，虽可解但增加复杂度。
• 忽略组合项：如 \( f(x) = e^x \cos x \)，应整体判断复根 \( 1+i \) 是否为特征根。
• 多项式次数不匹配：若 \( f(x) = x^2 \)，应设 \( y_p = Ax^2 + Bx + C \)，而非仅 \( Ax^2 \)。
• 算子法替代方案：可用逆算子法 \( y_p = \frac{1}{L(D)} f(x) \) 避免猜测，适合编程实现。

8. 编程视角：自动化特解构造算法框架

在科学计算系统（如SymPy、MATLAB）中，可通过以下伪代码实现：

def construct_particular_solution(f, char_roots):
    alpha, beta = extract_exponent_and_frequency(f)
    root = complex(alpha, beta)
    k = multiplicity(root, char_roots)
    
    base_form = guess_base_form(f)
    corrected_form = multiply_by_x_power(base_form, k)
    
    yp = substitute_and_solve(corrected_form, f)
    return yp

该逻辑可用于构建微分方程求解器的核心模块。

9. 拓展思考：高阶方程与系统耦合

对于高阶方程或线性微分方程组，相同原则适用。例如，在控制系统中，输入信号（自由项）若与系统自然模态（齐次解）共振，则输出响应需引入时间增长因子（如 \( t^n e^{st} \)），这正是乘以 \( x^k \) 的物理意义——反映共振放大效应。

在数值仿真中，若未正确处理此类情况，可能导致稳态误差或发散。

10. 实践建议与调试技巧

1. 先求齐次解，明确所有基本模态。
2. 将 \( f(x) \) 分解为独立项，逐项处理。
3. 对每项判断其对应特征根重数。
4. 使用表格法快速匹配形式。
5. 代入后若方程退化为0=常数，说明 \( k \) 不足，需提升幂次。
6. 利用计算机代数系统验证手算结果。
7. 注意线性无关性：特解不能与齐次解线性相关。
8. 若 \( f(x) \) 含多个频率，分别求特解再叠加。
9. 避免重复设定：如 \( e^x \) 和 \( xe^x \) 已在齐次解中，则特解必须从 \( x^2 e^x \) 起。
10. 记录典型错误案例，建立调试清单。