# 第六章 微分方程

242 阅读2分钟

§1§ 基础定义

  1. 什么是微分方程?
  2. 什么是微分方程的阶?
  3. 什么是微分方程的特解和通解?

1.1 概念

思考:

当函数y(x)y(x)满足什么条件时,dydx=y\frac{\mathrm{d}y}{\mathrm{d}x}=y.

1.2 阶

【例 6.1 改】

  1. y2+xy=xy^2+xy=x\underline{\qquad}方程.

  2. y2+xy=xy^2+xy'=x\underline{\qquad}阶微分方程.

  3. y2+xy=exy'^2+xy=e^x\underline{\qquad}阶微分方程.

  4. y2+xy=xy^{2}+xy''=x\underline{\qquad}阶微分方程.

  5. yyy=xyy'''y''y'=xy\underline{\qquad}阶微分方程.

1.3 特解、通解

【例 6.2】
下列选项中\underline{\qquad}是微分方程y=6x+2y''=6x+2的特解,\underline{\qquad}是该方程的通解.

A. y=x3+x2+x+Cy=x^{3}+x^{2}+x+C.

B. y=x3+x2+x+1y=x^{3}+x^{2}+x+1.

C. y=x3+x2+C1x+C2y=x^3+x^2+C_1x+C_2.

D. y=x3+C1x2+x+C2y=x^3+C_1x^2+x+C_2.

1.4 线性微分方程

y(n)+y(n1)++y+y=\Box y^{(n)}+\Box y^{(n-1)}+\cdots+\Box y'+\Box y=\Box

§2§ 一阶方程求解方法

  1. 分离变量法
  2. 换元法
  3. 一阶线性方程公式法

2.1 分离变量

【例 6.4 】【1+y2=C(x21)1+y^2=C(x^2-1)

求微分方程 (xy2+x)dx+(yx2y)dy=0(xy^2+x)\mathrm{d}x+(y-x^2y)\mathrm{d}y=0 的通解.

【例 6.5】【y=xexy=xe^x

求微分方程 y=1+xxyy'=\frac{1+x}xy 满足 y(1)=ey(1)=e 的特解.

2.2 换元

yx=u(x)\frac{y}{x}=u(x) ,则 dydx=u+xdudx\frac{\mathrm{d}y}{\mathrm{d}x}=u+x\frac{\mathrm{d}u}{\mathrm{d}x} .

注:y=uxy=ux,则y=u+xuy'=u+xu'

【例 6.6】【答案:eyx=lnx+C-e^{-\frac{y}{x}}=\ln x+C

解方程dydx=eyx+yx\frac{\mathrm{d}y}{\mathrm{d}x}=e^{\frac yx}+\frac yx.

【例 6.7】【答案:y=2x1+x2y=\frac{2x}{1+x^2}

x2y+xy=y2x^2y'+xy=y^2,y(1)=1y(1)=1的特解.

【补充】【答案:arctanyx12ln(1+y2x2)=lnx+C\arctan\frac{y}{x}-\frac{1}{2}\ln(1+\frac{y^2}{x^2})=\ln x+C

dydx=x+yxy\frac{\mathrm{d}y}{\mathrm{d}x}=\frac{x+y}{x-y}的通解.

2.3 其它换元

【补充】【答案:yarctan(x+y)=π4y-\arctan(x+y)=-\frac{\pi}{4}

y=1(x+y)2y'=\frac{1}{(x+y)^2}满足y(1)=0y(1)=0的解.

2.4 线性方程

  1. 一阶线性微分方程:
    dydx+P(x)y=Q(x)\frac{\mathrm{d}y}{\mathrm{d}x}+P(x)y=Q(x)

  2. 通解公式:
    y=eP(x)dx(Q(x)eP(x)dxdx+C)y=e^{-\int P(x)\mathrm{d}x}\left(\int Q(x)e^{\int P(x)\mathrm{d}x}\mathrm{d}x+C\right)

【例6.8】【答案:y=(1+x)2[23(1+x)32+C]y=(1+x)^2\left[\frac23(1+x)^{\frac32}+C\right]

求方程dydx2yx+1=(x+1)52\frac{\mathrm{d}y}{\mathrm{d}x}-\frac{2y}{x+1}=(x+1)^{\frac52}的通解.

【例 6.9】【答案:y=13xlnx19xy=\frac13x\ln x-\frac19x

求微分方程的特解xy+2y=xlnxxy'+2y=x\ln x,其中y(1)=19y(1)=-\frac19.

  1. yyxx可以互换

【补充】【答案:2x=y22x=y^2

求方程(y26x)y+2y=0(y^2-6x)y'+2y=0的特解,其中y(2)=2y(2)=2.

【补充】【答案:xy2=lny+Cxy^2=\ln y+C

求方程y3dx+(2xy21)dy=0y^3\mathrm{d}x+(2xy^2-1)\mathrm{d}y=0的通解.

2.5 伯努利方程-伯努利换元『数一』

一阶线性微分方程:
dydx+P(x)y=Q(x)yn,(n0,1)\frac{\mathrm{d}y}{\mathrm{d}x}+P(x)y=Q(x)y^n,(n\ne 0,1)

【补充】下列微分方程中属于伯努利方程的是
A. dydx3xy=e2x\frac{\mathrm{d}y}{\mathrm{d}x}-3xy=e^{-2x}.

B. y+3yx=ysinxy'+\frac{3y}{x}=y\sin x.

C. dydxycosx=y4(sinx+2cosx)\frac{\mathrm{d}y}{\mathrm{d}x}-y\cos x=y^4(\sin x+2\cos x).

D. 2xdy[yxy(1+sinx)]dx=02x\mathrm{d}y-\left[y-xy\left(1+\sin x\right)\right]\mathrm{d}x=0.

【补充】【答案:z=y2z=y^{-2},z=x1+Ce2xz=x-1+Ce^{-2x},xy2y2+Cy2e2x=1xy^2-y^2+Cy^2e^{-2x}=1
dydxy=12(12x)y3\frac{\mathrm{d}y}{\mathrm{d}x}-y=\frac12\left(1-2x\right)y^3.

【补充】【答案:z=y4z=y^{-4},z=32x+316+Ce8xz=-\frac{3}{2}x+\frac{3}{16}+Ce^{-8x},32xy4+316y4+Cy4e8x=1-\frac{3}{2}xy^4+\frac{3}{16}y^4+Cy^4e^{-8x}=1
dydx2y=3xy5\frac{\mathrm{d}y}{\mathrm{d}x}-2y=3xy^5.

§3§ 可降解的二阶方程『仅数一、数二』

3.1 二阶少yy型:

y=p(x)y'=p(x),则y=py''=p',y=py'=p.

【例 6.11】【答案:y=x3+3x+1y=x^3+3x+1

求微分方程(1+x2)y=2xy(1+x^{2})y''=2xy'满足初始条件yx=0=1,yx=0=3y|_{x=0}=1,y'|_{x=0}=3的特解.

【补充】【答案:y=1+2x1y=\sqrt{1+2x}-1

求微分方程y+y3=0y''+y'^3=0满足y(0)=0y(0)=0,y(0)=1y'(0)=1的特解.

3.2 二阶少xx型:

y=p(y)y'=p(y),则y=ppy''=pp',y=py'=p.

注:y=dydx=dpdx=dpdydydx=py=ppy''=\frac{\mathrm{d}y'}{\mathrm{d}x}=\frac{\mathrm{d}p}{\mathrm{d}x}=\frac{\mathrm{d}p}{\mathrm{d}y}\cdot\frac{\mathrm{d}y}{\mathrm{d}x}=p'y'=p'p

【例 6.12】【答案:y=C2eC1xy=C_2e^{C_1x}

求微分方程yyy2=0yy''-y'^2=0的通解.

【补充】【答案:y=tanxy=\tan x

求微分方程y=2yyy''=2yy'满足y(0)=0y(0)=0,y(0)=1y'(0)=1的特解.

【补充:此题可以用两种方法】【答案:y=ln(cos(x+C))+Cy=-\ln(\cos (x+C))+C

求微分方程y=1+y2y''=1+y'^2的通解.

3.3 降维之后是一阶方程

【补充:仅数一、二】【答案:y=13x3+C1x2+C2y=\frac{1}{3} x^3+C_1x^2+C_2

求微分方程xyy=x2xy''-y'=x^2的通解.

【补充:仅数一、二】【答案:y=C1ex12x2x+C2y=C_1e^x-\frac{1}{2}x^2-x+C_2

求微分方程y=y+xy''=y'+x的通解.

3.4 可降阶的高维方程

【例 6.10】【答案:y=18e2x+sinx+C1x2+C2x+C3y=\frac{1}{8}e^{2x}+\sin x+C_1x^2+C_2x+C_3

求微分方程y=e2xcosxy'''=e^{2x}-\cos x的通解.

§4§ 线性方程基础知识

  1. 什么是线性组合?
  2. 什么是常系数线性组合?
  3. 什么是线性相关、线性无关?
  4. 什么是线性微分方程?
  5. 齐次和非齐次是什么意思?

4.1 线性组合

f1+f2++fn\Box f_1+\Box f_2+\cdots+\Box f_n』被称为『f1f_1,f2f_2,\cdots,fnf_n』的线性组合,『\Box,\Box,\cdots,\Box』被称为组合系数.若『\Box,\Box,\cdots,\Box』均为常数,则被称为常系数线性组合.

例如:
1.『xy+exy+2y(sinx)y+lnxxy'''+e^xy''+2y-(\sin x)y+\ln x』是『 yy''' , yy'' , yy' , yy , 11 』的非常系数线性组合,组合系数为『 xx , exe^x , 22 , sinx-\sin x , lnx\ln x 』.

2.『y2y+3y+4y1y'''-2y''+3y+4y-1』是『 yy''' , yy'' , yy' , yy , 11 』的常系数线性组合,组合系数为『 11 , 2-2 , 33 , 44 , 1-1 』.

3.『4ex+3e2x+2e3x+e4x4e^x+3e^{2x}+2e^{3x}+e^{4x}』是『 exe^{x} , e2xe^{2x} , e3xe^{3x} , e4xe^{4x} 』的常系数线性组合,组合系数为『 11 , 22 , 33 , 44 』.

4.2 线性无关、线性相关

若『f1f_1,f2f_2,\cdots,fnf_n』中任何一个函数都不能表示为其它函数的常系数线性组合,则称『f1f_1,f2f_2,\cdots,fnf_n线性无关,否则称线性相关.

例如:

  1. sinx\sin x , exe^x , lnx\ln x , x2x^2 』线性无关.
  2. 11 , xx , x2x^2 , (x+1)2(x+1)^2 』线性相关.

4.3 线性微分方程

齐次非齐次
一阶y+y=0\Box y'+\Box y=0y+y=\Box y'+\Box y=\Box
二阶y+y+y=0\Box y''+\Box y'+\Box y=0y+y+y=\Box y''+\Box y'+\Box y=\Box
\cdots\cdots\cdots
n阶y(n)+y(n1)++y=0\Box y^{(n)}+\Box y^{(n-1)}+\cdots+\Box y=0y(n)+y(n1)++y=\Box y^{(n)}+\Box y^{(n-1)}+\cdots+\Box y=\Box

请判断下列方程的类型:

  1. y+xy=0y'+xy=0
  2. (sinx)y(cosx)y=0(\sin x)y'-(\cos x) y=0
  3. y=yy'=y
  4. x2y+xy+y=0x^2y''+xy'+y=0
  5. y+(sinx)y(cosx)y=0y''+(\sin x)y'-(\cos x)y=0
  6. y=xyy''=xy
  7. y=yxyy'''=y''-xy
  8. y+xy=xy'+xy=x
  9. (sinx)y(cosx)y=tanx(\sin x)y'-(\cos x) y=\tan x
  10. y=yxy'=y-x
  11. x2y+xy+y=2xx^2y''+xy'+y=2x
  12. y+(sinx)y(cosx)y=1y''+(\sin x)y'-(\cos x)y=1
  13. y=y+exy''=y+e^x
  14. y+y+y+y=sinxy'''+y''+y'+y=\sin x
  15. yy=yyy''=y'

【例6.3】下列选项中\underline{\qquad}是线性微分方程.
A. y2+xy=xy^2+xy=x
B. y+x(y)3=xy+x(y')^3=x
C. y+xy=exy'+xy=e^{x}
D. y2+xy=xy^{2}+xy''=x

4.4 解的结构

齐次非齐次
一阶y=Cy(x)y=Cy(x)y=Cy(x)+yy=Cy(x)+y^*
二阶y=C1y1(x)+C2y2(x)y=C_1y_1(x)+C_2y_2(x)y=C1y1(x)+C2y2(x)+yy=C_1y_1(x)+C_2y_2(x)+y^*
三阶y=C1y1(x)+C2y2(x)+C3y3(x)y=C_1y_1(x)+C_2y_2(x)+C_3y_3(x)y=C1y1(x)+C2y2(x)+C3y3(x)+yy=C_1y_1(x)+C_2y_2(x)+C_3y_3(x)+y^*
\cdots\cdots\cdots
n阶y=C1y1(x)+C2y2(x)++Cnyn(x)y=C_1y_1(x)+C_2y_2(x)+\cdots+C_ny_n(x)y=C1y1(x)+C2y2(x)++Cnyn(x)+yy=C_1y_1(x)+C_2y_2(x)+\cdots+C_ny_n(x)+y^*

通解

  1. 齐通=C1齐特+C2齐特+C3齐特+\text{齐通}=C_1\cdot\text{齐特}+C_2\cdot\text{齐特}+C_3\cdot\text{齐特}+\cdots.

  2. 非齐通=齐通+非齐特=C1齐特+C2齐特+C3齐特++非齐特\text{非齐通}=\text{齐通}+\text{非齐特}=C_1\cdot\text{齐特}+C_2\cdot\text{齐特}+C_3\cdot\text{齐特}+\cdots+\text{非齐特}.

这两个公式里的这些齐特必须线性无关.

齐特

  1. k齐特=齐特k\text{齐特}=\text{齐特}.

  2. 齐特±齐特=齐特\text{齐特}\pm\text{齐特}=\text{齐特}.

  3. k1齐特+k2齐特+k2齐特+=齐特k_1\cdot\text{齐特}+k_2\cdot\text{齐特}+k_2\cdot\text{齐特}+\cdots=\text{齐特}.

齐特的线性组合仍然是齐特.

非齐特

  1. 齐特+非齐特=非齐特\text{齐特}+\text{非齐特}=\text{非齐特}.

  2. 非齐特非齐特=齐特\text{非齐特}-\text{非齐特}=\text{齐特}.

  3. k非齐特非齐特,除非k=1k\cdot\text{非齐特}\ne\text{非齐特},\text{除非}k=1.

特定系数组合

k1+k2+=0,k_1+k_2+\cdots=0,

k1非齐特+k2非齐特+k2非齐特+=齐特.k_1\cdot\text{非齐特}+k_2\cdot\text{非齐特}+k_2\cdot\text{非齐特}+\cdots=\text{齐特}.
k1+k2+=1,k_1+k_2+\cdots=1,

k1非齐特+k2非齐特+k2非齐特+=非齐特.k_1\cdot\text{非齐特}+k_2\cdot\text{非齐特}+k_2\cdot\text{非齐特}+\cdots=\text{非齐特}.

非齐通的进阶表达

\qquadnn阶线性方程能找到n+1n+1个线性无关的非齐特,则非齐通可以表达为

非齐通=C1非齐特++Cn+1非齐特,\text{非齐通}=C_1\cdot\text{非齐特}+\cdots+C_{n+1}\text{非齐特},

其中

C1++Cn+1=1.C_1+\cdots+C_{n+1}=1.

\qquad但注意,并不是所有的

C1非齐特++Cn+1非齐特C_1\cdot\text{非齐特}+\cdots+C_{n+1}\text{非齐特}

都是非齐通.它必须要满足

C1++Cn+1=1C_1+\cdots+C_{n+1}=1

且这n1n-1个非奇特都线性无关.

例题

【例 6.13】
设非齐次线性微分方程y+P(x)y=Q(x)y'+P(x)y=Q(x)有两个不同的解y1(x)y_1(x)y2(x)y_2(x),CC 为任意常数,则该方程的通解是\underline{\qquad}.

A. C[y1(x)y2(x)]C[ y_{1}(x)-y_{2}(x)]

B. y1(x)+C[y1(x)y2(x)]y_{1}(x)+C[ y_{1}(x)-y_{2}(x)]

C. C[y1(x)+y2(x)]C[ y_{1}(x)+y_{2}(x)]

D. y1(x)+C[y1(x)+y2(x)]y_{1}(x)+C[ y_{1}(x)+y_{2}(x)]

【例 6.14】【答案:y=C1+C2x+C3x2y=C_1+C_2x+C_3x^2
已知y1=1,y2=x,y3=x2y_1=1,y_2=x,y_3=x^2是某二阶非齐次线性微分方程的三个解,则该方程的通解为\underline{\qquad\qquad}.

【补充】 设y1y_1,y2y_2为一阶非齐次线性微分方程y+P(x)y=Q(x)y^{\prime}+P(x)y=Q(x)的两个特解,若有常数λ\lambda,μ\mu,使λy1+μy2\lambda y_1+\mu y_2为该方程的解,λy1μy2\lambda y_1-\mu y_2为对应的齐次方程的解,则\underline{\qquad}.
A. λ=12\lambda=\frac{1}{2},μ=12\mu=\frac{1}{2}.
B. λ=12\lambda=-\frac{1}{2},μ=12\mu=-\frac{1}{2}.
C. λ=23\lambda=\frac{2}{3},μ=13\mu=\frac{1}{3}.
D. λ=23\lambda=\frac{2}{3},μ=23\mu=\frac{2}{3}.

【补充】 设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)的三个线性无关的解,则该方程的通解为\underline{\qquad}.
A. C1y1(x)+C2y2(x)+y3(x)C_1y_1(x)+C_2y_2(x)+y_3(x).
B. C1y1(x)+(12C1)y2(x)+C1y3(x)C_1y_1(x)+(1-2C_1)y_2(x)+C_1y_3(x).
C. (C1C2)y1(x)+C2y2(x)+y3(x)(C_1-C_2)y_1(x)+C_2y_2(x)+y_3(x).
D. C1y1(x)+C2y2(x)+C3y3(x)(C1+C2+C3=1)C_1y_1(x)+C_2y_2(x)+C_3y_3(x)(C_1+C_2+C_3=1).

§5§ 线性方程的计算

5.1 常系数齐次

实根、无重复

若:r1=0r_1=0,r2=1r_2=1,
则:y=C1e0x+C2e1xy=C_1e^{0\cdot x}+C_2e^{1\cdot x}

若:r1=1r_1=-1,r2=1r_2=1,r3=0r_3=0
则:y=C1e1x+C2e1x+C3e0xy=C_1e^{-1\cdot x}+C_2e^{1\cdot x}+C_3e^{0\cdot x}

若:r1=1r_1=1,r2=2r_2=2,r3=3r_3=3,r4=4r_4=4
则:y=C1e1x+C2e2x+C3e3x+C4e4xy=C_1e^{1\cdot x}+C_2e^{2\cdot x}+C_3e^{3\cdot x}+C_4e^{4\cdot x}

实根、有重复

若:r1,2=1r_{1,2}=1,
则:y=(C1+C2x)e1xy=(C_1+C_2x)e^{1\cdot x}

若:r1,2,3=0,r4,5=1,r6=1r_{1,2,3}=0,r_{4,5}=1,r_6=-1
则:y=(C1+C2x+C3x2)e0x+(C4+C5x)e1x+C6e1xy=(C_1+C_2x+C_3x^2)e^{0\cdot x}+(C_4+C_5x)e^{1\cdot x}+C_6e^{-1\cdot x}

复根、无重复

若:r1,2=3±2ir_{1,2}=3\pm2i
则:y=(C1cos2x+C2sin2x)e3xy=(C_1\cos 2x+C_2\sin 2x)e^{3\cdot x}

若:r1,2=2±4i,r3,4=ir_{1,2}=-2\pm 4i,r_{3,4}=i
则:y=(C1cos4x+C2sin4x)e2x+(C3cos(1x)+C4sin(1x))e0xy=(C_1\cos 4x+C_2\sin 4x)e^{-2\cdot x}+(C_3\cos (1\cdot x)+C_4\sin (1\cdot x))e^{0\cdot x}

复数、有重复(了解)

若:r1,2;3,4=2±3ir_{1,2;3,4}=2\pm 3i,
则:y=((C1+C2x)cos3x+(C3+C4x)sin3x)e2xy=((C_1+C_2x)\cos 3x+(C_3+C_4x)\sin 3x)e^{2\cdot x}

若:r1,2;3,4=2±3i,r5,6;7,8;9,10=ir_{1,2;3,4}=2\pm 3i,r_{5,6;7,8;9,10}=-i,
则:y=((C1+C2x)cos3x+(C3+C4x)sin3x)e2x+((C5+C6x+C7x2)cos(x)+(C8+C9x+C10)sin(x))e0xy=((C_1+C_2x)\cos 3x+(C_3+C_4x)\sin 3x)e^{2\cdot x}+((C_5+C_6x+C_7x^2)\cos (-x)+(C_8+C_9x+C_{10})\sin (-x))e^{0\cdot x}

综合

若:r1=1,r2,3,4=0,r5,6=3ir_1=1,r_{2,3,4}=0,r_{5,6}=3i,
则:y=C1e1x+(C2+C3x+C4x2)e0x+(C5cos3x+C6sin3x)e0xy=C_1e^{1\cdot x}+(C_2+C_3x+C_4x^2)e^{0\cdot x}+(C_5\cos 3x+C_6\sin 3x)e^{0\cdot x}

【例 6.15】【答案:y=C1ex+C2e6xy=C_1e^x+C_2e^{6x}
求微分方程y7y+6y=0y''-7y'+6y=0的通解.

【例 6.16】【答案:y=(C1+C2x)exy=(C_1+C_2x)e^x
求微分方程y2y+y=0y''-2y'+y=0的通解.

【例 6.17】【答案:(C1cos2x+C2sin2x)e3x(C_1\cos 2x+C_2\sin 2x)e^{3x}
求微分方程y6y+13y=0y''-6y'+13y=0的通解.

【例 6.18】【答案:C1+C2x+(C3cos2x+C4sin2x)exC_1+C_2x+(C_3\cos 2x+C_4\sin 2x)e^x
求方程y(4)2y+5y=0y^{(4)}-2y'''+5y''=0的通解.

5.2 常系数非齐次

多项式×\times指数

=(x2+2x+1)e2x\cdots=(x^2+2x+1)e^{2x}.
y=(Ax2+Bx+C)e2xxy^*=(Ax^2+Bx+C)\cdot e^{2x}\cdot x^{\Box}
其中,\Box方框里填22作为特征方程根的重数.

=xex\cdots=xe^{-x}.
y=(Ax+B)exxy^*=(Ax+B)\cdot e^{-x}\cdot x^{\Box}
其中,\Box方框里填1-1作为特征方程根的重数.

=x2+2x+1\cdots=x^2+2x+1.
y=(Ax2+Bx+C)e0xxy^*=(Ax^2+Bx+C)\cdot e^{0\cdot x}\cdot x^{\Box}
其中,\Box方框里填00作为特征方程根的重数.

=1\cdots=1.
y=Ae0xxy^*=A\cdot e^{0\cdot x}\cdot x^{\Box}
其中,\Box方框里填00作为特征方程根的重数.

多项式×\times指数×\times三角函数

=(x2+2x+1)e2xsin(3x)\cdots=(x^2+2x+1)e^{-2x}\sin(3x).
y=(Ax2+Bx+C)e2xsin(3x)x+(A~x2+B~x+C~)e2xcos(3x)xy^*=(Ax^2+Bx+C)\cdot e^{-2 x}\cdot\sin(3x)\cdot x^{\Box}+(\widetilde{A}x^2+\widetilde{B}x+\widetilde{C})\cdot e^{-2 x}\cdot\cos(3x)\cdot x^{\Box}
其中,\Box方框里填2±3i-2\pm3i作为特征方程根的重数.

=xcos(x)\cdots=x\cos(x).
y=(Ax+B)e0xcos(x)x+(A~x+B~)e0xsin(x)xy^*=(Ax+B)\cdot e^{0 x}\cdot\cos(x)\cdot x^{\Box}+(\widetilde{A}x+\widetilde{B})\cdot e^{0 x}\cdot\sin(x)\cdot x^{\Box}
其中,\Box方框里填0±i0\pm i作为特征方程根的重数.

多个类型的线性组合

=x2+ex+sinx\cdots=x^2+e^x+\sin x.
y=(Ax2+Bx+C)e0xx+Dexx+Ee0xsinxx+Fe0xcosxxy^*=(Ax^2+Bx+C)\cdot e^{0 x}\cdot x^{\Box}+D\cdot e^x\cdot x^{\Box}+E\cdot e^{0\cdot x}\cdot \sin x\cdot x^{\Box}+F\cdot e^{0\cdot x}\cdot \cos x\cdot x^{\Box}.
其中,第一个\Box方框里填00作为特征方程根的重数. 第二个\Box方框里填11作为特征方程根的重数.第三、四个\Box方框里填0±i0\pm i作为特征方程根的重数.

=sin2x\cdots=\sin^2 x.
sin2x\sin^2 x看作1212cos2x\frac{1}{2}-\frac{1}{2}\cos 2x即可.
y=Ae0xx+Be0xcos2xx+Ce0xsin2xxy^*=A\cdot e^{0 x}\cdot x^{\Box}+B\cdot e^{0\cdot x}\cdot \cos 2x\cdot x^{\Box}+C\cdot e^{0\cdot x}\cdot \sin 2x\cdot x^{\Box}.
其中,第一个\Box方框里填00作为特征方程根的重数. 第二、三个\Box方框里填0±2i0\pm 2i作为特征方程根的重数.

习题练习

【补充】

  1. y+2y+2y=exsinxy''+2y'+2y=e^{-x}\sin x的特解形式y=y^*=\underline{\hspace{5cm}}.

  2. y6y+8y=ex+e2xy''-6y'+8y=e^x+e^{2x}的特解形式y=y^*=\underline{\hspace{5cm}}.

  3. y4y+4y=x2+8e2xy''-4y'+4y=x^2+8e^{2x}的特解形式y=y^*=\underline{\hspace{5cm}}.

  4. y+2y+2y=2excos2x2y''+2y'+2y=2e^{-x}\cos^2{\frac{x}{2}}的特解形式y=y^*=\underline{\hspace{5cm}}.

系数求解

【例题】求y+2y+2y=2excos2x2y''+2y'+2y=2e^{-x}\cos^2{\frac{x}{2}}的通解y=y^*=\underline{\hspace{5cm}}.

  • 第一步 求特征方程,并求根
    此方程的特征方程为
r2+2r+2=0r^2+2r+2=0

解的此方程

r1,2=1±ir_{1,2}=-1\pm i
  • 第二步 写出此方程对应的齐次通解
    此方程的通解形式为
y=C1excosx+C2exsinx+y.y=C_1e^{-x}\cos x +C_2e^{-x}\sin x +y^*.

其中,yy^*为此方程的一个非齐次特解.

  • 第三步 根据方程的右端项设yy^*
    方程的右端项为
2excos2x2=ex+excosx.2e^{-x}\cos^2{\frac{x}{2}}=e^{-x}+e^{-x}\cos x.

则设

y=Aexx+Bexcosxx+Cexsinxx.y^*=A\cdot e^{-x}\cdot x^{\Box}+B\cdot e^{-x}\cdot \cos x\cdot x^{\Box}+C\cdot e^{-x}\cdot \sin x\cdot x^{\Box}.

其中,第一个\Box方框里填1-1作为特征方程根的重数.由于1-1不是特征方程的根,即填00. 第二、三个\Box方框里填1±i-1\pm i作为特征方程根的重数.1±i-1\pm i是特征方程的一重根,即填11.即

y=Aex+Bxexcosx+Cxexsinx.y^*=Ae^{-x}+Bxe^{-x}\cos x+Cxe^{-x}\sin x.
  • 第四步 将yy^*带入原方程中求解未知参数

2y=2Aex+2Bxexcosx+2Cxexsinx.2y^*=2Ae^{-x}+2Bxe^{-x}\cos x+2Cxe^{-x}\sin x.

2(y)=(2)Aex+2Bexcosx+(2B)xexcosx+(2B)xexsinx+2Cexsinx+(2C)xexsinx+2Cxexcosx.2(y^*)'=(-2)Ae^{-x}\\ \hspace{1.5cm}+2Be^{-x}\cos x+(-2B)xe^{-x}\cos x+(-2B)xe^{-x}\sin x\\ \hspace{1.5cm}+2Ce^{-x}\sin x+(-2C)xe^{-x}\sin x+2Cxe^{-x}\cos x.

(y)=Aex+(B)ex+(B)exsinx+(B)excosx+Bxexcosx+Bxexsinx+(B)exsinx+Bxexsinx+(B)xexcosx+(C)exsinx+Cexcosx+(C)exsinx+Cxxsinx+(C)xexcosx+Cexcosx+(C)xexcosx+(C)xexsinx.(y^*)''=Ae^{-x}\\ \hspace{1.5cm}+(-B)e^{-x}+(-B)e^{-x}\sin x\\ \hspace{1.5cm}+(-B)e^{-x}\cos x+Bxe^{-x}\cos x+Bxe^{-x}\sin x\\ \hspace{1.5cm}+(-B)e^{-x}\sin x+Bxe^{-x}\sin x+(-B)xe^{-x}\cos x\\ \hspace{1.5cm}+(-C)e^{-x}\sin x+Ce^{-x}\cos x\\ \hspace{1.5cm}+(-C)e^{-x}\sin x+Cx^{-x}\sin x+(-C)xe^{-x}\cos x\\ \hspace{1.5cm}+Ce^{-x}\cos x+(-C)xe^{-x}\cos x+(-C)xe^{-x}\sin x.

(y)+2(y)+2y=(AB)ex+(B+2C)excosx+(2B)=ex+excosx(y^*)''+2(y^*)'+2y^*\\ \hspace{1.5cm}=(A-B)e^{-x}+(B+2C)e^{-x}\cos x+(-2B)\\ \hspace{1.5cm}=e^{-x}+e^{-x}\cos x

{AB=1B+2C=12B=0{A=1B=0C=12\left\{\begin{matrix}A-B=1 \\B+2C=1\\-2B=0\end{matrix}\right.\Rightarrow \left\{\begin{matrix}A=1 \\B=0\\C=\frac{1}{2}\end{matrix}\right.

y=ex+12xexsinxy^*=e^{-x}+\frac{1}{2}xe^{-x}\sin x
  • 第五步 组合出非齐次通解
    此方程的解为
y=C1excosx+C2exsinx+ex+12xexsinx.y=C_1e^{-x}\cos x +C_2e^{-x}\sin x +e^{-x}+\frac{1}{2}xe^{-x}\sin x.

习题练习

【例 6.19】【答案:y=C1ex+C2e3x+13xy=C_1e^{-x}+C_2e^{3x}+\frac{1}{3}-x
求微分方程y2y3y=3x+1y''-2y'-3y=3x+1的通解.

【例 6.20】【答案:y=C1e2x+C2e3xx(12x+1)e2xy=C_1e^{2x}+C_2e^{3x}-x(\frac{1}{2}x+1)e^{2x}
求微分方程y5y+6y=xe2xy''-5y'+6y=xe^{2x}的通解.

【例 6.21】【答案:y=C1e2x+C2ex310cos2x+110sin2xy=C_1e^{-2x}+C_2e^x-\frac{3}{10}\cos 2x+\frac{1}{10}\sin 2x
求微分方程y+y2y=2cos2xy''+y'-2y=2\cos2x的通解.

【例6.22】微分方程y+y=x2+1+sinxy''+y=x^2+1+\sin x的特解形式可设为
A. y=ax2+bx+c+x(Asinx+Bcosx)y^*=ax^2+bx+c+x(A\sin x+B\cos x).
B. y=ax2+bx+c+Asinxy^*=ax^2+bx+c+A\sin x.
C. y=x(ax2+bx+c+Asinx+Bcosx)y^*=x(ax^2+bx+c+A\sin x+B\cos x).
D. y=ax2+bx+c+Acosxy^*=ax^2+bx+c+A\cos x.

5.3 非常系数一阶

  1. 一阶线性微分方程:
    dydx+P(x)y=Q(x)\frac{\mathrm{d}y}{\mathrm{d}x}+P(x)y=Q(x)

  2. 通解公式:
    y=eP(x)dx(Q(x)eP(x)dxdx+C)y=e^{-\int P(x)\mathrm{d}x}\left(\int Q(x)e^{\int P(x)\mathrm{d}x}\mathrm{d}x+C\right)

5.4 非常系数高阶-欧拉方程『仅数一』

【补充】下列微分方程中属于欧拉方程的是
A. d3ydx31x3d2ydx2+1x2dydx=1\frac{\mathrm{d}^3y}{\mathrm{d}x^3}-\frac1{x^3}\frac{\mathrm{d}^2y}{\mathrm{d}x^2}+\frac1{x^2}\frac{\mathrm{d}y}{\mathrm{d}x}=1.
B. 2y+y3yx=cosx2y''+y'-\frac{3y}x=\cos x.
C. y+2yx3yx2=5xsinxy''+\frac{2y'}x-\frac{3y}{x^2}=5x\sin x.
D. x2y4y+2xy=lnx12xx^2y''-4y'+2xy=\ln x-\frac12x.

【补充】 【y3y+2y=0y'''-3y'+2y=0,y=(C1+C2lnx)x+C3x2y=(C_1+C_2\ln x)x+C_3x^{-2}
x3d3ydx3+3x2d2ydx22xdydx+2y=0x^3\frac{\mathrm{d}^3y}{\mathrm{d}x^3}+3x^2\frac{\mathrm{d}^2y}{\mathrm{d}x^2}-2x\frac{\mathrm{d}y}{\mathrm{d}x}+2y=0.

【补充】 【y2y+4y=etsinty''-2y'+4y=e^t\sin t,y=x[C1sin(3lnx)+C2cos(3lnx)+12xsin(lnx)]y=x[C_1\sin(\sqrt{3}\ln x)+C_2\cos(\sqrt{3}\ln x)+\frac{1}{2}x\sin (\ln x)]
yyx+4yx2=sin(lnx)x.y''-\frac{y'}x+\frac{4y}{x^2}=\frac{\sin\left(\ln x\right)}x.

§6§ 差分方程『仅数三』

【补充】 【答案: yt=C2t+t2t1y_t=C\cdot 2^t+t\cdot 2^{t-1}
差分方程yt+12yt=2ty_{t+1}-2y_t=2^t的通解 yt=y_t=\underline{\hspace{2cm}}.

【补充】 【答案: yt=C+12t(t1)y_t=C+\frac{1}{2}t(t-1)
差分方程Δyt=t\Delta y_t=t的通解 yt=y_t=\underline{\hspace{2cm}}.

【补充】 【答案: yx=C2x5y_x=C\cdot 2^x-5
差分方程Δ2yxyx=5\Delta^2y_x-y_x=5的通解 yx=y_x=\underline{\hspace{2cm}}.

§7§ 微分方程应用

几何应用

【例题】设y(x)y(x)是区间(0,32)(0,\frac32)内的可导函数,且y(1)=0y(1)=0,点PP是曲线 L:y=y(x)L:y=y(x)上的任意一点,LL在点PP处的切线与yy轴相交于点(0,Yp)(0,Y_p),法线与 xx轴相交于点(Xp,0)(X_p,0),若Xp=YpX_p=Y_p,求LL上点的坐标(x,y)(x,y)满足的方程.
【解:】
LLP(x0,y0)P(x_0,y_0)处的切线方程为

y=y(x0)(xx0)+y0,y=y'(x_0)(x-x_0)+y_0,

所以

Yp=y0x0y(x0);Y_p=y_0-x_0y'(x_0);

对应的法线方程为

y=1y(x0)(xx0)+y0,y=-\frac1{y'(x_0)}(x-x_0)+y_0,

所以

Xp=x0+y0y(x0);X_p=x_0+y_0y'(x_0);

因此

y0x0y(x0)=x0+y0y(x0),y_0-x_0y'(x_0)=x_0+y_0y'(x_0),

由于点PP的任意性知

yxy=x+yy,y-xy'=x+yy',

dydx=yxy+x=yx1yx+1\frac{\mathrm{d}y}{\mathrm{d}x}=\frac{y-x}{y+x}=\frac{\frac yx-1}{\frac yx+1}

这是一个齐次方程,可令

u(x)=yx,u(x)=\frac yx,

最终求得方程的通解为

ln(u2+1)+2arctanu=2lnx+C,\ln(u^2+1)+2\arctan u=-2\ln|x|+C\quad,

再由y(1)=0y(1)=0

ln(x2+y2)+2arctanyx=0,x(0,32).\ln(x^2+y^2)+2\arctan\frac yx=0,x\in(0,\frac32).

【补充】【答案:y=Cx^3】

已知f(x)f(x)可导,且f(x)>0(x0)f^{\prime}(x)>0(x\geqslant0),曲线 y=f(x)y=f(x)过原点,点MM为曲线y=f(x)y=f(x)上任意一点,过点MM的切线与xx轴相交于点TT,过点MMMPMP垂直于xx轴于点PP,且曲线y=f(x)y=f(x)与直线MPMP以及xx轴所围成图形的面积与三角形MTPMTP的面积比恒为3:23:2,求曲线满足的方程.