Skip to content

第二章 一阶微分方程的初等积分法 ​

第二章 一阶微分方程的初等积分法 ​

来源:课堂讲义《第二章 一阶微分方程的初等解法 2.1/2.2/2.3/2.4》 笔记由课堂 PDF 整理,公式按数学含义重建;模糊处标 [?] 待校对。

§2.1 变量分离方程与变量变换 ​

一、变量分离方程 ​

定义 1:形如

dydx=f(x)φ(y)(2.1)\frac{dy}{dx} = f(x)\varphi(y) \tag{2.1}

的方程称为变量分离方程,其中 f(x),φ(y)f(x), \varphi(y) 分别是 x,yx, y 的连续函数。

求解方法 ​

当 φ(y)≠0\varphi(y) \neq 0 时,将 (2.1) 写成

dyφ(y)=f(x) dx(2.2)\frac{dy}{\varphi(y)} = f(x)\,dx \tag{2.2}

两边积分得

∫dyφ(y)=∫f(x) dx+c\int \frac{dy}{\varphi(y)} = \int f(x)\,dx + c

由 (2.2) 所确定的函数 Φ(x,y)=c\Phi(x, y) = c 即为 (2.1) 的解。

注:若存在 y0y_0 使 φ(y0)=0\varphi(y_0) = 0,则 y=y0y = y_0 也是方程的解,但它不包含在通解中,必须予以补上。

例 1:求微分方程 dydx=y(10−y)\dfrac{dy}{dx} = y(10 - y) 的所有解。

解:方程两边同除以 y(10−y)y(10-y),分离变量得

dyy(10−y)=dx\frac{dy}{y(10-y)} = dx

积分得

ln⁡∣y10−y∣=x+c1\ln\left\vert \frac{y}{10-y}\right\vert = x + c_1

从上式中解出 yy(将常数记为 c≠0c \neq 0):

y=101+ce−xy = \frac{10}{1 + ce^{-x}}

此外由 φ(y)=y(10−y)=0\varphi(y) = y(10-y) = 0 得 y=0y = 0 和 y=10y = 10 也是方程的解。

故方程的所有解为:y=101+ce−xy = \dfrac{10}{1+ce^{-x}}(cc 为任常数)以及 y=0y = 0、y=10y = 10。

例 2:求微分方程 dydx=23 yx\dfrac{dy}{dx} = \dfrac{2\sqrt{3}\,y}{x} 的通解 [?]。

解:分离变量后得 dyy=23x dx\dfrac{dy}{y} = \dfrac{2\sqrt{3}}{x}\,dx,两边积分得

23ln⁡∣y∣=ln⁡∣x∣+c12\sqrt{3}\ln\vert y\vert = \ln\vert x\vert + c_1

整理后得通解 y2=cxy^2 = cx(cc 为任常数)[?]。

校对点

例 2 原文公式较乱("232\sqrt{3}" 与 "ln⁡2\ln 2" 混排),此处按常规解读为 dydx=23 yx\dfrac{dy}{dx} = \dfrac{2\sqrt{3}\,y}{x};另有解 y=0y = 0 未包含在通解中,应补上。

例 3:求微分方程 dydx=p(x)y\dfrac{dy}{dx} = p(x)y 的通解,其中 p(x)p(x) 是 xx 的连续函数。

解:分离变量 dyy=p(x) dx\dfrac{dy}{y} = p(x)\,dx,两边积分 ln⁡∣y∣=∫p(x) dx+c1\ln\vert y\vert = \int p(x)\,dx + c_1,由对数定义

y=ec1e∫p(x) dx=ce∫p(x) dxy = e^{c_1}e^{\int p(x)\,dx} = ce^{\int p(x)\,dx}

允许 c=0c = 0 时也包括 y=0y = 0,故方程的通解为

y=ce∫p(x) dx(c 为任常数)y = ce^{\int p(x)\,dx} \quad (\text{$c$ 为任常数})

例 4:求初值问题 {dydx=y2cos⁡xy(0)=1\begin{cases} \dfrac{dy}{dx} = y^2\cos x \\ y(0) = 1 \end{cases} 的特解。

解:先求通解。当 y≠0y \neq 0 时分离变量 dyy2=cos⁡x dx\dfrac{dy}{y^2} = \cos x\,dx,两边积分得 −1y=sin⁡x+c-\dfrac{1}{y} = \sin x + c,通解为

y=−1sin⁡x+cy = -\frac{1}{\sin x + c}

此外 y=0y = 0 也是方程的解。代入初始条件 y(0)=1y(0) = 1 得 c=−1c = -1,故所求特解为

y=11−sin⁡xy = \frac{1}{1 - \sin x}

二、可化为变量分离方程的类型 ​

(I) 齐次方程 ​

形如

dydx=g(yx)(2.5)\frac{dy}{dx} = g\left(\frac{y}{x}\right) \tag{2.5}

的方程称为齐次方程,其中 g(u)g(u) 是 uu 的连续函数。

求解方法:作变量代换 u=yxu = \dfrac{y}{x}(即 y=xuy = xu),由 y′=u+xdudxy' = u + x\dfrac{du}{dx} 得

xdudx=g(u)−ux\frac{du}{dx} = g(u) - u

这是关于 uu 的变量分离方程,解出后再变量还原。

例 4:求解方程 xdydx=y+xyx\dfrac{dy}{dx} = y + \sqrt{xy}(x<0x < 0)。

解:方程变形为 dydx=yx+yx\dfrac{dy}{dx} = \dfrac{y}{x} + \sqrt{\dfrac{y}{x}},这是齐次方程。令 u=yxu = \dfrac{y}{x},代入得 xdudx=ux\dfrac{du}{dx} = \sqrt{u},分离变量 duu=dxx\dfrac{du}{\sqrt{u}} = \dfrac{dx}{x},积分得

2u=ln⁡∣x∣+c2\sqrt{u} = \ln\vert x\vert + c

变量还原得原方程的通解 2yx=ln⁡∣x∣+c2\sqrt{\dfrac{y}{x}} = \ln\vert x\vert + c(x<0x<0)[?]。

例 6:求初值问题 {(y2+x2) dx+xy dy=0y(1)=0\begin{cases} (y^2 + x^2)\,dx + xy\,dy = 0 \\ y(1) = 0 \end{cases} 的解。

解:方程变形为 dydx=−x2+y2xy=−1+(y/x)2y/x\dfrac{dy}{dx} = -\dfrac{x^2+y^2}{xy} = -\dfrac{1+(y/x)^2}{y/x},这是齐次方程。令 u=yxu = \dfrac{y}{x},代入得 xdudx=−1+u2u−u=−1+2u2ux\dfrac{du}{dx} = -\dfrac{1+u^2}{u} - u = -\dfrac{1+2u^2}{u}[?],分离变量并积分:

∫u1+u2 du=−∫dxx\int \frac{u}{1+u^2}\,du = -\int \frac{dx}{x}

得 ln⁡1+u2=−ln⁡∣x∣+c\ln\sqrt{1+u^2} = -\ln\vert x\vert + c,整理得 x2+y2=c\sqrt{x^2+y^2} = c[?]。由 y(1)=0y(1)=0 定出 c=1c=1,故初值问题的解为 x2+y2=1x^2+y^2=1(x>0x>0)[?]。

校对点

例 6 讲义原文中间步骤(ln⁡∣u∣+ln⁡1+u2\ln\vert u\vert +\ln\sqrt{1+u^2} 等)提取较乱,最终结论为 x2+y2=cx\sqrt{x^2+y^2}=cx,由 y(1)=0y(1)=0 得 c=1c=1,即 x2+y2=1x^2+y^2=1,请对照原讲义核实。

(II) 形如 dydx=a1x+b1y+c1a2x+b2y+c2\dfrac{dy}{dx} = \dfrac{a_1x+b_1y+c_1}{a_2x+b_2y+c_2} 的方程 ​

其中 ai,bi,cia_i, b_i, c_i 均为常数,可经过变量变换化为变量分离方程。分三种情况:

情形 1:c1=c2=0c_1 = c_2 = 0,方程为齐次方程

dydx=g(yx)\frac{dy}{dx} = g\left(\frac{y}{x}\right)

由 (I) 可化为变量分离方程。

情形 2:a1a2=b1b2=k\dfrac{a_1}{a_2} = \dfrac{b_1}{b_2} = k(即 ∣a1b1a2b2∣=0\begin{vmatrix} a_1 & b_1 \\ a_2 & b_2 \end{vmatrix} = 0),令 u=a2x+b2yu = a_2 x + b_2 y,则方程化为变量分离方程。

情形 3:c1,c2c_1, c_2 不同时为零,且 ∣a1b1a2b2∣≠0\begin{vmatrix} a_1 & b_1 \\ a_2 & b_2 \end{vmatrix} \neq 0。解方程组

{a1x+b1y+c1=0a2x+b2y+c2=0\begin{cases} a_1x + b_1y + c_1 = 0 \\ a_2x + b_2y + c_2 = 0 \end{cases}

得交点 (α,β)≠(0,0)(\alpha, \beta) \neq (0,0),作坐标变换

X=x−α,Y=y−βX = x - \alpha, \quad Y = y - \beta

则方程化为 dYdX=g(YX)\dfrac{dY}{dX} = g\left(\dfrac{Y}{X}\right) 的齐次方程,进而求解。

求解步骤:

  1. 解方程组得 x=α,y=βx = \alpha, y = \beta
  2. 作变换 X=x−α,Y=y−βX = x - \alpha, Y = y - \beta,方程化为 dYdX=g(YX)\dfrac{dY}{dX} = g\left(\dfrac{Y}{X}\right)
  3. 再经变换 u=YXu = \dfrac{Y}{X} 化为变量分离方程
  4. 求解
  5. 变量还原

例 7:求微分方程 dydx=x−y+1x+y−3\dfrac{dy}{dx} = \dfrac{x - y + 1}{x + y - 3} 的通解。

解:解方程组 {x−y+1=0x+y−3=0\begin{cases} x - y + 1 = 0 \\ x + y - 3 = 0 \end{cases} 得 x=1,y=2x = 1, y = 2。令 X=x−1,Y=y−2X = x - 1, Y = y - 2,代入得

dYdX=X−YX+Y\frac{dY}{dX} = \frac{X - Y}{X + Y}

令 u=YXu = \dfrac{Y}{X},分离变量 1+u21−2u−u2 du=dXX\dfrac{1+u^2}{1-2u-u^2}\,du = \dfrac{dX}{X}[?],积分得

arctan⁡u−12ln⁡(1−2u−u2)=ln⁡∣X∣+c\arctan u - \frac{1}{2}\ln(1-2u-u^2) = \ln\vert X\vert + c

变量还原并整理得原方程的通解:

arctan⁡y−2x−1−12ln⁡(1−2y−2x−1−(y−2)2(x−1)2)=ln⁡∣x−1∣+c\arctan\frac{y-2}{x-1} - \frac{1}{2}\ln\left(1 - 2\frac{y-2}{x-1} - \frac{(y-2)^2}{(x-1)^2}\right) = \ln\vert x-1\vert + c

校对点

例 7 讲义原文为 "12arctan⁡2u1−u2+ln⁡(1+u2)=ln⁡X+c\frac{1}{2}\arctan\frac{2u}{1-u^2} + \ln(1+u^2) = \ln X + c" 形式,与上述写法等价但细节请对照原讲义。

更一般的类型(可经适当变量变换化为变量分离方程):

  • dydx=f(ax+by+c)\dfrac{dy}{dx} = f(ax + by + c) —— 令 u=ax+by+cu = ax+by+c
  • xdydx=f(xy)x\dfrac{dy}{dx} = f(xy) —— 令 u=xyu = xy
  • xdydx=f(yx2)x\dfrac{dy}{dx} = f\left(\dfrac{y}{x^2}\right) —— 令 u=yx2u = \dfrac{y}{x^2}
  • M(x,y)(x dx+y dy)+N(x,y)(x dy−y dx)=0M(x,y)(x\,dx + y\,dy) + N(x,y)(x\,dy - y\,dx) = 0,其中 M,NM, N 为齐次函数(次数可以不相同)

例 8:求微分方程 (xy+y2) dx+(x2−xy) dy=0(xy + y^2)\,dx + (x^2 - xy)\,dy = 0 的通解 [?]。

解:令 u=xyu = xy,则 du=x dy+y dxdu = x\,dy + y\,dx,代入方程并整理得

(u2−u) dx+x2 du=0(u^2 - u)\,dx + x^2\,du = 0

分离变量 dxx=duu−u2\dfrac{dx}{x} = \dfrac{du}{u - u^2},积分得 ln⁡∣x∣=ln⁡∣uu−1∣+c\ln\vert x\vert = \ln\left\vert \dfrac{u}{u-1}\right\vert + c,变量还原得通解

xy−y=cxxy - y = cx

校对点

例 8 讲义最终结论为 ln⁡∣xy∣−xy=c\ln\vert xy\vert - xy = c(原文 "ln⁡c=xy−ln⁡xy\ln c = xy - \ln xy"),与这里得到的 xy−y=cxxy - y = cx 不一致,请对照原讲义核实系数。

作业(2.1):p25:1.(2、6、7);2.(1、4、6);4;5;7


§2.2 线性方程与常数变易法 ​

一、一阶线性微分方程 ​

一阶线性微分方程在 a(x)≠0a(x) \neq 0 的区间上可写成

dydx+P(x)y=Q(x)(1)\frac{dy}{dx} + P(x)y = Q(x) \tag{1}

其中 P(x),Q(x)P(x), Q(x) 在考虑的区间上是连续函数。

  • 若 Q(x)≡0Q(x) \equiv 0,得 y′+P(x)y=0y' + P(x)y = 0,称为一阶齐次线性方程
  • 若 Q(x)≢0Q(x) \not\equiv 0,称为一阶非齐次线性方程

常数变易法 ​

步骤:

  1. 解对应的齐次方程 y′+P(x)y=0y' + P(x)y = 0,分离变量得

y=ce−∫P(x) dxy = ce^{-\int P(x)\,dx}

  1. 常数变易:将常数 cc 变为待定函数 c(x)c(x),令

y=c(x)e−∫P(x) dxy = c(x)e^{-\int P(x)\,dx}

  1. 代入原方程 (1),得 c′(x)=Q(x)e∫P(x) dxc'(x) = Q(x)e^{\int P(x)\,dx},积分得

c(x)=∫Q(x)e∫P(x) dx dx+c~c(x) = \int Q(x)e^{\int P(x)\,dx}\,dx + \tilde{c}

故 (1) 的通解为

y=e−∫P(x) dx(∫Q(x)e∫P(x) dx dx+c~)(3)y = e^{-\int P(x)\,dx}\left(\int Q(x)e^{\int P(x)\,dx}\,dx + \tilde{c}\right) \tag{3}

注:求 (1) 的通解可直接用公式 (3)。

例 1:求方程 (1+x)dydx−ny=ex(1+x)n+1(1+x)\dfrac{dy}{dx} - ny = e^x(1+x)^{n+1} 的通解(nn 为常数)。

解:改写为 dydx−n1+xy=ex(1+x)n\dfrac{dy}{dx} - \dfrac{n}{1+x}y = e^x(1+x)^n。齐次方程通解 y=c(1+x)ny = c(1+x)^n。常数变易,令 y=c(x)(1+x)ny = c(x)(1+x)^n,代入得 c′(x)=exc'(x) = e^x,积分得 c(x)=ex+c~c(x) = e^x + \tilde{c},故通解为

y=(1+x)n(ex+c~)y = (1+x)^n(e^x + \tilde{c})

例 2:求方程 dydx=x2−yx\dfrac{dy}{dx} = x^2 - \dfrac{y}{x} 的通解 [?]。

解:改写为 dydx+1xy=x2\dfrac{dy}{dx} + \dfrac{1}{x}y = x^2,用公式 (3)(P=1x,Q=x2P = \dfrac{1}{x}, Q = x^2):

y=e−ln⁡x(∫x2eln⁡x dx+c)=1x(x44+c)y = e^{-\ln x}\left(\int x^2 e^{\ln x}\,dx + c\right) = \frac{1}{x}\left(\frac{x^4}{4} + c\right)

例 3:求初值问题 {dydx=4yx2+x2+1y(1)=1\begin{cases} \dfrac{dy}{dx} = \dfrac{4y}{x^2} + x^2 + 1 \\ y(1) = 1 \end{cases} 的解 [?]。

解:原方程即 y′−4x2y=x2+1y' - \dfrac{4}{x^2}y = x^2+1,用公式 (3) 得通解

y=e4/x(∫(x2+1)e−4/x dx+c~)y = e^{4/x}\left(\int (x^2+1)e^{-4/x}\,dx + \tilde{c}\right)

代入 y(1)=1y(1) = 1 定出常数,得初值问题的解[?]。

校对点

例 2、例 3 原文公式与数字提取混乱("1+4x2+3x31+4x^2+3x^3" 等),此处按通解公式常规推导;讲义例 3 最终答案为 y=x2+3x+2x2ln⁡xy = x^2 + \frac{3}{x} + 2x^2\ln x[?],请对照原讲义核实。

练习:

  1. (1+x2)y′+xy=e−x(1+x^2)y' + xy = e^{-x} → 通解 y=e−x+c1+x2y = \dfrac{e^{-x} + c}{\sqrt{1+x^2}}[?]
  2. y′+2xy=4xy' + 2xy = 4x → 通解 y=2+ce−x2y = 2 + ce^{-x^2}
  3. dydx+y=xy3\dfrac{dy}{dx} + y = xy^3 → 改写为 dxdy−xy=y2\dfrac{dx}{dy} - \dfrac{x}{y} = y^2[?],通解 x=cy−y2x = cy - y^2[?]

二、伯努利(Bernoulli)方程 ​

形如

dydx+P(x)y=Q(x)yn(n≠0,1)\frac{dy}{dx} + P(x)y = Q(x)y^n \quad (n \neq 0, 1)

的方程称为伯努利方程,其中 P(x),Q(x)P(x), Q(x) 为 xx 的连续函数。

解法:引入变量变换 z=y1−nz = y^{1-n},方程化为

dzdx+(1−n)P(x)z=(1−n)Q(x)\frac{dz}{dx} + (1-n)P(x)z = (1-n)Q(x)

求上述线性方程的通解,再变量还原。

例 4:求方程 dydx+1xy=2y2\dfrac{dy}{dx} + \dfrac{1}{x}y = 2y^2 的通解。

解:这是 n=2n=2 的伯努利方程。令 z=y−1z = y^{-1},方程化为 dzdx−1xz=−2\dfrac{dz}{dx} - \dfrac{1}{x}z = -2,解得 z=cx+2z = cx + 2,故通解为

y=1cx+2y = \frac{1}{cx + 2}

例 5:求方程 2xy′+y=y2ln⁡x2xy' + y = y^2\ln x 的通解[?]。

解:改写为 y′+12xy=ln⁡x2xy2y' + \dfrac{1}{2x}y = \dfrac{\ln x}{2x}y^2,这是 n=2n=2 的伯努利方程。令 z=y−1z = y^{-1},得 z′−12xz=−ln⁡x2xz' - \dfrac{1}{2x}z = -\dfrac{\ln x}{2x},解得

z=x(c−12(ln⁡x)2)z = \sqrt{x}\left(c - \frac{1}{2}(\ln x)^2\right)

故通解 y=1x(c−12ln⁡2x)y = \dfrac{1}{\sqrt{x}\left(c - \frac{1}{2}\ln^2 x\right)}。

三、线性微分方程的应用举例:R-L 串联电路 ​

由电感 LL、电阻 RR 和电源组成的串联电路,电动势 EE 为常数。开关 KK 合上后(t=0t=0 时刻,初始电流 I(0)=0I(0) = 0),求电流强度 I(t)I(t)。

由 Kirchhoff 第二定律(闭合回路电压代数和为零):电感压降 LdIdtL\dfrac{dI}{dt},电阻压降 RIRI,故

LdIdt+RI=E⇒dIdt+RLI=ELL\frac{dI}{dt} + RI = E \quad \Rightarrow \quad \frac{dI}{dt} + \frac{R}{L}I = \frac{E}{L}

解线性方程得通解 I(t)=ER+ce−Rt/LI(t) = \dfrac{E}{R} + ce^{-Rt/L},由 I(0)=0I(0) = 0 得 c=−ERc = -\dfrac{E}{R},故

I(t)=ER(1−e−Rt/L)I(t) = \frac{E}{R}\left(1 - e^{-Rt/L}\right)

作业(2.2):p30 1.(2、4、8、12);2;p31 4;5.(2);7.(2)


§2.3 恰当方程与积分因子 ​

一、恰当方程的定义及条件 ​

设 u(x,y)u(x,y) 是连续可微的函数,其全微分为 du=∂u∂xdx+∂u∂ydydu = \dfrac{\partial u}{\partial x}dx + \dfrac{\partial u}{\partial y}dy。如果恰好有方程

∂u(x,y)∂xdx+∂u(x,y)∂ydy=0\frac{\partial u(x,y)}{\partial x}dx + \frac{\partial u(x,y)}{\partial y}dy = 0

就可以马上写出它的隐式解 u(x,y)=cu(x,y) = c。

定义 1:若有函数 u(x,y)u(x,y) 使得 du(x,y)=M(x,y)dx+N(x,y)dydu(x,y) = M(x,y)dx + N(x,y)dy,则称微分方程

M(x,y)dx+N(x,y)dy=0(1)M(x,y)dx + N(x,y)dy = 0 \tag{1}

是恰当方程。此时 (1) 的通解为 u(x,y)=cu(x,y) = c。

例:xdy+ydx=0xdy + ydx = 0 是恰当方程(=d(xy)= d(xy));(3x2+2xy)dy+(y3+x2)dx=0(3x^2 + 2xy)dy + (y^3 + x^2)dx = 0 是恰当方程(=d(x2y+y3)= d(x^2y + y^3)[?])。

需考虑的问题:

  1. 方程 (1) 是否为恰当方程?
  2. 若是恰当方程,怎样求解?
  3. 若不是恰当方程,有无可能转化为恰当方程求解?

方程为恰当方程的充要条件 ​

定理 1:设函数 M(x,y),N(x,y)M(x,y), N(x,y) 在矩形区域 RR 内连续且有连续的一阶偏导数,则方程 M(x,y)dx+N(x,y)dy=0M(x,y)dx + N(x,y)dy = 0 为恰当方程的充要条件是

∂M(x,y)∂y=∂N(x,y)∂x(2)\frac{\partial M(x,y)}{\partial y} = \frac{\partial N(x,y)}{\partial x} \tag{2}

证明:

  • 必要性:设 (1) 是恰当方程,存在 u(x,y)u(x,y) 使 du=Mdx+Ndydu = Mdx + Ndy,则 ∂u∂x=M,∂u∂y=N\dfrac{\partial u}{\partial x} = M, \dfrac{\partial u}{\partial y} = N,从而 ∂2u∂y∂x=∂M∂y\dfrac{\partial^2 u}{\partial y\partial x} = \dfrac{\partial M}{\partial y},∂2u∂x∂y=∂N∂x\dfrac{\partial^2 u}{\partial x\partial y} = \dfrac{\partial N}{\partial x},由混合偏导数连续相等得 ∂M∂y=∂N∂x\dfrac{\partial M}{\partial y} = \dfrac{\partial N}{\partial x}。
  • 充分性:若 ∂M∂y=∂N∂x\dfrac{\partial M}{\partial y} = \dfrac{\partial N}{\partial x},构造 u(x,y)=∫M(x,y) dx+φ(y)u(x,y) = \int M(x,y)\,dx + \varphi(y),其中 φ(y)\varphi(y) 由 ∂u∂y=N\dfrac{\partial u}{\partial y} = N 确定:φ′(y)=N−∂∂y∫M dx\varphi'(y) = N - \dfrac{\partial}{\partial y}\int M\,dx,可证右端与 xx 无关,积分即得。

通解公式:若 (1) 为恰当方程,则其通解为

∫M(x,y) dx+∫(N−∂∂y∫M dx)dy=c\int M(x,y)\,dx + \int\left(N - \frac{\partial}{\partial y}\int M\,dx\right)dy = c

二、恰当方程的求解 ​

关键是求原函数 u(x,y)u(x,y) 的问题。三种方法:

1. 不定积分法 ​

步骤:① 判断是否恰当方程;② 设 u(x,y)=∫M dx+φ(y)u(x,y) = \int M\,dx + \varphi(y);③ 由 ∂u∂y=N\dfrac{\partial u}{\partial y} = N 求 φ(y)\varphi(y)。

例 1:验证方程 (ex+y) dx+(x−2sin⁡y) dy=0(e^x + y)\,dx + (x - 2\sin y)\,dy = 0 是恰当方程并求通解。

解:M=ex+y,N=x−2sin⁡yM = e^x + y, N = x - 2\sin y,∂M∂y=1=∂N∂x\dfrac{\partial M}{\partial y} = 1 = \dfrac{\partial N}{\partial x},故是恰当方程。

设 u=∫(ex+y)dx=ex+yx+φ(y)u = \int(e^x + y)dx = e^x + yx + \varphi(y),由 ∂u∂y=x+φ′(y)=x−2sin⁡y\dfrac{\partial u}{\partial y} = x + \varphi'(y) = x - 2\sin y 得 φ′(y)=−2sin⁡y\varphi'(y) = -2\sin y,积分 φ(y)=2cos⁡y\varphi(y) = 2\cos y。故通解

ex+xy+2cos⁡y=ce^x + xy + 2\cos y = c

2. 线积分法 ​

由曲线积分与路径无关定理,取 (x0,y0)∈R(x_0, y_0) \in R:

u(x,y)=∫x0xM(x,y0) dx+∫y0yN(x,y) dyu(x,y) = \int_{x_0}^{x} M(x, y_0)\,dx + \int_{y_0}^{y} N(x, y)\,dy

通解为 ∫x0xM(x,y0) dx+∫y0yN(x,y) dy=c\int_{x_0}^{x} M(x,y_0)\,dx + \int_{y_0}^{y} N(x,y)\,dy = c。

例 2:求解方程 (2xcos⁡y+ey) dx+(ey−x2sin⁡y) dy=0(2x\cos y + e^y)\,dx + (e^y - x^2\sin y)\,dy = 0。

解:M=2xcos⁡y+ey,N=ey−x2sin⁡yM = 2x\cos y + e^y, N = e^y - x^2\sin y,∂M∂y=−2xsin⁡y+ey=∂N∂x\dfrac{\partial M}{\partial y} = -2x\sin y + e^y = \dfrac{\partial N}{\partial x},是恰当方程。取 (x0,y0)=(0,0)(x_0,y_0) = (0,0):

u=∫0x2x dx+∫0y(ey−x2sin⁡y) dy=x2+ey−1+x2cos⁡y−x2u = \int_0^x 2x\,dx + \int_0^y (e^y - x^2\sin y)\,dy = x^2 + e^y - 1 + x^2\cos y - x^2

整理得通解

x2cos⁡y+ey=cx^2\cos y + e^y = c

3. 分组凑微法 ​

采用"分项组合"的方法,把本身已构成全微分的项分出来,再把剩余项凑成全微分。应熟记简单二元函数的全微分:

xdy+ydx=d(xy),xdy−ydxy2=d(xy),ydx−xdyx2=d(yx)xdy + ydx = d(xy), \quad \frac{xdy - ydx}{y^2} = d\left(\frac{x}{y}\right), \quad \frac{ydx - xdy}{x^2} = d\left(\frac{y}{x}\right)

xdy−ydxx2+y2=d(arctan⁡yx),xdy−ydxxy=d(ln⁡∣yx∣),xdy+ydxx−y=d(ln⁡x+yx−y)[?]\frac{xdy - ydx}{x^2 + y^2} = d\left(\arctan\frac{y}{x}\right), \quad \frac{xdy - ydx}{xy} = d\left(\ln\left\vert \frac{y}{x}\right\vert \right), \quad \frac{xdy + ydx}{x - y} = d\left(\ln\frac{x+y}{x-y}\right)[?]

例 3:求方程 (3x2+6xy2) dx+(6x2y+4y3) dy=0(3x^2 + 6xy^2)\,dx + (6x^2y + 4y^3)\,dy = 0 的通解。

解:∂M∂y=12xy=∂N∂x\dfrac{\partial M}{\partial y} = 12xy = \dfrac{\partial N}{\partial x},是恰当方程。分项组合:

(3x2 dx+4y3 dy)+(6xy2 dx+6x2y dy)=d(x3)+d(y4)+d(3x2y2)=0(3x^2\,dx + 4y^3\,dy) + (6xy^2\,dx + 6x^2y\,dy) = d(x^3) + d(y^4) + d(3x^2y^2) = 0

故通解为 x3+y4+3x2y2=cx^3 + y^4 + 3x^2y^2 = c。

例 4:验证方程 (cos⁡x−xsin⁡x+xy2) dx+(x2y−y) dy=0(\cos x - x\sin x + xy^2)\,dx + (x^2y - y)\,dy = 0[?] 是恰当方程,并求满足 y(0)=2y(0) = 2 的解。

解:分项组合得 d(xcos⁡x)+d(x2y22)−d(y22)=0d(x\cos x) + d\left(\dfrac{x^2y^2}{2}\right) - d\left(\dfrac{y^2}{2}\right) = 0[?],通解为

xcos⁡x+x2y22−y22=cx\cos x + \frac{x^2y^2}{2} - \frac{y^2}{2} = c

由 y(0)=2y(0) = 2 得 c=−2c = -2,故 xcos⁡x+x2y2−y22=−2x\cos x + \dfrac{x^2y^2 - y^2}{2} = -2,即 x2y2−y2+2xcos⁡x+4=0x^2y^2 - y^2 + 2x\cos x + 4 = 0。

校对点

例 4 讲义原文为 "(cos⁡x−xsin⁡x) dx(\cos x - x\sin x)\,dx" 与 "xy2dx+(x2y−y)dyxy^2 dx + (x^2y - y)dy" 组合成 d(xcos⁡x)+12d(x2y2−y2)d(x\cos x) + \dfrac{1}{2}d(x^2y^2 - y^2),最终 xcos⁡x+x2y2−y22=cx\cos x + \dfrac{x^2y^2 - y^2}{2} = c,c=−2c = -2。若 MM 的 xy2xy^2 系数不同请以原讲义为准。

三、积分因子 ​

对非恰当方程,乘上一个因子后可变为恰当方程。

例:一阶线性方程 dy−(P(x)y+Q(x))dx=0dy - (P(x)y + Q(x))dx = 0 不是恰当方程,乘以 e−∫P(x)dxe^{-\int P(x)dx} 后变为恰当方程。

1. 定义 ​

如果存在连续可微函数 μ(x,y)≠0\mu(x,y) \neq 0,使得

μ(x,y)M(x,y) dx+μ(x,y)N(x,y) dy=0\mu(x,y)M(x,y)\,dx + \mu(x,y)N(x,y)\,dy = 0

为恰当方程,则称 μ(x,y)\mu(x,y) 是方程 (1) 的一个积分因子。

例 5:验证 μ=xy\mu = xy 是方程 (3xy+4y2) dx+(2x2+3xy) dy=0(3xy + 4y^2)\,dx + (2x^2 + 3xy)\,dy = 0 的积分因子并求通解。

解:乘以 μ=xy\mu = xy 后 M′=3x2y2+4xy3,N′=2x3y+3x2y2M' = 3x^2y^2 + 4xy^3, N' = 2x^3y + 3x^2y^2,∂M′∂y=6x2y+12xy2=∂N′∂x\dfrac{\partial M'}{\partial y} = 6x^2y + 12xy^2 = \dfrac{\partial N'}{\partial x},是恰当方程。分项组合得 d(x3y2)+d(x2y3)=0d(x^3y^2) + d(x^2y^3) = 0,通解

x3y2+x2y3=cx^3y^2 + x^2y^3 = c

2. 积分因子的确定 ​

μ(x,y)\mu(x,y) 是方程 (1) 的积分因子的充要条件是

μ(∂M∂y−∂N∂x)=N∂μ∂x−M∂μ∂y\mu\left(\frac{\partial M}{\partial y} - \frac{\partial N}{\partial x}\right) = N\frac{\partial \mu}{\partial x} - M\frac{\partial \mu}{\partial y}

这本身是偏微分方程,一般比原方程更困难,但提供了寻找特殊形式积分因子的途径。

仅依赖于 xx 的积分因子:方程有仅依赖于 xx 的积分因子 μ=μ(x)\mu = \mu(x) 的充要条件是

ψ(x)=∂M∂y−∂N∂xN\psi(x) = \frac{\dfrac{\partial M}{\partial y} - \dfrac{\partial N}{\partial x}}{N}

只是 xx 的函数(与 yy 无关),此时积分因子为

μ(x)=e∫ψ(x) dx(10)\mu(x) = e^{\int \psi(x)\,dx} \tag{10}

仅依赖于 yy 的积分因子:同理,充要条件是

φ(y)=∂N∂x−∂M∂yM\varphi(y) = \frac{\dfrac{\partial N}{\partial x} - \dfrac{\partial M}{\partial y}}{M}

只是 yy 的函数,此时积分因子为 μ(y)=e∫φ(y) dy\mu(y) = e^{\int \varphi(y)\,dy}。

例 6:求方程 (2y+2y2) dx+(y+ex) dy=0(2y + 2y^2)\,dx + (y + e^x)\,dy = 0[?] 的通解。

解:M=2y+2y2,N=y+exM = 2y + 2y^2, N = y + e^x,∂M∂y=2+4y\dfrac{\partial M}{\partial y} = 2 + 4y,∂N∂x=ex\dfrac{\partial N}{\partial x} = e^x,不是恰当方程。计算

∂M∂y−∂N∂xN=2+4y−exy+ex\frac{\dfrac{\partial M}{\partial y} - \dfrac{\partial N}{\partial x}}{N} = \frac{2 + 4y - e^x}{y + e^x}

不满足仅依赖 xx[?];换用仅依赖 yy 的判别:

∂N∂x−∂M∂yM=ex−2−4y2y+2y2\frac{\dfrac{\partial N}{\partial x} - \dfrac{\partial M}{\partial y}}{M} = \frac{e^x - 2 - 4y}{2y + 2y^2}

亦非仅依赖 yy[?],此题积分因子需另寻(讲义结论 μ=ex\mu = e^x,通解 2xyex+2y2ex=c2xye^x + 2y^2e^x = c[?])。

校对点

例 6 讲义原文:M=2y+2y2M = 2y+2y^2,N=y+exN = y + e^x,∂M∂y−∂N∂x=(2+4y)−ex\dfrac{\partial M}{\partial y} - \dfrac{\partial N}{\partial x} = (2+4y) - e^x,讲义写为 2+exN\dfrac{2+e^x}{N} 形式有笔误;最终取 μ=ex\mu = e^x,通解 2xyex+2y2ex=c2xye^x + 2y^2e^x = c。请对照原讲义。

例 7:求解方程 dydx=−x+x2+y2y\dfrac{dy}{dx} = -\dfrac{x + \sqrt{x^2 + y^2}}{y}(y>0y > 0)[?]。

解:改写为 x dx+y dy+x2+y2 dx=0x\,dx + y\,dy + \sqrt{x^2+y^2}\,dx = 0,即 12d(x2+y2)+x2+y2 dx=0\dfrac{1}{2}d(x^2+y^2) + \sqrt{x^2+y^2}\,dx = 0,取积分因子 μ=1x2+y2\mu = \dfrac{1}{\sqrt{x^2+y^2}},得

d(x2+y2)+dx=0d\left(\sqrt{x^2+y^2}\right) + dx = 0

故通解为 x2+y2+x=c\sqrt{x^2+y^2} + x = c。

例 8:求解方程 ydx−(x+y) dy=0ydx - (x + y)\,dy = 0 的通解(四种方法)。

解:

  • 方法 1(积分因子):∂M∂y−∂N∂x=1−(−1)=2\dfrac{\partial M}{\partial y} - \dfrac{\partial N}{\partial x} = 1 - (-1) = 2,∂N∂x−∂M∂yM=−2y\dfrac{\dfrac{\partial N}{\partial x} - \dfrac{\partial M}{\partial y}}{M} = -\dfrac{2}{y} 仅依赖 yy,μ=e∫−2/y dy=1y2\mu = e^{\int -2/y\,dy} = \dfrac{1}{y^2}。乘以 μ\mu 得 y dx−x dyy2−dyy=0\dfrac{y\,dx - x\,dy}{y^2} - \dfrac{dy}{y} = 0,即 d(xy)−d(ln⁡y)=0d\left(\dfrac{x}{y}\right) - d(\ln y) = 0,通解 xy−ln⁡y=c\dfrac{x}{y} - \ln y = c,即 x=yln⁡y+cyx = y\ln y + cy。
  • 方法 2(凑微分):ydx−xdy=y dyydx - xdy = y\,dy,ydx−xdyy2=dyy\dfrac{ydx - xdy}{y^2} = \dfrac{dy}{y},得 d(xy)=d(ln⁡y)d\left(\dfrac{x}{y}\right) = d(\ln y),通解 x=yln⁡y+cyx = y\ln y + cy。
  • 方法 3(齐次方程):dydx=yx+y\dfrac{dy}{dx} = \dfrac{y}{x + y} 是齐次方程,令 u=yxu = \dfrac{y}{x},解得 ln⁡∣x∣=ln⁡∣u1−u∣+c\ln\vert x\vert = \ln\left\vert \dfrac{u}{1-u}\right\vert + c[?],变量还原得 x=yln⁡y+cyx = y\ln y + cy。
  • 方法 4(线性方程):改写为 dxdy−1yx=1\dfrac{dx}{dy} - \dfrac{1}{y}x = 1,以 yy 为自变量、xx 为未知函数的线性方程,公式 (3) 得 x=yln⁡y+cyx = y\ln y + cy。

校对点

例 8 讲义最终通解为 xy+ln⁡y=c\dfrac{x}{y} + \ln y = c(即 x+yln⁡y=cyx + y\ln y = cy)[?],与上面 x=yln⁡y+cyx = y\ln y + cy 符号不同,请以原讲义为准。

作业(2.3):p41 1.(1)(5);2.(1)(4)(7)(9);3、4、5、6、9


§2.4 一阶隐方程与参数表示 ​

一阶隐式方程 F(x,y,y′)=0F(x, y, y') = 0,当 y′y' 未能解出或解出形式复杂时,采用引进参数的办法使其变为导数已解出的方程类型。主要研究四种类型:

  1. y=f(x,y′)y = f(x, y')
  2. x=f(y,y′)x = f(y, y')
  3. F(x,y′)=0F(x, y') = 0(不显含 yy)
  4. F(y,y′)=0F(y, y') = 0(不显含 xx)

参数形式解的定义:如果存在定义在 (α,β)(\alpha, \beta) 上的函数 φ(t),ψ(t)\varphi(t), \psi(t) 使 x=φ(t),y=ψ(t)x = \varphi(t), y = \psi(t),且 F(φ(t),ψ(t),φ′(t)/ψ′(t))=0F(\varphi(t), \psi(t), \varphi'(t)/\psi'(t)) = 0[?],则称其为方程的参数形式解。

一、可解出 y′y' 的方程:y=f(x,y′)y = f(x, y') ​

设 y′=py' = p,方程变为 y=f(x,p)y = f(x, p),两边对 xx 求导,以 p′=dpdxp' = \dfrac{dp}{dx} 代入:

p=∂f∂x+∂f∂pdpdxp = \frac{\partial f}{\partial x} + \frac{\partial f}{\partial p}\frac{dp}{dx}

这是关于变量 x,px, p 的一阶微分方程。求出其通解后回代:

  • 若通解为 p=φ(x,c)p = \varphi(x, c),代入 y=f(x,p)y = f(x, p) 即得通解 y=f(x,φ(x,c))y = f(x, \varphi(x,c))
  • 若通解为 x=ψ(p,c)x = \psi(p, c),得参数形式通解 {x=ψ(p,c)y=f(ψ(p,c),p)\begin{cases} x = \psi(p, c) \\ y = f(\psi(p,c), p) \end{cases}
  • 若通解为 Φ(x,p,c)=0\Phi(x, p, c) = 0,得参数形式通解 {Φ(x,p,c)=0y=f(x,p)\begin{cases} \Phi(x, p, c) = 0 \\ y = f(x, p) \end{cases}

例 1:求解方程 y=(dydx)2−2xdydx+2x2y = \left(\dfrac{dy}{dx}\right)^2 - 2x\dfrac{dy}{dx} + 2x^2[?]。

解:令 p=y′p = y',则 y=p2−2xp+2x2y = p^2 - 2xp + 2x^2,两边对 xx 求导:

p=2pdpdx−2p−2xdpdx+4xp = 2p\frac{dp}{dx} - 2p - 2x\frac{dp}{dx} + 4x

整理得 (p−2x)(dpdx−1)=0(p - 2x)\left(\dfrac{dp}{dx} - 1\right) = 0[?],解得 p=x+cp = x + c 或 p=2xp = 2x。

  • p=x+cp = x + c 代入得通解 y=c2+cx+x2y = c^2 + cx + x^2[?]
  • p=2xp = 2x 代入得一个解 y=2x2y = 2x^2

曲线 y=2x2y = 2x^2 是积分曲线族的包络(通解中 c=2xc=2x 时的极限情形),在微分方程中称 y=2x2y = 2x^2 为原方程的奇解。

校对点

例 1 讲义通解为 y=c2+2cx+x2y = c^2 + 2cx + x^2[?](原文 "12p2+...\frac{1}{2}p^2 + ..." 混乱),奇解 y=2x2y = 2x^2 明确。请对照原讲义核实系数。

二、可解出 xx 的方程:x=f(y,y′)x = f(y, y') ​

设 y′=py' = p,方程变为 x=f(y,p)x = f(y, p),两边对 yy 求导,以 dxdy=1p\dfrac{dx}{dy} = \dfrac{1}{p} 代入:

1p=∂f∂y+∂f∂pdpdy\frac{1}{p} = \frac{\partial f}{\partial y} + \frac{\partial f}{\partial p}\frac{dp}{dy}

得到关于变量 y,py, p 的一阶微分方程,求出通解 Φ(y,p,c)=0\Phi(y, p, c) = 0 后得参数形式通解 {Φ(y,p,c)=0x=f(y,p)\begin{cases} \Phi(y,p,c) = 0 \\ x = f(y, p) \end{cases}。

例 2:求解方程 x(dydx)3−ydydx+2=0x\left(\dfrac{dy}{dx}\right)^3 - y\dfrac{dy}{dx} + 2 = 0[?]。

解:设 p=y′p = y',则 x=yp−2p3x = \dfrac{yp - 2}{p^3}(p≠0p \neq 0),两边对 yy 求导得关于 y,py, p 的方程,解得 yp4=c+2pyp^4 = c + 2p[?],参数形式通解为

{x=p4c+2p3−2pp4[?]y=c+2pp4[?]\begin{cases} x = \dfrac{p^4 c + 2p^3 - 2p}{p^4}[?] \\ y = \dfrac{c + 2p}{p^4}[?] \end{cases}

另有解 y=0y = 0。

校对点

例 2 讲义最终参数通解:x=2p3−cp4x = \dfrac{2p^3 - c}{p^4},y=2p+cp4y = \dfrac{2p + c}{p^4}(p≠0p \neq 0 为参数),请以原讲义为准。

三、不显含 yy 的方程:F(x,y′)=0F(x, y') = 0 ​

设 y′=py' = p,方程 F(x,p)=0F(x, p) = 0 表示 xpxp 平面上的一条曲线。若能找到参数表示 x=φ(t),p=ψ(t)x = \varphi(t), p = \psi(t)(满足 F(φ(t),ψ(t))=0F(\varphi(t), \psi(t)) = 0),则由 dy=p dx=ψ(t)φ′(t) dtdy = p\,dx = \psi(t)\varphi'(t)\,dt 积分得

y=∫ψ(t)φ′(t) dt+cy = \int \psi(t)\varphi'(t)\,dt + c

参数形式通解为 {x=φ(t)y=∫ψ(t)φ′(t) dt+c\begin{cases} x = \varphi(t) \\ y = \int \psi(t)\varphi'(t)\,dt + c \end{cases}。关键一步(也是最困难一步)是找曲线的参数表示。

例 3:求解方程 y′=x1+(y′)2y' = x\sqrt{1 + (y')^2}[?]。

解:设 p=y′p = y',则 p=x1+p2p = x\sqrt{1+p^2},即 x=p1+p2x = \dfrac{p}{\sqrt{1+p^2}}。取参数 tt:令 p=tan⁡tp = \tan t,则 x=sin⁡tx = \sin t。由 dy=p dx=tan⁡tcos⁡t dt=sin⁡t dtdy = p\,dx = \tan t \cos t\,dt = \sin t\,dt,积分得 y=−cos⁡t+cy = -\cos t + c。参数形式通解

{x=sin⁡ty=c−cos⁡t\begin{cases} x = \sin t \\ y = c - \cos t \end{cases}

消去参数得 (x)2+(y−c)2=1(x)^2 + (y - c)^2 = 1,即单位圆族。

四、不显含 xx 的方程:F(y,y′)=0F(y, y') = 0 ​

设 y′=py' = p,方程 F(y,p)=0F(y, p) = 0 表示 ypyp 平面上的一条曲线,参数表示 y=φ(t),p=ψ(t)y = \varphi(t), p = \psi(t)。由 dx=dyp=φ′(t)ψ(t) dtdx = \dfrac{dy}{p} = \dfrac{\varphi'(t)}{\psi(t)}\,dt 积分得

x=∫φ′(t)ψ(t) dt+cx = \int \frac{\varphi'(t)}{\psi(t)}\,dt + c

参数形式通解为 {y=φ(t)x=∫φ′(t)ψ(t) dt+c\begin{cases} y = \varphi(t) \\ x = \int \dfrac{\varphi'(t)}{\psi(t)}\,dt + c \end{cases}。

例 4:求解方程 y=(y′)2−2y′+1y = (y')^2 - 2y' + 1[?]。

解:设 p=y′p = y',令 y=t2−2t+1=(t−1)2y = t^2 - 2t + 1 = (t-1)^2,p=tp = t[?]。由 dx=dyp=2(t−1)t dt=2(1−1t)dtdx = \dfrac{dy}{p} = \dfrac{2(t-1)}{t}\,dt = 2\left(1 - \dfrac{1}{t}\right)dt,积分 x=2t−2ln⁡∣t∣+cx = 2t - 2\ln\vert t\vert + c。参数形式通解

{y=(t−1)2x=2t−2ln⁡∣t∣+c\begin{cases} y = (t-1)^2 \\ x = 2t - 2\ln\vert t\vert + c \end{cases}

此外 y=±1y = \pm 1 也是方程的解[?]。

五、利用变量代换的微分方程积分法 ​

当 F(x,y,y′)=0F(x, y, y') = 0 中 x,y,y′x, y, y' 都不易解出,或虽能解出但积分较复杂时,可先作适当的变量代换再求解。如何选择变量代换没有一定规律,需要大量练习积累经验。

例 5:求解方程 y(y′)2cos⁡x+y′cos⁡x−xy′−x=0y(y')^2\cos x + y'\cos x - xy' - x = 0[?](令 u=sin⁡y,v=sin⁡xu = \sin y, v = \sin x)。

解:令 u=sin⁡y,v=sin⁡xu = \sin y, v = \sin x,则 du=cos⁡y dy,dv=cos⁡x dxdu = \cos y\,dy, dv = \cos x\,dx,y′=cos⁡xcos⁡ydvduy' = \dfrac{\cos x}{\cos y}\dfrac{dv}{du}[?]。代入原方程化为关于 u,vu, v 的克莱洛(Clairaut)方程 u=vdudv+(dudv)2u = v\dfrac{du}{dv} + \left(\dfrac{du}{dv}\right)^2,其通解 u=cv+c2u = cv + c^2,奇解 u=−v24u = -\dfrac{v^2}{4}。变量还原得

sin⁡y=csin⁡x+c2,sin⁡y=−sin⁡2x4\sin y = c\sin x + c^2, \qquad \sin y = -\frac{\sin^2 x}{4}

例 6:求方程 y(y′)2−2xy′+y=0y(y')^2 - 2xy' + y = 0 的通解[?](令 u=y2,v=x2u = y^2, v = x^2)。

解:令 u=y2,v=x2u = y^2, v = x^2,代入化为克莱洛方程,通解与奇解分别为[?]:

y=cx−c22[?],x2=2p2+2p2(1+p2)2[?]y = cx - \frac{c^2}{2}[?], \qquad x^2 = 2p^2 + \frac{2p^2}{(1+p^2)^2}[?]

校对点

例 5、例 6 讲义参数代换细节提取严重混乱(克莱洛方程通解 u=cv+c2u = cv + c^2、u=−v24u = -\frac{v^2}{4} 等),此处仅给出结构框架,建议对照原讲义第 26–27 页核实。

作业(2.4):p49-50:1.(15、19、20、25、29、30);2、3、4;5(只做证明部分)


本章小结 ​

  • 一阶方程初等解法体系:变量分离 → 齐次方程 → 线性方程(常数变易法)→ 伯努利方程 → 恰当方程(积分因子)→ 隐方程(参数法)
  • 能初等解法的微分方程很有限,如 Riccati 方程 dydx=P(x)y2+Q(x)y+R(x)\dfrac{dy}{dx} = P(x)y^2 + Q(x)y + R(x) 一般没有初等解法

相关笔记 ​

  • ODE常微分方程
  • 第一章 绪论

用 VitePress 搭建 · 托管在 GitHub Pages