初等积分法
Last updated
Last updated
从形式上,对于一阶常微分方程,我们希望找到连续函数 $P$ 和 $Q$
定义——恰当方程
若存在连续函数 $\phi(x,y)$ 满足
则称 $P\text{d}x+Q\text{d}y$ 为恰当的,此时
定理:由 $\phi=C$ 确定的函数 $y=u(x)$ 或 $x=u(y)$ 满足原方程
证明:设 $\phi=C$ 确定 $y=u(x)$
从而 $y=u(x)$ 是方程 $P\text{d}x+Q\text{d}y=0$ 的解
反之,若 $y=u(x)$ 是原方程的解,也一定有 $\phi(x,u(x))=C$
例:
进而有一系列问题
如何判断 $P\text{d}x+Q\text{d}y=0$ 为恰当的
$P\text{d}x+Q\text{d}y=0$ 恰当 $\Leftrightarrow$ $\exists$ 连续函数 $\phi$,$\phi'_x=P, \phi'_y=Q$
一般情况很难判断,但如果我们有更强的条件,我们就可以有以下定理
定理:$P,Q\in C^1(D)$,$D\subset R^2$ 是单连通区域,则 $P\text{d}x+Q\text{d}y=0$ 恰当等价于
证明:
充分性:设 $P\text{d}x+Q\text{d}y=0$ 恰当,则 $\exists$ 连续函数 $\phi$,$\phi'_x=P, \phi'_y=Q, (x,y)\in D$,由于 $\phi\in C^2(D)$
必要性:设我们已经有 ${\partial P}/{\partial y}={\partial Q}/{\partial x}, (x,y)\in D$,取定 $p_0(x_0,y_0)\in D$,则 $\forall p(x,y)\in D$,可以定义与路径无关的积分
它显然满足 $\text{d}\phi=P\text{d}x+Q\text{d}y$
具体地,该积分可以构造在一个与坐标轴平行的矩形区域上,即
或
取决于先对 $x$ 或 $y$ 积分
如果是恰当的,如何求 $\phi$
例:
$\partial P/\partial y=e^x+2y=\partial Q/\partial x$,从而原方程是恰当的
构造 $\phi(x,y)=ye^x+2e^x+y^2x+f(y)$
则不妨直接另 $f(y)\equiv0$
如果不是恰当的,如何求 $\phi$
定义:称 $P\text{d}x+Q\text{d}y=0$ 为分离变量型,若
一般这样的方程不是恰当的,但是若 $Q_1(x)P_2(y)\neq0$,我们有
显然是恰当的,因 $\phi$ 的二阶偏导数为零,通积分为
若 $Q_1(x_i)=P_2(y_j)=0$,则 $x\equiv x_i, y\equiv y_j$ 都是解
例:
此时的解为
$x\equiv0$ 和 $y\equiv\pm1$ 也是解
例:
若 $y\neq0$,极易解得
为此得到通积分
外加特解 $y=0$
此时在平面上几乎所有点 $P_0$,在局部有且仅有一条积分曲线(局部唯一),而例外是 $y=0$,在其上每一点的局部范围内都有无穷多条积分曲线(局部不唯一)
例——跟踪问题:设某 $A$ 从原点出发,沿 $x$ 轴正向前进,同时 $B$ 从 $(0,b)$ 开始跟踪,且与 $A$ 保持等距 $b$,求 $B$ 的光滑运动轨迹
其中
因此
因为
最终得到
例:
令 $y=ux$,则
若采用极坐标,解即为 $r=Ce^\theta$,事实上方程化为
$q(x)\equiv0,y\neq0$ 时
$y=0$ 显然也是一个解,则最终
$q(x)\neq0$
从而
通解
线性方程解的性质
$q(x)\equiv0$ 时(齐次线性方程),解恒为零或恒不为零
所有解都为整体的(大范围的),即 $p(x)$ 和 $q(x)$ 的定义域上该解恒有意义
齐次方程任两个解的和还是解,任意解乘以常数后还是解
方程的任意一个解与对应的齐次方程的通解的和构成方程的通解
方程的初值问题的解存在且唯一
例:
例:设 $f\in C^1[0,+\infty),a(x)\in C[0,\infty),a(x)\ge c_0>0$,已知
求证 $\lim_{x\to+\infty}f(x)=0$
设 $g(x)=f'(x)+a(x)f(x)\in C[0,+\infty)$,则
则
例:若 $p(x)$ 和 $q(x)$ 均为 $2\pi$ 周期的,且
有唯一的有界解,求证这个有界解为 $2\pi$ 周期的
设 $y(x)$ 是这个有界解,则 $y(x+2\pi)$ 仍然为有界解,由有界解的唯一性,$y(x)=y(x+2\pi)$
令 $y=u(x)$,则
例:
$c_1=c_2=0$,则有
$a_1b_2-a_2b_1=0$,则有
$a_1b_2-a_2b_1\neq0$,设 $\alpha,\beta$ 是
的解,则有
令 $X=x-\alpha,Y=y-\beta$,有
化为了齐次方程
事实上,用类似的方法,形如
的方程我们也可以解
若 $n>0$,则 $y\equiv0$ 是解,下设 $y\neq0$
令 $z=y^{1-n}$,则
是一阶线性微分方程
若 $r\equiv0$,转为 Bernoulli 方程,否则一般难以直接解出
定理:若已知一个解 $y=\varphi(x)$,则可以求通解
证明:令 $y=u(x)+\varphi(x)$,有
化为 Bernoulli 方程
定理:$P(x)\equiv a,Q(x)\equiv0,r(x)=bx^m$ 时能化为可积分求解的方程,其中 $a,b,m$ 是常数,$m=0,-2,-4k/(2k\pm1)$,$k=1,2,\cdots$
证明:不妨设 $a=1$,将方程记为 $Tm$,$T_0$ 显然是分离变量的,对于 $T{-2}$,作变换 $z=xy$
同样是变量分离的
一般情形下,记 $X_k^\pm$ 对应 $m=-4k/(2k\pm1)$ 的方程,对 $X_k^+$ 作变换
再作变换
即将方程化为了一个 $X^+_{k-1}$ 的方程;如此重复 $k$ 次即可化为 $m=0$ 的方程并求解
对于 $X^-_{k}$,情况完全类似
定义:
若非零函数 $\mu(x,y)$ 满足 $\mu P\text dx+\mu Q\text dy=0$ 是恰当的,则称 $\mu$ 为原方程的一个积分因子
定理:若 $P,Q$ 为 $C^1$ 的,则 $C^1$ 函数 $\mu$ 为积分因子,当且仅当
即解一个一般的一阶常微分方程的难度相当于解一个一阶线性偏微分方程,一般而言解不出来
一些特殊情形
例:可分离变量的方程
取 $\mu=Q_1^{-1}(x)P_2^{-1}(y)$
例:一阶线性常微分方程
取
定理:方程 $P(x,y)\text dx+Q(x,y)\text dy=0$ 有一个只依赖于 $x$ 的 $C^1$ 积分因子,当且仅当
只依赖于 $x$;有一个只依赖于 $y$ 的 $C^1$ 积分因子,当且仅当
只依赖于 $y$,并且
证明:$\mu$ 为只依赖于 $x$ 的积分因子等价于
另一情况的证明完全类似
例:回到上例中一阶线性常微分方程的情形
积分因子为
定理:方程有形如 $\mu=f[\phi(x,y)]$ 的积分因子,当且仅当
证明:
例:
与 $y$ 无关,则 $\mu=e^{3x}$ 是个积分因子
定理:若 $\mu(x,y)$ 是一个积分因子,使得
则 $\mu(x,y)f[\Phi(x,y)]$ 也是积分因子,其中 $f$ 是任意连续函数
证明:
分组求积分因子法
假设方程左端可以分成两组,即
其中第一组和第二组各有积分因子 $\mu_1,\mu_2$
从而如果能适当选取 $f_1$ 和 $f_2$,使得 $\mu=\mu_1f_1(\Phi_1)=\mu_1f_1(\Phi_1)$,则 $\mu$ 就是方程的一个积分因子
例:
第一部分积分因子 $\mu_1=x^{-3}$,且
第二部分积分因子$\mu_2=y^{-2}$,且 $\Phi_2=x=C$
从而需要寻找可微函数 $g_1$ 和 $g_2$,使得
只需取 $g_1(xy)=(xy)^{-2},g_2(x)=x^{-5}$,得到原方程的积分因子
进而得到全微分
定理:设 $\mu_1,\mu_2$ 都为积分因子,若 $\mu_1/\mu_2$ 不为常数,则 $\mu_1/\mu_2$ 是方程的通解
证明:设 $\mu_1P\text dx+\mu_2Q\text dy=\text d\phi_1, \mu_2P\text dx+\mu_2Q\text dy=\text d\phi_2$
设 $\omega=P\text dx+Q\text dy$
例:当 $P,Q$ 为 $n$ 次齐次函数的时候,求积分因子
令 $y=ux$,
有积分因子
从而原方程有积分因子
也可直接证明它为积分因子,由齐次性
从而
同理
定义:给定一族曲线 $\phi(x,y,c)=0$ ,找一族曲线 $\psi(x,y,k)=0$,使得两组曲线之间夹角为等值;若夹角为直角,则称第二族曲线是第一族曲线的正交轨线族
代入 $c=c(x,y)$
设与上述曲线族夹角为 $\alpha$ 的曲线斜率为 $y'$,则若 $\alpha\neq\pi/2$
若 $\alpha=\pi/2$
例:求与抛物线族 $y=cx^2$ 正交的曲线族
其正交曲线族满足的方程为
例:求与曲线 $y=x\ln(cx)$ 夹角为 $-\pi/4$ 的曲线族
其正交曲线族满足的方程为
令 $y=ux$,则当 $y\neq-x$ 时
同时 $y=-x$ 也是一个解
见习题 2-3 的第 6 题
记 $N(t)$ 为 $t$ 时刻病人数(足够大),有
令 $\Delta t\to 0$,有
$r$ 被称为传染率,若数据显示一周之内病人变为原来十倍,则可以估计 $r$
然而由于疾病只能传播给健康人,随着健康人数的减少,传染率应该会下降,假设它是线性的
食肉动物吃食草动物,食草动物吃草,可以用一组方程表现两种动物的种群大小是如何相互制约的
设 $x$ 是捕食者数量,$y$ 是被捕食者数量
最简单的假设:捕食者增长率与被捕食者数目线性正相关,被捕食者增长率与捕食者数目线性负相关
消去时间参量有
首先 $F$ 是一个凸函数,因为
从而 $\forall K\in(F(x_0,y_0),+\infty)$,$F(x,y)=K$ 是一个闭凸曲线,绕着 $(x_0,y_0)$ 旋转
很难求出 $x(t)$ 和 $y(t)$ 的具体表达式,但是可以求平均值,假设周期为 $T$
两边积分
同理
若对两种动物都进行捕杀,方程组变为
平均值变为
即捕食者数目减少,被捕食者数目增加
从 $(0,0)$ 至 $(x,y)$
历史上,这是约翰·伯努利的证明
任意一点处满足“光的折射定律”——入射角与出射角的正弦值之比为常数,该常数为两种介质中的光速之比,当物体纵坐标为 $y$ 时,有
这里 $\sin\theta$ 满足
令 $\tan \theta=\sqrt{\frac{y}{C_1-y}}$,则
令 $C=C_1/2, \phi=2\theta$,则
这是摆线的方程