非齐次方程特解求解：待定系数法与特征根共振修正

1. 从一道“猜答案”的题说起为什么要单独学非齐次方程很多人学微分方程时被卡住的第一个坎不是齐次方程而是非齐次方程。齐次方程好歹有套固定的特征方程流程代入、求根、写通解三步走完清清楚楚。但一遇到非齐次等式右边多了个函数很多人就懵了——这多出来的东西到底怎么处理先说结论非齐次方程的通解永远等于“对应齐次方程的通解”加上“一个特解”。用公式写就是 y y_h y*。这个结构是整章的地基后面所有方法都是围绕“怎么找到那个 y*”展开的。常系数非齐次线性微分方程的标准形式是y p y q y f(x)其中 p、q 是常数f(x) 是非齐次项也叫右端项。这门课的核心问题只有一个已知 f(x) 的具体形式怎么把特解 y* 猜出来。为什么说“猜”因为常系数线性微分方程有一个特别好的性质——右端项是某种类型的函数时特解往往也是同一类型。于是数学家们总结出了一套“待定系数法”先根据右端项的形式设出特解的基本形状再把特解代回方程去解那几个未知系数。整个过程看起来像是猜实际上每一步都有规则可循猜得有章法。这篇文章只讲两种最常见的右端项形式多项式型 f(x) e^(λx) P_m(x)以及三角型 f(x) e^(λx) [P_m(x) cos(ωx) Q_n(x) sin(ωx)]。掌握了这两种应付考试和工程里的绝大多数问题都够用了。2. 核心思路待定系数法的底层逻辑与适用边界2.1 为什么可以“设出特解的形式”很多人第一次接触待定系数法心里会有个疑问凭什么说右端项是 e^(λx) 乘以多项式特解就一定是同样形式这背后的原因得从“微分算子”的角度理解才不至于变成背公式。把求导看成一种运算记作 D则 D(e^(λx) x^k) e^(λx)(λ x^k k x^(k-1))。你会发现e^(λx) 乘以任何多项式经过求导后仍然保持 e^(λx) 乘以另一个多项式的形式。而常系数线性微分方程左边的每项都是 y 的各阶导数乘以常数所以整个方程左边作用在 e^(λx) Q_m(x) 上的结果仍然是 e^(λx) 乘一个多项式。关键结论如果让特解 y* e^(λx) Q_m(x)那么代入方程后两边都含有公因子 e^(λx)可以约掉。剩下的就是两个多项式之间的恒等式通过对比系数就能把 Q_m(x) 里的每个系数确定下来。这个逻辑的妙处在于它把“微分方程求解”转化成了“多项式恒等对比”机械化程度非常高。哪怕你对微分方程本身不太有感觉只要会做多项式展开和合并同类项就能算出正确的特解。2.2 特解“乘 x^k”的本质共振修正这里有一个所有初学者都会踩的点什么时候特解需要额外乘以 x 的幂次规则一句话概括设出特解形式后如果发现特解中的每一项都与对应齐次方程的通解线性相关也就是被“特征根”撞上了就必须通解乘 x^kk 取使两者不再线性相关的最小正整数。打个现实里的比方你在荡秋千如果你推秋千的节奏恰好和秋千自然振动的频率完全一致秋千就会越荡越高这就是共振。微分方程里的“共振”也一样——外力 f(x) 的形式如果恰好和齐次方程本身的某个解“同频”特解就不能简单地沿用原形式必须乘以 x^k 来调整“震幅的多项式次数”才能保证方程左边能产生那个右端项。这里的 k 取多少取决于 λ或 λ ± iω作为特征根的重数λ 不是特征根k 0不乘λ 是单特征根重数 1k 1乘一个 xλ 是二重特征根重数 2k 2乘 x²。为了把这个判断做对你必须先把齐次方程的特征方程解出来。所以做题顺序永远是先求齐次通解再判断 λ 是不是特征根最后才设特解。跳过第一步直接设特解是错误率最高的一种习惯。2.3 使用条件什么时候“猜”得出来待定系数法虽然方便却有限制。它只适用于常系数线性方程而且右端项必须是以下类型的有限组合多项式、指数函数、正弦余弦函数以及它们的乘积。也就是说 f(x) 必须是“指数多项式三角型”的有限线性组合。凡是右端项出现分式、对数、绝对值、分段函数、含任意函数形式的项待定系数法就失效了。这时候得换武器比如常数变易法。常数变易法适用面更广但积分往往比较痛苦。在工程应用中能用待定系数法的题目占大多数所以优先掌握它是性价比最高的路径。3. 第一种形式拆解f(x) e^(λx) P_m(x)3.1 特解设法的完整表格右端项是 e^(λx) 乘以 m 次多项式 P_m(x)这在物理和工程中极其常见比如强迫振动中的指数衰减驱动力、电路中的指数激励等。拿到这种题直接按表设特解条件特解形式λ 不是特征根y* e^(λx) Q_m(x)λ 是单特征根y* x e^(λx) Q_m(x)λ 是二重特征根y* x² e^(λx) Q_m(x)这里的 Q_m(x) 是系数待定的 m 次多项式写成 a_m x^m a_(m-1) x^(m-1) ... a_0有多少个系数要确定就有多少个未知数。每个未知数通过代入后对比同类项系数列出一个线性方程最后解方程组得到。一个容易被忽略的细节即使原方程右端的 P_m(x) 缺了某些低次项比如只有 x² 项没有一次项和常数项你设 Q_m(x) 时也必须把完整的多项式形式写出来——所有次数从 m 到 0 的项都要带。原因是求导和约化过程中低次项会反复出现你只设高次项的话代入方程后对比系数时会发现低次项的方程对不上结果永远解不出来。3.2 实操案例λ 非特征根的情形先看一个最简单的。解方程y - 3y 2y x² 1第一步求齐次方程 y - 3y 2y 0 的通解。特征方程是 r² - 3r 2 0因式分解得 (r-1)(r-2) 0特征根 r₁ 1r₂ 2。齐次通解为y_h C₁ e^x C₂ e^(2x)第二步判断 λ。右端项是多项式 x² 1可以理解为 e^(0·x) 乘以二次多项式所以 λ 0。特征根是 1 和 20 不是特征根于是特解设为y* A x² B x C注意虽然右端没有一次项Q(x) 依然要保留 Bx 项和常数项 C。第三步代入求系数。逐项求导y* 2A x By* 2A代入原方程(2A) - 3(2A x B) 2(A x² B x C) x² 1整理得2A x² (2B - 6A) x (2A - 3B 2C) x² 1对比系数得到方程组x² 系数2A 1所以 A 1/2x 系数2B - 6A 0代入 A 1/2 得 2B - 3 0B 3/2常数项2A - 3B 2C 1代入 A 1/2B 3/2得 1 - 9/2 2C 1即 2C 9/2C 9/4。特解为 y* (1/2)x² (3/2)x 9/4原方程通解y C₁ e^x C₂ e^(2x) (1/2)x² (3/2)x 9/4这题虽然简单但系数对比的过程完整展示了方法的全貌。建议你在纸上亲自算一遍尤其是代入整理那一步慢一点没关系算多了自然就快。3.3 实操案例λ 是特征根的情形现在把难度升一级。解方程y - 2y y e^x齐次方程 y - 2y y 0 的特征方程是 r² - 2r 1 0即 (r-1)² 0特征根 r 1 是二重根。齐次通解y_h (C₁ C₂ x) e^x这里有个关键判断右端是 e^x对应 λ 1恰好是二重特征根。正是最需要警惕的情况——如果不乘 x直接设 y* A e^x代入后左边会得 0怎么都凑不出右端的 e^x。根据规则二重根要乘 x²所以特解设为y* A x² e^x可能有人会问为什么不乘 x 而乘 x²验证一下就知道。逐项求导y* A(2x x²) e^xy* A(2 4x x²) e^x代入 y - 2y y每一项都含有 e^x提取公因子后得到A[(2 4x x²) - 2(2x x²) x²] e^x A[2 4x x² - 4x - 2x² x²] e^x 2A e^x方程要求这个结果等于 e^x所以 2A 1A 1/2。通解为y (C₁ C₂ x) e^x (1/2) x² e^x注意特解里 x² 的来源因为特征根的重数是 2虚设多项式 Q(x) 的次数本应为 0常数乘以 x² 后次数升到 2。如果把特解写成 (A x² B x C)e^x 再代入B 和 C 最终都会被解成 0所以直接从 A x² e^x 设起省去无效计算。这个案例其实透露了一个通用经验一旦判定 λ 是特征根特解中的多项式次数会整体抬高。二重根抬高 2 次单根抬高 1 次。抬高后原来低次项对应的系数在代入后多半会自行消失但为了严谨和避免初学者出错初期还是建议写完整形式熟练后再跳过那些注定为 0 的项。4. 第二种形式拆解三角型右端项4.1 特解设法的完整表格与原理第二类常见右端项是含三角函数的f(x) e^(λx) [P_m(x) cos(ωx) Q_n(x) sin(ωx)]这类问题的物理背景非常直白——强迫振动中外力本身是简谐激励可能还带指数衰减或增长因子。求解时依然用待定系数法但形式上要复杂一些因为三角函数通过求导会互相转化sin 会变成 coscos 会变成 sin所以特解必须同时包含正余弦两项。特解设法的规则条件特解形式λ iω 不是特征根y* e^(λx)[S_N(x) cos(ωx) T_N(x) sin(ωx)]λ iω 是特征根单根y* x e^(λx)[S_N(x) cos(ωx) T_N(x) sin(ωx)]这里 N 取 m 和 n 中的较大值。S_N(x) 和 T_N(x) 都是次数为 N 的待定多项式各项系数均为待定未知数。有些教材只写了右端是 cos 或 sin 单项的情形但实际做题时右端往往是二者同时出现。哪怕右端只有 A cos(ωx) 一项特解也必须设成 C₁ cos(ωx) C₂ sin(ωx) 两项。原因在于求导会改变函数类型如果你只设 cos 形式代入方程后 sin 项的系数将无法匹配整体无解。这是这部分最容易丢掉的特解形式完整性问题。判断是否与特征根相撞看的是复数 λ iω严格来说不是单独看 λ 或 ω。比如 λ 0ω 2那就看 i2 是否是特征方程 r² pr q 0 的根。只有特征方程存在纯虚根且虚部恰好等于 ω 时才需要乘 x。这个判断方式很多人记混值得反复强调。4.2 实操案例右端只有 cos 的情形解方程y 4y cos(2x)齐次方程特征方程是 r² 4 0特征根 r ±2i。齐次通解y_h C₁ cos(2x) C₂ sin(2x)注意右端 cos(2x) 对应 λ 0ω 2那么 λ iω 2i恰好是特征根。也就是说外力频率和系统固有频率相同纯粹是共振情形。特解必须乘 xy* x[A cos(2x) B sin(2x)]逐项求导y* A cos(2x) B sin(2x) x[-2A sin(2x) 2B cos(2x)]y* -4A sin(2x) 4B cos(2x) x[-4A cos(2x) - 4B sin(2x)]代入 y 4y第一部分的非 x 项为 -4A sin(2x) 4B cos(2x) 第二部分的 x 项为 x[-4A cos(2x) - 4B sin(2x)] 加上 4y* 得 4x[A cos(2x) B sin(2x)]。 x 的项正好抵消。剩下的-4A sin(2x) 4B cos(2x) cos(2x)对比系数cos(2x) 系数4B 1B 1/4sin(2x) 系数-4A 0A 0。特解 y* (1/4) x sin(2x)通解y C₁ cos(2x) C₂ sin(2x) (1/4) x sin(2x)注意几个细节特解中 A 最终为 0意味着 cos 分量在调整后没有出现但这不意味着一开始可以不设 A。因为求导过程中B sin(2x) 项的导数会产生 cos 项为了平衡右端的 cos(2x)A 才被解为 0。如果你只设 xB sin(2x)结果不可能对得上。4.3 实操案例带指数的三角型再看一个混合型y - y e^x cos(2x)齐次方程 y - y 0 的特征方程是 r² - 1 0特征根 r ±1。齐次通解y_h C₁ e^x C₂ e^(-x)右端项是 e^x cos(2x)对应 λ 1ω 2λ iω 1 2i。这个复数不是特征根特征根只有 ±1所以特解形式直接设为y* e^x[A cos(2x) B sin(2x)]这里右端没有 sin 项但特解必须保留 sin 项理由同上。逐项求导。先记住一个恒等式如果 u e^x cos(2x)则 u e^x(cos(2x) - 2 sin(2x))u e^x(-3 cos(2x) - 4 sin(2x))。同理若 v e^x sin(2x)则 v e^x(sin(2x) 2 cos(2x))v e^x(-3 sin(2x) 4 cos(2x))。于是y* e^x[(A 2B) cos(2x) (B - 2A) sin(2x)]y* e^x[(-3A 4B) cos(2x) (-4A - 3B) sin(2x)]代入 y - y公因子 e^x 提取后[(-3A 4B) - A] cos(2x) [(-4A - 3B) - B] sin(2x) cos(2x)整理(-4A 4B) cos(2x) (-4A - 4B) sin(2x) cos(2x)对比系数cos(2x)-4A 4B 1sin(2x)-4A - 4B 0联立解得 A -1/8B 1/8。特解 y* e^x[-1/8 cos(2x) 1/8 sin(2x)]通解y C₁ e^x C₂ e^(-x) (1/8) e^x[sin(2x) - cos(2x)]这题的运算量明显上来了。初学时建议先把 e^x cos(2x) 和 e^x sin(2x) 各自的一阶二阶导数单独写在草稿上再代入不要直接在原式里硬展开那样极易出错。熟练之后可以尝试用“指数换元法”——把 cos 和 sin 换成复指数形式求解后再取实部计算会对称一些但前期不推荐。5. 实战综合叠加原理与多形式右端项5.1 右端项是多块相加时的处理策略实际题目里右端项经常不是单一形式而是多个不同类型函数相加比如y 3y 2y x e^(-x)很多新手会试图直接设一个同时包含 x 和 e^(-x) 的特解这种做法不是不可以但特解形式容易设错尤其当两个子项对应的“撞根”情况不一致时很容易漏乘或重乘 x。最稳妥的做法是调用叠加原理如果 f(x) f₁(x) f₂(x)那么分别求方程 y 3y 2y f₁(x) 和 y 3y 2y f₂(x) 的特解 y₁* 和 y₂*原方程的一个特解就是 y₁* y₂*。这套处理方式和线性代数里“非齐次线性方程组特解分解”是同一个思想核心依据就是微分算子本身是线性的L(y₁* y₂*) L(y₁*) L(y₂*) f₁ f₂ f所以分批处理思路清晰每一小步只需要判断一个子项与特征根的关系出错概率大幅降低。5.2 多块右端项的完整解题演示解方程y 3y 2y x e^(-x)齐次方程特征方程 r² 3r 2 0因式分解后 (r1)(r2) 0特征根 r₁ -1r₂ -2。齐次通解y_h C₁ e^(-x) C₂ e^(-2x)先处理 f₁(x) x对应 λ 0。0 不是特征根设 y₁* A x B。代入y₁* Ay₁* 00 3A 2(Ax B) x整理得 2Ax (3A 2B) x。对比系数2A 1A 1/23A 2B 0代得 3/2 2B 0B -3/4。所以 y₁* (1/2)x - 3/4。再处理 f₂(x) e^(-x)对应 λ -1。这里关键点来了λ -1 恰好是特征方程的单根所以特解必须乘 x。设 y₂* C x e^(-x)。求导y₂* C e^(-x)(1 - x)y₂* C e^(-x)(x - 2)代入 y 3y 2y提取 e^(-x)C[(x - 2) 3(1 - x) 2x] C[(x - 2 3 - 3x 2x)] C[1]方程右边 e^(-x) 的系数是 1所以 C 1。y₂* x e^(-x)。完了吗没有。注意 y₂* 实际上乘 x 后是 x e^(-x)而齐次通解里有 C₁ e^(-x)x e^(-x) 与 C₂ e^(-2x) 线性无关也没有撞到二重根所以不需要再乘更高次幂。两个子特解加起来y* (1/2)x - 3/4 x e^(-x)最终通解y C₁ e^(-x) C₂ e^(-2x) (1/2)x - 3/4 x e^(-x)这个例子值得反复体会。因为它展示了叠加原理下你只需要分别处理每个子项各自判断“撞不撞根”。第二步如果不乘 x 直接设 y₂* C e^(-x)代入左边得到 0右边 e^(-x) 永远凑不出来这是非齐次方程里最经典的错误。5.3 特解系数解不出来先检查这三件事在实际做题时一旦对比系数后方程出现矛盾比如求出 0 1不要急着怀疑“方法失效”按顺序排查特解形式是否设全了比如三角型只设了 cos 忘了 sin多项式漏了低次项。是否漏乘了 xλ 或 λ iω 是特征根却没乘 x必然导致代入后无论如何调节系数都无法匹配。有没有化简错误求导后合并同类项时指数函数和三角函数的系数特别容易算串建议每一步合并同类项后都回看一眼别急着往方程里代。这三条自查顺序几乎能定位掉九成以上的错误。剩下的一成往往是计算失误无他重算一遍就好。6. 常见错误与排查技巧实录6.1 问题速查表下面这张表是我学生时代和后来辅导别人时常常见到的错误汇总整理出来方便你对照检查常见错误现象解决办法右端多项式缺项特设多项式也跟着缺项对比系数时发现方程对不上特设多项式必须写完整m 次就要带着 m 到 0 全部项λ 是特征根却忘乘 x代入后左边为 0无法匹配右端解齐次特征根后先列出来再与 λ 比对λ iω 与特征根的关系判断错在共振情形漏乘 x先判断 r² pr q 0 是否有根等于 λ ± iω三角型特设只写 cos 或只写 sin对比另一项系数时始终矛盾特设必须同时包含 cos 和 sin即使右端只有一种求导过程不化简就代入计算量爆炸极易出错分步求导先合并同类项再代入把齐次通解的任意常数代入到特解最终结果“看似正确”但通解被破坏特解系数完全由方程确定不含任意常数每一个问题背后都对应不少“惨案”。比如漏乘 x 的错误在考试中属于典型的概念性失误一旦出现整道题分基本拿不到。所以做题之前先花 30 秒把特征根列清楚再动笔这点时间非常值。6.2 一个隐蔽性很强的细节λ 是特征根时特解多项式次数的抬升很多人记住了“单根乘 x二重根乘 x²”但忽略了乘完之后待定多项式的次数应当如何变化。我们回头看看 3.3 那题原本右端 e^x 对应的多项式次数 m 0设 Q₀(x) A。但因为二重根特解变成 x²·Ae^x。如果机械理解成“把 Q 的次数抬升”那 Q 的次数就从 0 变到 2即特解含 A x² e^x正好与正确结果一致。但如果你把“乘 x²”和“多项式次数抬升”两件事叠加起来处理比如右端本来是 e^x·xm 1且 λ 是二重根正确的特解应该设成 x² e^x (A x B) e^x(A x³ B x²)。有人会误设成 x² e^x (A x³ B x² C x D) 或者漏掉低次项这两种都是不稳定的。稳妥的理解方式永远是先写下基本形式 e^x Q_m(x)然后整体乘以 x^k其中 k 是特征根重数。乘完之后展开整理再取出未知系数不要先抬升多项式的次数再去乘。6.3 用初始条件确定任意常数的完整收尾做题时还经常要求满足初始条件。比如方程y 4y cos(2x)y(0) 1y(0) 0我们前面已经得到通解y C₁ cos(2x) C₂ sin(2x) (1/4) x sin(2x)代入 y(0) 1C₁ 1。求导y -2C₁ sin(2x) 2C₂ cos(2x) (1/4) sin(2x) (1/2) x cos(2x)代入 y(0) 02C₂ 0C₂ 0。所以满足初值条件的解是y cos(2x) (1/4) x sin(2x)注意这里出现了一个很有意义的细节特解中 x sin(2x) 这一项在 t 增大时会呈现“幅度线性增长”的振荡这就是共振现象在数学上的体现。系统以固有频率受到激励振幅随时间不断增大工程上要尽量避免否则结构可能被破坏。用初始条件确定常数时要注意的是必须在“齐次通解 特解”这个完整式子上代入而不是只在齐次通解里代入。因为特解不含任意常数本身也参与了初值的匹配。我见过不少人把特解丢在一边只用齐次通解去满足初值结果整条曲线从一开始就偏离正确解后面全错。7. 学习建议与两个高频问题的深度追问7.1 先练手还是先记表我的建议很多初学者看到特解表格第一反应是“背下来就好了”。真上了考场给的右端项往往和表格里的形式略有出入一紧张就不知道套哪个。我的建议是别先死记先做三道完整例题把设特解、代入、对比系数、解方程这套流程完整走一遍。有了手感之后再看表格那些条件是不是特征根、要不要乘 x、多项式次数多高就自然有了画面感记忆成本大幅降低。练习时特意选三个类型λ 不是特征根的、λ 是单根的、λ 是二重根的。把三题并排比较你会发现特解形式从 e^x 多项式到 x e^x 多项式再到 x² e^x 多项式变化规律一目了然。这就是“做三道题胜过背三页表”的学习逻辑。三角型同理分别练非共振和共振各一题就足够应付绝大多数考查场景了。7.2 学了待定系数法还需要常数变易法吗要。待定系数法虽然快但是对右端项形式有严格限制。反过来常数变易法可以处理任意右端项代价是计算过程包含积分且每一步都要带着齐次方程的线性无关解做运算。我的建议是优先掌握待定系数法因为考试和基础应用里绝大多数是非齐次项简单的情形用待定系数法一步到位、不容易出错。但必须理解常数变易法的推导思路至少知道它怎么由齐次通解中的常数“变易”而来。万一右端项里出现 e^(1/x)、ln x、分式这类非标准函数你还有退路。7.3 特征根是复数时特解该怎么设特征根是复数时齐次通解里的指数函数会转化成三角函数这在物理上对应振荡行为。非齐次项如果是三角函数实际上就是在判断 λ ± iω 是否等于特征根 α ± iβ。比较时需同时比较实部和虚部只有实部、虚部都相等才算“撞根”。举个例子y y sin(x)这里的 λ 0ω 1λ iω i。特征根是 ±i因此发生共振特解需要乘 x。如果你只是看“右端有三角函数”就直接设 y* A cos x B sin x代入后左边同样什么都凑不出来因为每一项都被齐次通解吸收了。这又一次说明解齐次方程这个前置动作永远不能省。8. 写在最后从形式到本质的一点体会我教这门课、自己复习备考的时候反复体会到一个道理微分方程求解的难点往往不在计算而在“判断”。你得先判断方程类型判断特征根判断共振条件判断特解形式——每一步判断都是一次对基本概念的检验。真正动手算的步骤反而机械只要细心就不会错。待定系数法的本质是用有限的函数类去匹配方程的结构。它之所以有效是因为常系数线性微分方程对某些函数类具有“封闭性”e^(λx)、多项式、正余弦函数它们经过线性组合和求导后还在这几类函数里。抓住这个本质你就不会在题目变化时发慌因为你始终知道特解长什么样、为什么要长这样。如果你正在备考或刚学到这一章我的建议就一句话不要赶进度把每一道例题都算到“结果代回去验证成立”为止。按我上面的思路一类一类地过先多项式指数型、再三角型、再做叠加型最后用初值条件收尾这一章的框架就能稳稳立起来。之后无论遇到什么组合、什么变化你都会发现核心方法始终是同一个设形式、代入、对比、解系数。这个小技巧送给你。