ARTICLE DETAIL

资讯详情

深耕郑州网站建设与运营推广的一线实战洞察。

常系数非齐次线性微分方程求解:待定系数法特解与共振判断全解析

常系数非齐次线性微分方程求解:待定系数法特解与共振判断全解析 一拿到形如 y ay by f(x) 的常系数非齐次线性微分方程很多人第一反应是赶紧套待定系数法设特解。这个方向是对的但真正动手前有一步比求导更重要——先看清楚 f(x) 长什么样。待定系数法能处理的非齐次项其实就两种标准形式一种是多项式乘指数型 f(x) P_n(x)e^{αx}另一种是带三角函数的振荡型 f(x) e^{αx}[P_n(x)cos βx Q_m(x)sin βx]。这两个形式覆盖了高数、考研和工科专业课里九成以上的题目剩下的一成才需要动用常数变易法。我教过很多学生发现大家记特解表都能背下来真正的分水岭出现在三个判断题上α 到底撞不撞特征根cos 和 sin 为什么要成对设乘的 x^k 里的 k 到底取几这篇文章就把这两种常见形式完整拆开从判断流程到完整例题再到我改作业时反复见到的五个错误一次性讲透。适合正在学微分方程的本科生也适合考研数学基础阶段或者准备补考的同学。1. 先搭框架非齐次方程的通解为什么不是“猜”出来的1.1 通解 齐次通解 非齐次特解常系数非齐次线性微分方程的标准写法是y ay by f(x)其中 a、b 是常数f(x) 不等于零等于零就是齐次方程。为什么书上总说“先解对应的齐次方程”再去求一个特解因为线性算子 L[y] y ay by 满足叠加性如果 L[y_c] 0L[y_p] f(x)那么 L[y_c y_p] 0 f(x) f(x)所以 y_c y_p 一定是原方程的解。反过来看如果 y₁ 和 y₂ 都是同一个非齐次方程的解那么它们的差满足 L[y₁ - y₂] f(x) - f(x) 0也就是说差是齐次方程的解。这句话的潜台词是任意两个特解之间只差一个齐次通解所以非齐次方程的通解一定是“齐次通解 随便一个特解”不多不少。这里有个很重要的认知通解里两个任意常数只藏在齐次通解 y_c 里特解 y_p 是一个没有任何自由度的具体函数。很多学生最后写答案时会把特解里的系数 A、B 错当常数保留然后和 y_c 的 C₁、C₂ 混在一起这属于框架没搭好后面我会专门说。1.2 齐次部分的三句话复习既然万变不离齐次通解特征方程必须先过关。对 y ay by 0令 y e^{rx}代入后得到特征方程 r² ar b 0按判别式分三种情况特征根情况齐次通解 y_c两个不等实根 r₁ ≠ r₂y_c C₁e^{r₁x} C₂e^{r₂x}二重实根 r₁ r₂ ry_c (C₁ C₂x)e^{rx}共轭复根 α ± iβy_c e^{αx}(C₁cos βx C₂sin βx)我上课时反复强调先写 y_c再想 y_p。因为判断“共振”时需要对比特征根而共振正是这两类形式里最容易出错的地方。1.3 为什么待定系数法只认这两类 f(x)待定系数法的核心逻辑是“猜一个和 f(x) 同形状的函数”。凭什么敢这么猜因为常系数线性微分方程求导不会改变函数的“家族”多项式乘 e^{αx} 求导后还是多项式乘 e^{αx}cos、sin 求导后只在二者之间互相转换。这类函数在求导运算下是封闭的所以设出带未知系数的同形状函数代回方程比较系数一定能解出这些未知系数。反过来如果 f(x) 是 ln x、tan x、1/x 这种求导后“变种”的函数再强行设同形状函数就行不通了因为求导之后会冒出新的函数形态系数根本没法比。这时候就要换常数变易法。所以你拿到的题目如果非齐次项是这两种标准形式的有限和待定系数法是首选不是的话果断换方法别硬猜。2. 第一种形式f(x)P_n(x)e^{αx}设特解的关键是“三个 k 的判定”2.1 判定流程特征根与 α 的关系第一种标准形式是 f(x) P_n(x)e^{αx}其中 P_n(x) 是 n 次多项式。注意两个退化情形也包含在这里α 0 时退化成纯多项式 f(x) P_n(x)n 0 时退化成纯指数 f(x) A e^{αx}。设特解的标准格式是y_p x^k Q_n(x) e^{αx}其中 Q_n(x) 是系数待定的 n 次多项式必须把从 x^n 到常数项的所有 n1 个系数全部设出来。k 的取值由 α 和特征根的关系决定条件k 取值α 不是特征根k 0α 是单特征根k 1α 是二重特征根k 2比如后面要讲的例题 y - 3y 2y 6x e^x特征根是 1 和 2α 1 恰是单根所以 k 1特解必须设成 y_p x(Ax B)e^x而不是 (Ax B)e^x。2.2 完整例题y - 3y 2y 6x e^x先解齐次特征方程 r² - 3r 2 (r - 1)(r - 2) 0所以 y_c C₁e^x C₂e^{2x}。再看非齐次项 f(x) 6x e^x属于第一种形式P₁(x) 6x 是一次多项式α 1。因为 1 是单特征根所以设y_p x(Ax B)e^x (Ax² Bx)e^x求导要仔细建议分步来y_p (2Ax B)e^x (Ax² Bx)e^x [Ax² (2A B)x B]e^xy_p [2Ax (2A B)]e^x [Ax² (2A B)x B]e^x [Ax² (4A B)x (2A 2B)]e^x代回原方程注意先合并每一项 e^x 前的多项式L[y_p] [Ax² (4A B)x (2A 2B)]e^x - 3[Ax² (2A B)x B]e^x 2(Ax² Bx)e^xx² 的系数A - 3A 2A 0自动抵消了这就是共振的特征。x 的系数4A B - 6A - 3B 2B -2A。常数项2A 2B - 3B 2A - B。所以 L[y_p] (-2Ax 2A - B)e^x令它等于 6x e^x比较系数得-2A 62A - B 0解得 A -3B -6。于是 y_p (-3x² - 6x)e^x通解为y C₁e^x C₂e^{2x} (-3x² - 6x)e^x这道题的完整流程必须走一遍因为它同时包含了“判断共振”“设含 x 的特解”“求导合并比较系数”三个核心动作。2.3 两个退化情形纯指数和纯多项式退化情形一纯指数 f(x) 3e^x方程为 y - y 3e^x。特征根是 ±1α 1 是单根所以设 y_p Ax e^x。这里很多人会设 y_p Ae^x代进去会发现 L[y_p] 0永远凑不出 3e^x原因就是 Ae^x 本来就在齐次通解里A 无论取什么都无法产生非零输出。按正确设法y_p A(x 1)e^xy_p A(x 2)e^x于是y_p - y_p A(x 2)e^x - Ax e^x 2Ae^x令 2A 3得 A 3/2。通解为 y C₁e^x C₂e^{-x} (3/2)x e^x。退化情形二纯多项式 f(x) x²方程为 y - y x²。这时候 α 0。特征根是 ±10 不是特征根所以 k 0直接设 y_p Ax² Bx C。有个初学者常犯的错误看到右边最高次是 x²设的时候只写 Ax²漏掉低次项。实际上 Q_n(x) 必须从 x^n 一直到常数项全部写齐因为你求导后低次项会产生新的低次项必须靠设出来的系数去平衡。y_p 2A代入得2A - (Ax² Bx C) -Ax² - Bx (2A - C)令它等于 x²比较系数-A 1 → A -1B 02A - C 0 → C -2。所以 y_p -x² - 2通解为 y C₁e^x C₂e^{-x} - x² - 2。2.4 为什么要乘 x^k一个公式看懂很多同学背下了表格但不知道 x^k 是哪来的。这里给一个能彻底看明白的推导。对 y_p Q(x)e^{αx}其中 Q(x) 是多项式求导规律是y_p (Q αQ)e^{αx}y_p (Q 2αQ α²Q)e^{αx}所以L[y_p] [Q (2α a)Q (α² aα b)Q]e^{αx}注意 α² aα b 就是特征多项式在 α 处的取值 P(α)。如果 deg Q n当 α 不是特征根时P(α) ≠ 0L[y_p] 里保留 n 次项正好能匹配 f(x) 的 n 次多项式。当 α 是单根时P(α) 0最高项被吃掉只剩下 Q (2α a)Q 提供的 n-1 次项。这时候如果还设 Q 为 n 次多项式输出最多是 n-1 次永远配不上 n 次。解决办法就是让 Q 的次数变成 n1等价于在原猜测前面乘一个 x。当 α 是二重根时连 2α a 也等于 0只剩下 Q 提供的 n-2 次项所以必须乘 x²。这就是 x^k 的由来它不是魔法也不是老师拍脑袋定的规则而是求导降次之后的必然补偿。搞懂这一点以后遇到再诡异的共振情形你都能自己推出来。3. 第二种形式f(x)e^{αx}[P_n(x)cos βx Q_m(x)sin βx]两个必须记住的规矩3.1 设特解的表格口诀第二种标准形式是f(x) e^{αx}[P_n(x)cos βx Q_m(x)sin βx]其中 P_n 和 Q_m 分别是 n 次和 m 次多项式。设特解的格式是y_p x^k e^{αx}[R_l(x)cos βx S_l(x)sin βx]这里有两个规矩缺一个都会翻车。规矩一l 取 n 和 m 的较大值也就是 l max(n, m)而且 R_l 和 S_l 都必须设成完整的 l 次多项式。规矩二判断 k 的参照对象不是 α而是复数 s α iβ。只要 s 是特征方程的根k 就取 1。对二阶常系数方程来说如果 α ± iβ 是共轭复根那 s 只能算单根所以 k 最多取到 1不会出现 k 2 的情形。这跟在第一种形式里 α 撞了单特征根就乘一个 x 是同一个道理——cos βx 和 sin βx 本质上是复指数 e^{(αiβ)x} 的实部和虚部。条件k 取值α iβ 不是特征根k 0α iβ 是特征根二阶时必为单根k 1两个特例顺手记一下α 0 时 e^{αx} 消失变成纯三角函数型β 0 时 sin 项消失、cos 项变成 1直接退化成第一种形式。3.2 完整例题y 4y sin 2x这道题是“共振”的经典教材题。齐次方程 y 4y 0 的特征根是 ±2i所以 y_c C₁cos 2x C₂sin 2x。也就是说sin 2x 本身就在齐次通解里躺着。如果你不管共振直接设 y_p A cos 2x B sin 2x代进去会发现 L[y_p] 0A、B 解不出来。这就是共振的信号必须乘 x。因为 α 0β 2α iβ 2i 是特征根所以设y_p x(A cos 2x B sin 2x)求导过程放这里你以后可以对照检查y_p A cos 2x B sin 2x x(-2A sin 2x 2B cos 2x) (A 2Bx)cos 2x (B - 2Ax)sin 2xy_p 2B cos 2x - 2(A 2Bx)sin 2x - 2A sin 2x 2(B - 2Ax)cos 2x (4B - 4Ax)cos 2x (-4A - 4Bx)sin 2x代回 y 4yL[y_p] (4B - 4Ax)cos 2x (-4A - 4Bx)sin 2x 4x(A cos 2x B sin 2x) 4B cos 2x - 4A sin 2x注意 x 的项全部抵消了只剩下两个纯三角函数。令它等于 sin 2x比较系数得4B 0-4A 1于是 A -1/4B 0。y_p -(x/4)cos 2x通解y C₁cos 2x C₂sin 2x - (x/4)cos 2x这道题在物理里的对应场景就是无阻尼系统受到与固有频率相同的外力激励振幅随时间线性增长那个 x 因子就是共振累积效应的数学体现。3.3 完整例题y - 2y 2y e^x cos x这道题的特点是特征根本身带有实部对应的齐次解本身就含 e^x 因子和右边的 e^x cos x 撞了个满怀。先看齐次特征方程 r² - 2r 2 0判别式 Δ 4 - 8 -4根为 r 1 ± i。所以 y_c e^x(C₁cos x C₂sin x)。右边的 f(x) e^x cos x 属于第二种形式α 1β 1n m 0。因为 α iβ 1 i 正是特征根所以 k 1设y_p x e^x(A cos x B sin x)这里我强烈推荐一个技巧先把 e^x 提出来设 y_p e^x w(x)其中 w(x) x(A cos x B sin x)。这样能用现成的求导公式如果 y e^x w那么y e^x(w w) y e^x(w 2w w)代进原方程L[y_p] y - 2y 2y e^x[(w 2w w) - 2(w w) 2w] e^x(w w)这个化简结果非常漂亮比硬算 y_p 省一半功夫。现在只要算 w ww Ax cos x Bx sin x w (A Bx)cos x (B - Ax)sin x w (2B - Ax)cos x (-2A - Bx)sin x于是w w (2B - Ax Ax)cos x (-2A - Bx Bx)sin x 2B cos x - 2A sin x所以 L[y_p] e^x(2B cos x - 2A sin x)令它等于 e^x cos x得2B 1-2A 0解得 B 1/2A 0。于是y_p (x/2)e^x sin x通解y e^x(C₁cos x C₂sin x) (x/2)e^x sin x也可以合并写成 y e^x[C₁cos x (C₂ x/2)sin x]但要注意 C₂ 和 x/2 不能真的合并成一个常数因为 x/2 不是常数最好保留分开的形式避免后续代初值条件时出错。3.4 最省事的计算技巧先提 e^{αx}上面那道题展示了“先提 e^{αx}”的巨大优势这里把它推广成通用公式。对一般的方程 y ay by f(x)如果设 y_p e^{αx}u(x)则y e^{αx}(u αu) y e^{αx}(u 2αu α²u)于是L[y_p] e^{αx}[u (2α a)u (α² aα b)u]只要 f(x) 也含 e^{αx} 因子两边同时约掉 e^{αx}剩下的只是一个关于 u 的、没有指数因子的方程全是多项式或三角函数的加减乘除计算量直接降一个档次。我给学生讲这两种形式时第一遍要求按最朴素的办法硬算让他们体会结构第二遍就要求用这个提公因子法提速。熟练之后解第二种形式基本不会在求导环节丢分。4. 五个高频错误每个都是真实改作业改出来的4.1 漏乘 x^k这个错误出现的频率最高。典型场景是 y 4y sin 2x 这种“右边已经在齐次通解里”的题还有 y - y 3e^x 这种“指数撞特征根”的题。漏乘 x^k 的直接后果是代回方程后左边恒等于 0或者出现根本无法匹配的多项式。判断方法就一句话先写齐次通解再对着 f(x) 看只要 f(x) 里的某个函数形式已经在 y_c 里出现过e^{αx} 撞根、sin/cos 撞复根、多项式里 0 撞根就必须乘 x。撞重根就乘 x²。4.2 三角函数只设一半典型错误右边是 cos x特解只设 y_p A cos x不写 B sin x。为什么必须成对设因为对 cos x 求导会产出 sin x对 sin x 求导又会变回 cos x代入二阶方程后左右两边必然同时出现 sin 和 cos只有两个系数都能调才能同时消掉两边同名的三角函数。右边的 sin 项就算系数为零左边求导过程中也会冒出来所以 sin 前面的多项式照样要设。4.3 多项式次数取错第二种形式里R_l 和 S_l 的次数统一取 max(n, m)不是各取各的。比如右边是 x cos x sin x也就是 n 1、m 0特解必须设成y_p (Ax B)cos x (Cx D)sin x而不是设成 (Ax B)cos x C sin x。原因是求导后 x cos x 会产生 x sin x 项你就需要 Cx D 这种完整的一次式去匹配。反过来如果只取较高次数而有顾虑记住宁可多设一个低次项也绝不能少设。多设的系数最后会算出来是 0不影响结果少设了就直接解不出来。4.4 忘了叠加原理当 f(x) 是多种标准类型的和时比如 f(x) x sin 2x正确的做法是把非齐次项拆开每一部分单独设特解然后把所有特解加起来。这个原理在线性系统里就是“各自响应叠加”。有人会把 x 和 sin 2x 强行塞进同一个猜测式里不是不行而是容易乱。更稳的流程是先拆开、分别待定、最后相加。叠出来的 y_p 可能又杂又长但这是通解的正确结构。4.5 比较系数和最终通解的细节错误比较系数之前必须先把同类项合并干净。把导数算完之后e^x 前面的多项式、cos 前面的多项式、sin 前面的多项式都要各自合并彻底再逐项对比。我见过太多人看着像“两行一样的式子”就急着抄答案结果一个符号错整道题全崩。另外一个细节是通解的写法。特解里求出来的 A、B 都是具体数字不能再保留成未知常数任意常数只能来自 y_c。写最终答案时如果 f(x) x sin 2x 这种多部分叠加特解是好几段加起来的检查一下每一段是否都不含 C₁、C₂。如果给了初值条件代初值时必须代入完整的 y也就是 y_c y_p不能只代入 y_c。5. 综合题拆解叠加原理 二重根共振同时出现怎么办5.1 拆解流程综合题的套路其实不会超出前面讲的范围只是把多个考点叠在一起。拿到题目按这四步走求特征根写出齐次通解 y_c。把 f(x) 拆成若干个标准形式之和。对每一部分分别判断 k设出各自的 y_{p,i}。代入原方程求系数可以逐个代入也可以合在一起代入最后把所有特解相加。关键在于第三步每个部分独立判断共振不能因为某一部分撞根就让所有部分都乘 x也不能因为某一部分不撞根就让所有部分都不乘 x。这里最容易出错的就是“惯性思维”。5.2 完整例题y - 2y y e^x x这道题把前面所有要点串起来了建议自己先做再看解答。第一步齐次。特征方程 r² - 2r 1 (r - 1)² 0r 1 是二重根所以y_c (C₁ C₂x)e^x第二步拆 f(x)。f(x) e^x x两项分开看。对 e^xα 1而 1 是二重特征根所以 k 2。注意这里不能只乘一个 x必须设y_{p1} Ax²e^x如果只设 Ax e^x代进去会发现左边仍然部分为 0这就是特征根重数没有尊重到位。对 x这是第一种形式的退化情形α 0而 0 不是特征根所以 k 0设y_{p2} Bx C第三步分别代入。先看 y_{p1}y_{p1} (2Ax Ax²)e^x y_{p1} (2A 4Ax Ax²)e^x代入L[y_{p1}] y_{p1} - 2y_{p1} y_{p1} [(2A 4Ax Ax²) - 2(2Ax Ax²) Ax²]e^x 2Ae^x所以 2A 1A 1/2。再看 y_{p2} Bx Cy By 0L[y_{p2}] 0 - 2B Bx C Bx (C - 2B)令它等于 x得 B 1C - 2B 0 → C 2。第四步合并。通解为y (C₁ C₂x)e^x (1/2)x²e^x x 2这道题每次讲完都有学生倒吸一口气原来二重根和 e^x 撞在一起要乘 x²而 x 尽管也作为多项式出现却因为 0 不是根而完全不需要乘 x。同一个方程里的两个非齐次项待遇天差地别这正是这类题的精髓。5.3 自查三件事算完一道题别急着合上习题集花三十秒做三个检查。第一特解里有没有残留 A、B、C 等待定字母有就说明没算完。第二把 y_p 大概代回原方程扫一眼确保代入后两边首项符号对得上。不求精确但至少不出现“左边是多项式右边是 e^x”这种离谱错位。第三对比一下 y_p 和 y_c 的每一项确认不是简单重复。如果 y_p 里的某个函数形式和 y_c 完全一样却没有乘 x那基本可以断定共振漏判了。这三个检查用熟了以后你的正确率会比闷头刷题高很多。最后再分享一个我带学生时的习惯性做法准备一张卡片正面写“先写 y_c再判共振”背面写两种形式的特解表格。每次做题前先看卡片做完题再回看卡片自查一遍。其实微分方程这门课计算量大是一回事真正的门槛在于有没有一个固定的“动手前检查清单”。等你形成肌肉记忆看到 f(x) 自动就会分类、自动想到 k 取几这两种标准形式就再也不会成为你的扣分项了。
返回列表