二阶常系数非齐次线性微分方程的一般形式是
y′′+py′+qy=f(x)(1) y''+py'+qy=f(x) \tag{1} y′′+py′+qy=f(x)(1)
其中p,q是常数。

一般而言,二阶常系数非齐次线性微分方程的通解,归结为求对应的齐次方程。
y′′+py′+qy=0(2) y''+py'+qy=0 \tag{2} y′′+py′+qy=0(2)
的通解和非齐次方程(1)本身的一个特解。由于二阶常系数齐次线性微分方程的通解的求法之前已得到解决,所以这里只需讨论求二阶常系数非齐次线性微分方程的一个特解y∗y^*y∗的方法。

这里只讨论当方程(1)中的f(x)f(x)f(x)取两种常见形式时求y∗y^*y∗的方法。这种方法的特点是不用积分就可求出y∗y^*y∗来,它叫做待定系数法。f(x)f(x)f(x)的两种形式是

(1)f(x)=eλxPm(x)f(x)=e^{\lambda x}P_m(x)f(x)=eλxPm​(x),其中λ\lambdaλ是常数,Pm(x)P_m(x)Pm​(x)是x的一个m次多项式:
pm(x)=a0xm+a1xm−1+⋅⋅⋅+am−1x+am p_m(x)=a_0x^m+a_1x^{m-1}+···+a_{m-1}x+a_m pm​(x)=a0​xm+a1​xm−1+⋅⋅⋅+am−1​x+am​
(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],其中λ,ω\lambda,\omegaλ,ω是常数,ω≠0,Pl(x)、Qn(x)\omega\neq 0,P_l(x)、Q_n(x)ω​=0,Pl​(x)、Qn​(x)分别是xxx的lll次、n次多项式,且仅有一个可为零。

一、f(x)=eλxPm(x)f(x)=e^{\lambda x}P_m(x)f(x)=eλxPm​(x)型

方程(1)的特解y∗y^*y∗是使(1)成为恒等式的函数。怎样的函数能使(1)成为恒等式呢?因为(1)式右端f(x)f(x)f(x)是多项式Pm(x)P_m(x)Pm​(x)与指数函数eλxe^{\lambda x}eλx的乘积,而多项式与指数函数乘积的导数仍然是多项式与指数函数的乘积,因此,我们推测y∗=R(x)eλxy^*=R(x)e^{\lambda x}y∗=R(x)eλx(其中R(x)R(x)R(x)是某个多项式)可能是方程(1)的特解。

把y∗,y∗′,y∗′′y^*,y^*{'},y^*{''}y∗,y∗′,y∗′′代入方程(1),然后考虑能否选取适当的多项式R(x)R(x)R(x),使y∗=R(x)eλxy^*=R(x)e^{\lambda}xy∗=R(x)eλx满足方程(1)。为此,将
y∗=R(x)eλxy∗′=eλx[λR(x)+R′(x)]y∗′′=eλx[λ2R(x)+2λR′(x)+R′′(x)] y^*=R(x)e^{\lambda x} \\ y^*{'}=e^{\lambda x}[\lambda R(x)+R'(x)] \\ y^*{''}=e^{\lambda x}[\lambda^2 R(x)+2\lambda R'(x)+R''(x)] y∗=R(x)eλxy∗′=eλx[λR(x)+R′(x)]y∗′′=eλx[λ2R(x)+2λR′(x)+R′′(x)]
代入方程(1)并消去eλxe^{\lambda x}eλx,得
R′′(x)+(2λ+p)R′(x)+(λ2+pλ+q)R(x)=Pm(x)(3) R''(x)+(2\lambda+p)R'(x)+(\lambda^2+p\lambda+q)R(x)=P_m(x) \tag{3} R′′(x)+(2λ+p)R′(x)+(λ2+pλ+q)R(x)=Pm​(x)(3)
(i)如果λ\lambdaλ不是(2)式的特征方程r2+pr+q=0r^2+pr+q=0r2+pr+q=0的根,即λ2+pλ+q≠0\lambda^2+p\lambda+q\neq 0λ2+pλ+q​=0,由于Pm(x)P_m(x)Pm​(x)是一个m次多项式,要使(3)的两端恒等,那么可令R(x)R(x)R(x)为另一个m次多项式Rm(x)R_m(x)Rm​(x):
Rm(x)=b0xm+b1xm−1+⋅⋅⋅+bm−1x+bm R_m(x)=b_0x^m+b_1x^{m-1}+···+b_{m-1}x+b_m Rm​(x)=b0​xm+b1​xm−1+⋅⋅⋅+bm−1​x+bm​
代入(3)式,比较等式两端x同次幂的系数,就得到以b0,b1,⋅⋅⋅,bmb_0,b_1,···,b_mb0​,b1​,⋅⋅⋅,bm​作为未知数的m+1个方程的联立方程组。从而可以定出这些bi(i=0,1,⋅⋅⋅,m)b_i(i=0,1,···,m)bi​(i=0,1,⋅⋅⋅,m),并得到所求的特解y∗=Rm(x)eλxy^*=R_m(x)e^{\lambda x}y∗=Rm​(x)eλx。

(ii)如果λ\lambdaλ是特征方程r2+pr+q=0r^2+pr+q=0r2+pr+q=0的单根,即λ2+pλ+q=0\lambda^2+p\lambda+q=0λ2+pλ+q=0,但2λ+p≠02\lambda+p\neq 02λ+p​=0,要使(3)的两端恒等,那么R′(x)R'(x)R′(x)必须是m次多项式。此时可令
R(x)=xRm(x) R(x)=xR_m(x) R(x)=xRm​(x)
并且可用同样的方法来确定Rm(x)R_m(x)Rm​(x)的系数bi(i=0,1,2,⋅⋅⋅,m)b_i(i=0,1,2,···,m)bi​(i=0,1,2,⋅⋅⋅,m)

(iii)如果λ\lambdaλ是特征方程r2+pr+q=0r^2+pr+q=0r2+pr+q=0的重根,即λ2+pλ+q=0\lambda^2+p\lambda+q=0λ2+pλ+q=0,且2λ+p=02\lambda+p=02λ+p=0,要使(3)的两端恒等,那么R′′(x)R''(x)R′′(x)必须是m次多项式。此时可令
R(x)=x2Rm(x) R(x)=x^2R_m(x) R(x)=x2Rm​(x)
并用同样的方法来确定Rm(x)R_m(x)Rm​(x)中的系数。

综上论述,有如下结论:

如果f(x)=eλxPm(x)f(x)=e^{\lambda x}P_m(x)f(x)=eλxPm​(x),那么二阶常系数非齐次线性微分方程(1)具有形如
y∗=xkRm(x)eλx(4) y^*=x^kR_m(x)e^{\lambda x} \tag{4} y∗=xkRm​(x)eλx(4)
的特解,其中Rm(x)R_m(x)Rm​(x)是与Pm(x)P_m(x)Pm​(x)同次(m次)的多项式,而k按λ\lambdaλ不是特征方程的根、是特征方程的单根或是特征方程的重根依次为0、1或2.

上式结论可推广到n阶常系数非齐次线性微分方程,但要注意(4)式中k是特征方程含根λ\lambdaλ的重复次数(即若λ\lambdaλ不是特征方程的根,则k取为0;若λ\lambdaλ是特征方程的s重根,则k取为s)。

例1:求微分方程y′′−2y′−3y=3x+1y''-2y'-3y=3x+1y′′−2y′−3y=3x+1的一个特解。

解:这是二阶常系数非齐次性微分方程,且函数f(x)f(x)f(x)是eλxPm(x)e^{\lambda x}P_m(x)eλxPm​(x)型(其中λ=0,Pm(x)=3x+1\lambda=0,P_m(x)=3x+1λ=0,Pm​(x)=3x+1)

与所给方程对应的齐次方程为
y′′−y′=3y=0 y''-y'=3y=0 y′′−y′=3y=0
它的特征方程为
r2−2r−3=0 r^2-2r-3=0 r2−2r−3=0
由于这里λ=0\lambda=0λ=0不是特征方程的根,所以应设特解为
y∗=b0x+b1 y^*=b_0x+b_1 y∗=b0​x+b1​
把它代入所给方程,得
−3b0x−2b0−3b1=3x+1 -3b_0x-2b_0-3b_1=3x+1 −3b0​x−2b0​−3b1​=3x+1
比较两端x同次幂的系数,得
{−3b0=3,−2b0−3b1=1 \begin{cases} -3b_0=3,\\ -2b_0-3b_1=1 \end{cases} {−3b0​=3,−2b0​−3b1​=1​
由此求得b0=−1,b1=13b_0=-1,b_1=\frac{1}{3}b0​=−1,b1​=31​。于是求得一个特解为
y∗=−x+13 y^*=-x+\frac{1}{3} y∗=−x+31​

二、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]型

应用欧拉公式
cos θ=eiθ+e−iθ2, sin θ=eiθ−e−iθ2i cos\,\theta=\frac{e^{i\theta}+e^{-i\theta}}{2},\, sin\,\theta=\frac{e^{i\theta}-e^{-i\theta}}{2i} cosθ=2eiθ+e−iθ​,sinθ=2ieiθ−e−iθ​
把f(x)f(x)f(x)表示成复变指数函数的形式,有
f(x)=eλx[P1cos wx+Qnsin wx]=P(x)e(λ+wi)x+P‾(x)e(λ−wi)x f(x)=e^{\lambda x}[P_1cos\, wx+Q_nsin\,wx] \\ =P(x)e^{(\lambda+wi)x}+\overline P(x)e^{(\lambda-wi)x} f(x)=eλx[P1​coswx+Qn​sinwx]=P(x)e(λ+wi)x+P(x)e(λ−wi)x
其中
P(x)=P12−Qn2i,P‾(x)=P12+Qn2i P(x)=\frac{P_1}{2}-\frac{Q_n}{2}i,\overline P(x)=\frac{P_1}{2}+\frac{Q_n}{2}i P(x)=2P1​​−2Qn​​i,P(x)=2P1​​+2Qn​​i
是互成共轭的m次多项式(即它们对应项的系数是共轭复数),而m=max[l,n]m=max[l,n]m=max[l,n]。

应用上一目的结果,对于f(x)f(x)f(x)中的第一项P(x)e(λ+ωi)xP(x)e^{(\lambda+\omega i)x}P(x)e(λ+ωi)x,可求出一个m次多项式Rm(x)R_m(x)Rm​(x),使得y1∗=xkRme(λ+ωi)xy_1^*=x^kR_me^{(\lambda+\omega i)x}y1∗​=xkRm​e(λ+ωi)x为方程
y′′+py′+qy=P(x)e(λ+wi)x y''+py'+qy=P(x)e^{(\lambda+wi)x} y′′+py′+qy=P(x)e(λ+wi)x
的特解,其中k按λ+wi\lambda+wiλ+wi不是特征方程的根或是特征方程的单根依次取0或1.由于f(x)f(x)f(x)的第二项P‾(x)e(λ−wi)x\overline P(x)e^{(\lambda-wi)x}P(x)e(λ−wi)x与第一项P(x)⋅e(λ+wi)xP(x)·e^{(\lambda+wi)x}P(x)⋅e(λ+wi)x成共轭,所以与y1∗y_1^*y1∗​成共轭的函数y2∗=xkR‾me(λ−wi)xy_2^*=x^k\overline R_me^{(\lambda-wi)x}y2∗​=xkRm​e(λ−wi)x必然是方程
y′′+py′+qy=P‾(x)e(λ−wi)x y''+py'+qy=\overline P(x)e^{(\lambda-wi)x} y′′+py′+qy=P(x)e(λ−wi)x
的特解,这里R‾m\overline R_mRm​表示与RmR_mRm​成共轭的m次多项式。于是,根据叠加原理知,方程(1)具有形如
y∗=xkRme(λ+wi)x+xkR‾me(λ−wi)x y^*=x^kR_me^{(\lambda+wi)x}+x^k\overline R_me^{(\lambda-wi)x} y∗=xkRm​e(λ+wi)x+xkRm​e(λ−wi)x
的特解。上式可写为
y∗=xkeλx[Rmewxi+R‾me−wxi]=xkeλx[Rm(cos wx+isin wx)+R‾m(cos wx−isin wx)] y^*=x^ke^{\lambda x}[R_me^{wxi}+\overline R_me^{-wxi}] \\ = x^ke^{\lambda x}[R_m(cos\,wx+isin\,wx)+\overline R_m(cos\,wx-isin\,wx)] y∗=xkeλx[Rm​ewxi+Rm​e−wxi]=xkeλx[Rm​(coswx+isinwx)+Rm​(coswx−isinwx)]
由于括号内的两项是互成共轭的,相加后既无虚部,所以可以写成实函数的形式
y∗=xkeλx[Rm(1)(x)cos wx+Rm(2)(x)sin wx] y^*=x^ke^{\lambda x}[R_m^{(1)}(x)cos\,wx+R_m^{(2)}(x)sin\,wx] y∗=xkeλx[Rm(1)​(x)coswx+Rm(2)​(x)sinwx]
综上所述,有如下结论:

如果f(x)=eλx[P1(x)cos wx+Qn(x)sin wx]f(x)=e^{\lambda x}[P_1(x)cos\,wx+Q_n(x)sin\,wx]f(x)=eλx[P1​(x)coswx+Qn​(x)sinwx],则二阶常系数非齐次线性微分方程(1)的特解可设为
y∗=xkeλx[Rm(1)(x)cos wx+Rm(2)(x)sin wx](5) y^*=x^ke^{\lambda x}[R^{(1)}_m(x)cos\,wx+R_m^{(2)}(x)sin\,wx] \tag{5} y∗=xkeλx[Rm(1)​(x)coswx+Rm(2)​(x)sinwx](5)
其中Rm(1)(x)、Rm(2)(x)R_m^{(1)}(x)、R_m^{(2)}(x)Rm(1)​(x)、Rm(2)​(x)是m次多项式,m=max{l,n}m=max\{l,n\}m=max{l,n},而k按λ+wi\lambda+wiλ+wi(或λ−wi\lambda-wiλ−wi)不是特征方程的根、或是特征方程的单根依次取0或1.

上述结论可推广到n阶常系数非齐次线性微分方程,但要注意(5)式中的k是特征方程中含根λ+wi\lambda+wiλ+wi(或λ−wi\lambda-wiλ−wi)的重复次数。

例2:求微分方程y′′+y=xcos 2xy''+y=xcos\,2xy′′+y=xcos2x的一个特解。

解:所给方程是二阶常系数非齐次线性方程,且f(x)f(x)f(x)属于eλx[Pl(x)cos wx+Qn(x)sin wx]e^{\lambda x}[P_l(x)cos\,wx+Q_n(x)sin\,wx]eλx[Pl​(x)coswx+Qn​(x)sinwx]型(其中λ=0,w=2,Pl(x)=x,Qn(x)=0\lambda=0,w=2,P_l(x)=x,Q_n(x)=0λ=0,w=2,Pl​(x)=x,Qn​(x)=0)。

与所给方程对应的齐次方程为
y′′+y=0 y''+y=0 y′′+y=0
它的特征方程为
r2+1=0 r^2+1=0 r2+1=0
由于这里λ+wi=2i\lambda+wi=2iλ+wi=2i不是特征方程的根,所以应设特解为
y∗=(ax+b)cos 2x+(cx+d)sin 2x y^*=(ax+b)cos\,2x+(cx+d)sin\,2x y∗=(ax+b)cos2x+(cx+d)sin2x
把它代入所给方程,得
(−3ax−3b+4c)cos 2x−(3cx+3d+4a)sin 2x=xcos 2x (-3ax-3b+4c)cos\,2x-(3cx+3d+4a)sin\,2x=xcos\,2x (−3ax−3b+4c)cos2x−(3cx+3d+4a)sin2x=xcos2x
比较两端同类项的系数,得
{−3a=1,−3b+4c=0,−3c=0,−3d−4a=0 \begin{cases} -3a=1,\\ -3b+4c=0,\\ -3c=0,\\ -3d-4a=0 \end{cases} ⎩⎪⎪⎪⎨⎪⎪⎪⎧​−3a=1,−3b+4c=0,−3c=0,−3d−4a=0​
由此解得
a=13,b=0,c=0,d=49 a=\frac{1}{3},b=0,c=0,d=\frac{4}{9} a=31​,b=0,c=0,d=94​
于是求得一个特解为
y∗=−13xcos 2x+49sin 2x y^*=-\frac{1}{3}xcos\,2x+\frac{4}{9}sin\,2x y∗=−31​xcos2x+94​sin2x

Logo

AtomGit 是由开放原子开源基金会联合 CSDN 等生态伙伴共同推出的新一代开源与人工智能协作平台。平台坚持“开放、中立、公益”的理念,把代码托管、模型共享、数据集托管、智能体开发体验和算力服务整合在一起,为开发者提供从开发、训练到部署的一站式体验。

更多推荐