503 字
3 分钟
6.8 常系数非齐次微分方程的解法

§\S6.8 常系数非齐次线性微分方程#

待定系数法#

方程通解 = 齐次通解 + 非齐次特解

1. 情况一#

形如 f(x)=eλxPm(x)f(x)=e^{\lambda x}P_m(x) 的式子,设 y=eλxQ(x)y^*=e^{\lambda x}Q(x),则

Q(x)+(p+2λ)Q(x)+(λ2+pλ+q)Q(x)=Pm(x)Q''(x)+(p+2\lambda)Q'(x)+(\lambda ^2+p\lambda +q)Q(x)=P_m(x)
  1. λ2+pλ+q0\lambda^2+p\lambda+q\ne 0λ\lambda 不是特征根,Q(x)=Qm(x)Q(x)=Q_m(x)
  2. λ2+pλ+q=0\lambda^2+p\lambda+q=0p+2λ0p+2\lambda\ne 0λ\lambda 为单根,Q(x)=Qm+1(x)[xQm(x)]Q(x)=Q_{m+1}(x)\quad [xQ_m(x)]
  3. λ2+pλ+q=0\lambda^2+p\lambda+q=0p+2λ=0p+2\lambda=0λ\lambda 为重根,Q(x)=Qm+2(x)[x2Qm(x)]Q(x)=Q_{m+2}(x)\quad [x^2Q_m(x)]

f(x)=eλxPm(x)f(x)=e^{\lambda x}P_m(x),设 y=eλxxkQm(x)y^*=e^{\lambda x}x^kQ_m(x),其中 Qm(x)Q_m(x) 的最高次幂与 Pm(x)P_m(x) 相同,且

k={0,λ不是特征根1,λ是单根2,λ是重根k=\begin{cases} 0, \lambda 不是特征根\\ 1, \lambda 是单根\\ 2, \lambda 是重根 \end{cases}

例题#

  1. y7y+12y=xy''-7y'+12y=x

    解:特征方程 r27r+12=0r^2-7r+12=0r1=3r_1=3r2=4r_2=4

    齐次通解 y^=C1e3x+C2e4x\hat{y}=C_1e^{3x}+C_2e^{4x}

    λ=0\lambda=0 不是特征根,设非齐次特解为 y=e0xx0(ax+b)=ax+by^*=e^{0\cdot x}x^0\cdot(ax+b)=ax+b

    07a+12(ax+b)=x{a=112b=71440-7a+12(ax+b)=x \Rightarrow \begin{cases} a=\dfrac{1}{12}\\ b=\dfrac{7}{144} \end{cases}

    y=C1e3x+C2e4x+112x+7144\therefore y=C_1e^{3x}+C_2e^{4x}+\dfrac{1}{12}x+\dfrac{7}{144}

  2. y5y+6y=xe2xy''-5y'+6y=xe^{2x}

    解:r25r+6=0r^2-5r+6=0r1=2r_1=2r2=3r_2=3

    齐次通解 y^=C1e2x+C2e3x\hat{y}=C_1e^{2x}+C_2e^{3x}

    λ=2\lambda=2 是单根,设 y=e2xx1(ax+b)=e2x(ax2+bx)y^*=e^{2x}x^1(ax+b)=e^{2x}(ax^2+bx)

    (y)=e2x(4ax2+4bx+4ax+2b+4ax+2b+2a)(y^*)''=e^{2x}(4ax^2+4bx+4ax+2b+4ax+2b+2a)

    • x2x^2 系数:4a52a+6a=04a-5\cdot 2a+6a=0,自动成立
    • x1x^1 系数:4b+8a5(2b+2a)+6b=14b+8a-5(2b+2a)+6b=1
    • x0x^0 系数:4b+2a5b+60=04b+2a-5b+6\cdot 0=0

    解得 a=12a=-\dfrac{1}{2}b=1b=-1

    y=C1e2x+C2e3xe2x(12x2x)\therefore y=C_1e^{2x}+C_2e^{3x}-e^{2x}(-\dfrac{1}{2}x^2-x)

  3. y+4y+4y=e2xy''+4y'+4y=e^{-2x}

    解:r2+4r+4=0r^2+4r+4=0r1=r2=2r_1=r_2=-2

    y^=C1e2x+C2xe2x\hat{y}=C_1e^{-2x}+C_2xe^{-2x}

    λ=2\lambda=-2 是重根,设 y=e2xx2ay^*=e^{-2x}x^2\cdot a

    (y)=e2x(4ax28ax+2a)(y^*)''=e^{-2x}(4ax^2-8ax+2a)

    • x2x^2 系数:自动为0
    • x1x^1 系数:8a+8a+40=0-8a+8a+4\cdot 0=0,自动成立
    • x0x^0 系数:2a=12a=1

    解得 a=12a=\dfrac{1}{2}

    y=C1e2x+C2xe2x+12x2e2x\therefore y=C_1e^{-2x}+C_2xe^{-2x}+\dfrac{1}{2}x^2e^{-2x}

2. 情况二#

形如 f(x)=eαx[Am(x)cosβx+Bn(x)sinβx]f(x)=e^{\alpha x}[A_m(x)\cos \beta x+B_n(x)\sin \beta x] 的式子,设 y=eαxxk[Ql(x)cosβx+Rl(x)sinβx]y^*=e^{\alpha x}x^k[Q_l(x)\cos \beta x+R_l(x)\sin \beta x]l=max{m,n}l=\max\{m,n\},其中

{k=0α+iβ不是特征根k=1α+iβ是特征根\begin{cases} k=0 \quad \alpha+i\beta 不是特征根\\ k=1 \quad \alpha+i\beta 是特征根 \end{cases}

例题#

  1. y+y=4sinxy''+y=4\sin x

    解:r2+1=0r^2+1=0r1=ir_1=ir2=ir_2=-i

    y^=C1cosx+C2sinx\hat{y}=C_1\cos x+C_2\sin xλ=0+i=i\lambda =0+i=i 是特征根

    y=x(acosx+bsinx)y^*=x(a\cos x+b\sin x)

    (y)=2(asinx+bcosx)+(acosxbsinx)(y^*)''=2(-a\sin x+b\cos x)+(-a\cos x-b\sin x)

    a=2a=-2b=0b=0

    y=C1cosx+C2sinx2xcosx\therefore y=C_1\cos x+C_2\sin x-2x \cos x

  2. y+y=3sin2xy''+y=3\sin 2x

    解:r2+1=0r^2+1=0r1=ir_1=ir2=ir_2=-i

    λ=0\lambda=0 不是特征根,设 y=acos2x+bsin2xy^*=a\cos 2x+b\sin 2x

    (y)=4acos2x4bsin2x(y^*)''=-4a\cos 2x-4b\sin 2x

    {3a=03b=3{a=0b=1\begin{cases} -3a=0\\ -3b=3 \end{cases}\Rightarrow \begin{cases} a=0\\ b=-1 \end{cases}

    y=C1cosx+C2sinxsin2x\therefore y=C_1\cos x+C_2\sin x-\sin 2x

  3. y+y=exsinxy''+y=e^x\sin x

    解:r2+1=0r^2+1=0r1=ir_1=ir2=ir_2=-i

    λ=1+i\lambda=1+i 不是特征根

    y=exx0(acosx+bsinx)=ex(acosx+bsinx)y^*=e^x\cdot x^0(a\cos x+b\sin x)=e^x(a\cos x+b\sin x)

    (y)=ex(2asinx+2bcosx)(y^*)''=e^x(-2a\sin x+2b\cos x)

    {b2a=1a+2b=0{a=25b=15\begin{cases} b-2a=1\\ a+2b=0 \end{cases} \Rightarrow \begin{cases} a=-\dfrac{2}{5}\\ b=\dfrac{1}{5} \end{cases}

    y=C1cosx+C2sinx+ex(25cosx+15sinx)\therefore y=C_1\cos x+C_2\sin x+e^x(-\dfrac{2}{5}\cos x+\dfrac{1}{5}\sin x)

  4. 已知 y+ay+by=Cexy''+ay'+by=Ce^xy=ex(1+xe2x)y^*=e^{-x}(1+xe^{2x}),求该方程的通解

    • 填空解法:y=exxkay^*=e^x\cdot x^k\cdot a,令1和-1是两个单根

      (r+1)(r1)=r21=0(r+1)(r-1)=r^2-1=0a=0a=0b=1b=-1

      yy=Cexy''-y=Ce^x

      2ex+xexxex=Cex2e^x+xe^x-xe^x=Ce^xC=2C=2

      y=C1ex+C2ex+xex\therefore y=C_1e^x+C_2e^{-x}+xe^x [注意 yy^* 的形式]

    • 大题解法:将 yy^* 代回原方程,比较系数确定a、b、c

  5. y2y+2y=x2+2excos2x2y''-2y'+2y=x^2+2e^x\cos^2{\dfrac{x}{2}}

    • 利用叠加原理,y=y1(f1(x)=0的特解)+y2(f2(x)=0的特解)y^*=y_1^*(f_1(x)=0的特解)+y_2^*(f_2(x)=0的特解)

    解:r1=1+ir_1=1+ir2=1ir_2=1-iy^=ex(C1cosx+C2sinx)\hat{y}=e^x(C_1\cos x+C_2\sin x)

    注意等式右方可以拆成 x2+excosx+exx^2+e^x\cos x+e^x

    y1=e0xx0(ax2+bx+c)y1=12(x+1)2(f(x)=x2)y2=exx1(acosx+bsinx)y2=x2exsinx(f(x)=excosx)y3=exx0ay3=ex(f(x)=ex)\begin{align*} y_1^*=e^{0x}\cdot x^0(ax^2+bx+c)&\rightarrow y_1^*=\dfrac{1}{2}(x+1)^2\quad(f(x)=x^2)\\ y_2^*=e^x\cdot x^1(a'\cos x+b'\sin x)&\rightarrow y_2^*=\dfrac{x}{2}e^x\sin x\quad(f(x)=e^x\cos x)\\ y_3^*=e^x\cdot x^0 \cdot a''&\rightarrow y_3^*=e^x \quad(f(x)=e^x) \end{align*}

    y=y^+y=ex(C1cosx+C2sinx)+12(x+1)2+ex+x2exsinx\therefore y=\hat{y}+y^*=e^x(C_1\cos x+C_2\sin x)+\dfrac{1}{2}(x+1)^2+e^x+\dfrac{x}{2}e^x\sin x

3. 推广#

y(n)+P1y(n1)++Pn1y+Pny=f(x)y^{(n)}+P_1y^{(n-1)}+\cdots+P_{n-1}y'+P_ny=f(x)
  1. 形如 f(x)=eλxPm(x)f(x)=e^{\lambda x}P_m(x),设 y=eλxxkQm(x)y^*=e^{\lambda x}x^kQ_m(x),其中 kkxx 的重数有关
  2. 形如 f(x)=eαx[Am(x)cosβx+Bn(x)sinβx]f(x)=e^{\alpha x}[A_m(x)\cos \beta x+B_n(x)\sin \beta x],设 y=eαxxk[Ql(x)cosβx+Rl(x)sinβx](l=max{m,n})y^*=e^{\alpha x}x^k[Q_l(x)\cos \beta x+R_l(x)\sin \beta x]\quad(l=\max\{m,n\}),其中 kkα+iβ\alpha + i\beta 的重数有关

例题#

y(4)+y=x+ex+3sinxy^{(4)}+y''=x+e^x+3\sin x

解:r4+r2=r2(r2+1)=0r^4+r^2=r^2(r^2+1)=0r1=r2=0r_1=r_2=0r3=ir_3=ir4=ir_4=-i

y^=C1+C2x+C3cosx+C4sinx\hat{y}=C_1+C_2x+C_3\cos x+C_4\sin x

y1=e0xx2(ax+b)y1=16x3(y(4)+y=x)y2=exx0cy2=12ex(y(4)+y=ex)y3=e0xx0(dcosx+esinx)y3=32xcosx(y(4)+y=3sinx)\begin{align*} y_1^*=e^{0x}\cdot x^2(ax+b)&\rightarrow y_1^*=\dfrac{1}{6}x^3\quad(y^{(4)}+y''=x)\\ y_2^*=e^x\cdot x^0\cdot c&\rightarrow y_2^*=\dfrac{1}{2}e^x\quad(y^{(4)}+y''=e^x)\\ y_3^*=e^{0x}\cdot x^0 \cdot (d\cos x+e\sin x)&\rightarrow y_3^*=\dfrac{3}{2}x\cos x\quad(y^{(4)}+y''=3\sin x) \end{align*}

y=y^+y=C1+C2x+C3cosx+C4sinx+16x3+12ex+32xcosx\therefore y=\hat{y}+y^*=C_1+C_2x+C_3\cos x+C_4\sin x+\dfrac{1}{6}x^3+\dfrac{1}{2}e^x+\dfrac{3}{2}x\cos x

6.8 常系数非齐次微分方程的解法
https://gitee.com/jason_ren/advanced-math-note
作者
Jason Ren
发布于
2025-12-03
许可协议
CC BY-SA 4.0

部分信息可能已经过时

封面
Music
Artist
封面
Music
Artist
0:00 / 0:00