WJF87二阶线性非齐方程的解法

二阶常系数线性非齐次微分方程的解法高等数学·第八章第七节｜待定系数法与解的叠加原理Contents课程目录常微分方程求解方法的系统梳理与典型问题解析01方程基本形式与理论回顾02待定系数法核心原理03三类特解形式详解04应用举例与典型题解05总结与拓展思考CHAPTER01方程基本形式与理论回顾从齐次到非齐次，建立完整的理论认知框架HigherMathematics·ODE二阶常系数线性微分方程的一般形式二阶常系数线性微分方程的标准形式为y''+py'+qy=f(x)，其中p、q为常数。当f(x)≡0时为齐次方程，当f(x)≢0时为非齐次方程。齐次方程的通解是非齐次方程通解的基础组成部分。StandardFormy''+py'+qy=f(x)01标准形式y''+py'+qy=f(x)中，p和q均为实常数，f(x)为已知的连续函数，称为非齐次项或自由项。p,q∈ℝ02当f(x)≡0时退化为齐次方程y''+py'+qy=0，其特征方程为r²+pr+q=0，通过判别式Δ=p²−4q确定特征根类型。Δ=p²−4q03非齐次方程的通解结构为y=Y(x)+y*(x)，其中Y(x)是对应齐次方程的通解，y*(x)是非齐次方程的一个特解。y=Y(x)+y*(x)SECOND-ORDERODE齐次方程通解的三种形式回顾齐次方程y″+py′+qy=0的通解由特征方程r²+pr+q=0的根决定。根据判别式Δ的符号，特征根分为两个不等实根、两个相等实根和一对共轭复根三种情况，分别对应指数型、多项式乘指数型和三角函数乘指数型三种通解形式。Δ>0：两个不等实根Y=C₁er₁x+C₂er₂x特征根r₁≠r₂，两个线性无关的解均为纯指数函数，通解由两项指数叠加构成。r₁≠r₂Δ=0：两个相等实根Y=(C₁+C₂x)erx特征根r₁=r₂=r，第二个线性无关解需乘以x因子，构成多项式与指数的乘积形式。r₁=r₂Δ<0：一对共轭复根Y=eαx(C₁cosβx+C₂sinβx)特征根r=α±βi，通解体现振荡衰减或增长的物理特征，三角函数与指数相乘。α±βiSecond-OrderLinearODE非齐次方程的通解结构定理非齐次方程y''+py'+qy=f(x)的通解可分解为y=Y(x)+y*(x)，其中Y(x)是对应齐次方程的通解，y*(x)是一个特解。01定理表述：若Y(x)是y''+py'+qy=0的通解，y*(x)是y''+py'+qy=f(x)的一个特解，则y=Y+y*就是非齐次方程的通解。该结构揭示了非齐次方程解的叠加本质，将复杂问题分解为两个相对简单的子问题。02证明思路：将y=Y+y*代入非齐次方程左端，利用线性性质展开，Y满足齐次方程贡献为零，y*贡献恰好等于f(x)。验证过程充分运用了微分方程的线性叠加原理，体现了数学结构的内在和谐性。03解的结构含义：通解中Y(x)含两个任意常数C₁、C₂，体现解的自由度；y*(x)是确定的，由f(x)的形式决定。齐次解刻画解空间的维度特征，特解则反映非齐次项对系统的特定驱动作用。04求解策略转化：非齐次方程求解=求齐次通解（已掌握）+求非齐次特解（本节核心），关键在于根据f(x)的形式高效求特解。待定系数法与常数变易法是两类经典方法，适用于不同形式的非齐次激励函数。非齐次微分方程待定系数法适用的三类非齐次项形式三类非齐次项可统一归纳为f(x)=e^(αx)[Pₘ(x)cosβx+Qₙ(x)sinβx]的综合形式，其中α、β为实数，Pₘ和Qₙ为多项式。第一类：多项式f(x)=Pₘ(x)非齐次项为m次多项式，如f(x)=3x²+2x+1，可视为α=0、β=0的综合形式特例。α=0·β=0第二类：指数×多项式f(x)=e^(αx)Pₘ(x)，如f(x)=e^(2x)(x+1)，物理中常见的强迫振动激励形式。强迫振动第三类：指数×三角函数f(x)=e^(αx)[Pₘ(x)cosβx+Qₙ(x)sinβx]，如f(x)=eˣ·sin2x，对应交流电路周期激励。周期激励统一归纳：综合形式三类均表示为f(x)=e^(αx)[Pₘcosβx+Qₙsinβx]的特殊情况，α=0时退化为多项式或三角函数。统一形式Chapter02待定系数法核心原理从代入验证到形式确定，掌握特解的系统求解方法MethodofUndeterminedCoefficients待定系数法的基本思想待定系数法的核心思想是"以形猜解"：根据非齐次项f(x)的函数类型，预设特解y*具有相似的函数形式但含未知系数，将预设形式代入原方程后通过比较等式两端同类项的系数，建立关于未知系数的代数方程组，从而确定特解。从基本逻辑、操作步骤、方法优势与关键难点四个维度，理解这一经典代数化求解策略。01基本逻辑f(x)的形式决定特解y*的"模板"，只需确定模板中的未知系数，将微分方程问题转化为代数方程问题f→y*02操作步骤根据f(x)写出预设形式→代入原方程→比较同类项系数→解代数方程组得待定系数3步03方法优势避免积分运算的复杂性，将求特解过程完全代数化，计算规范可操作，适合工程快速求解代数化04关键难点正确预设特解形式是核心，当预设形式与齐次通解"冲突"时，需乘以x的适当幂次修正×xk常系数非齐次线性微分方程f(x)为多项式时的特解形式推导当f(x)=Pₘ(x)时，设特解为多项式代入方程y''+py'+qy=Pₘ(x)。由于求导会降低多项式次数，需要根据特征方程的根与0的关系来确定特解多项式的次数和是否需要乘以x的幂次。01代入与降次设特解y*=Q(x)代入方程后，qy*保持原次数，py*'降低1次，y*''降低2次。通过比较等式两边多项式次数，建立系数对应关系求解。SUBSTITUTE02q≠0·零非特征根最高次项由qy*决定，特解与Pₘ(x)同次，设为Qₘ(x)=a₀xᵐ+a₁xᵐ⁻¹+…+aₘ。此时0不是特征方程的根，无需乘以x的幂次。Qₘ(x)03q=0,p≠0·单根最高次项由py*'决定，特解需乘以x，设为y*=xQₘ(x)。此时0是特征方程的单根，乘以x保证特解形式的有效性。xQₘ(x)04q=0,p=0·重根最高次项由y*''决定，特解需乘以x²，设为y*=x²Qₘ(x)。此时0是特征方程的重根，乘以x²确保特解与齐次解线性无关。x²Qₘ(x)MethodofUndeterminedCoefficients待定系数法的核心判定规则核心规则取决于参数α与特征方程根的关系：不是特征根直接设，单根乘x，重根乘x²。特解形式确定规则表特征根条件修正因子特解预设形式α(或α+iβ)不是特征根x⁰=1直接按f(x)类型设同形式特解α(或α+iβ)是特征单根x¹=x在同形式特解前乘以xα(或α+iβ)是特征重根x²在同形式特解前乘以x²修正因子xk中k等于α（或α+iβ）作为特征根的重数（0、1或2）SUPERPOSITIONPRINCIPLE解的叠加原理非齐次项分解为子项之和时，各子特解之和即为原方程特解，可分别处理各组成部分。叠加原理表述若yᵢ*是y''+py'+qy=fᵢ(x)的特解，则Σyᵢ*是y''+py'+qy=Σfᵢ(x)的特解，这是线性方程独有的优良性质。Σyᵢ*应用场景当f(x)=3x²+e²ˣ+sin3x这样的混合形式时，可分别对多项式、指数和三角函数部分各求一个特解再相加。3x²+e²ˣ+sin3x注意事项各子特解的预设形式需独立确定，每个子特解都要单独判断是否需要乘以x的修正因子，不能遗漏。×xᵏ修正与齐次通解的关系最终通解仍为y=Y(x)+Σyᵢ*，叠加原理只影响特解部分，齐次通解Y(x)不受分解影响。Y(x)+Σyᵢ*AdvancedMethod复数方法与辅助方程技巧利用欧拉公式e^(iβx)=cosβx+isinβx，将含三角函数的非齐次项转化为复指数形式，构造辅助方程求解后取实部或虚部即得原方程特解。01转化原理当f(x)=Pₘ(x)cosβx时，利用cosβx=Re[e^(iβx)]，构造辅助方程z''+pz'+qz=Pₘ(x)e^(iβx)，求得其特解z*后取实部即为所求。Re[z*]→特解02虚部取法当f(x)=Pₙ(x)sinβx时，利用sinβx=Im[e^(iβx)]，同样构造辅助方程，求得特解z*后取虚部即为所求。Im[z*]→特解03计算优势复指数函数的求导规则简单——每求一次导乘以iβ，避免了三角函数反复使用积化和差公式的繁琐计算。d/dx→×iβ04适用前提原方程系数p、q必须为实数，这样实部和虚部才能独立满足各自方程，保证取实部或虚部后仍为原方程的解。p,q∈ℝChapter03三类特解形式详解多项式型、指数多项式型、三角函数型的系统求解二阶常系数非齐次线性微分方程第一类特解：f(x)=Pₘ(x)·eαx对于y″+py′+qy=Pₘ(x)·eαx，设特解y*=xk·Qₘ(x)·eαx，k取决于α是否为特征根CASE01α不是特征根特解设为y*=Qₘ(x)·eαx，Qₘ(x)含m+1个待定系数，代入方程比较系数后可逐一确定k=0CASE02α是特征单根特解设为y*=x·Qₘ(x)·eαx，乘以x是因为eαx已是齐次解的一部分，直接设会与齐次通解线性相关k=1CASE03α是特征重根特解设为y*=x²·Qₘ(x)·eαx，此时eαx和xeαx都是齐次解，需乘以x²才能保证线性无关k=2常系数非齐次线性微分方程第一类特解求解示例以y''−5y'+6y=xe2x为例，特征根r₁=2、r₂=3，α=2是特征单根，故k=1，特解设为y*=x(ax+b)e2x。代入方程比较系数后确定a=−1/2、b=−1，最终特解为y*=(−x²/2−x)e2x。01求特征根，判断根型解特征方程r²−5r+6=0得r₁=2,r₂=3。判断α=2是特征单根，确定修正因子k=1。02设定特解形式设y*=x(ax+b)e2x=(ax²+bx)e2x，其中Pₘ(x)=x为一次多项式，Qₘ(x)=ax+b也为一次。03代入方程，比较系数计算y*'和y*''，代入原方程，消去公因子e2x后整理得−2ax+(2a−b)=x，比较系数得a=−1/2,b=−1。04写出完整通解特解y*=(−x²/2−x)e2x，完整通解为y=C₁e2x+C₂e3x+(−x²/2−x)e2x。Non-Homogeneous·Trigonometric第二类特解：含三角函数的非齐次项对于f(x)=e^(αx)[Pₗ(x)cosβx+Pₙ(x)sinβx]，设特解y*=xᵏ·e^(αx)[Rₘ⁽¹⁾(x)cosβx+Rₘ⁽²⁾(x)sinβx]，其中m=max(l,n)，k取决于α+iβ是否为特征根。即使f(x)只含cos或sin，特解中cos和sin两项都必须设置。01特解形式y*=xᵏ·e^(αx)[Rₘ⁽¹⁾(x)cosβx+Rₘ⁽²⁾(x)sinβx]，两个多项式的次数均为m=max(l,n)02k的判定看α+iβ是否为特征方程的根：不是特征根则k=0，是特征根则k=1（复根成对出现，不可能为重根）03易错提醒即使f(x)只含cos没有sin，特解中仍需同时设cos和sin两项，因为求导会使cos和sin互相转化04系数个数两个m次多项式各有m+1个系数，共2(m+1)个待定系数，需比较cosβx和sinβx各项系数分别建立方程常系数非齐次线性微分方程第二类特解求解示例以y''+y=cosx为例，特征根为±i，而α+iβ=i是特征根，故k=1。设特解y*=x(acosx+bsinx)，代入后比较系数得a=0、b=1/2，特解为y*=(x/2)sinx。01特征分析：r²+1=0的根为r=±i，非齐次项cosx对应α=0、β=1，α+iβ=i是特征根，确定k=102设特解形式：y*=x(a·cosx+b·sinx)，计算y*'=acosx+bsinx+x(−asinx+bcosx)，再求y*''03代入比较系数：y*''+y*=cosx化简得2bcosx−2asinx=cosx，比较得2b=1→b=1/2，−2a=0→a=004最终结论：特解y*=(x/2)sinx，通解y=C₁cosx+C₂sinx+(x/2)sinx，特解中含xsinx项体现了共振现象常系数线性微分方程·特解方法第三类特解：纯三角函数型f(x)=Acosβx或Asinβx纯三角函数型是第二类中α=0、多项式为零次的特殊情况。设特解y*=xᵏ(acosβx+bsinβx)，k取决于±βi是否为特征根。当±βi是特征根时出现共振，特解中出现x因子，振幅随时间线性增长。基本形式f(x)=Acosβx或Asinβx时，设特解y*=xᵏ(acosβx+bsinβx)，a、b为待定系数y*=xᵏ(a·c+b·s)k的判定简化看±βi是否为特征根，将−β²代入特征方程检验：若−β²不是特征根则k=0，是则k=1k∈{0,1}共振现象当βi恰为特征根时（如弹簧系统固有频率等于外力频率），特解含x因子，物理上表现为振幅不断增大的共振振幅∝x实用技巧：复数方法构造z''+pz'+qz=Ae^(iβx)，求得z*后取实部或虚部，避免直接处理三角函数运算z*→Re/ImReferenceTable三类特解形式完整总结将待定系数法的三类特解形式统一汇总，根据非齐次项f(x)的类型和参数与特征根的关系，可快速确定特解的预设形式。特解形式速查表f(x)的类型非特征根(k=0)是特征根(k=1或2)Pₘ(x)·e^(αx)Qₘ(x)·e^(αx)xᵏ·Qₘ(x)·e^(αx)（k=1单根，k=2重根）e^(αx)[Pₗcosβx+Pₙsinβx]e^(αx)[Rₘ⁽¹⁾cosβx+Rₘ⁽²⁾sinβx]x·e^(αx)[Rₘ⁽¹⁾cosβx+Rₘ⁽²⁾sinβx]Acosβx或Asinβxacosβx+bsinβxx(acosβx+bsinβx)Pₘ(x)（纯多项式）Qₘ(x)（α=0非特征根）xᵏ·Qₘ(x)（α=0为特征根，k=1或2）m=max(l,n)，Qₘ、Rₘ⁽¹⁾、Rₘ⁽²⁾均为m次待定多项式；k等于相应参数作为特征根的重数CHAPTER04应用举例与典型题解从基础到综合，通过典型例题掌握完整求解流程待定系数法·二阶常系数非齐次例1：纯多项式型f(x)=3x+1求解y''-2y'-3y=3x+1。特征根为r₁=-1、r₂=3，α=0不是特征根，设特解y*=ax+b，代入比较系数得a=-1、b=1/3，特解为y*=-x+1/3，通解为y=C₁e⁻ˣ+C₂e³ˣ-x+1/3。01求特征根，判定k值特征方程r²-2r-3=0分解为(r+1)(r-3)=0，得r₁=-1,r₂=3；α=0不是特征根，k=002设特解形式并求导设特解y*=ax+b，则y*'=a,y*''=0，代入原方程得0-2a-3(ax+b)=3x+103比较系数，解出待定常数整理得-3ax+(-2a-3b)=3x+1，比较系数：-3a=3→a=-1；-2(-1)-3b=1→b=1/304合成最终通解齐次通解Y=C₁e⁻ˣ+C₂e³ˣ，最终通解y=C₁e⁻ˣ+C₂e³ˣ-x+1/3SECOND-ORDERODE例2：α为特征重根的情况求解y''-4y'+4y=(1+x)e^(2x)。特征方程r²-4r+4=0有重根r=2，α=2是二重特征根，故k=2。设特解y*=x²(ax+b)e^(2x)，代入化简后得a=1/6、b=1/2。01确定特征根与重数特征方程r²-4r+4=(r-2)²=0，r=2是二重特征根；α=2是特征重根，确定k=202设定特解形式设特解y*=x²(ax+b)e^(2x)=(ax³+bx²)e^(2x)，Pₘ(x)=1+x为一次多项式，Qₘ(x)=ax+b也为一次03代入化简利用公式简化：对y=xᵏQ(x)e^(αx)代入后可利用特征方程性质消去多项项，整理得6ax+2b=1+x04比较系数求通解比较系数：6a=1→a=1/6，2b=1→b=1/2；特解y*=x²(x/6+1/2)e^(2x)，通解y=(C₁+C₂x)e^(2x)+x²(x/6+1/2)e^(2x)HIGHERMATHEMATICS·二阶常系数例3：复数方法求解含三角函数的特解求解y''+4y=sin2x。特征根为±2i，α+iβ=2i是特征根，k=1。构造辅助方程z''+4z=e^(2ix)，设z*=axe^(2ix)，代入求得a=-i/4，取虚部得特解y*=-(x/4)cos2x。01特征根判定特征方程r²+4=0得r=±2i；非齐次项sin2x对应α=0、β=2，α+iβ=2i是特征根，k=102构造辅助方程构造辅助方程z''+4z=e^(2ix)，设z*=axe^(2ix)（因2i是特征单根，乘以x），代入得2a·2ie^(2ix)=e^(2ix)03求解复数系数解得a=1/(4i)=-i/4，故z*=(-i/4)xe^(2ix)=(-ix/4)(cos2x+isin2x)=(x/4)sin2x-(ix/4)cos2x04取虚部得特解取虚部（因sin2x=Im[e^(2ix)]）：y*=Im[z*]=-(x/4)cos2x，通解y=C₁cos2x+C₂sin2x-(x/4)cos2xEXAMPLE04例4：叠加原理处理复合非齐次项当f(x)=f₁(x)+f₂(x)时，分别求出对应f₁和f₂的特解y₁*和y₂*，则y₁*+y₂*即为原方程的特解。叠加原理将复合非齐次项分解为简单项分别处理，是解决复杂问题的重要策略。01题设分析y''−3y'+2y=ex+cosx，f(x)由指数项ex和三角项cosx两部分组成，需分别处理02处理指数项对y''−3y'+2y=ex，特征根r=1,2，α=1是单根，设y₁*=axex，代入得a=−1，y₁*=−xex03处理三角项对y''−3y'+2y=cosx，α+iβ=i不是特征根，设y₂*=a·cosx+b·sinx，代入比较系数确定a、b04叠加得通解原方程特解y*=y₁*+y₂*=−xex+a·cosx+b·sinx，通解y=C₁ex+C₂e2x+y*Methodology例5：综合应用与解题策略综合题求解的标准流程：第一步解特征方程确定特征根，第二步分析f(x)类型和参数，第三步查表确定特解形式，第四步代入比较系数，第五步写出完整通解。严格按流程操作可避免遗漏和错误。Phase1·标准化求解五步法01解特征方程r²+pr+q=0得特征根r₁、r₂，明确齐次通解Y(x)的形式02识别f(x)类型（多项式/指数/三角/组合），提取关键参数α、β、多项式次数mPhase2·关键检验与最终求解03判断α（或α+iβ）是否为特征根及重数，确定修正因子k，写出完整特解预设形式04计算y*′、y*″，代入原方程，按同类项分组比较系数，解代数方程组确定所有待定系数05组合最终结果：y=Y(x)+y*(x)，如有初值条件则代入确定C₁、C₂MISTAKESTOAVOID常见错误与注意事项待定系数法求解中的典型错误包括：忘记检验特征根、漏设三角函数项、比较系数遗漏、混淆特解与通解忘记检验特征根未检验α或α+iβ是否为特征根，直接设特解导致与齐次通解线性相关，结果恒为零无法匹配f(x)。ROOTCHECK漏设三角函数项f(x)只含cosβx时漏设sinβx项。求导使cos和sin互相转化，两项必须同时设置。SIN+COS比较系数遗漏代入方程后遗漏某些次数的同类项，或混淆了e^(αx)因子前后的系数对应关系。MATCHALL混淆特解与通解将特解y*误认为最终答案，忘记加上齐次通解Y(x)。y=Y+y*才是完整结果。y=Y+y*InitialValueProblem初值问题的完整求解流程初值问题要求在通解的基础上，利用初始条件确定任意常数，从而得到唯一的确定解，对应给定初始状态的系统响应。01求解流程：先求通解y=C₁y₁(x)+C₂y₂(x)+y*(x)，再对通解求导得y'(x)，代入x₀处的两个初始条件。02建立方程组：将y(x₀)=y₀和y'(x₀)=y₀'分别代入，得到关于C₁、C₂的二元一次方程组，求解即可。03物理意义：齐次通解对应系统自由振动（由初始条件激发），非齐次特解对应强迫振动（由外力驱动），叠加为完整响应。04验算建议：将最终结果代入原方程验证等式成立，同时验证初始条件是否满足，双重检验可有效避免计算错误。CHAPTER05总结与拓展思考梳理核心方法体系，拓展思维深度Summary知识体系与核心方