常系数非齐次线性微分方程(两种常见形式)
二阶常系数非齐次线性微分方程的一般形式是
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)=a0xm+a1xm−1+⋅⋅⋅+am−1x+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)=b0xm+b1xm−1+⋅⋅⋅+bm−1x+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∗=b0x+b1
把它代入所给方程,得
−3b0x−2b0−3b1=3x+1
-3b_0x-2b_0-3b_1=3x+1
−3b0x−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[P1coswx+Qnsinwx]=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−2Qni,P(x)=2P1+2Qni
是互成共轭的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∗=xkRme(λ+ω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∗=xkRme(λ−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∗=xkRme(λ+wi)x+xkRme(λ−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[Rmewxi+Rme−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∗=−31xcos2x+94sin2x
AtomGit 是由开放原子开源基金会联合 CSDN 等生态伙伴共同推出的新一代开源与人工智能协作平台。平台坚持“开放、中立、公益”的理念,把代码托管、模型共享、数据集托管、智能体开发体验和算力服务整合在一起,为开发者提供从开发、训练到部署的一站式体验。
更多推荐



所有评论(0)