【高等数学基础进阶】常微分方程-part2

文章目录

    • 常考题型与典型例题
      • 微分方程求解
      • 综合题
      • 应用题

常考题型与典型例题

微分方程求解

例7:微分方程 x y ′ + y ( ln ⁡ x − ln ⁡ y ) = 0 xy'+y(\ln x-\ln y)=0 xy+y(lnxlny)=0满足条件 y ( 1 ) = e 3 y(1)=e^{3} y(1)=e3的解为 y = ( ) y=() y=()

原式两边同除 x x x
y ′ = y x ln ⁡ y x y'= \frac{y}{x}\ln \frac{y}{x} y=xylnxy
u = y x u= \frac{y}{x} u=xy
u ′ = u ( ln ⁡ u − 1 ) x ln ⁡ ∣ ln ⁡ u − 1 ∣ = ln ⁡ x + C ln ⁡ u − 1 = C x ln ⁡ y x − 1 = C x \begin{aligned} u'&= \frac{u(\ln u-1)}{x}\\ \ln \left|\ln u-1\right|&=\ln x+C\\ \ln u-1&=Cx\\ \ln \frac{y}{x}-1&=Cx\\ \end{aligned} ulnlnu1lnu1lnxy1=xu(lnu1)=lnx+C=Cx=Cx
y ( 1 ) = e 3 ⇒ C = 2 y(1)=e^{3}\Rightarrow C=2 y(1)=e3C=2,则
ln ⁡ y x − 1 = 2 x y = x e 2 x + 1 \begin{aligned} \ln \frac{y}{x}-1&=2x\\ y&=xe^{2x+1} \end{aligned} lnxy1y=2x=xe2x+1

例8:微分方程 y d x + ( x − 3 y 2 ) d y = 0 ydx+(x-3y^{2})dy=0 ydx+(x3y2)dy=0满足条件 y ∣ x = 1 = 1 y \Big|_{x=1}^{}=1 y x=1=1的解为 y = ( ) y=() y=()

本题可以用全微分方程的方式,但现在暂时未涉及到相关知识,该种方法会在以后补充

观察本题,现在对于 y y y不满足三类的任一一种,一般有两种常用思路:考虑 x , y x,y x,y对调;做变量代换

y d x + ( x − 3 y 2 ) d y = 0 d x d y + x y = 3 y x = e − ∫ 1 y d y [ ∫ 3 y e ∫ 1 y d y d y + C ] = 1 y [ y 3 + C ] 代入 x = 1 , y = 1 ⇒ C = 0 ⇒ x = y 2 ⇒ y = ± x 由于 x = 1 , y = 1 y = x \begin{aligned} ydx+(x-3y^{2})dy&=0\\ \frac{dx}{dy}+ \frac{x}{y}&=3y\\ x&=e^{-\int_{}^{} \frac{1}{y}dy}\left[\int_{}^{}3ye^{\int_{}^{}\frac{1}{y}dy} dy+C\right]\\ &= \frac{1}{y}[y^{3}+C]\\ 代入x=1,y=1\Rightarrow C&=0\\\Rightarrow x&=y^{2}\\ \Rightarrow y&=\pm \sqrt{x}\\ &由于x=1,y=1\\ y&=\sqrt{x} \end{aligned} ydx+(x3y2)dydydx+yxx代入x=1,y=1Cxyy=0=3y=ey1dy[3yey1dydy+C]=y1[y3+C]=0=y2=±x 由于x=1,y=1=x

例8:微分方程 y ′ ′ + 2 y ′ + 3 y = 0 y''+2y'+3y=0 y′′+2y+3y=0的通解为 y = ( ) y=() y=()

特征方程 r 2 + 2 r + 3 = 0 ⇒ r 1 , 2 = − 1 ± 2 i y = e − 1 ( C 1 cos ⁡ 2 x + C 2 sin ⁡ 2 x ) \begin{aligned} 特征方程\quad r^{2}+2r+3&=0\Rightarrow r_{1,2}=-1\pm \sqrt{2}i\\ y&=e^{-1}(C_{1}\cos \sqrt{2}x+C_{2}\sin \sqrt{2}x) \end{aligned} 特征方程r2+2r+3y=0r1,2=1±2 i=e1(C1cos2 x+C2sin2 x)

例9:微分方程 y ′ ′ − 4 y ′ + 8 y = e 2 x ( 1 + cos ⁡ 2 x ) y''-4y'+8y=e^{2x}(1+\cos 2x) y′′4y+8y=e2x(1+cos2x)的特解可设为 y ∗ = ( ) y ^{*}=() y=()

根据定理4

定理4:如果 y 1 ∗ ( x ) , y 2 ∗ y^{*}_{1}(x),y^{*}_{2} y1(x),y2分别是方程
y ′ ′ + p ( x ) y ′ + q ( x ) y = f 1 ( x ) y ′ ′ + p ( x ) y ′ + q ( x ) y = f 2 ( x ) \begin{aligned}y''+p(x)y'+q(x)y=f_{1}(x)\\y''+p(x)y'+q(x)y=f_{2}(x)\end{aligned} y′′+p(x)y+q(x)y=f1(x)y′′+p(x)y+q(x)y=f2(x)
的特解,则
y 1 ∗ ( x ) + y 2 ∗ y^{*}_{1}(x)+y^{*}_{2} y1(x)+y2
是方程 y ′ ′ + p ( x ) y ′ + q ( x ) y = f 1 ( x ) + f 2 ( x ) y''+p(x)y'+q(x)y=f_{1}(x)+f_{2}(x) y′′+p(x)y+q(x)y=f1(x)+f2(x)的一个特解

原微分方程可分解为
y ′ ′ − 4 y ′ + 8 y = e 2 x y ′ ′ − 4 y ′ + 8 y = e 2 x cos ⁡ 2 x \begin{aligned} y''-4y'+8y&=e^{2x}\\ y''-4y'+8y&=e^{2x}\cos 2x \end{aligned} y′′4y+8yy′′4y+8y=e2x=e2xcos2x
因此
特征方程 r 2 − 4 r + 8 = 0 r 1 , 2 = 4 + ± − 16 2 = 2 ± 2 i 对于 y ′ ′ − 4 y ′ + 8 y = e 2 x y ∗ = C 1 e 2 x 对于 y ′ ′ − 4 y ′ + 8 y = e 2 x cos ⁡ 2 x y ∗ = x e 2 x [ C 2 cos ⁡ 2 x + C 3 sin ⁡ 2 x ] y ∗ = C 1 e 2 x + x e 2 x [ C 2 cos ⁡ 2 x + C 3 sin ⁡ 2 x ] \begin{aligned} 特征方程\quad r^{2}-4r+8&=0\\ r_{1,2}&=\frac{4+\pm \sqrt{-16}}{2}=2 \pm 2i\\ 对于y''-4y'+8y&=e^{2x}\\ y^{*}&=C_{1}e^{2x}\\ 对于y''-4y'+8y&=e^{2x}\cos 2x \\ y^{*}&=xe^{2x}[C_{2}\cos 2x+C_{3}\sin 2x]\\ y^{*}&=C_{1}e^{2x}+xe^{2x}[C_{2}\cos 2x+C_{3}\sin 2x] \end{aligned} 特征方程r24r+8r1,2对于y′′4y+8yy对于y′′4y+8yyy=0=24+±16 =2±2i=e2x=C1e2x=e2xcos2x=xe2x[C2cos2x+C3sin2x]=C1e2x+xe2x[C2cos2x+C3sin2x]

例10:设 y = 1 2 e 2 x + ( x − 1 3 ) e x y= \frac{1}{2}e^{2x}+(x- \frac{1}{3})e^{x} y=21e2x+(x31)ex是二阶常系数非齐次线性方程微分方程 y ′ ′ + a y ′ + b y = c e x y''+ay'+by=ce^{x} y′′+ay+by=cex的一个特解,求 a , b , c a,b,c a,b,c

y = 1 2 e 2 x + ( x − 1 3 ) e x y= \frac{1}{2}e^{2x}+(x- \frac{1}{3})e^{x} y=21e2x+(x31)ex y ′ ′ + a y ′ + b y = c e x y''+ay'+by=ce^{x} y′′+ay+by=cex的一个特解,根据定理2

定理2:如果 y ∗ y^{*} y是非齐次方程的一个特解, y 1 ( x ) y_{1}(x) y1(x) y 2 ( x ) y_{2}(x) y2(x)是齐次方程的两个线性无关的特解,则
y = C 1 y 1 ( x ) + C 2 y 2 ( x ) + y ∗ ( x ) y=C_{1}y_{1}(x)+C_{2}y_{2}(x)+y^{*}(x) y=C1y1(x)+C2y2(x)+y(x)
是非齐次方程的通解
y = 1 2 e 2 x + ( x − 1 3 ) e x y= \frac{1}{2}e^{2x}+(x- \frac{1}{3})e^{x} y=21e2x+(x31)ex有一个非齐次的特解,两个齐次的特解,因此需要进行判断
判断前,建议先分成三项
以下 是判断过程,可以省略

y = 1 2 e 2 x + x e x − 1 3 e x y=\frac{1}{2}e^{2x}+xe^{x}- \frac{1}{3}e^{x} y=21e2x+xex31ex
如果 1 2 e 2 x \frac{1}{2}e^{2x} 21e2x是非齐次方程的一个特解,不管系数,代入 y ′ ′ + a y ′ + b y = c e x y''+ay'+by=ce^{x} y′′+ay+by=cex
A e 2 x = B e x Ae^{2x}=Be^{x} Ae2x=Bex
显然不成立,因此 1 2 e 2 x \frac{1}{2}e^{2x} 21e2x是齐次方程的一个特解
如果 x e x xe^{x} xex是非齐次方程的一个特解,显然代入不好判断,考虑其实齐次方程的一个特解的可能性,如果 x e x xe^{x} xex是齐次方程的一个特解,那么 1 1 1一定是对应特征方程的二重根,由于之前已经确定齐次方程的一个特解为 1 2 e 2 x \frac{1}{2}e^{2x} 21e2x,因此,其不是齐次方程的一个特解, x e x xe^{x} xex是非齐次方程的一个特解
剩余的 − 1 3 e x - \frac{1}{3}e^{x} 31ex一定是齐次方程的一个特解

可知 y 1 = e 2 x , y 2 = e x y_{1}=e^{2x},y_{2}=e^{x} y1=e2x,y2=ex是齐次方程的两个线性无关的解, y ∗ x e x y^{*}xe^{x} yxex是非齐次方程的一个解,因此齐次方程的特征方程为
( r − 1 ) ( r − 2 ) = 0 ⇒ a = − 3 , b = 2 (r-1)(r-2)=0\Rightarrow a=-3,b=2 (r1)(r2)=0a=3,b=2
y = x e x y=xe^{x} y=xex代入 y ′ ′ − 3 y ′ + 2 y = c e x y''-3y'+2y=ce^{x} y′′3y+2y=cex,可得 c = − 1 c=-1 c=1

综合题

例11:设 y = f ( x ) y=f(x) y=f(x)是微分方程 y ′ ′ − y ′ − e sin ⁡ x = 0 y''-y'-e^{\sin x}=0 y′′yesinx=0的解,且 f ′ ( x 0 ) = 0 f'(x_{0})=0 f(x0)=0,说明 f ( x ) f(x) f(x) x 0 x_{0} x0处取得极小值

x 0 x_{0} x0代入微分方程
f ′ ′ ( x 0 ) = e sin ⁡ x 0 > 0 f''(x_{0})=e^{\sin x_{0}}>0 f′′(x0)=esinx0>0
因此 f ( x ) f(x) f(x) x 0 x_{0} x0处取得极小值

有关微分方程的综合题有时需不要把解求出来

例12:设 y = y ( x ) y=y(x) y=y(x)是二阶常系数微分方程 y ′ ′ + p y ′ + q y = e 3 x y''+py'+qy=e^{3x} y′′+py+qy=e3x满足初始条件 y ( 0 ) = y ′ ( 0 ) = 0 y(0)=y'(0)=0 y(0)=y(0)=0的特解,则当 x → 0 x \to 0 x0时,函数 ln ⁡ ( 1 + x 2 ) y ( x ) \frac{\ln (1+x^{2})}{y(x)} y(x)ln(1+x2)的极限为()

y ′ ′ + p y ′ + q y = e 3 x y''+py'+qy=e^{3x} y′′+py+qy=e3x y ′ ′ ( x ) y''(x) y′′(x)连续且 y ′ ′ ( 0 ) = 1 y''(0)=1 y′′(0)=1

由于 p y ′ , q y , e 3 x py',qy,e^{3x} py,qy,e3x连续,因此 y ′ ′ y'' y′′连续

lim ⁡ x → 0 ln ⁡ ( 1 + x 2 ) y ( x ) = lim ⁡ x → 0 x 2 y ( x ) = lim ⁡ x → 0 2 x y ′ ( x ) = lim ⁡ x → 0 2 y ′ ′ ( x ) = 2 y ′ ′ ( 0 ) = 2 \begin{aligned} \lim\limits_{x \to 0}\frac{\ln (1+x^{2})}{y(x)}&=\lim\limits_{x\to0} \frac{x^{2}}{y(x)}=\lim\limits_{x\to0} \frac{2x}{y'(x)}=\lim\limits_{x\to0} \frac{2}{y''(x)}=\frac{2}{y''(0)}=2 \end{aligned} x0limy(x)ln(1+x2)=x0limy(x)x2=x0limy(x)2x=x0limy′′(x)2=y′′(0)2=2

例13:已知连续函数 f ( x ) f(x) f(x)满足条件 f ( x ) = ∫ 0 3 x f ( t 3 ) d t + e 2 x f(x)=\int_{0}^{3x}f(\frac{t}{3})dt+e^{2x} f(x)=03xf(3t)dt+e2x,求 f ( x ) f(x) f(x)

未知函数在积分里面叫做积分方程,解积分方程的一般方法两边同时求导,化为微分方程

两边同时求导
f ′ ( x ) = 3 f ( x ) + 2 e 2 x f ′ ( x ) − 3 f ( x ) = 2 e 2 x f ( x ) = e 3 x ( − 2 e − x + C ) = C e 3 x − 2 e 2 x (1) \begin{aligned} f'(x)&=3f(x)+2e^{2x}\\ f'(x)-3f(x)&=2e^{2x}\\ f(x)&=e^{3x}(-2e^{-x}+C)=Ce^{3x}-2e^{2x}\quad \text{(1)} \end{aligned} f(x)f(x)3f(x)f(x)=3f(x)+2e2x=2e2x=e3x(2ex+C)=Ce3x2e2x(1)
x = 0 x=0 x=0时,代入题目条件
f ( 0 ) = ∫ 0 0 f ( t 3 ) d t + e 2 x f ( 0 ) = 1 \begin{aligned} f(0)&=\int_{0}^{0}f(\frac{t}{3})dt+e^{2x}\\ f(0)&=1 \end{aligned} f(0)f(0)=00f(3t)dt+e2x=1
f ( 0 ) = 1 f(0)=1 f(0)=1代入 ( 1 ) (1) (1),得 C = 2 C=2 C=2,因此
f ( x ) = 3 e 3 x − 2 e 2 x f(x)=3e^{3x}-2e^{2x} f(x)=3e3x2e2x

例14:设函数 f ( x ) f(x) f(x)连续,且满足 ∫ 0 x f ( x − t ) d t = ∫ 0 x ( x − t ) f ( t ) d t + e − x − 1 \begin{aligned} \int_{0}^{x}f(x-t)dt=\int_{0}^{x}(x-t)f(t)dt+e^{-x}-1\end{aligned} 0xf(xt)dt=0x(xt)f(t)dt+ex1,求 f ( x ) f(x) f(x)的表达式

u = x − t u=x-t u=xt,则 ∫ 0 x f ( x − t ) d t = ∫ 0 x f ( u ) d u \begin{aligned} \int_{0}^{x}f(x-t)dt=\int_{0}^{x}f(u)du\end{aligned} 0xf(xt)dt=0xf(u)du
∫ 0 x f ( u ) d u = x ∫ 0 x f ( t ) d t − ∫ 0 x t f ( t ) d t + e − x − 1 两端求导 f ( x ) = ∫ 0 x f ( t ) d t − e − x 令 x = 0 , f ( 0 ) = − 1 f ′ ( x ) − f ( x ) = e − x f ( x ) = C e x − e − x 2 代入 f ( 0 ) = − 1 , 得 C = − 1 2 f ( x ) = − e x + e − x 2 \begin{aligned} \int_{0}^{x}f(u)du&=x \int_{0}^{x}f(t)dt-\int_{0}^{x}tf(t)dt+e^{-x}-1\\ &两端求导\\ f(x)&=\int_{0}^{x}f(t)dt-e^{-x}\quad 令x=0,f(0)=-1\\ f'(x)-f(x)&=e^{-x}\\ f(x)&=Ce^{x}-\frac{e^{-x}}{2}\\ &代入f(0)=-1,得C=- \frac{1}{2}\\ f(x)&=- \frac{e^{x}+e^{-x}}{2} \end{aligned} 0xf(u)duf(x)f(x)f(x)f(x)f(x)=x0xf(t)dt0xtf(t)dt+ex1两端求导=0xf(t)dtexx=0,f(0)=1=ex=Cex2ex代入f(0)=1,C=21=2ex+ex

应用题

例15:设函数 f ( x ) f(x) f(x)在定义域 I I I上的导数大于零,若对于任意的 x 0 ∈ I x_{0}\in I x0I,曲线 y = f ( x ) y=f(x) y=f(x)在点 ( x 0 , f ( x 0 ) ) (x_{0},f(x_{0})) (x0,f(x0))处的切线与直线 x = x 0 x=x_{0} x=x0 x x x轴所围成区域的面积恒为 4 4 4,且 f ( 0 ) = 2 f(0)=2 f(0)=2,求 f ( x ) f(x) f(x)的表达式

曲线 y = f ( x ) y=f(x) y=f(x)在点 ( x 0 , f ( x 0 ) ) (x_{0},f(x_{0})) (x0,f(x0))处的切线方程为
y − f ( x 0 ) = f ′ ( x 0 ) ( x − x 0 ) 令 y = 0 x = x 0 − f ( x 0 ) f ′ ( x 0 ) \begin{aligned} y-f(x_{0})&=f'(x_{0})(x-x_{0})\\ &令y=0\\ x&=x_{0}-\frac{f(x_{0})}{f'(x_{0})} \end{aligned} yf(x0)x=f(x0)(xx0)y=0=x0f(x0)f(x0)
因此该面积为
S = 1 2 ∣ f ( x 0 ) ∣ f ′ ( x 0 ) ⋅ ∣ f ( x 0 ) ∣ = 4 记 y = f ( x 0 ) , 则 1 2 y 2 = 4 y ′ − 8 y = x + C 代入 y ( 0 ) = 2 C = − 4 带回上式 y = 8 4 − x \begin{aligned} S&= \frac{1}{2}\frac{\left|f(x_{0})\right|}{f'(x_{0})}\cdot |f(x_{0})|=4\\ &记y=f(x_{0}),则\\ \frac{1}{2}y^{2}&=4y'\\- \frac{8}{y}&=x+C\\ &代入y(0)=2\\ C&=-4\quad 带回上式\\ y&=\frac{8}{4-x} \end{aligned} S21y2y8Cy=21f(x0)f(x0)f(x0)=4y=f(x0),=4y=x+C代入y(0)=2=4带回上式=4x8

例16:设非负函数 y = y ( x ) ( x ≥ 0 ) y=y(x)(x \geq 0) y=y(x)(x0)满足微分方程 x y ′ ′ − y ′ + 2 = 0 xy''-y'+2=0 xy′′y+2=0。当曲线 y = y ( x ) y=y(x) y=y(x)过原点时,其与直线 x = 1 x=1 x=1 y = 0 y=0 y=0围成的平面区域 D D D的面积为 2 2 2,求 D D D y y y轴旋转所得旋转体的体积

y ′ = p y'=p y=p,则 y ′ ′ = p ′ y''=p' y′′=p,代入微分方程得
p ′ − 1 x p = − 2 x ( x > 0 ) y ′ = p = e ∫ 1 x d x ( ∫ − 1 x e ∫ − 1 x d x d x + C 1 ) = 2 + C 1 x y = 2 x + 1 2 C 1 x 2 + C 2 \begin{aligned} p'- \frac{1}{x}p&=- \frac{2}{x}\quad (x>0)\\ y'=p&=e^{\int_{}^{} \frac{1}{x}dx}(\int_{}^{} - \frac{1}{x}e^{\int_{}^{}- \frac{1}{x}dx}dx +C_{1})\\ &=2+C_{1}x\\ y&=2x+ \frac{1}{2}C_{1}x^{2}+C_{2} \end{aligned} px1py=py=x2(x>0)=ex1dx(x1ex1dxdx+C1)=2+C1x=2x+21C1x2+C2
由方程 x y ′ ′ − y ′ + 2 = 0 ( x ≥ 0 ) xy''-y'+2=0\quad (x \geq 0) xy′′y+2=0(x0),可知 y y y x = 0 x=0 x=0处连续。由 y ( 0 ) = 0 y(0)=0 y(0)=0,有
lim ⁡ x → 0 + y = 0 ⇒ C 2 = 0 ⇒ y = 2 x + 1 2 C 1 x 2 \lim\limits_{x\to0^{+}}y=0 \Rightarrow C_{2}=0\Rightarrow y=2x+ \frac{1}{2}C_{1}x^{2} x0+limy=0C2=0y=2x+21C1x2

由于上面出现了 1 x \begin{aligned} \frac{1}{x}\end{aligned} x1,因此此处不能直接把 x = 0 x=0 x=0代入方程
由于一阶导数在 x = 0 x=0 x=0处存在,因此 y y y x = 0 x=0 x=0处连续,可以用逼近的方式求

由与直线 x = 1 x=1 x=1 y = 0 y=0 y=0围成的平面区域 D D D的面积为 2 2 2
2 = ∫ 0 1 ( 2 x + 1 2 C 1 x 2 ) d x = 1 + 1 6 C 1 ⇒ C 1 = 6 2=\int_{0}^{1}\left(2x+ \frac{1}{2}C_{1}x^{2}\right)dx=1+ \frac{1}{6}C_{1}\Rightarrow C_{1}=6 2=01(2x+21C1x2)dx=1+61C1C1=6

y = 2 x + 3 x 2 y=2x+3x^{2} y=2x+3x2
旋转体体积为
V = 2 π ∫ 0 1 x y ( x ) d x = 17 π 6 V=2 \pi \int_{0}^{1}xy(x)dx= \frac{17\pi}{6} V=2π01xy(x)dx=617π

一般应用考的较多的就是定积分的应用,和微分方程的应用

你可能感兴趣的:(高等数学,学习,高等数学)