考研数二17题｜常系数线性微分方程通解与稳定性深度精讲

待定系数法：非齐次特解的快速构造

对于非齐次方程 y'' + py' + qy = f(x)，当 f(x) 为以下形式时，可采用待定系数法：

• f(x) = Pₘ(x)（m次多项式）：设特解 y = x^k Qₘ(x)e^{0·x}，其中 k 为 0（若 r=0 不是特征根）、1（r=0 是单根）、2（r=0 是重根）。
• f(x) = e^{λx}Pₘ(x)：设 y = x^k Qₘ(x)e^{λx}，k 根据 λ 是否为特征根确定。
• f(x) = e^{λx}[Pₘ(x)cosωx + Qₙ(x)sinωx]：设 y = x^k e^{λx}[Rₛ(x)cosωx + Sₛ(x)sinωx]，s = max{m,n}，k 根据 λ±ωi 是否为特征根确定。

在2021考研数二17题中，非齐次项常为 e^{-x} 或 sinx 形式。例如当 f(x) = e^{-x} 时，若 -1 不是特征根，则设 y = Ae^{-x}；若 -1 是单根，则设 y = Axe^{-x}。关键在于准确判断 λ 与特征根的关系，这是高频易错点。

典型错误：忽略 x^k 因子。例如当 f(x)=e^{rx} 且 r 是特征根时，仍设 y=Ae^{rx}，代入后必得 0=非零值，导致矛盾。必须补充 x 因子以保证线性无关性。

常数变易法：普适性更强的特解求解

当 f(x) 形式复杂（如 ln x、tan x）或待定系数法失效时，常数变易法是通用解法。其核心思想是将齐次通解中的常数 C₁、C₂ 替换为待定函数 u₁(x)、u₂(x)：

其中 y₁、y₂ 是齐次方程的两个线性无关解。通过约束条件 u₁'y₁ + u₂'y₂ = 0，可得：

其中 W = y₁y₂'
- y₁'y₂ 是朗斯基行列式。积分后得 u₁、u₂，进而得到特解。

在2021数二17题的变式中，若非齐次项为 e^{-x} 且 -1 是特征根，待定系数法需设 y = Axe^{-x}，但代入后可能仍遇困难（如系数为零），此时常数变易法更可靠。例如对 y'' + 2y' + y = e^{-x}，特征根 r=-1（重根），设 y = u₁(x)e^{-x} + u₂(x)xe^{-x}，经计算可得 u₁ = x³/6, u₂ = -x²/2，故 y = x³e^{-x}/6。

注意：常数变易法计算量较大，但逻辑清晰，是验证待定系数法结果的有力工具。

稳定性判定：解的渐近行为分析

对于线性系统，解的稳定性完全由特征方程的根决定：

• 渐近稳定：所有特征根实部 < 0 → y→0 (x→∞)
• 不稳定：至少一个特征根实部 > 0 → |y|→∞
• 临界稳定：特征根实部 ≤ 0，且存在实部=0的单根 → 有界但不趋于零（如纯振荡）

在2021考研数二17题中，题目常隐含稳定性要求。例如若特征方程为 r² + 2r + 5 = 0，根为 r = -1±2i，实部均为 -1 < 0，故系统渐近稳定，解随时间衰减至零。

物理意义：在力学中，稳定性对应振动系统的阻尼特性。实部负值表示存在阻尼，振幅逐渐减小；实部正值表示能量输入，振幅发散；实部为零表示无阻尼自由振动。

常见误区：仅关注特征根的虚部（振荡），忽略实部（衰减/增长）。稳定性仅由实部决定！

步骤2：构造非齐次特解

2.1 分析非齐次项：f(x) = e^{-x}

2.2 判断特征根：特征根为 -1±2i，其中 -1 是实部，但 -1 本身不是特征根（特征根是复数）

2.3 设特解形式：因 -1 不是特征根，设 y = Ae^{-x}

2.4 代入原方程求 A：
y' = -Ae^{-x}, y'' = Ae^{-x}
代入：Ae^{-x} + 2(-Ae^{-x}) + 5(Ae^{-x}) = e^{-x}
(A
- 2A + 5A)e^{-x} = e^{-x} ⇒ 4A = 1 ⇒ A = 1/4

2.5 特解结果：y = (1/4)e^{-x}

关键技巧：计算系数时，合并同类项要仔细。此处 A
- 2A + 5A = 4A，易算错为 3A 或 5A。

变式处理：若 f(x) = e^{-x}cos2x（与齐次解振荡部分相同），则需设 y = x e^{-x}(Acos2x + Bsin2x)，补充 x 因子以保证线性无关性。

步骤3：合并通解

3.1 写出通解表达式：
y = y_h + y = e^{-x}(C₁cos2x + C₂sin2x) + (1/4)e^{-x}

3.2 简化通解（可选）：
y = e^{-x}[C₁cos2x + C₂sin2x + 1/4]

3.3 应用初始条件（若题目给出）：
例如 y(0)=1, y'(0)=0：
y(0) = C₁ + 1/4 = 1 ⇒ C₁ = 3/4
y' = -e^{-x}[C₁cos2x + C₂sin2x + 1/4] + e^{-x}[-2C₁sin2x + 2C₂cos2x]
y'(0) = -C₁
- 1/4 + 2C₂ = 0 ⇒ -3/4
- 1/4 + 2C₂ = 0 ⇒ C₂ = 1/2

3.4 写出特解：
y = e^{-x}[(3/4)cos2x + (1/2)sin2x + 1/4]

关键技巧：初始条件代入时，注意 y' 的乘积法则求导。可分步计算：y = e^{-x}·u(x)，则 y' = -e^{-x}u + e^{-x}u'。

常见错误：忘记 y 也需满足原方程。验证：将 y = (1/4)e^{-x} 代入原方程左边，应等于 e^{-x}。

步骤4：稳定性分析

4.1 提取特征根：r = -1 ± 2i

4.2 判断实部：实部 α = -1 < 0

4.3 得出结论：所有特征根实部为负，系统渐近稳定，即当 x→∞ 时，y→0

4.4 物理解释：解中 e^{-x} 因子主导长期行为，振荡项被指数衰减压制，故系统最终趋于平衡状态。

拓展应用：在控制系统中，稳定性决定系统是否收敛。本题解的稳定性意味着微小扰动不会导致系统发散，具有工程实用价值。

验证技巧：取 x=10，计算 y(10) ≈ e^{-10}[...] ≈ 0，直观体现稳定性。

例2：阻尼振动模型应用

题目：质量-弹簧-阻尼系统满足方程 my'' + cy' + ky = 0，其中 m=1kg, c=2N·s/m, k=5N/m。求位移 y(t)，并分析系统稳定性。

解析：
方程：y'' + 2y' + 5y = 0
特征根 r = -1±2i
通解：y = e^{-t}(C₁cos2t + C₂sin2t)
稳定性：实部 -1<0 → 渐近稳定，振动幅度随时间衰减至零。

物理意义：系统存在阻尼（c>0），初始扰动后，振动物体最终停在平衡位置。若 c=0（无阻尼），则 r=±√5i，解为纯振荡，临界稳定；若 c<0（负阻尼），则不稳定，振幅发散。

工程启示：设计减震系统时，需确保 c>0 且足够大以实现快速衰减，但过大阻尼会延长响应时间。本题对应欠阻尼情形（c²<4mk），是工程中最常见的工作状态。

例3：RLC电路瞬态响应

题目：RLC串联电路满足 Lq'' + Rq' + q/C = E(t)，其中 L=0.1H, R=2Ω, C=0.2F, E(t)=10e^{-t}。求电荷 q(t) 的通解及稳态行为。

解析：
方程：0.1q'' + 2q' + 5q = 10e^{-t}
化为标准型：q'' + 20q' + 50q = 100e^{-t}
齐次解：特征根 r²+20r+50=0 → r=-10±5√2
y_h = C₁e^{(-10+5√2)t} + C₂e^{(-10-5√2)t}
特解：设 y = Ae^{-t}，代入得 A=100/(1-20+50)=100/31
通解：q = y_h + 100e^{-t}/31
稳定性：两特征根均为负（-10±7.07≈-2.93,-17.07）→ 渐近稳定

工程意义：电荷 q(t) 随时间衰减至稳态值 100/31≈3.23C。电路从暂态过渡到稳态，符合物理实际。

与2021考研数二17题的关联：电路分析是微分方程的重要应用，稳定性分析决定电路是否收敛，避免振荡失控。

偏微分方程基础：从常微分到偏微分

典型方程：
• 热传导方程：∂u/∂t = k∇²u
• 波动方程：∂²u/∂t² = c²∇²u
• 拉普拉斯方程：∇²u = 0

解法思路：分离变量法
设 u(x,t) = X(x)T(t)，代入方程分离为常微分方程组，再结合边界条件求解。

实例：一维热传导 u_t = ku_{xx}, u(0,t)=u(L,t)=0, u(x,0)=f(x)
解得：u(x,t) = Σ B_n sin(nπx/L) e^{-k(nπ/L)²t}
其中 B_n = (2/L)∫₀ᴸ f(x)sin(nπx/L)dx

与2021考研数二17题的关联：偏微分方程的分离变量法最终归结为常微分方程的求解，特征方程法仍是核心工具。

数值解法：欧拉法与龙格-库塔法

欧拉法：
y_{n+1} = y_n + hf(x_n, y_n)
一阶精度，简单但误差累积快。

四阶龙格-库塔法（RK4）：
k₁ = hf(x_n, y_n)
k₂ = hf(x_n+h/2, y_n+k₁/2)
k₃ = hf(x_n+h/2, y_n+k₂/2)
k₄ = hf(x_n+h, y_n+k₃)
y_{n+1} = y_n + (k₁+2k₂+2k₃+k₄)/6
四阶精度，工程中广泛使用。

应用示例：对 y' = -2y + e^{-x}, y(0)=1，取 h=0.1，计算 y(0.1)：
k₁ = 0.1×(-2×1 + 1) = -0.1
k₂ = 0.1×(-2×(1-0.05) + e^{-0.05}) ≈ 0.1×(-1.9 + 0.9512) ≈ -0.09488
...（略）
y(0.1) ≈ 1 + ( -0.1 + 2×(-0.09488) + 2×(-0.09497) + (-0.09048) )/6 ≈ 0.9048

与2021考研数二17题的关联：数值解法适用于无法解析求解的方程，是理论解的重要补充。理解数值方法需以解析解为基础。