第七章 微分方程(一)

在许多问题中,往往不能直接找出所需要的函数关系,但是根据问题所提供的情况,有时可以列出含有要找的函数及其导数的关系式。这样的关系式就是所谓的微分方程。微分方程建立以后,对它进行研究,找出未知函数来,这就是解微分方程。本章主要介绍微分方程的一些基本概念和几种常用的微分方程的解法。——高等数学同济版

目录

习题7-1 微分方程的基本概念

  本节主要介绍了微分方程的基本概念。

7.一个半球体形状的雪堆,其体积融化率与半球面面积 A A A成正比,比例系数 k > 0 k>0 k>0。假设在融化过程中雪堆始终保持半球体状,已知半径为 r 0 r_0 r0的雪堆在开始的3小时内,融化了其体积的 7 8 \frac{7}{8} 87,问雪堆全部融化需要多少小时?

  设雪堆在时刻 t t t的体积为 V = 2 3 π r 3 V=\cfrac{2}{3}\pi r^3 V=32πr3,侧面积 S = 2 π r 2 S=2\pi r^2 S=2πr2。由题设知
d V d t = 2 π r 2 d r d t = − k S = − 2 π k r 2 . \cfrac{\mathrm{d}V}{\mathrm{d}t}=2\pi r^2\cfrac{\mathrm{d}r}{\mathrm{d}t}=-kS=-2\pi kr^2. dtdV=2πr2dtdr=kS=2πkr2.
  于是
d r d t = − k . \cfrac{\mathrm{d}r}{\mathrm{d}t}=-k. dtdr=k.
  积分得
r = − k t + C . r=-kt+C. r=kt+C.
  由 r = r 0 r=r_0 r=r0,得 C = r 0 , r = r 0 − k t C=r_0,r=r_0-kt C=r0,r=r0kt。又 V ∣ t = 3 = 1 8 V ∣ t = 0 V\biggm\vert_{t=3}=\cfrac{1}{8}V\biggm\vert_{t=0} Vt=3=81Vt=0,即 2 3 π ( r 0 − 3 k ) 3 = 1 8 ⋅ 2 3 π r 0 3 \cfrac{2}{3}\pi(r_0-3k)^3=\cfrac{1}{8}\cdot\cfrac{2}{3}\pi r_0^3 32π(r03k)3=8132πr03,得 k = 1 6 r 0 k=\cfrac{1}{6}r_0 k=61r0,从而
r = r 0 − 1 6 r 0 t . r=r_0-\cfrac{1}{6}r_0t. r=r061r0t.
  因雪堆全部融化时, r = 0 r=0 r=0,故得 t = 6 t=6 t=6,即雪堆全部融化需6小时。(这道题主要难点在体积融化率即体积融化速率的理解上

习题7-2 可分离变量的微分方程

  本节主要介绍可分离变量的微分方程的解法。

习题7-3 齐次方程

  本节主要介绍齐次方程的解法。

2. 求下列齐次方程满足所给初始条件的的特解:

(1) ( y 2 − 3 x 2 ) d y + 2 x y x = 0 , y ∣ x = 0 = 1 ; (y^2-3x^2)\mathrm{d}y+2xy\mathrm{x}=0,y\biggm\vert_{x=0}=1; (y23x2)dy+2xyx=0,yx=0=1;

  原方程可写成 1 − 3 x 2 y 2 + 2 x y d x d y = 0 1-3\cfrac{x^2}{y^2}+2\cfrac{x}{y}\cfrac{\mathrm{d}x}{\mathrm{d}y}=0 13y2x2+2yxdydx=0。令 u = x y u=\cfrac{x}{y} u=yx,即 x = y u x=yu x=yu,有 d x d y = u + y d u d y \cfrac{\mathrm{d}x}{\mathrm{d}y}=u+y\cfrac{\mathrm{d}u}{\mathrm{d}y} dydx=u+ydydu,则原方程成为
1 − 3 u 2 + 2 u ( u + y d u d y ) = 0. 1-3u^2+2u\left(u+y\cfrac{\mathrm{d}u}{\mathrm{d}y}\right)=0. 13u2+2u(u+ydydu)=0.
  分离变量得
2 u u 2 − 1 d u = d y y . \cfrac{2u}{u^2-1}\mathrm{d}u=\cfrac{\mathrm{d}y}{y}. u212udu=ydy.
  积分得
ln ⁡ ∣ u 2 − 1 ∣ = ln ⁡ ∣ y ∣ + ln ⁡ C 1 . \ln|u^2-1|=\ln|y|+\ln C_1. lnu21=lny+lnC1.
  即
u 2 − 1 = C y . u^2-1=Cy. u21=Cy.
  代入 y = x y y=\cfrac{x}{y} y=yx并整理,得通解 x 2 − y 2 = C y 3 x^2-y^2=Cy^3 x2y2=Cy3。由初值条件 x = 0 x=0 x=0 y = 1 y=1 y=1,得 C = − 1 C=-1 C=1。于是所求特解为
y 3 = y 2 − x 2 . y^3=y^2-x^2. y3=y2x2.

习题7-4 一阶线性微分方程

  本节主要介绍了一阶线性微分方程的求解方法。

7.用适当的变量代换将下列方程化为可分离变量的方程,然后求出通解:

(3) x y ′ + y = y ( ln ⁡ x + ln ⁡ y ) ; xy'+y=y(\ln x+\ln y); xy+y=y(lnx+lny);

  令 u = x y u=xy u=xy,则 u ′ = y + x y ′ u'=y+xy' u=y+xy,且原方程变为 u ′ = u x ln ⁡ u u'=\cfrac{u}{x}\ln u u=xulnu,即
d u u ln ⁡ u = d x x . \cfrac{\mathrm{d}u}{u\ln u}=\cfrac{\mathrm{d}x}{x}. ulnudu=xdx.
  积分得 ln ⁡ ∣ ln ⁡ u ∣ = ln ⁡ ∣ x ∣ + ln ⁡ C 1 \ln|\ln u|=\ln|x|+\ln C_1 lnlnu=lnx+lnC1,即 u = e C x u=e^{Cx} u=eCx。代入 u = x y u=xy u=xy,得原方程的通解 x y = e C x xy=e^{Cx} xy=eCx,即 y = e C x x y=\cfrac{e^{Cx}}{x} y=xeCx。(这道题利用变量换元时对所换的元的导数可以进行整体代换

(4) y ′ = y 2 + 2 ( sin ⁡ x − 1 ) y + sin ⁡ 2 x − 2 sin ⁡ x − cos ⁡ x + 1 ; y'=y^2+2(\sin x-1)y+\sin^2x-2\sin x-\cos x+1; y=y2+2(sinx1)y+sin2x2sinxcosx+1;

  将原方程写成 y ′ = ( y + sin ⁡ x − 1 ) 2 − cos ⁡ x y'=(y+\sin x-1)^2-\cos x y=(y+sinx1)2cosx,令 u = y + sin ⁡ x − 1 u=y+\sin x-1 u=y+sinx1,则 u ′ = y ′ + cos ⁡ x u'=y'+\cos x u=y+cosx,且原方程变为 u ′ = u 2 u'=u^2 u=u2,即 d u u 2 = d x \cfrac{\mathrm{d}u}{u^2}=\mathrm{d}x u2du=dx
  积分得 − 1 u = x + C -\cfrac{1}{u}=x+C u1=x+C,即 u = − 1 x + C u=-\cfrac{1}{x+C} u=x+C1。代入 u = y + sin ⁡ x − 1 u=y+\sin x-1 u=y+sinx1,得原方程的通解
y = 1 − sin ⁡ x − 1 x + C . y=1-\sin x-\cfrac{1}{x+C}. y=1sinxx+C1.
这道题说明换元时可以先配方再换元

(5) y ( x y + 1 ) d x + x ( 1 + x y + x 2 y 2 ) d y = 0. y(xy+1)\mathrm{d}x+x(1+xy+x^2y^2)\mathrm{d}y=0. y(xy+1)dx+x(1+xy+x2y2)dy=0.

  原方程改写成 x y ( x y + 1 ) + x 2 ( 1 + x y + x 2 y 2 ) d y d x = 0 xy(xy+1)+x^2(1+xy+x^2y^2)\cfrac{\mathrm{d}y}{\mathrm{d}x}=0 xy(xy+1)+x2(1+xy+x2y2)dxdy=0。令 u = x y u=xy u=xy,即 y = u x y=\cfrac{u}{x} y=xu,则 d y d x = x d u d x − u x 2 \cfrac{\mathrm{d}y}{\mathrm{d}x}=\cfrac{x\cfrac{\mathrm{d}u}{\mathrm{d}x}-u}{x^2} dxdy=x2xdxduu,且原方程变为
u ( u + 1 ) + ( 1 + u + u 2 ) ( x d u d x − u ) = 0. u(u+1)+(1+u+u^2)\left(x\cfrac{\mathrm{d}u}{\mathrm{d}x}-u\right)=0. u(u+1)+(1+u+u2)(xdxduu)=0.
  整理并分离变量,得 1 + u + u 2 u 3 d u = d x x \cfrac{1+u+u^2}{u^3}\mathrm{d}u=\cfrac{\mathrm{d}x}{x} u31+u+u2du=xdx。积分得 − 1 2 u 2 − 1 u + ln ⁡ ∣ u ∣ = ln ⁡ ∣ x ∣ + C 1 -\cfrac{1}{2u^2}-\cfrac{1}{u}+\ln|u|=\ln|x|+C_1 2u21u1+lnu=lnx+C1,代入 u = x y u=xy u=xy,并整理,得原方程的通解为
2 x 2 y 2 ln ⁡ ∣ y ∣ − 2 x y − 1 = C x 2 y 2 ( C = 2 C 1 ) . 2x^2y^2\ln|y|-2xy-1=Cx^2y^2(C=2C_1). 2x2y2lny2xy1=Cx2y2(C=2C1).
这道题说明可以在换元前可以先乘以某一相同变量凑整再换元

习题7-5 可降阶的高阶微分方程

从这一节其我们将讨论二阶及二阶以上的微分方程,即所谓高阶微分方程。——高等数学同济版

  本节主要介绍了三种高阶微分方程的降阶方法。

1.求下列各微分方程的通解:

(8) y 3 y ′ ′ − 1 = 0 ; y^3y''-1=0; y3y1=0;

  令 y ′ = p y'=p y=p,则 y ′ ′ = p d p d y y''=p\cfrac{\mathrm{d}p}{\mathrm{d}y} y=pdydp,且原方程化为 y 3 p d p d y − 1 = 0 y^3p\cfrac{\mathrm{d}p}{\mathrm{d}y}-1=0 y3pdydp1=0。分离变量,得
p d p = 1 y 3 d y . p\mathrm{d}p=\cfrac{1}{y^3}\mathrm{d}y. pdp=y31dy.
  积分得 p 2 = − 1 y 2 + C 1 p^2=-\cfrac{1}{y^2}+C_1 p2=y21+C1,故
y ′ = p = ± C 1 − 1 y 2 = ± 1 ∣ y ∣ C 1 y 2 − 1 . y'=p=\pm\sqrt{C_1-\cfrac{1}{y^2}}=\pm\cfrac{1}{|y|}\sqrt{C_1y^2-1}. y=p=±C1y21 =±y1C1y21 .
  分离变量,得
∣ y ∣ d y C 1 y 2 − 1 = ± d x . \cfrac{|y|\mathrm{d}y}{\sqrt{C_1y^2-1}}=\pm\mathrm{d}x. C1y21 ydy=±dx.
  由于 ∣ y ∣ = y s g n ( y ) |y|=y\mathrm{sgn}(y) y=ysgn(y),故上式两端积分,
s g n ( y ) ∫ y d y C 1 y 2 − 1 = ± ∫ d x , s g n ( y ) C 1 y 2 − 1 = ± C 1 x + C 2 . \mathrm{sgn}(y)\displaystyle\int\cfrac{y\mathrm{d}y}{\sqrt{C_1y^2-1}}=\pm\displaystyle\int\mathrm{d}x,\\ \mathrm{sgn}(y)\sqrt{C_1y^2-1}=\pm C_1x+C_2. sgn(y)C1y21 ydy=±dx,sgn(y)C1y21 =±C1x+C2.
  两边平方,得 C 1 y 2 − 1 = ( C 1 x + C 2 ) 2 C_1y^2-1=(C_1x+C_2)^2 C1y21=(C1x+C2)2。(这道题主要利用绝对值的函数表示来求解

(9) y ′ ′ = 1 y ; y''=\cfrac{1}{\sqrt{y}}; y=y 1;

  方程 y ′ ′ = 1 y y''=\cfrac{1}{\sqrt{y}} y=y 1属于 y ′ ′ = f ( x ) y''=f(x) y=f(x)型方程,除了设 y ′ = p y'=p y=p y ′ ′ = p d p d y y''=p\cfrac{\mathrm{d}p}{\mathrm{d}y} y=pdydp来降阶求解外,还可以用如下方法求解:
  在 y ′ ′ = f ( y ) y''=f(y) y=f(y)的两端乘以 2 y ′ 2y' 2y,得
2 y ′ y ′ ′ = 2 f ( y ) y ′ . 2y'y''=2f(y)y'. 2yy=2f(y)y.
  即 ( y ′ 2 ) ′ = 2 f ( y ) y ′ (y'^2)'=2f(y)y' (y2)=2f(y)y。若 F ( y ) F(y) F(y) f ( y ) f(y) f(y)的原函数,则有
( y ′ 2 ) ′ = 2 [ F ( y ) ] ′ . (y'^2)'=2[F(y)]'. (y2)=2[F(y)].
  积分得到降阶的方程 y ′ 2 = 2 F ( y ) + C 1 y'^2=2F(y)+C_1 y2=2F(y)+C1
  本小题按上述方法求解:用 2 y ′ 2y' 2y乘方程 y ′ ′ = 1 y y''=\cfrac{1}{\sqrt{y}} y=y 1的两端,得
2 y ′ y ′ ′ = 2 y ′ y . 2y'y''=\cfrac{2y'}{\sqrt{y}}. 2yy=y 2y.
  即 ( y ′ 2 ) ′ = ( 4 y ) ′ (y'^2)'=(4\sqrt{y})' (y2)=(4y ),故 y ′ 2 = 4 y + C 1 ′ y'^2=4\sqrt{y}+C_1' y2=4y +C1,有 y ′ = ± 2 y + C 1 ( C 1 = C 1 ′ 4 ) y'=\pm2\sqrt{\sqrt{y}+C_1}\left(C_1=\cfrac{C_1'}{4}\right) y=±2y +C1 (C1=4C1)
  分离变量,得
d x = ± d y 2 y + C 1 . \mathrm{d}x=\pm\cfrac{\mathrm{d}y}{2\sqrt{\sqrt{y}+C_1}}. dx=±2y +C1 dy.
  积分,得
x = ± ∫ d ( y ) 2 2 y + C 1 = ± ∫ y d y 2 y + C 1 = ± ∫ ( y + C 1 ) − C 1 2 y + C 1 d ( y ) = ± [ ∫ y + C 1 d ( y + C 1 ) − C 1 ∫ 1 y + C 1 d ( y + C 1 ) ] = ± [ 2 3 ( y + C 1 ) 3 2 − 2 C 1 ( y + C 1 ) 1 2 ] + C 2 . \begin{aligned} x&=\pm\displaystyle\int\cfrac{\mathrm{d}(\sqrt{y})^2}{2\sqrt{\sqrt{y}+C_1}}=\pm\displaystyle\int\cfrac{\sqrt{y}\mathrm{d}\sqrt{y}}{2\sqrt{\sqrt{y}+C_1}}=\pm\displaystyle\int\cfrac{(\sqrt{y}+C_1)-C_1}{2\sqrt{\sqrt{y}+C_1}}\mathrm{d}(\sqrt{y})\\ &=\pm\left[\displaystyle\int\sqrt{\sqrt{y}+C_1}\mathrm{d}(\sqrt{y}+C_1)-C_1\displaystyle\int\cfrac{1}{\sqrt{\sqrt{y}+C_1}}\mathrm{d}(\sqrt{y}+C_1)\right]\\ &=\pm\left[\cfrac{2}{3}(\sqrt{y}+C_1)^{\frac{3}{2}}-2C_1(\sqrt{y}+C_1)^{\frac{1}{2}}\right]+C_2. \end{aligned} x=±2y +C1 d(y )2=±2y +C1 y dy =±2y +C1 (y +C1)C1d(y )=±[y +C1 d(y +C1)C1y +C1 1d(y +C1)]=±[32(y +C1)232C1(y +C1)21]+C2.
这道题主要利用凑整的方法求解

(10) y ′ ′ = ( y ′ ) 3 + y ′ . y''=(y')^3+y'. y=(y)3+y.

  令 y ′ = p y'=p y=p,则 y ′ ′ = p d p d y y''=p\cfrac{\mathrm{d}p}{\mathrm{d}y} y=pdydp,原方程化为 p d p d y = p 3 + p p\cfrac{\mathrm{d}p}{\mathrm{d}y}=p^3+p pdydp=p3+p,即
p [ d p d y − ( 1 + p 2 ) ] = 0. p\left[\cfrac{\mathrm{d}p}{\mathrm{d}y}-(1+p^2)\right]=0. p[dydp(1+p2)]=0.
  若 p ≡ 0 p\equiv0 p0,则 y ≡ C y\equiv C yC y ≡ C y\equiv C yC是原方程的解,但不是通解。
  若 p ≢ 0 p\not\equiv0 p0,由于 p p p的连续性,必在 x x x的某区间有 p ≠ 0 p\neq0 p=0。于是
d p d y − ( 1 + p 2 ) = 0. \cfrac{\mathrm{d}p}{\mathrm{d}y}-(1+p^2)=0. dydp(1+p2)=0.
  分离变量,得 d p 1 + p 2 = d y \cfrac{\mathrm{d}p}{1+p^2}=\mathrm{d}y 1+p2dp=dy,积分得 arctan ⁡ p = y − C \arctan p=y-C arctanp=yC,即 p = tan ⁡ ( y − C 1 ) p=\tan(y-C_1) p=tan(yC1),亦即 cot ⁡ ( y − C 1 ) d y = d x \cot(y-C_1)\mathrm{d}y=\mathrm{d}x cot(yC1)dy=dx。积分得 ln ⁡ sin ⁡ ( y − C 1 ) = x + ln ⁡ C 2 \ln\sin(y-C_1)=x+\ln C_2 lnsin(yC1)=x+lnC2。即 sin ⁡ ( y − C 1 ) = C 2 e x \sin(y-C_1)=C_2e^x sin(yC1)=C2ex,也可写成 y = arcsin ⁡ ( C 2 e x ) + C 1 y=\arcsin(C_2e^x)+C_1 y=arcsin(C2ex)+C1
  由于当 C 2 = 0 C_2=0 C2=0时, y = C 1 y=C_1 y=C1,故前面所得的解 y ≡ y\equiv y也包含在这个通解内。(这道题主要利用了降阶的方法

2.求下列各微分方程满足所给初值条件的特解:

(1) y 3 y ′ ′ + 1 = 0 , y ∣ x = 1 = 1 , y ′ ∣ x = 1 = 0 ; y^3y''+1=0,y\biggm\vert_{x=1}=1,y'\biggm\vert_{x=1}=0; y3y+1=0,yx=1=1,yx=1=0;

  将原方程写成 y ′ ′ + 1 y 3 = 0 y''+\cfrac{1}{y^3}=0 y+y31=0,两端乘以 2 y ′ 2y' 2y,得
2 y ′ y ′ ′ + 2 y ′ y 3 = 0. 2y'y''+\cfrac{2y'}{y^3}=0. 2yy+y32y=0.
  即 ( y ′ 2 − 2 y ′ y 3 ) ′ = 0 \left(y'^2-\cfrac{2y'}{y^3}\right)'=0 (y2y32y)=0,由此得
y ′ 2 − 1 y 2 = C 1 . y'^2-\cfrac{1}{y^2}=C_1. y2y21=C1.
  代入初值条件: y = 1 y=1 y=1 y ′ = 0 y'=0 y=0,得 C 1 = − 1 C_1=-1 C1=1,故有
y ′ 2 = 1 y 2 − 1 = 1 − y 2 y 2 , y ′ = ± 1 − y 2 y . y'^2=\cfrac{1}{y^2}-1=\cfrac{1-y^2}{y^2},\\ y'=\pm\cfrac{\sqrt{1-y^2}}{y}. y2=y211=y21y2,y=±y1y2 .
  分离变量,得
y d y 1 − y 2 = ± d x . \cfrac{y\mathrm{d}y}{\sqrt{1-y^2}}=\pm\mathrm{d}x. 1y2 ydy=±dx.
  积分得 − 1 − y 2 = ± x + C 2 -\sqrt{1-y^2}=\pm x+C_2 1y2 =±x+C2。代入初值条件: x = 1 x=1 x=1 y = 1 y=1 y=1,得 C = ∓ 1 C=\mp1 C=1。于是有
− 1 − y 2 = ± ( x − 1 ) . -\sqrt{1-y^2}=\pm(x-1). 1y2 =±(x1).
  两边平方,得 x 2 + y 2 = 2 x x^2+y^2=2x x2+y2=2x。由于在点 x = 1 x=1 x=1处, y = 1 y=1 y=1,故在 x = 1 x=1 x=1的某邻域内 y > 0 y>0 y>0,因而特解可以表示为
y = 2 x − x 2 . y=\sqrt{2x-x^2}. y=2xx2 .
这道题判断正负主要根据条件的点所给的正负

(4) y ′ ′ = e 2 y ; y ∣ x = 0 = y ′ ∣ x = 0 = 0 ; y''=e^{2y};y\biggm\vert_{x=0}=y'\biggm\vert_{x=0}=0; y=e2y;yx=0=yx=0=0;

  在原方程两端同乘以 2 y ′ 2y' 2y,得 2 y ′ y ′ ′ = 2 y ′ e 2 y 2y'y''=2y'e^{2y} 2yy=2ye2y,即 ( y ′ 2 ) ′ = ( e 2 y ) ′ (y'^2)'=(e^{2y})' (y2)=(e2y),积分得
y ′ 2 = e 2 y + C 1 . y'^2=e^{2y}+C_1. y2=e2y+C1.
  代入初值条件: x = 0 x=0 x=0 y = y ′ = 0 y=y'=0 y=y=0,得 C 1 = − 1 C_1=-1 C1=1,从而有
y ′ = ± e 2 y − 1 . y'=\pm\sqrt{e^{2y}-1}. y=±e2y1 .
  分离变量后积分
∫ d y e 2 y − 1 = ± ∫ d x . \displaystyle\int\cfrac{\mathrm{d}y}{\sqrt{e^{2y}-1}}=\pm\displaystyle\int\mathrm{d}x. e2y1 dy=±dx.
  即
∫ d ( e − y ) e 2 y − 1 = ∓ ∫ d x . \displaystyle\int\cfrac{\mathrm{d}(e^{-y})}{\sqrt{e^{2y}-1}}=\mp\displaystyle\int\mathrm{d}x. e2y1 d(ey)=dx.
  得 arcsin ⁡ ( e − y ) = ∓ x + C 2 \arcsin(e^{-y})=\mp x+C_2 arcsin(ey)=x+C2。代入初值条件: x = 0 x=0 x=0 y = 0 y=0 y=0,得 C 2 = π 2 C_2=\cfrac{\pi}{2} C2=2π。于是得特解
e − y = sin ⁡ ( π 2 ± x ) = cos ⁡ x . e^{-y}=\sin\left(\cfrac{\pi}{2}\pm x\right)=\cos x. ey=sin(2π±x)=cosx.
  即 y = − ln ⁡ cos ⁡ x = ln ⁡ sec ⁡ x y=-\ln\cos x=\ln\sec x y=lncosx=lnsecx。(这道题利用三角替换公式求解

(5) y ′ ′ = 3 y , y ∣ x = 0 = 1 , y ′ ∣ x = 0 = 2 ; y''=3\sqrt{y},y\biggm\vert_{x=0}=1,y'\biggm\vert_{x=0}=2; y=3y ,yx=0=1,yx=0=2;

  在原方程两端同乘以 2 y ′ 2y' 2y,得 2 y ′ y ′ ′ = 6 y ′ y 2y'y''=6y'\sqrt{y} 2yy=6yy ,即 ( y ′ 2 ) ′ = ( 4 y 3 2 ) ′ (y'^2)'=(4y^{\frac{3}{2}})' (y2)=(4y23),积分得 y ′ 2 = 4 y 3 2 y'^2=4y^{\frac{3}{2}} y2=4y23。代入初值条件 x = 0 x=0 x=0 y = 1 y=1 y=1 y ′ = 2 y'=2 y=2,得 C 1 = 0 C_1=0 C1=0,从而有 y ′ = ± 2 y 3 4 y'=\pm2y^{\frac{3}{4}} y=±2y43。并由于 y ′ ∣ x = 0 = 2 y'\biggm\vert_{x=0}=2 yx=0=2,故取 y ′ = 2 y 3 4 y'=2y^{\frac{3}{4}} y=2y43。分离变量后积分 ∫ d y y 3 4 = 2 ∫ d x \displaystyle\int\cfrac{\mathrm{d}y}{y^{\frac{3}{4}}}=2\displaystyle\int\mathrm{d}x y43dy=2dx 4 y 1 4 = 2 x + C 2 4y^{\frac{1}{4}}=2x+C_2 4y41=2x+C2。代入初值条件: x = 0 x=0 x=0 y = 1 y=1 y=1,得 C 2 = 4 C_2=4 C2=4,于是得特解
y = ( 1 4 + 1 ) 4 . y=\left(\cfrac{1}{4}+1\right)^4. y=(41+1)4.
这道题采用凑整的方法

(6) y ′ ′ + ( y ′ ) 2 = 1 , y ∣ x = 0 = 0 , y ′ ∣ x = 0 = 0 ; y''+(y')^2=1,y\biggm\vert_{x=0}=0,y'\biggm\vert_{x=0}=0; y+(y)2=1,yx=0=0,yx=0=0;

  令 y ′ = p y'=p y=p,则 y ′ ′ = p d p d y y''=p\cfrac{\mathrm{d}p}{\mathrm{d}y} y=pdydp,原方程变为 p d p d y + p 2 = 1 p\cfrac{\mathrm{d}p}{\mathrm{d}y}+p^2=1 pdydp+p2=1。分离变量,得
p d p 1 − p 2 = d y . \cfrac{p\mathrm{d}p}{1-p^2}=\mathrm{d}y. 1p2pdp=dy.
  由初值条件: y = 0 y=0 y=0 p = 0 p=0 p=0,积分
∫ 0 p p d p 1 − p 2 = ∫ 0 y d y . \displaystyle\int^p_0\cfrac{p\mathrm{d}p}{1-p^2}=\displaystyle\int^y_0\mathrm{d}y. 0p1p2pdp=0ydy.
  得
− 1 2 ln ⁡ ( 1 − p 2 ) = y . -\cfrac{1}{2}\ln(1-p^2)=y. 21ln(1p2)=y.
  即 p = ± 1 − e − 2 y p=\pm\sqrt{1-e^{-2y}} p=±1e2y 。又分离变量,得
d y 1 − e − 2 y = ± d x . \cfrac{\mathrm{d}y}{\sqrt{1-e^{-2y}}}=\pm\mathrm{d}x. 1e2y dy=±dx.
  由初值条件: x = 0 x=0 x=0 y = 0 y=0 y=0,积分
∫ 0 y d y 1 − e − 2 y = ± ∫ 0 x d x , ∫ 0 y d ( e y ) 1 − e − 2 y = ± ∫ 0 x d x . \displaystyle\int^y_0\cfrac{\mathrm{d}y}{\sqrt{1-e^{-2y}}}=\pm\displaystyle\int^x_0\mathrm{d}x,\\ \displaystyle\int^y_0\cfrac{\mathrm{d}(e^y)}{\sqrt{1-e^{-2y}}}=\pm\displaystyle\int^x_0\mathrm{d}x. 0y1e2y dy=±0xdx,0y1e2y d(ey)=±0xdx.
  得
ln ⁡ ( e y + e 2 y − 1 ) = ± x . \ln(e^y+\sqrt{e^{2y}-1})=\pm x. ln(ey+e2y1 )=±x.
  即
e y = e x + e − x 2 . e^y=\cfrac{e^x+e^{-x}}{2}. ey=2ex+ex.
  或写成
y = ln ⁡ e x + e − x 2 . y=\ln\cfrac{e^x+e^{-x}}{2}. y=ln2ex+ex.
这道题利用积分结果一致求解

4.设有一质量为 m m m的物体,在空中由静止开始下落,如果空气阻力为 R = c v R=cv R=cv(其中 c c c为常数, v v v为物体运动的速度),试求物体下落的距离 s s s与时间 t t t的函数关系。

  根据牛顿第二定律,有关系式
m d 2 s d t 2 = m g − c d s d t . m\cfrac{\mathrm{d}^2s}{\mathrm{d}t^2}=mg-c\cfrac{\mathrm{d}s}{\mathrm{d}t}. mdt2d2s=mgcdtds.
  并依据题设条件,得初值问题
d 2 s d t 2 = g − c m d s d t , s ∣ t = 0 = 0 , d s d t ∣ t = 0 = 0. \cfrac{\mathrm{d}^2s}{\mathrm{d}t^2}=g-\cfrac{c}{m}\cfrac{\mathrm{d}s}{\mathrm{d}t},\quad s\biggm\vert_{t=0}=0,\quad\cfrac{\mathrm{d}s}{\mathrm{d}t}\biggm\vert_{t=0}=0. dt2d2s=gmcdtds,st=0=0,dtdst=0=0.
  令 d s d t = v \cfrac{\mathrm{d}s}{\mathrm{d}t}=v dtds=v,方程成为 d v d t = g − c m v \cfrac{\mathrm{d}v}{\mathrm{d}t}=g-\cfrac{c}{m}v dtdv=gmcv,分离变量后积分
∫ d v g − c m v = ∫ d t . \displaystyle\int\cfrac{\mathrm{d}v}{g-\cfrac{c}{m}v}=\displaystyle\int\mathrm{d}t. gmcvdv=dt.
  得
ln ⁡ ( g − c m v ) = − c m t + C 1 . \ln\left(g-\cfrac{c}{m}v\right)=-\cfrac{c}{m}t+C_1. ln(gmcv)=mct+C1.
  代入初值条件 v ∣ t = 0 = 0 v|_{t=0}=0 vt=0=0,得 C 1 = ln ⁡ g C_1=\ln g C1=lng。于是有
v = d s d t = m g c ( 1 − e − c m t ) . v=\cfrac{\mathrm{d}s}{\mathrm{d}t}=\cfrac{mg}{c}(1-e^{-\frac{c}{m}t}). v=dtds=cmg(1emct).
  积分得
s = m g c ( t + m c e − c m t ) + C 2 . s=\cfrac{mg}{c}\left(t+\cfrac{m}{c}e^{-\frac{c}{m}t}\right)+C_2. s=cmg(t+cmemct)+C2.
  代入初值条件 s ∣ t = 0 = 0 s|_{t=0}=0 st=0=0,得 C 2 = m 2 g c 2 C_2=\cfrac{m^2g}{c^2} C2=c2m2g。故所求特解(即下落的距离与时间的关系)为
s = m g c ( t + m c e − c m t − m c ) = m g c t + m 2 g c 2 ( e − c m t − 1 ) . \begin{aligned} s&=\cfrac{mg}{c}\left(t+\cfrac{m}{c}e^{-\frac{c}{m}t}-\cfrac{m}{c}\right)\\ &=\cfrac{mg}{c}t+\cfrac{m^2g}{c^2}(e^{-\frac{c}{m}t}-1). \end{aligned} s=cmg(t+cmemctcm)=cmgt+c2m2g(emct1).
这道题主要利用物理公式求解

习题7-6 高阶线性微分方程

本节和以下两节,我们将讨论在实际问题中应用较多的所谓高阶线性微分方程。——高等数学同济版

  本节主要介绍了二阶线性微分方程为主的高阶线性微分方程。

5.已知 y 1 ( x ) = e x y_1(x)=e^x y1(x)=ex是齐次线性方程 ( 2 x − 1 ) y ′ ′ − ( 2 x + 1 ) y ′ + 2 y = 0 (2x-1)y''-(2x+1)y'+2y=0 (2x1)y(2x+1)y+2y=0的一个解,求此方程的通解。

  设 y 2 ( x ) = y 1 u = e x u y_2(x)=y_1u=e^xu y2(x)=y1u=exu是方程的解,则 y 2 ′ = e x ( u + u ′ ) y'_2=e^x(u+u') y2=ex(u+u) y 2 ′ ′ = e x ( u + 2 u ′ + u ′ ′ ) y''_2=e^x(u+2u'+u'') y2=ex(u+2u+u),代入方程并整理,得
e x [ ( 2 x − 1 ) u ′ ′ + ( 2 x − 3 ) u ′ ] = 0. e^x[(2x-1)u''+(2x-3)u']=0. ex[(2x1)u+(2x3)u]=0.
  即
( 2 x − 1 ) u ′ ′ + ( 2 x − 3 ) u ′ = 0. (2x-1)u''+(2x-3)u'=0. (2x1)u+(2x3)u=0.
  令 u ′ = p u'=p u=p,则 u ′ ′ = p ′ u''=p' u=p,且上式成为
( 2 x − 1 ) p ′ + ( 2 x − 3 ) p = 0. (2x-1)p'+(2x-3)p=0. (2x1)p+(2x3)p=0.
  分离变量后积分
∫ d p p = − ∫ 2 x − 3 2 x − 1 d x \displaystyle\int\cfrac{\mathrm{d}p}{p}=-\displaystyle\int\cfrac{2x-3}{2x-1}\mathrm{d}x pdp=2x12x3dx
  得
ln ⁡ ∣ p ∣ = − x + ln ⁡ ∣ 2 x − 1 ∣ + ln ⁡ C . \ln|p|=-x+\ln|2x-1|+\ln C. lnp=x+ln2x1+lnC.
  取 C = 1 C=1 C=1,即 p = ( 2 x − 1 ) e − x p=(2x-1)e^{-x} p=(2x1)ex。再积分得
u = ∫ ( 2 x − 1 ) e − x d x = − [ ( 2 x − 1 ) e − x + 2 e − x + C 0 ] . u=\displaystyle\int(2x-1)e^{-x}\mathrm{d}x=-[(2x-1)e^{-x}+2e^{-x}+C_0]. u=(2x1)exdx=[(2x1)ex+2ex+C0].
  取 C 0 = 0 C_0=0 C0=0,即 u = − ( 2 x + 1 ) e − x u=-(2x+1)e^{-x} u=(2x+1)ex,故
y 2 = e x u = − ( 2 x + 1 ) . y_2=e^xu=-(2x+1). y2=exu=(2x+1).
   y 2 y_2 y2 y 1 y_1 y1线性无关,故原方程的通解为
y = C 1 ( 2 x + 1 ) + C 2 e x . y=C_1(2x+1)+C_2e^x. y=C1(2x+1)+C2ex.
这道题利用常数变易法求解

6.已知 y 1 ( x ) = x y_1(x)=x y1(x)=x是齐次线性方程 x 2 y ′ ′ − 2 x y ′ + 2 y = 0 x^2y''-2xy'+2y=0 x2y2xy+2y=0的一个解,求非齐次线性方程 x 2 y ′ ′ − 2 x y ′ + 2 y = 2 x 3 x^2y''-2xy'+2y=2x^3 x2y2xy+2y=2x3的通解。

  设 y 2 = y 1 u = x u y_2=y_1u=xu y2=y1u=xu是非齐次线性方程的解,则 y 2 ′ = u + x u ′ y'_2=u+xu' y2=u+xu y 2 ′ ′ = 2 u ′ + x u ′ ′ y''_2=2u'+xu'' y2=2u+xu,代入方程并整理,得
u ′ ′ = 0. u''=0. u=0.
  不妨取 u = x u=x u=x,则 y 2 = y 1 u = x 2 y_2=y_1u=x^2 y2=y1u=x2,且 y 2 y_2 y2 y 1 y_1 y1线性无关。
  将非齐次方程化为标准型
y ′ ′ − 2 x y ′ + 2 x 2 y = 2 x . y''-\cfrac{2}{x}y'+\cfrac{2}{x^2}y=2x. yx2y+x22y=2x.
  则它的通解为
y = C 1 y 1 + C 2 y 2 − y 1 ∫ y 2 f W d x + y 2 ∫ y 1 f W d x . y=C_1y_1+C_2y_2-y_1\displaystyle\int\cfrac{y_2f}{W}\mathrm{d}x+y_2\displaystyle\int\cfrac{y_1f}{W}\mathrm{d}x. y=C1y1+C2y2y1Wy2fdx+y2Wy1fdx.
  其中 f = 2 x , W = ∣ y 1 y 2 y 1 ′ y 2 ′ ∣ = ∣ x x 2 1 2 x ∣ = x 2 . f=2x,W=\begin{vmatrix}y_1&y_2\\y'_1&y'_2\end{vmatrix}=\begin{vmatrix}x&x^2\\1&2x\end{vmatrix}=x^2. f=2x,W=y1y1y2y2=x1x22x=x2.
  故
y = C 1 x + C 2 x 2 − x ∫ 2 x 3 x 2 d x + x 2 ∫ 2 x 2 x 2 d x = C 1 x + C 2 x 2 + x 3 . \begin{aligned} y&=C_1x+C_2x^2-x\displaystyle\int\cfrac{2x^3}{x^2}\mathrm{d}x+x^2\displaystyle\int\cfrac{2x^2}{x^2}\mathrm{d}x\\ &=C_1x+C_2x^2+x^3. \end{aligned} y=C1x+C2x2xx22x3dx+x2x22x2dx=C1x+C2x2+x3.
这道题主要考察在高阶导数为0的情况

7.已知齐次线性方程 y ′ ′ + y = 0 y''+y=0 y+y=0的通解为 Y ( x ) = C 1 cos ⁡ x + C 2 sin ⁡ x Y(x)=C_1\cos x+C_2\sin x Y(x)=C1cosx+C2sinx,求非齐次线性方程 y ′ ′ + y = sec ⁡ x y''+y=\sec x y+y=secx的通解。

  由题设知, y 1 = cos ⁡ x y_1=\cos x y1=cosx y 2 = sin ⁡ x y_2=\sin x y2=sinx都是齐次方程 y ′ ′ + y = 0 y''+y=0 y+y=0的解,且 y 1 y_1 y1 y 2 y_2 y2线性无关,则非齐次方程 y ′ ′ + y = sec ⁡ x y''+y=\sec x y+y=secx的通解为
y = C 1 cos ⁡ x + C 2 sin ⁡ x − y 1 ∫ y 2 f W d x + y 2 ∫ y 1 f W d x . y=C_1\cos x+C_2\sin x-y_1\displaystyle\int\cfrac{y_2f}{W}\mathrm{d}x+y_2\displaystyle\int\cfrac{y_1f}{W}\mathrm{d}x. y=C1cosx+C2sinxy1Wy2fdx+y2Wy1fdx.
  其中
f = sec ⁡ x , W = ∣ y 1 y 2 y 1 ′ y 2 ′ ∣ = ∣ cos ⁡ x sin ⁡ x − sin ⁡ x cos ⁡ x ∣ = 1. f=\sec x,\quad W=\begin{vmatrix}y_1&y_2\\y'_1&y'_2\end{vmatrix}=\begin{vmatrix}\cos x&\sin x\\-\sin x&\cos x\end{vmatrix}=1. f=secx,W=y1y1y2y2=cosxsinxsinxcosx=1.
  故
y = C 1 cos ⁡ x + C 2 sin ⁡ x − cos ⁡ x ∫ sin ⁡ x cos ⁡ x d x + sin ⁡ x ∫ cos ⁡ x cos ⁡ x d x = C 1 cos ⁡ x + C 2 sin ⁡ x − cos ⁡ x ln ⁡ ∣ cos ⁡ x ∣ + x sin ⁡ x . \begin{aligned} y&=C_1\cos x+C_2\sin x-\cos x\displaystyle\int\cfrac{\sin x}{\cos x}\mathrm{d}x+\sin x\displaystyle\int\cfrac{\cos x}{\cos x}\mathrm{d}x\\ &=C_1\cos x+C_2\sin x-\cos x\ln|\cos x|+x\sin x. \end{aligned} y=C1cosx+C2sinxcosxcosxsinxdx+sinxcosxcosxdx=C1cosx+C2sinxcosxlncosx+xsinx.
这道题利用已知公式求解

8.已知齐次线性方程 x 2 y ′ ′ − x y ′ + y = 0 x^2y''-xy'+y=0 x2yxy+y=0的通解为 Y ( x ) = C 1 x + C 2 x ln ⁡ ∣ x ∣ Y(x)=C_1x+C_2x\ln|x| Y(x)=C1x+C2xlnx,求线性非齐次方程 x 2 y ′ ′ − x y ′ + y = x x^2y''-xy'+y=x x2yxy+y=x的通解。

  由题设知 y 1 = x y_1=x y1=x y 2 = x ln ⁡ ∣ x ∣ y_2=x\ln|x| y2=xlnx都是齐次方程的解, y 1 y_1 y1 y 2 y_2 y2显然线性无关的。将非齐次方程化为标准型 y ′ ′ − 1 x y ′ + 1 x 2 y = 1 x y''-\cfrac{1}{x}y'+\cfrac{1}{x^2}y=\cfrac{1}{x} yx1y+x21y=x1,则方程的通解为
y = C 1 y 1 + C 2 y 2 − y 1 ∫ y 2 f W d x + y 2 ∫ y 1 f W d x . y=C_1y_1+C_2y_2-y_1\displaystyle\int\cfrac{y_2f}{W}\mathrm{d}x+y_2\displaystyle\int\cfrac{y_1f}{W}\mathrm{d}x. y=C1y1+C2y2y1Wy2fdx+y2Wy1fdx.
  其中
f = 1 x , W = ∣ y 1 y 2 y 1 ′ y 2 ′ ∣ = ∣ x x ln ⁡ ∣ x ∣ 1 ln ⁡ ∣ x ∣ + 1 ∣ = 1. f=\cfrac{1}{x},\quad W=\begin{vmatrix}y_1&y_2\\y'_1&y'_2\end{vmatrix}=\begin{vmatrix}x&x\ln|x|\\1&\ln|x|+1\end{vmatrix}=1. f=x1,W=y1y1y2y2=x1xlnxlnx+1=1.
  因
∫ y 2 f W d x = ∫ ln ⁡ ∣ x ∣ x d x = 1 2 ln ⁡ 2 ∣ x ∣ , ∫ y 1 f W d x = ∫ 1 x d x = ln ⁡ ∣ x ∣ . \begin{aligned} \displaystyle\int\cfrac{y_2f}{W}\mathrm{d}x&=\displaystyle\int\cfrac{\ln|x|}{x}\mathrm{d}x=\cfrac{1}{2}\ln^2|x|,\\ \displaystyle\int\cfrac{y_1f}{W}\mathrm{d}x&=\displaystyle\int\cfrac{1}{x}\mathrm{d}x=\ln|x|. \end{aligned} Wy2fdxWy1fdx=xlnxdx=21ln2x,=x1dx=lnx.
  故非齐次方程的通解为
y = C 1 x + C 2 x ln ⁡ ∣ x ∣ − x 1 2 ln ⁡ 2 ∣ x ∣ + x ln ⁡ ∣ x ∣ ⋅ ln ⁡ ∣ x ∣ = C 1 x + C 2 x ln ⁡ ∣ x ∣ + x 2 ln ⁡ 2 ∣ x ∣ . \begin{aligned} y&=C_1x+C_2x\ln|x|-x\cfrac{1}{2}\ln^2|x|+x\ln|x|\cdot\ln|x|\\ &=C_1x+C_2x\ln|x|+\cfrac{x}{2}\ln^2|x|. \end{aligned} y=C1x+C2xlnxx21ln2x+xlnxlnx=C1x+C2xlnx+2xln2x.
这道题利用行列式和已知公式求解

习题7-7 常系数齐次线性微分方程

  这节介绍了常系数齐次线性微分方程的解法即特征方程法。

习题7-8 常系数非齐次线性微分方程

本节着重讨论二阶常系数非齐次线性微分方程的解法,并对 n n n阶方程的解法作必要的说明。——高等数学同济版

  这节主要介绍了常系数非齐次线性微分方程的解法。

6.设函数 φ ( x ) \varphi(x) φ(x)连续,且满足 φ ( x ) = e x + ∫ 0 x t φ ( t ) d t − x ∫ 0 x φ ( t ) d t , \varphi(x)=e^x+\displaystyle\int^x_0t\varphi(t)\mathrm{d}t-x\displaystyle\int^x_0\varphi(t)\mathrm{d}t, φ(x)=ex+0xtφ(t)dtx0xφ(t)dt, φ ( x ) \varphi(x) φ(x)

  由所给方程可得 φ ( 0 ) = 1 \varphi(0)=1 φ(0)=1,在该方程两端对 x x x求导,得
φ ′ ( x ) = e x + x φ ( x ) − ∫ 0 x φ ( t ) d t − x φ ( x ) . \varphi'(x)=e^x+x\varphi(x)-\displaystyle\int^x_0\varphi(t)\mathrm{d}t-x\varphi(x). φ(x)=ex+xφ(x)0xφ(t)dtxφ(x).
  即
φ ′ ( x ) = e x + − ∫ 0 x φ ( t ) d t . \varphi'(x)=e^x+-\displaystyle\int^x_0\varphi(t)\mathrm{d}t. φ(x)=ex+0xφ(t)dt.
  可见 φ ′ ( 0 ) = 1 \varphi'(0)=1 φ(0)=1
  又在方程 φ ′ ( x ) = e x + − ∫ 0 x φ ( t ) d t \varphi'(x)=e^x+-\displaystyle\int^x_0\varphi(t)\mathrm{d}t φ(x)=ex+0xφ(t)dt的两端对 x x x求导,得
φ ′ ′ ( x ) = e x − φ ( x ) . \varphi''(x)=e^x-\varphi(x). φ(x)=exφ(x).
  若记 φ ( x ) = 0 \varphi(x)=0 φ(x)=0,则有初值问题
{ y ′ ′ + y = e x , y ∣ x = 0 = 1 , y ′ ∣ x = 0 = 1. \begin{cases} y''+y=e^x,\\ y|_{x=0}=1,\quad y'|_{x=0}=1. \end{cases} {y+y=ex,yx=0=1,yx=0=1.
  上述非齐次方程对应的齐次方程的特征方程为 r 2 + 1 = 0 r^2+1=0 r2+1=0,解得 r 1 , 2 = ± i r_{1,2}=\pm i r1,2=±i,而 f ( x ) = e x f(x)=e^x f(x)=ex λ = 1 \lambda=1 λ=1不是方程的特征根,故令 y ∗ = A e x y^*=Ae^x y=Aex是上述方程的特解,代入上式并消去 e x e^x ex,得 A = 1 2 A=\cfrac{1}{2} A=21,于是上式有通解
y = C 1 cos ⁡ x + C 2 sin ⁡ x + 1 2 e x . y=C_1\cos x+C_2\sin x+\cfrac{1}{2}e^x. y=C1cosx+C2sinx+21ex.
  且有
y ′ = − C 1 sin ⁡ x + C 2 cos ⁡ x + 1 2 e x . y'=-C_1\sin x+C_2\cos x+\cfrac{1}{2}e^x. y=C1sinx+C2cosx+21ex.
  代入初值条件 x = 0 x=0 x=0 y = 1 y=1 y=1 y ′ = 1 y'=1 y=1,有
{ C 1 + 1 2 = 1 ′ C 2 + 1 2 = 1. \begin{cases} C_1+\cfrac{1}{2}=1'\\ C_2+\cfrac{1}{2}=1. \end{cases} C1+21=1C2+21=1.
  即 C 1 = C 2 = 1 2 C_1=C_2=\cfrac{1}{2} C1=C2=21。于是得
y = φ ( x ) = 1 2 ( cos ⁡ x + sin ⁡ x + e x ) . y=\varphi(x)=\cfrac{1}{2}(\cos x+\sin x+e^x). y=φ(x)=21(cosx+sinx+ex).
这道题主要利用了常系数微分方程的求解

习题7-9 欧拉方程

  这节主要介绍了欧拉方程的求解方法。

习题7-10 常系数线性微分方程组解法举例

  本节主要介绍了常系数线性微分方程组的两种解法。

写在最后

  如果觉得文章不错就点个赞吧。另外,如果有不同的观点,欢迎留言或私信。
   欢迎非商业转载,转载请注明出处。
  本文未完,总习题七见《第七章 微分方程(二)》,传送门在这里。

  • 1
    点赞
  • 7
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值