例7:微分方程 x y ′ + y ( ln x − ln y ) = 0 xy'+y(\ln x-\ln y)=0 xy′+y(lnx−lny)=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} u′ln∣lnu−1∣lnu−1lnxy−1=xu(lnu−1)=lnx+C=Cx=Cx
由 y ( 1 ) = e 3 ⇒ C = 2 y(1)=e^{3}\Rightarrow C=2 y(1)=e3⇒C=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} lnxy−1y=2x=xe2x+1
例8:微分方程 y d x + ( x − 3 y 2 ) d y = 0 ydx+(x-3y^{2})dy=0 ydx+(x−3y2)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+(x−3y2)dydydx+yxx代入x=1,y=1⇒C⇒x⇒yy=0=3y=e−∫y1dy[∫3ye∫y1dydy+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=0⇒r1,2=−1±2i=e−1(C1cos2x+C2sin2x)
例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} 特征方程r2−4r+8r1,2对于y′′−4y′+8yy∗对于y′′−4y′+8yy∗y∗=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+(x−31)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+(x−31)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+(x−31)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+xex−31ex
如果 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} y∗xex是非齐次方程的一个解,因此齐次方程的特征方程为
( r − 1 ) ( r − 2 ) = 0 ⇒ a = − 3 , b = 2 (r-1)(r-2)=0\Rightarrow a=-3,b=2 (r−1)(r−2)=0⇒a=−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′′−y′−esinx=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 x→0时,函数 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} x→0limy(x)ln(1+x2)=x→0limy(x)x2=x→0limy′(x)2x=x→0limy′′(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(−2e−x+C)=Ce3x−2e2x(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)=3e3x−2e2x
例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(x−t)dt=∫0x(x−t)f(t)dt+e−x−1,求 f ( x ) f(x) f(x)的表达式
令 u = x − t u=x-t u=x−t,则 ∫ 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(x−t)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)=x∫0xf(t)dt−∫0xtf(t)dt+e−x−1两端求导=∫0xf(t)dt−e−x令x=0,f(0)=−1=e−x=Cex−2e−x代入f(0)=−1,得C=−21=−2ex+e−x
例15:设函数 f ( x ) f(x) f(x)在定义域 I I I上的导数大于零,若对于任意的 x 0 ∈ I x_{0}\in I x0∈I,曲线 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} y−f(x0)x=f′(x0)(x−x0)令y=0=x0−f′(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} S21y2−y8Cy=21f′(x0)∣f(x0)∣⋅∣f(x0)∣=4记y=f(x0),则=4y′=x+C代入y(0)=2=−4带回上式=4−x8
例16:设非负函数 y = y ( x ) ( x ≥ 0 ) y=y(x)(x \geq 0) y=y(x)(x≥0)满足微分方程 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} p′−x1py′=py=−x2(x>0)=e∫x1dx(∫−x1e∫−x1dxdx+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(x≥0),可知 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} x→0+limy=0⇒C2=0⇒y=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+61C1⇒C1=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π
一般应用考的较多的就是定积分的应用,和微分方程的应用