Skip to content

第四节 高阶线性微分方程

Tables of Content

I. 线性微分方程解的性质和结构

1.1 解的性质与结构

  1. 二阶齐次线性微分方程

    形式y+P(x)y+Q(x)y=0\displaystyle y'' + P(x)y' + Q(x)y = 0

    解的性质:若y1y_1y2y_2为方程的两个解,则C1y1+C2y2C_1 y_1 + C_2 y_2也是此方程的一个解

    通解的结构:若y1y_1y2y_2为方程的两个线性无关[1]的解,则C1y1+C2y2C_1 y_1 + C_2 y_2是此方程的通解

  2. 二阶非齐次线性微分方程

    形式y+P(x)y+Q(x)y=f(x)y'' + P(x)y' + Q(x)y = f(x)

    解的性质

    • y1y_1y2y_2为方程的两个解,则y1y2y_1 - y_2是其所对应的齐次微分方程的一个解
    • η\eta为方程的一个解、而ξ\xi是其所对应的齐次微分方程的一个解,则η+ξ\eta + \xi也是此方程的一个解

    通解的结构:若y(x)y^*(x)y+P(x)y+Q(x)y=f(x)y'' + P(x) y' + Q(x) y = f(x)的一个特解、Y(x)Y(x)y+P(x)y+Q(x)y=0y'' + P(x) y' + Q(x) y = 0的通解,则y+P(x)y+Q(x)y=f(x)y'' + P(x) y' + Q(x) y = f(x)的通解为:

1.2 解的叠加性原理

  • 若:

    1. y1(x)y_1(x)y+P(x)y+Q(x)y=f1(x)y'' + P(x) y' + Q(x) y = f_1(x)的一个解
    2. y2(x)y_2(x)y+P(x)y+Q(x)y=f2(x)y'' + P(x) y' + Q(x) y = f_2(x)的一个解
    3. …………
    4. yn(x)y_n(x)y+P(x)y+Q(x)y=fn(x)y'' + P(x) y' + Q(x) y = f_n(x)的一个解
  • 则:

    是方程:

    的一个解

Proof
  • 例1:设y1(x)y_1(x)y2(x)y_2(x)y3(x)y_3(x)是方程y+p(x)y+q(x)y=f(x)y'' + p(x) y' + q(x) y = f(x)的三个线性无关的解,若C1C_1C2C_2是任意常数,则该方程的通解是_____.

    • A:C1y1(x)+C2y2(x)+y3(x)C_1 y_1(x) + C_2 y_2(x) + y_3(x)
    • B:C1y1(x)+C2y2(x)(C1+C2)y3(x)C_1 y_1(x) + C_2 y_2(x) - (C_1 + C_2) y_3(x)
    • C:C1y1(x)+C2y2(x)(1C1C2)y3(x)C_1 y_1(x) + C_2 y_2(x) - (1 - C_1 - C_2) y_3(x)
    • D:C1y1(x)+C2y2(x)+(1C1C2)y3(x)C_1 y_1(x) + C_2 y_2(x) + (1 - C_1 - C_2) y_3(x)
    Answer

    根据解的叠加性原理,对于选项A,其作为通解所对应的方程是:

    对于选项B,其作为通解所对应的方程是:

    对于选项C,其作为通解所对应的方程是:

    对于选项D,其作为通解所对应的方程是:

    故选D

II. 高阶常系数齐次线性微分方程的通解

2.1 二阶常系数齐次线性微分方程

其形式为:

其通解的求法为:

  1. 写出特征方程:λ2+pλ+q=0\lambda^2 + p \lambda + q = 0

  2. 求出特征方程的两个根λ1\lambda_1λ2\lambda_2

  3. 根据两个根的不同形式,依照下表得出通解:

    根的形式通解
    两个不相等的实根λ1λ2\lambda_1 \neq \lambda_2y=C1eλ1x+C2eλ2xy = C_1 e^{\lambda_1 x} + C_2 e^{\lambda_2 x}
    两个相等的实根λ1=λ2=r\lambda_1 = \lambda_2 = ry=(C1+C2x)erxy = (C_1 + C_2 x) e^{rx}
    一对共轭复根λ1,2=α±βi\lambda_{1,2} = \alpha \pm \beta iy=eαx(C1cosβx+C2sinβx)y = e^{\alpha x} (C_1 \cos \beta x + C_2 \sin \beta x)
Proof

注意到:

如这般的函数,其导数与自身仅差常数倍,因此可以设其解形如:

于是将这个形式代入原方程可得:

由此便将一个微分方程转换为普通一元二次代数方程,根据根的判别式:

可知存在三种情况:

  • Δ>0\Delta \gt 0时,根为λ1λ2\lambda_1 \neq \lambda_2

    于是得出两个特解:

    根据二阶齐次线性微分方程的解的性质可知,其通解为:

  • Δ=0\Delta = 0时,根为λ1=λ2=r\lambda_1 = \lambda_2 = r,且r=p2r = -\dfrac{p}{2}p2=4qp^2 = 4q

    此时仅可得出一个特解:

    又因为所要求的两个特解不能是线性相关的,所以设两个特解之比为一个函数:

    然后分别求出y2(x)y_2'(x)

    以及y2(x)y_2''(x)

    再将其代入原方程:

    于是将u(x)u(x)代入即可得到y2(x)=xerxy_2(x) = x e^{rx},然后通过这两个解便可得到通解:

  • Δ<0\Delta \lt 0时,根为λ1,2=α±βi\lambda_{1, 2} = \alpha \pm \beta i

    此时依旧可以得到两个特解,并通过欧拉公式[2]进行化简:

    又因为齐次微分方程的解的线性组合依旧是它的解,所以可以设:

    于是可以得到通解:

  • 例2:求解下列微分方程:

    1. yy6y=0\displaystyle y'' - y' - 6y = 0

      Answer
    2. y2y+y=0\displaystyle y'' - 2y' + y = 0

      Answer
    3. y4y+5y=0\displaystyle y'' - 4y' + 5y = 0

      Answer
    4. y+2y=0\displaystyle y'' + 2y = 0

      Answer

2.2 nn阶常系数齐次线性微分方程

其形式为:

其特征方程为:

根据根的不同情况得到对应的通解:

  1. 对于其中每一个不重复的实根λ\lambda,在通解中所对应的项为:

  2. 对于其中每一组kk个重复的实根λ1=λ2=...=λk=r\lambda_1 = \lambda_2 = ... = \lambda_k = r,在通解中所对应的项为:

  3. 对于其中每一对不重复的复根λ1,2=α±βi\lambda_{1, 2} = \alpha \pm \beta i,在通解中所对应的项为:

  4. 对于其中每一组kk对重复的复根λ1,2=λ3,4=...=λ2k1,2k=α±βi\lambda_{1,2} = \lambda_{3,4} = ... = \lambda_{2k - 1, 2k} = \alpha \pm \beta i,在通解中所对应的项为:

  • 例3:求解微分方程y(4)2y+5y=0\displaystyle y^{(4)} - 2y''' + 5y'' = 0.

    Answer
  • 例4:在下列微分方程中,以y=C1ex+C2cos2x+C3sin2xy = C_1 e^x + C_2 \cos 2x + C_3 \sin 2x为通解的是:

    • A:y+y4y4y=0y''' + y'' - 4y' - 4y = 0
    • B:y+y+4y+4y=0y''' + y'' + 4y' + 4y = 0
    • C:yy4y+4y=0y''' - y'' - 4y' + 4y = 0
    • D:yy+4y4y=0y''' - y'' + 4y' - 4y = 0
    Answer

III. 二阶常系数非齐次线性微分方程

其形式为:

写出它的特征方程:

求出方程的两个根λ1,λ2\lambda_1, \lambda_2,然后根据f(x)f(x)不同形式得到特解:

  1. f(x)=eλxPn(x)\displaystyle f(x) = e^{\lambda x} P_n(x)Pn(x)P_n(x)为关于xxnn次多项式

    对于特征方程λ2+pλ+q=0\lambda^2 + p\lambda + q = 0

    • λ\lambda不是特征方程的根:

    • λ\lambda是特征方程的单根:

    • λ\lambda是特征方程的二重根:

    Tip

    其中Rn(x)R_n(x)是一个与Pn(x)P_n(x)同阶的待定系数多项式,设其为:

  2. f(x)=eαx[Pl(x)cosβx+Qn(x)sinβx]\displaystyle f(x) = e^{\alpha x} \left[ P_l(x) \cos \beta x + Q_n(x) \sin \beta x \right]Pl(x)P_l(x)Qn(x)Q_n(x)分别是关于xxll次和nn次多项式

    • λ1,2=α±βi\lambda_{1,2} = \alpha \pm \beta i是特征方程的根:

    • λ1,2=α±βi\lambda_{1,2} = \alpha \pm \beta i不是特征方程的根:

    Tip

    其中m=max{l,n}m = \max \{ l, n \}Rm(1)(x)R_m^{(1)}(x)Rm(2)(x)R_m^{(2)}(x)是关于xxmm次多项式,分别设其为:

    以及:

最后再将此特解与对应的齐次微分方程y+py+qy=0y'' + py' + qy = 0的通解相加即可得到通解

  • 例5:求下列微分方程的特解:

    1. y2y3y=3x+1\displaystyle y'' - 2y' - 3y = 3x + 1

      Answer

      首先根据f(x)f(x)的形式以及特征方程的根求出特解的形式:

      然后将y=A0+A1xy^* = A_0 + A_1 x代入原微分方程:

      可知:

      于是:

      所以最终可得特解为:

    2. yy=2ex\displaystyle y'' - y' = 2e^x

      Answer

      首先根据f(x)f(x)的形式以及特征方程的根求出特解的形式:

      然后将y=A0xexy^* = A_0 x e^x代入原微分方程以求得特解:

    3. y4y+3y=xex\displaystyle y'' - 4y' + 3y = xe^x

      Answer

      首先根据f(x)f(x)的形式以及特征方程的根求出特解的形式:

      然后将y=A0xex+A1x2exy^* = A_0 x e^x + A_1 x^2 e^x代入原微分方程:

      可知:

      于是:

      所以最终可得特解为:

    4. yy=excos2x\displaystyle y'' - y = e^x \cos 2x

      Answer

      首先根据f(x)f(x)的形式以及特征方程的根求出特解的形式:

      然后将y=ex[Acos2x+Bsin2x]y^* = e^x \left[ A \cos 2x + B \sin 2x \right]代入原微分方程:

      可知:

      于是:

      所以最终可得特解为:

  • 例6:求微分方程y+y=x2+1+sinxy'' + y = x^2 + 1 + \sin x的通解

    Answer

    初步来看,这个微分方程的f(x)f(x)部分不好处理,所以我们可以将其拆成两个函数的和再分别求出特解,于是两个特解相加便可得到原微分方程的特解,最后再求出齐次形式下的通解,便可得到最终结果

    首先写出其特征方程并求根:

    • f(x)=x2+1f(x) = x^2 + 1时的特解:

      根据f(x)f(x)的形式,设特解为:

      代入可得:

      解得:

      于是可以得到特解为:

    • f(x)=sinxf(x) = \sin x时的特解:

      根据f(x)f(x)的形式,设特解为:

      代入可得:

      解得:

      于是可以得到特解为:

    • f(x)=0f(x) = 0时的通解:

      根据特征方程的根,可以得到通解为:

    所以将三个解相加,即可得到原方程的通解:

  • 例6:求解二阶微分方程y+2y=sin3xy'' + 2y = \sin 3x当满足初始条件y(0)=1y(0) = 1y(0)=1y'(0) = -1时的解

    Answer

    首先写出特征方程并求出其根:

    于是设其特解为:

    代入原方程可得:

    解得:

    于是可求得特解为:

    然后再根据根的形式得到齐次的通解:

    将两个解相加即可得到原微分方程的通解:

    再将两个初始条件代入可得到方程组:

    解得:

    所以最终可求得满足初始条件的解为:

IV. 欧拉方程

形如:

的微分方程,称之为欧拉方程;作变换x=etx = e^t,将自变量从xx换成tt,并采用记号DD表示对tt的求导运算,则有:

将其代入欧拉方程便可得到一个以tt为自变量的常系数线性微分方程,求出其解后再把t=lnxt = \ln x回代,即可得到原方程的解

Proof

设:

则:

其中:

由此可以得到结论:

  • 例7:求欧拉方程x2y+3xy+2y=5sinlnxx^2 y'' + 3x y' + 2y = 5 \sin \ln x的通解

    Answer

    x=etx = e^t,于是有:

    这是一个二阶常系数非齐次线性微分方程,找出其特征方程并求根:

    于是其对应的齐次微分方程的通解为:

    设其特解为:

    将其代入可得:

    于是有:

    解得:

    则可得到特解为:

    相加即可得到通解为:

    最后将t=lnxt = \ln x回代即可得到此欧拉方程的通解:


  1. 线性无关:y1(x)y_1(x)​除以y2(x)y_2(x)​并非恒等于一个常数CC↩︎

  2. 欧拉公式:eix=cosx+isinxe^{ix} = \cos x + i \sin x ↩︎