D78常系数非齐次线性微分方程

常系数非齐次线性微分方程高等数学·微分方程章节|第七节Contents本节内容概览常系数非齐次线性微分方程的求解方法与理论基础01基本概念与解的结构02f(x)=e^(λx)Pm(x)型特解求法03三角函数型f(x)特解求法04n阶方程推广与实际应用05方法总结与综合练习CHAPTER01基本概念与解的结构明确方程形式，理解通解的组成结构Definition·标准形式二阶常系数非齐次线性微分方程二阶常系数非齐次线性微分方程的标准形式为y''+py'+qy=f(x)，其中p、q为常数，f(x)为非零已知函数。01方程一般形式y''+py'+qy=f(x)，p、q为常数，f(x)为非零连续函数，称为自由项或驱动项。f(x)≠002常系数含义y及其各阶导数的系数均为常数，而非常量函数，具有特殊求解方法，可用特征方程法处理。CONSTANT03非齐次判定f(x)≠0为非齐次；f(x)≡0则退化为对应齐次方程y''+py'+qy=0，两者解结构密切相关。NON-HOMO04求解目标在齐次通解基础上寻找特解y*，两者相加得非齐次方程通解，这是叠加原理的核心应用。y=Y+y*常系数非齐次线性微分方程非齐次方程解的结构定理非齐次方程通解=齐次通解Y+特解y*，求解由此转化为先求齐次通解、再求一个特解。通解结构通解公式非齐次方程通解y=Y+y*，Y为齐次方程y''+py'+qy=0的通解，y*为非齐次方程的任一特解y=Y+y*通解结构解的性质齐次通解Y含两个任意常数C₁、C₂，体现解的完整性；特解y*不含任意常数，是某一具体函数C₁·C₂定理意义化繁为简将复杂问题分解：齐次通解已在上节解决，本节只需专注于求特解y*的方法求特解y*定理意义特解不唯一只要找到满足方程的任何一个特解即可，不同特解之间相差一个齐次解y₁*−y₂*=YMethodofUndeterminedCoefficients待定系数法的基本思想待定系数法是一种无需积分即可求特解的代数方法，通过推测特解结构、代入方程、比较系数三步完成求解。方法核心根据f(x)的形式预先假设特解y*的结构形式，其中包含若干待确定的常数系数。假设结构求解步骤将假设的特解y*及其各阶导数代入原方程，整理后比较等式两边同类项的系数。代入比较建立方程组通过系数比较得到关于待定系数的线性方程组，解方程组即可确定所有待定系数。线性求解适用范围f(x)为多项式、指数函数、正弦/余弦函数，以及这些函数的乘积形式。三类函数NonhomogeneousEquations两类常见的f(x)形式本节重点讨论两种典型f(x)形式的特解求法，涵盖了大多数工程应用问题。第一类：指数多项式型f(x)=eλx·Pm(x)其中λ为常数，Pm(x)为m次多项式SpecialCases01当λ=0时退化为纯多项式型f(x)=Pm(x)02当Pm(x)为常数时退化为纯指数型第二类：指数三角型f(x)=eλx[Pl(x)cosωx+Qn(x)sinωx]其中λ、ω为常数且ω≠0SpecialCases01当λ=0时为纯三角函数型02多项式为常数时简化为eλx(Acosωx+Bsinωx)Chapter02f(x)=eλxPm(x)型特解求法指数函数与多项式乘积型右端项的特解推导与分情况讨论Methodology特解形式的推测依据由于多项式与指数函数乘积的导数仍为同类形式，因此当f(x)=eλxPm(x)时，可推测特解y*=R(x)eλx，其中R(x)为待定多项式。01结构特征f(x)=eλx·Pm(x)是指数函数与m次多项式的乘积，这是推测特解形式的出发点。EXPONENTIAL×POLYNOMIAL02导数性质多项式与指数函数乘积的导数仍为"多项式×指数函数"形式：[Q(x)eλx]′=eλx[Q′(x)+λQ(x)]CLOSUREPROPERTY03特解推测设y*=R(x)·eλx，其中R(x)为待定多项式，代入方程后左端也将呈现相同结构。y*=R(x)·eλx04核心问题确定R(x)的次数（是否与Pm(x)同次）以及各项系数，这需要分情况讨论。DEGREE&COEFFICIENTSDERIVATION代入推导过程将y*=R(x)eλx代入原方程并消去eλx，得R''+(2λ+p)R'+(λ²+pλ+q)R=Pm(x)，其中λ²+pλ+q恰为特征方程在r=λ处的值。01一阶导数y*'=[R(x)eλx]'=eλx[λR(x)+R'(x)]02二阶导数y*''=eλx[λ²R(x)+2λR'(x)+R''(x)]03代入方程将y*、y*'、y*''代入y''+py'+qy=eλxPm(x)，两边消去eλx≠004化简结果R''(x)+(2λ+p)R'(x)+(λ²+pλ+q)R(x)=Pm(x)，这是确定R(x)的关键方程CONSTANT-COEFFICIENTNON-HOMOGENEOUSLINEARODE情况一：λ不是特征方程的根当λ不是特征方程r²+pr+q=0的根时，λ²+pλ+q≠0，此时R(x)与Pm(x)同次，设R(x)=Rm(x)，代入化简后的方程，通过比较系数确定m+1个待定系数，特解为y*=Rm(x)e^(λx)。条件与结论CONDITIONλ²+pλ+q≠0，即λ不是特征方程的根CONCLUSIONy*=Rm(x)·eλx，其中Rm(x)是与Pm(x)同次的m次多项式λ²+pλ+q≠0求解方法STEP1·SETUP设Rm(x)=b₀xᵐ+b₁xᵐ⁻¹+…+bm₋₁x+bm，共m+1个待定系数STEP2·SOLVE代入R''+(2λ+p)R'+(λ²+pλ+q)R=Pm(x)，比较x各次幂系数，得m+1个方程的方程组m+1待定系数EXAMPLE·CASE1情况一例题解析通过求解y''-5y'+6y=(x+1)e^(4x)的特解，演示情况一的完整流程：确认λ=4不是特征根，设R(x)=ax+b，代入比较系数求得a=1/2,b=-3/4，最终特解y*=(x/2-3/4)e^(4x)。01识别参数f(x)=(x+1)e4x，故λ=4，Pm(x)=x+1为一次多项式（m=1）λ=4,m=102检验特征根特征方程r²−5r+6=0的根为r₁=2、r₂=3，λ=4不是特征根，属于情况一r₁=2,r₂=303设特解形式y*=(ax+b)e4x，其中a、b为待定系数，R(x)=ax+b与Pm(x)同次R(x)=ax+b04代入求解计算y*'、y*''后代入原方程，比较系数得a=1/2、b=−3/4y*=(x/2−3/4)e4xCHAPTER·SECOND-ORDERODE情况二：λ是特征方程的单根当λ是特征方程的单根时，λ²+pλ+q=0但2λ+p≠0，此时R(x)必须提高一次，设R(x)=x·Rm(x)。这是因为R项系数为零后，若R仍为m次则左端最高次为m-1，无法匹配右端Pm(x)的次数，需要乘以x来补偿。条件分析01条件：λ²+pλ+q=0（λ是特征根）且2λ+p≠0（非重根，是单根）特征方程有单根λ，意味着λ在特征多项式中只出现一次，不是重根情形。02方程退化：R''+(2λ+p)R'+0·R=Pm(x)，R项消失，左端次数降低由于R的系数2λ+p≠0但R本身的系数为零，导致方程结构发生变化，需要调整特解形式。λ²+pλ+q=0(单根条件)特解修正01设R(x)=x·Rm(x)，在原Rm(x)前乘以x，使R的次数提高一次乘以x是为了补偿R项消失导致的次数不足，确保方程两边多项式次数能够匹配。02特解形式：y*=x·Rm(x)·eλx，其中Rm(x)仍为m次多项式确定系数的方法与情况一相同，通过比较系数法求解Rm(x)的各项系数。y*=x·Rm(x)·eλx(修正后的特解)CASESTUDY·情况二情况二例题解析求解y''−5y'+6y=xe2x的特解：λ=2是特征方程的单根，设y*=x(ax+b)e2x，解得a=½,b=−1，特解为y*=x(x/2−1)e2x。STEP01参数识别f(x)=xe2x，故λ=2，Pm(x)=x（m=1），需判断λ与特征根的关系。λ=2STEP02特征根判断r²−5r+6=0解得r₁=2、r₂=3，λ=2是单根，属于情况二。r₁=2STEP03设特解形式y*=x·(ax+b)·e2x，比情况一多了x因子，R(x)为二次多项式。x·R(x)STEP04代入求解代入化简，利用λ是特征根的性质，解得a=½、b=−1。a=½,b=−1CASEIII·特征方程重根情况三：λ是特征方程的重根当λ是特征方程的重根时，λ²+pλ+q=0且2λ+p=0，原方程退化为R''=Pm(x)。为使R''为m次多项式，R必须是m+2次，因此设R(x)=x²·Rm(x)，特解y*=x²·Rm(x)·eλx。这种情况在物理上对应共振现象。条件与退化λ²+pλ+q=0,2λ+p=0λ²+pλ+q=0且2λ+p=0，即λ是二重特征根方程完全退化为R''=Pm(x)，特征根项与一阶导项系数均为零R''=Pm(x)完全退化形式特解形式y*=x²·Rm(x)·eλx设R(x)=x²·Rm(x)，需在Rm(x)前乘以x²才能匹配多项式次数物理意义：共振导致振幅随时间增长，x²项对应二次增长共振物理对应现象EXAMPLE·CASEIII情况三例题解析求解y''-4y'+4y=e^(2x)的特解：特征方程(r-2)²=0有二重根r=2，而λ=2正是重根，属于情况三。设y*=ax²e^(2x)，代入后利用重根性质，大量项消去，最终解得a=1/2，特解为y*=(x²/2)e^(2x)。此结果体现了共振现象：外激励频率与系统固有频率相同时振幅持续增长。01参数识别f(x)=e2x，故λ=2，Pm(x)=1（0次多项式，m=0）。这是确定特解形式的第一步，需要准确识别外激励函数的参数特征。02特征根判断r²-4r+4=(r-2)²=0，r=2是二重根，λ=2为重根，属于情况三。当外激励指数与特征重根相等时，必须采用修正的特解形式。03设特解形式y*=a·x²·e2x，在常数前乘以x²（注意不是x）。由于是二重根情况，修正因子为x的k次方，k=2。04代入求解计算各阶导数代入，利用重根性质化简，得2a=1，解得a=1/2。最终特解y*=(x²/2)e^(2x)，完成了情况三的完整求解。常系数非齐次线性微分方程统一结论与k值判定规则对于f(x)=e^(λx)Pm(x)型，特解统一形式为y*=xᵏ·Rm(x)·e^(λx)，解题核心是先确定k值。k值判定规则表k值等于λ作为特征根的重数，决定了特解中x因子的幂次λ与特征根的关系k值特解形式λ不是特征根（λ²+pλ+q≠0）0y*=Rm(x)·eλxλ是特征方程的单根（λ²+pλ+q=0,2λ+p≠0）1y*=x·Rm(x)·eλxλ是特征方程的重根（λ²+pλ+q=0,2λ+p=0）2y*=x²·Rm(x)·eλxStandardProcedure第一类f(x)求解标准流程求解f(x)=e^(λx)Pm(x)型特解的标准流程为六步法：求特征根→识别参数→确定k值→设特解形式→代入比较系数→解方程组。掌握这一流程可以系统化地解决所有此类问题，避免遗漏关键步骤。前期准备步骤1–301求特征根解特征方程r²+pr+q=0，求出特征根r₁、r₂02识别参数从f(x)中识别λ值和多项式Pm(x)的次数m03确定k值判断λ是特征根的几重根，确定k值（0、1或2）求解执行步骤4–604设特解形式设Rm(x)为完整m次多项式，写出y*=xᵏRm(x)eλx05代入方程计算y*'和y*''，代入原方程，两边消去eλx06比较系数比较x各次幂系数，建立m+1元线性方程组，求解Rm(x)系数COMPREHENSIVEEXAMPLE综合例题：求通解求y''+y'-2y=(3x+2)eˣ的通解。先求齐次通解Y=C₁eˣ+C₂e⁻²ˣ；再由λ=1是单根确定k=1，设y*=x(ax+b)eˣ，代入求得a=1、b=-1；最终通解y=C₁eˣ+C₂e⁻²ˣ+x(x-1)eˣ。STEP01求齐次通解特征方程r²+r−2=0解得r₁=1、r₂=−2齐次通解Y=C₁eˣ+C₂e⁻²ˣSTEP02确定k值f(x)=(3x+2)eˣ中λ=1是特征单根故k=1，Pₘ(x)=3x+2为一次多项式STEP03设特解形式y*=x·(ax+b)·eˣ其中R₁(x)=ax+b为一次多项式STEP04代入求解计算导数代入，比较系数得a=1、b=−1特解y*=x(x−1)eˣFINAL通解y=C₁eˣ+C₂e⁻²ˣ+x(x−1)eˣChapter03三角函数型f(x)特解求法利用欧拉公式将三角函数转化为复指数形式求解SecondType·f(x)第二类f(x)的形式与特征第二类f(x)含指数·三角·多项式三类因子。λ=0退化为纯三角函数型，求解策略是利用欧拉公式将三角函数转化为复指数形式。一般形式f(x)=eλx[Pl(x)cosωx+Qn(x)sinωx]，其中λ、ω为常数，ω≠0Pl(x)为l次多项式，Qn(x)为n次多项式，令m=max(l,n)为最高次数m=max(l,n)常见特例λ=0：f(x)=Pl(x)cosωx+Qn(x)sinωx，纯三角函数型多项式为常数：f(x)=eλx(Acosωx+Bsinωx)，最简振荡激励欧拉公式→复指数METHODII·EXTENSION欧拉公式的应用欧拉公式e^(iθ)=cosθ+isinθ建立了三角函数与复指数的联系。由此可将含三角函数的f(x)转化为复指数形式求解。01欧拉公式e^(iθ)=cosθ+i·sinθ，连接指数函数与三角函数的核心公式，被誉为"数学中最美的公式"当θ=π时，得到e^(iπ)+1=0，将五个基本数学常数统一COREIDENTITY02三角函数指数化cosωx=(e^(iωx)+e^(−iωx))/2，sinωx=(e^(iωx)−e^(−iωx))/(2i)将三角函数表示为复指数的线性组合，便于微分运算EXPONENTIALFORM03转化示例e^(λx)cosωx=½[e^((λ+iω)x)+e^((λ−iω)x)]，变为两个复指数之和同理e^(λx)sinωx=(1/2i)[e^((λ+iω)x)−e^((λ−iω)x)]DECOMPOSITION04求解思路f(x)分解为复指数形式后，每个分量按第一类方法求特解，再利用叠加原理合成最终取实部或虚部即得原方程的特解SUPERPOSITIONMethodology第二类f(x)特解的统一形式对于f(x)=eλx[Pl(x)cosωx+Qn(x)sinωx]，特解形式为y*=xᵏ·eλx[Rm⁽¹⁾(x)cosωx+Rm⁽²⁾(x)sinωx]，其中m=max(l,n)，k取决于复数λ+iω是否为特征根。特解形式y*=xk·eλx·[Rm(1)(x)cosωx+Rm(2)(x)sinωx]m=max(l,n)为多项式最高次数，Rm(1)和Rm(2)各含m+1个待定系数m+1系数k值判定判断复数λ+iω是否为特征方程r²+pr+q=0的根（注意是复数，非实数λ）λ+iω不是特征根时k=0；是特征根时k=1（二阶方程最多为单根）k=0或1ParticularSolution·kValuek值判定的关键：复数特征根第二类f(x)的k值判定依据是复数λ+iω是否为特征根，而非实数λ。这是因为欧拉公式转化后，指数参数变为复数λ±iω。对于实系数二次特征方程，复数根必然共轭成对出现，故λ+iω是根则λ-iω也是根，此时k=1。核心区别第一类看实数λ是否为特征根；第二类看复数λ+iω是否为特征根λvsλ+iω原因解析欧拉公式转化后，f(x)对应e^((λ+iω)x)和e^((λ-iω)x)，需按指数参数判定Euler转化共轭性质实系数方程若有复数根，必为共轭对α±βi，λ+iω是根则λ-iω也是α±βi常见错误误用实数λ判断k值，正确做法是计算(λ+iω)²+p(λ+iω)+q是否为零代入验算例题·第二类f(x)第二类f(x)例题解析特解：特征根r=±2i与λ+iω=2i重合，k=1共振。设y*=x(acos2x+bsin2x)，解得y*=(x/4)sin2x，x因子体现振幅线性增长。01参数识别f(x)=cos2x，故λ=0、ω=2、l=n=0（多项式为常数），m=max(0,0)=002判断k值特征方程r²+4=0的根为r=±2i，λ+iω=2i是特征根，故k=1（共振）03设特解形式y*=x·(a·cos2x+b·sin2x)，其中R₀⁽¹⁾=a，R₀⁽²⁾=b为待定常数04代入求解计算导数代入，比较cos2x和sin2x系数得a=0、b=1/4，特解y*=(x/4)sin2xEXAMPLE·特解求解非共振情况例题求y''+y'+y=e^x·cosx的特解：特征根r=(-1±i√3)/2，复数λ+iω=1+i不是特征根，故k=0（非共振）。设y*=e^x(acosx+bsinx)，无需x因子。代入求解可得a=3/5、b=-1/5，特解y*=e^x(3cosx-sinx)/5。EQUATIONy''+y'+y=ex·cosxf(x)=ex·cosx01参数识别由f(x)=ex·cosx识别：λ=1、ω=1、l=n=0、m=0λ=1ω=102判断k值特征根r=(−1±i√3)/2，λ+iω=1+i不是特征根，故k=0（非共振）k=0非共振03设特解形式y*=ex(a·cosx+b·sinx)，无需乘以x因子y*=ex(acosx+bsinx)04代入求解计算各阶导数代入原方程，比较系数得a=3/5、b=−1/5y*=ex(3cosx−sinx)/5CHAPTER04n阶方程推广与实际应用从二阶推广到n阶，探索微分方程在物理工程中的应用特解结构n阶方程的推广结论n阶常系数非齐次线性微分方程的特解形式与二阶相同：y*=xk·Rm(x)·eλx。关键区别在于k的取值：k等于λ作为n阶特征方程根的重数，可以是0到n之间的整数。EQUATIONn阶方程形式y(n)+p₁y(n−1)+···+pn−1y′+pny=f(x)其中p₁,…,pn为常数p₁…pₙ=constCHARACTERISTIC特征方程rn+p₁rn−1+···+pn=0n次代数方程，可能有n个根（含重根）deg=n01k的取值范围k=λ作为特征根的重数，取值范围为0,1,2,…,n；二阶方程时k最多为2，n阶时可达n。0→n02高重根示例若λ是四阶方程的三重根，则k=3，特解需乘以x³；推广至n阶时，k值可能大于2。k=3Application应用一：受迫振动问题弹簧振子受迫振动的运动方程mx''+cx'+kx=F₀cosωt是典型的二阶常系数非齐次线性微分方程。当外力频率ω接近系统固有频率ω₀=√(k/m)时发生共振（对应k=1），振幅随时间线性增长。这一现象解释了桥梁共振、机械振动等工程问题。01运动方程mx''+cx'+kx=F₀cosωt，m为质量，c为阻尼系数，k为弹簧刚度，F₀cosωt为周期外力02固有频率ω₀=√(k/m)，系统自由振动时的频率，由系统本身参数决定03共振条件当ω≈ω₀且阻尼c很小时，系统发生共振，对应λ+iω=iω₀为特征根的情况04工程启示设计机械结构时需避开共振频率，如建筑抗震设计、桥梁防共振措施Application·DifferentialEquations应用二：RLC串联电路RLC串联电路的电荷方程Lq''+Rq'+q/C=E(t)与机械振动方程具有相同数学结构。当外加交流电源频率ω接近电路固有频率ω₀=1/√(LC)时发生谐振，电流幅值显著增大。收音机调台正是利用谐振原理，通过调节电容匹配电台频率。电路方程核心方程Lq''+Rq'+q/C=E(t)q为电荷量，L为电感，R为电阻，C为电容。该二阶常微分方程完整描述了RLC电路中电荷随时间演化的动力学行为。力学类比电感L↔质量m（惯性储能），电阻R↔阻尼c（能量耗散），电容倒数1/C↔弹簧刚度k（弹性恢复）。这种对应关系揭示了电磁系统与机械系统深层的统一性。Lq''+Rq'+q/C=E(t)谐振现象固有频率ω₀=1/√(LC)电路自由振荡的特征频率，由电感L与电容C共同决定。当外加电源频率等于ω₀时，电路发生谐振，阻抗最小，电流达到峰值。收音机选台调节可变电容C改变ω₀，使其与目标电台频率ω精确匹配，实现谐振选台。频率选择性由品质因数Q=ω₀L/R决定，Q值越高选频特性越尖锐。ω₀=1/√(LC)RESONANCECHAPTER05方法总结与综合练习系统梳理两类求解方法，识别常见错误，巩固解题能力COMPARISON两类f(x)求解方法对比两类f(x)的求解核心差异在于：第一类特解含一个待定多项式Rm(x)，k值由实数λ与特征根的关系决定；第二类特解含两个待定多项式，k值由复数λ+iω与特征根的关系决定。两类f(x)求解要点对比表对比项

TYPEIeλxPm(x)

TYPEIIeλx[Plcosωx+Qnsinωx]特解形式y*=xk·Rm(x)·eλxy*=xk·eλx[Rm(1)cosωx+Rm(2)sinωx]待定多项式数1个：Rm(x)2个：Rm(1)(x)和Rm(2)(x)多项式次数m次（与Pm同次）m=max(l,n)次k值判定依据实数λ是否为特征根复数λ+iω是否为特征根k的可能取值0、1、2（二阶）0、1（二阶最多单根）两类方法的本质区别在于特解结构和k值判定依据的不同COMMONMISTAKES常见错误与注意事项求解常系数非齐次线性微分方程时的四大常见错误：设特解时漏项、k值判定依据混淆、忘记乘xk因子、混淆特解与通解。设特解时的错误漏项错误—Rm(x)必须设为完整的m次多项式如m=2时应设ax²+bx+c，不能只设ax²。漏设低次项会导致方程无法求解，必须保留所有次数项。漏乘xk—当λ是特征根时忘记乘以xkk值由λ的重数决定：单根k=1，重根k=2。漏乘会导致特解形式错误，代入后无法满足原方程。k值判定错误依据混淆—第一类看实数λ，第二类看复数λ+iω第一类非齐次项f(x)=eλxPm(x)用λ判定；第二类f(x)=eλx[Pcosωx+Qsinωx]用λ+iω判定。两类标准不可混用。概念混淆—特解y*不含任意常数C₁、C₂特解是满足非齐次方程的某个具体解，通解y=Y+y*才含任意常数。求特解时若保留C₁、C₂会导致结果错误。SuperpositionPrinciple叠加原理：f(x)为多项之和当f(x)=f₁(x)+f₂(x)+...+fₖ(x)时，可分别求出对应各fᵢ(x)的特解yᵢ*，则y*=y₁*+y₂*+...+yₖ*为原方程的特解。这一叠加原理使得复杂的f(x)可以分解为若干简单项分别求解，是处理复合右端项的有效策略。叠加原理若f(x)=f₁(x)+f₂(x)，y₁*和y₂*分别是对应f₁和f₂的特解，则y₁*+y₂*是原方程的特解。y₁*+y₂*应用示例f(x)=3eˣ+cos2