*本文略去了很多证明,只记录结论
*文中的微分方程均指代二阶常系数非线性微分方程
二阶常系数非齐次线性微分方程的形式为:ay′′+by′+cy=f(x)ay'' + by' + cy = f(x)ay′′+by′+cy=f(x)
微分方程的通解 = 对应的二阶常系数齐次线性微分方程通解 + 自身的一个特解
简单记为:通解 = 齐次通解 + 特解。
二阶常系数齐次线性微分方程通解的解法:二阶常系数齐次线性微分方程的通解
下面只需要解出微分方程的特解即可
对应微分方程:ay′′+by′+cy=f(x)ay'' + by' + cy = f(x)ay′′+by′+cy=f(x)右式f(x)f(x)f(x)有两种形式:
①f(x)=eλxPm(x)f(x) = e^{\lambda x}P_m(x)f(x)=eλxPm(x)型
此时微分方程对应的特解为:y∗=xkRm(x)eλxy^* = x^kR_m(x)e^{\lambda x}y∗=xkRm(x)eλx
其中:
Rm(x)R_m(x)Rm(x)是于pm(x)p_m(x)pm(x)同次(mmm次)的多项式(例如:b1x2+b1x+b2b_1x^2 +b_1x + b_2b1x2+b1x+b2)kkk按λ\lambdaλ不是特征方程的根、是单根、是重根依次取0、1、2
②f(x)=eλx[Pl(x)cosωx+Qn(x)sinωx]f(x) = e^{\lambda x}[P_l(x)\cos \omega x + Q_n(x)\sin \omega x]f(x)=eλx[Pl(x)cosωx+Qn(x)sinωx]型
此时微分方程对应的特解为:y∗=xkeλx[Rm(1)cosωx+Rm(2)sinωx]y^* = x^ke^{\lambda x}[R_m^{(1)}\cos \omega x + R_m^{(2)}\sin \omega x]y∗=xkeλx[Rm(1)cosωx+Rm(2)sinωx]
其中:
Rm(1)和Rm(2)R_m^{(1)}和R_m^{(2)}Rm(1)和Rm(2)是同型mmm次多项式,m=max(l,n)m = \max{(l ,n)}m=max(l,n)(即m = l、n中的最大值;所谓同型即例如:Rm(1)=ax+b,Rm(2)=cx+dR_m^{(1)} = ax + b, R_m^{(2)} = cx + dRm(1)=ax+b,Rm(2)=cx+d)kkk按 λ+ωi\lambda + \omega iλ+ωi 不是特征方程的根、是单根依次取0、1
得到这个不完全的特解后根据需要求出其不同阶的导数然后带入微分方程,即可解出特解中的系数,到这里,就得到了微分方程的完整特解,于齐次通解相加即的微分方程的通解。
例:
求微分方程 2y′′+y′−y=2ex2y'' + y' - y = 2e^x2y′′+y′−y=2ex 的通解
解:
微分方程对应的齐次微分方程的特征方程为 2r2+r−1=02r^2 + r - 1 = 02r2+r−1=0
可得通解:y=C1e−x+C2e12xy = C_1e^{-x} + C_2e^{\frac{1}{2}x}y=C1e−x+C2e21x
微分方程的右式f(x)=2exf(x) = 2e^xf(x)=2ex满足f(x)=eλxPm(x)f(x) = e^{\lambda x}P_m(x)f(x)=eλxPm(x)型,且λ=1,m=0\lambda = 1, m = 0λ=1,m=0,
所以,设特解为:y∗=aexy^* = ae^xy∗=aex
所以y∗=aexy^* = ae^xy∗=aex、y∗′=aexy^{*'} = ae^xy∗′=aex、y∗′′=aexy^{*''} = ae^xy∗′′=aex
带入微分方程左式得:2aex+aex−aex=2ex2ae^x + ae^x - ae^x = 2e^x2aex+aex−aex=2ex
得:a=1a = 1a=1
所以特解为:y∗=exy^* = e^xy∗=ex
微分方程的通解为:y=C1e−x+C2e12x+exy = C_1e^{-x} + C_2e^{\frac{1}{2}x} + e^xy=C1e−x+C2e21x+ex