在许多问题中,往往不能直接找出所需要的函数关系,但是根据问题所提供的情况,有时可以列出含有要找的函数及其导数的关系式。这样的关系式就是所谓的
微分方程
。微分方程建立以后,对它进行研究,找出未知函数来,这就是解微分方程
。本章主要介绍微分方程的一些基本概念和几种常用的微分方程的解法。——高等数学同济版
本节主要介绍了微分方程的基本概念。
解 设雪堆在时刻 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=r0−kt。又 V ∣ t = 3 = 1 8 V ∣ t = 0 V\biggm\vert_{t=3}=\cfrac{1}{8}V\biggm\vert_{t=0} V∣∣∣∣t=3=81V∣∣∣∣t=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π(r0−3k)3=81⋅32π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=r0−61r0t.
因雪堆全部融化时, r = 0 r=0 r=0,故得 t = 6 t=6 t=6,即雪堆全部融化需6小时。(这道题主要难点在体积融化率即体积融化速率的理解上)
本节主要介绍可分离变量的微分方程的解法。
本节主要介绍齐次方程的解法。
解 原方程可写成 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 1−3y2x2+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. 1−3u2+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}. u2−12udu=ydy.
积分得
ln ∣ u 2 − 1 ∣ = ln ∣ y ∣ + ln C 1 . \ln|u^2-1|=\ln|y|+\ln C_1. ln∣u2−1∣=ln∣y∣+lnC1.
即
u 2 − 1 = C y . u^2-1=Cy. u2−1=Cy.
代入 y = x y y=\cfrac{x}{y} y=yx并整理,得通解 x 2 − y 2 = C y 3 x^2-y^2=Cy^3 x2−y2=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=y2−x2.
本节主要介绍了一阶线性微分方程的求解方法。
解 令 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 ln∣lnu∣=ln∣x∣+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。(这道题利用变量换元时对所换的元的导数可以进行整体代换)
解 将原方程写成 y ′ = ( y + sin x − 1 ) 2 − cos x y'=(y+\sin x-1)^2-\cos x y′=(y+sinx−1)2−cosx,令 u = y + sin x − 1 u=y+\sin x-1 u=y+sinx−1,则 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+sinx−1,得原方程的通解
y = 1 − sin x − 1 x + C . y=1-\sin x-\cfrac{1}{x+C}. y=1−sinx−x+C1.
(这道题说明换元时可以先配方再换元)
解 原方程改写成 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=x2xdxdu−u,且原方程变为
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)(xdxdu−u)=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 −2u21−u1+ln∣u∣=ln∣x∣+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). 2x2y2ln∣y∣−2xy−1=Cx2y2(C=2C1).
(这道题说明可以在换元前可以先乘以某一相同变量凑整再换元)
从这一节其我们将讨论二阶及二阶以上的微分方程,即所谓
高阶微分方程
。——高等数学同济版
本节主要介绍了三种高阶微分方程的降阶方法。
解 令 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 y3pdydp−1=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=±C1−y21=±∣y∣1C1y2−1.
分离变量,得
∣ y ∣ d y C 1 y 2 − 1 = ± d x . \cfrac{|y|\mathrm{d}y}{\sqrt{C_1y^2-1}}=\pm\mathrm{d}x. C1y2−1∣y∣dy=±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)∫C1y2−1ydy=±∫dx,sgn(y)C1y2−1=±C1x+C2.
两边平方,得 C 1 y 2 − 1 = ( C 1 x + C 2 ) 2 C_1y^2-1=(C_1x+C_2)^2 C1y2−1=(C1x+C2)2。(这道题主要利用绝对值的函数表示来求解)
解 方程 y ′ ′ = 1 y y''=\cfrac{1}{\sqrt{y}} y′′=y1属于 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'. 2y′y′′=2f(y)y′.
即 ( y ′ 2 ) ′ = 2 f ( y ) y ′ (y'^2)'=2f(y)y' (y′2)′=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)]'. (y′2)′=2[F(y)]′.
积分得到降阶的方程 y ′ 2 = 2 F ( y ) + C 1 y'^2=2F(y)+C_1 y′2=2F(y)+C1。
本小题按上述方法求解:用 2 y ′ 2y' 2y′乘方程 y ′ ′ = 1 y y''=\cfrac{1}{\sqrt{y}} y′′=y1的两端,得
2 y ′ y ′ ′ = 2 y ′ y . 2y'y''=\cfrac{2y'}{\sqrt{y}}. 2y′y′′=y2y′.
即 ( y ′ 2 ) ′ = ( 4 y ) ′ (y'^2)'=(4\sqrt{y})' (y′2)′=(4y)′,故 y ′ 2 = 4 y + C 1 ′ y'^2=4\sqrt{y}+C_1' y′2=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+C1dy.
积分,得
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+C1d(y)2=±∫2y+C1ydy=±∫2y+C1(y+C1)−C1d(y)=±[∫y+C1d(y+C1)−C1∫y+C11d(y+C1)]=±[32(y+C1)23−2C1(y+C1)21]+C2.
(这道题主要利用凑整的方法求解)
解 令 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 p≡0,则 y ≡ C y\equiv C y≡C。 y ≡ C y\equiv C y≡C是原方程的解,但不是通解。
若 p ≢ 0 p\not\equiv0 p≡0,由于 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=y−C,即 p = tan ( y − C 1 ) p=\tan(y-C_1) p=tan(y−C1),亦即 cot ( y − C 1 ) d y = d x \cot(y-C_1)\mathrm{d}y=\mathrm{d}x cot(y−C1)dy=dx。积分得 ln sin ( y − C 1 ) = x + ln C 2 \ln\sin(y-C_1)=x+\ln C_2 lnsin(y−C1)=x+lnC2。即 sin ( y − C 1 ) = C 2 e x \sin(y-C_1)=C_2e^x sin(y−C1)=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≡也包含在这个通解内。(这道题主要利用了降阶的方法)
解 将原方程写成 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. 2y′y′′+y32y′=0.
即 ( y ′ 2 − 2 y ′ y 3 ) ′ = 0 \left(y'^2-\cfrac{2y'}{y^3}\right)'=0 (y′2−y32y′)′=0,由此得
y ′ 2 − 1 y 2 = C 1 . y'^2-\cfrac{1}{y^2}=C_1. y′2−y21=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}. y′2=y21−1=y21−y2,y′=±y1−y2.
分离变量,得
y d y 1 − y 2 = ± d x . \cfrac{y\mathrm{d}y}{\sqrt{1-y^2}}=\pm\mathrm{d}x. 1−y2ydy=±dx.
积分得 − 1 − y 2 = ± x + C 2 -\sqrt{1-y^2}=\pm x+C_2 −1−y2=±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). −1−y2=±(x−1).
两边平方,得 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=2x−x2.
(这道题判断正负主要根据条件的点所给的正负)
解 在原方程两端同乘以 2 y ′ 2y' 2y′,得 2 y ′ y ′ ′ = 2 y ′ e 2 y 2y'y''=2y'e^{2y} 2y′y′′=2y′e2y,即 ( y ′ 2 ) ′ = ( e 2 y ) ′ (y'^2)'=(e^{2y})' (y′2)′=(e2y)′,积分得
y ′ 2 = e 2 y + C 1 . y'^2=e^{2y}+C_1. y′2=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′=±e2y−1.
分离变量后积分
∫ d y e 2 y − 1 = ± ∫ d x . \displaystyle\int\cfrac{\mathrm{d}y}{\sqrt{e^{2y}-1}}=\pm\displaystyle\int\mathrm{d}x. ∫e2y−1dy=±∫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. ∫e2y−1d(e−y)=∓∫dx.
得 arcsin ( e − y ) = ∓ x + C 2 \arcsin(e^{-y})=\mp x+C_2 arcsin(e−y)=∓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. e−y=sin(2π±x)=cosx.
即 y = − ln cos x = ln sec x y=-\ln\cos x=\ln\sec x y=−lncosx=lnsecx。(这道题利用三角替换公式求解)
解 在原方程两端同乘以 2 y ′ 2y' 2y′,得 2 y ′ y ′ ′ = 6 y ′ y 2y'y''=6y'\sqrt{y} 2y′y′′=6y′y,即 ( y ′ 2 ) ′ = ( 4 y 3 2 ) ′ (y'^2)'=(4y^{\frac{3}{2}})' (y′2)′=(4y23)′,积分得 y ′ 2 = 4 y 3 2 y'^2=4y^{\frac{3}{2}} y′2=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 y′∣∣∣∣x=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=2∫dx得 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.
(这道题采用凑整的方法)
解 令 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. 1−p2pdp=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. ∫0p1−p2pdp=∫0ydy.
得
− 1 2 ln ( 1 − p 2 ) = y . -\cfrac{1}{2}\ln(1-p^2)=y. −21ln(1−p2)=y.
即 p = ± 1 − e − 2 y p=\pm\sqrt{1-e^{-2y}} p=±1−e−2y。又分离变量,得
d y 1 − e − 2 y = ± d x . \cfrac{\mathrm{d}y}{\sqrt{1-e^{-2y}}}=\pm\mathrm{d}x. 1−e−2ydy=±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. ∫0y1−e−2ydy=±∫0xdx,∫0y1−e−2yd(ey)=±∫0xdx.
得
ln ( e y + e 2 y − 1 ) = ± x . \ln(e^y+\sqrt{e^{2y}-1})=\pm x. ln(ey+e2y−1)=±x.
即
e y = e x + e − x 2 . e^y=\cfrac{e^x+e^{-x}}{2}. ey=2ex+e−x.
或写成
y = ln e x + e − x 2 . y=\ln\cfrac{e^x+e^{-x}}{2}. y=ln2ex+e−x.
(这道题利用积分结果一致求解)
解 根据牛顿第二定律,有关系式
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=mg−cdtds.
并依据题设条件,得初值问题
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=g−mcdtds,s∣∣∣∣t=0=0,dtds∣∣∣∣t=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=g−mcv,分离变量后积分
∫ d v g − c m v = ∫ d t . \displaystyle\int\cfrac{\mathrm{d}v}{g-\cfrac{c}{m}v}=\displaystyle\int\mathrm{d}t. ∫g−mcvdv=∫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(g−mcv)=−mct+C1.
代入初值条件 v ∣ t = 0 = 0 v|_{t=0}=0 v∣t=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(1−e−mct).
积分得
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+cme−mct)+C2.
代入初值条件 s ∣ t = 0 = 0 s|_{t=0}=0 s∣t=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+cme−mct−cm)=cmgt+c2m2g(e−mct−1).
(这道题主要利用物理公式求解)
本节和以下两节,我们将讨论在实际问题中应用较多的所谓高阶线性微分方程。——高等数学同济版
本节主要介绍了二阶线性微分方程为主的高阶线性微分方程。
解 设 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[(2x−1)u′′+(2x−3)u′]=0.
即
( 2 x − 1 ) u ′ ′ + ( 2 x − 3 ) u ′ = 0. (2x-1)u''+(2x-3)u'=0. (2x−1)u′′+(2x−3)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. (2x−1)p′+(2x−3)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=−∫2x−12x−3dx
得
ln ∣ p ∣ = − x + ln ∣ 2 x − 1 ∣ + ln C . \ln|p|=-x+\ln|2x-1|+\ln C. ln∣p∣=−x+ln∣2x−1∣+lnC.
取 C = 1 C=1 C=1,即 p = ( 2 x − 1 ) e − x p=(2x-1)e^{-x} p=(2x−1)e−x。再积分得
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=∫(2x−1)e−xdx=−[(2x−1)e−x+2e−x+C0].
取 C 0 = 0 C_0=0 C0=0,即 u = − ( 2 x + 1 ) e − x u=-(2x+1)e^{-x} u=−(2x+1)e−x,故
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.
(这道题利用常数变易法求解)
解 设 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. y′′−x2y′+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+C2y2−y1∫Wy2fdx+y2∫Wy1fdx.
其中 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=∣∣∣∣y1y1′y2y2′∣∣∣∣=∣∣∣∣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+C2x2−x∫x22x3dx+x2∫x22x2dx=C1x+C2x2+x3.
(这道题主要考察在高阶导数为0的情况)
解 由题设知, 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+C2sinx−y1∫Wy2fdx+y2∫Wy1fdx.
其中
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=∣∣∣∣y1y1′y2y2′∣∣∣∣=∣∣∣∣cosx−sinxsinxcosx∣∣∣∣=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+C2sinx−cosx∫cosxsinxdx+sinx∫cosxcosxdx=C1cosx+C2sinx−cosxln∣cosx∣+xsinx.
(这道题利用已知公式求解)
解 由题设知 y 1 = x y_1=x y1=x与 y 2 = x ln ∣ x ∣ y_2=x\ln|x| y2=xln∣x∣都是齐次方程的解, 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} y′′−x1y′+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+C2y2−y1∫Wy2fdx+y2∫Wy1fdx.
其中
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=∣∣∣∣y1y1′y2y2′∣∣∣∣=∣∣∣∣x1xln∣x∣ln∣x∣+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} ∫Wy2fdx∫Wy1fdx=∫xln∣x∣dx=21ln2∣x∣,=∫x1dx=ln∣x∣.
故非齐次方程的通解为
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