f(x)=eλxP(x)(λ可以是复数,Pm(x)是m次多项式)

一、f(x)=e λxP m(x)(λ可以是复数,P m(x)是m次多项式)

根据函数f(x)的特点,推测方程具有形如y=Q(x)eλx的特解(其中Q(x)是某个待定的多项式).由y=Q(x)eλx,有

y*′=[Q′(x)+λQ(x)]eλx, y*″=[Q″(x)+2λQ′(x)+λ2 Q(x)]eλx

代入原方程,整理得

[Q″(x)+(2λ+p)Q′(x)+(λ2+pλ+q)Q(x)]eλx=Pm(x)eλx

消去eλx,得

由特征多项式φ(r)=r2+pr+q,知φ′(r)=2r+p.上式可简记为

注意到,方程(3)是关于多项式Q(x)的微分方程,只要设出Q(x)的一般形式,可利用比较系数的方法求出其各项系数.

例1 解微分方程y″-7y′+12y=x.

解 齐次方程为φ(r)=r2-7r+12=0,解得r1=3,r2=4.

对应齐次方程的通解为

=C1e3x+C2e4x,其中C1,C2为任意常数.

设原方程特解为y=Q(x),将φ(0)=12,φ′(0)=-7,Pm(x)=x代入方程(2)得

Q″(x)-7Q′(x)+12Q(x)=x,

比较等式两边次数,Q(x)应为一次多项式.设Q(x)=ax+b,则

-7a+12(ax+b)=x,

比较系数得

解得a=

因此原方程的一个特解为

故原方程的通解为

依照上面分析,也可以直接写出Q(x)的具体待定形式.事实上,由(2)式知:

若λ不是特征根,则λ2+pλ+q≠0,故Q(x)应为m次多项式,可设特解y=Qm(x)eλx=(a0xm+a1xm-1+…+am),eλx(a0,a1,…,am为特定系数);

若λ是特征根并且为单根,则λ2+pλ+q=0,但2λ+p≠0,若Q′(x)应为m次多项式,可设特解y=xQm(x)eλx=x(a0xm+a1xm1+…+am)eλx(a0,a1…,am为待定系数);

若λ是特征根,且为二重根,则λ2+pλ+q=0,且2λ+p=0,故Q″(x)应为m次多项式,可设特解y=x2 Qm(x)eλx=x2(a0xm+a1xm1+…+am)eλx(a0,a1,…,am为待定系数).

综上,非齐次线性方程特解设法为考虑方程

其对应的齐次方程为

特征方程为

于是可设方程(4)的特解为

其中Qm(x)=a0xm+a1xm+…+am,(a0,a1,…,am为待定系数).

k取值由λ是特征方程(6)的k重根决定,k=0,1,2.

接下来,用比较系数方法可以求出a0,a1,…,am从而写出方程(4)的特解.

例2 解微分方程y″+4y′+4y=.

解 特征方程φ(r)=r2+4r+4=0,解得r=-2(重根).对应的齐次微分方程的通解为

其中C1,C2为任意常数.

由于-2是二重根,可设原方程特解为

y=e-2x·x2·A,

于是

(y)′=(2xA-2x2 A)e-2x, (y)″=(2A-8xA+4x2 A)e-2x

代入原方程得

因此原方程的一个特解为

故原方程的通解为