常微分方程(4)

奇解

一阶隐式微分方程

本节讨论一阶隐式方程 \[ F(x, y,\frac{\mathrm{d}y}{\mathrm{d}x})=0 \tag{1} \] 的几个特殊解法,隐式指的是 \(\frac{\mathrm{d}y}{\mathrm{d}x}\)没有预先表示为 \((x, y)\)的显函数.

微分法

设从微分方程(1)中可显式解出未知函数 \[ y=f(x, p) \tag{2} \] 其中 \(p=\frac{\mathrm{d}y}{\mathrm{d}x}\). 设函数 \(f(x,p)\)\((x,p)\)是连续可微的. 则由方程(2)对 \(x\)进行微分,得 \[ p=f'_x(x,p)+f'_p(x,p)\frac{\mathrm{d}p}{\mathrm{d}x} \]\[ [f'_(x,p)-p]\mathrm{d}x+f'_p(x,p)\mathrm{d}p=0 \tag{3} \] 这是一个关于变量 \(x\)\(p\)的一阶显式微分方程.

求得(3)的通解 \(p=u(x,C)\),得到(2)的通解 \[ y=f(x,u(x,C)) \] 其中 \(C\)是一个任意常数;另外若(3)有特解 \(p=w(x)\)则方程(2)有相应的特解 \[ y=f(x,w(x)) \] 另一方面,若(3)的通解可写成 \(x=v(p,C)\)的形式,则(2)的通解可写成 \[ \begin{cases} x=v(p,C) \\ y=f(v(p,C),p) \end{cases} \] 这里 \(p\)视作一个参变量;同样如果(3)有特解 \(x=z(p)\),则(2)有相应的特解 \[ \begin{cases} x=z(p) \\ y=f(z(p),) \end{cases} \]

参数法

设微分方程不明显包含自变量,即 \[ F(y,p)=0 \quad \biggl(p=\frac{\mathrm{d}y}{\mathrm{d}x}\biggr) \tag{4} \] 作为变元 \(y\)\(p\)之间的联系,方程(4)在 \((y,p)\)平面上一般表示若干条曲线. 设 \[ y=g(t),\quad p=h(t) \tag{5} \] 是其中一条. 称(5)为(4)的一个参数表示. 一般来说 \(g(t)\)\(g'(t)\)\(h(t)\)都是参数 \(t\)的连续函数,并且 \(h(t) \neq 0\). 根据上述微分方程的参数表示,我们有 \[ \mathrm{d}x=\frac{1}{p}\mathrm{d}y=\frac{g'(t)}{h(t)}\mathrm{d}t \] 再利用积分,可得 \[ x=\int_{}^{} \frac{g'(t)}{h(t)} \mathrm{d}t+C \] 因此,微分方程(4)有通解 \[ x=\int_{}^{} \frac{g'(t)}{h(t)} \mathrm{d}t+C, \quad y=g(t) \tag{6} \]

一般而言,一阶隐式微分方程 \[ F(x, y,p)=0 \quad \biggl( p=\frac{\mathrm{d}y}{\mathrm{d}x} \biggr) \tag{7} \]\((x, y, p)\)空间表示曲面. 设它的参数表达式为 \[ x=f(u,v), \quad y=g(u,v), \quad p=h(u,v) \] 这里 \(u\)\(v\) 是两个参数. 则 \(\mathrm{d}y=p\mathrm{d}x\)可写成如下形式 \[ M(u,v)\mathrm{d}u+N(u,v)\mathrm{d}v=0 \tag{8} \] 其中 \[ \begin{cases} M(u,v)=g'_u(u,v)-h(u,v)f'_u(u,v) \\ N(u,v)=g'_v(u,v)-h(u,v)f'_v(u,v) \end{cases} \] 若能求得一阶显式微分方程(8)的通解 \[ v=Q(u,C) \tag{9} \] 则微分方程(7)有通解 \[ x=f(u,Q(u,C)), \quad y=g(u,Q(u,C)) \] 其中 \(u\)是参变量,而 \(C\)是一个积分常数;若其有特解 \(v=S(u)\),则 \[ x=f(u,S(u)),\quad y=g(u,S(u)) \] 是微分方程(7)的特解.

奇解

【定义】设一阶微分方程 \[ F(x, y, \frac{\mathrm{d}y}{\mathrm{d}x})=0 \tag{10} \] 有一特解 \[ \Gamma: y=\varphi(x) \quad(x\in J) \] 如果对每一点 \(Q\in \Gamma\),在 \(Q\)任何领域内方程(10)有一个不同于 \(\Gamma\)的解在 \(Q\)点与 \(\Gamma\)相切,则称 \(\Gamma\)是微分方程(10)的奇解.

定理4.1(奇解存在的必要条件):设函数 \(F(x, y,p)\)\((x, y, p)\in G\)是连续的,而且对 \(y\)\(p\)有连续的偏导数 \(F'_y\)\(F'_p\),若函数 \(y=\varphi(x),(x\in J)\)是微分方程(10)的一个奇解,并且 \[ (x,\varphi(x)\varphi'(x))\in G \quad(x\in J) \] 则奇解 \(y=\varphi(x)\)满足一个称之为p-判别式的联立方程 \[ F(x, y, p)=0, \quad F'_p(x, y, p)=0 \quad(p=\frac{\mathrm{d}y}{\mathrm{d}x}) \tag{11} \] 设从(11)中消去 \(p\)得到方程 \[ \Delta(x, y)=0 \tag{12} \] 则称由此所决定的曲线为方程(10)的P-判别曲线. 因此,微分方程(10)的奇解是一条p-判别曲线.

:用反证法,假设 \(\exists x_0\in J\),使得 \[ F'_p(x_0, y_0, p_0) \neq 0 \] 利用隐函数定理,\(F\)关于 \(p\)的偏导数不为零保证在 \((x_0,y_0)\)附近唯一确定函数 \[ \frac{\mathrm{d}y}{\mathrm{d}x}=f(x, y) \] 推得 \(f(x, y)\)\((x_0, y_0)\)附近关于 \(y\)有连续的偏导数,再用Picard定理知在 \((x_0,y_0)\)附近有唯一解,从而与奇解的假设矛盾.

注意,由p-判别式确定的函数 \(y=\psi(x)\)不一定是相应微分方程的解;即使是解,也不一定是奇解.

上面的定理要求出通解才能验证奇解与否,下面的定理在某种条件下克服了这一困难.

定理4.2:设函数 \(F(x, y, p)\)\((x, y,p)\in G\)是二阶连续可微的. 又设由微分方程(10)的p-判别式 \[ F(x, y, p)=0, \quad F'_p(x, y, p)=0 \] 消去\(p\)后得到的函数 \(y=\psi(x)(x\in J)\)是微分方程(10)的解(注意上面定理的注).而且设条件 \[ F'_y(x,\psi(x),\psi'(x)) \neq 0,\quad F''_{pp}(x,\psi(x),\psi'(x)) \neq 0 \tag{13} \] 以及 \[ F'_p(x,\psi(x),\psi'(x))=0 \tag{14} \]\(x\in J\)成立. 则 \(y=\psi(x)\)是微分方程(10)的奇解.

:放到了文章最后.

我们说明该定理中(13)(14)都是不可去的.

考虑方程 \[ \biggl( \frac{\mathrm{d}y}{\mathrm{d}x} \biggr)^{2}-y^{2}=0 \] 得到 \(y=0\),此时有 \(F'_y(x,0,0)=0\),而原方程的通解为 \(y=x\exp {\pm x}\).

考虑方程 \[ \sin (y \frac{\mathrm{d}y}{\mathrm{d}x})=y \] 得到 \(y=0\),此时有 \(F''_{pp}(x,0,0)=0\)\(y=0\)不是原方程的奇解.

考虑方程 \[ y=2x+y'-\frac{1}{3}(y')^3 \] 此时 \(F'_p(x,\psi(x),\psi'(x))=0\)

另一个类似于2.2中的结果:设连续函数 \(E(y)\)满足条件: \[ \begin{cases} E(y)=0, \quad y=0\\ E(y)\neq 0, \quad 0<y\leqslant 1 \end{cases} \]\(y=0\)是微分方程 \[ \frac{\mathrm{d}y}{\mathrm{d}x}=E(y) \] 的奇解当且仅当瑕积分 \[ \int_{0}^{1} \frac{\mathrm{d}y}{E(y)} \] 收敛. \(\Rightarrow\):若 \(y=0\)为奇解,则存在局部解使得 \(y\neq 0\),故有 \(\frac{\mathrm{d}y}{E}=\mathrm{d}x \Rightarrow \int_{0}^{1} \frac{\mathrm{d}y}{E}\)收敛.

\(\Leftarrow\):$ _{0}^{1} $收敛,可令 \(x(y)=x_0+\int_{0}^{y} \frac{\mathrm{d}y}{E}\)是异于 \(y=0\)的解,由于 \(y'(0)=0\),得到 \(y=0\)是奇解.

包络

这一节想利用有关曲线族的包络的概念来阐明奇解与通解之间的联系,以及讨论寻求奇解的方法.

设单参数 \(C\)的曲线族 \[ K(C): \qquad V(x,y,C)=0 \tag{15} \] 其中函数 \(V(x, y, C)\)\((x, y, C)\in D\)是连续可微的.

【定义】设在平面上有一条连续可微的曲线 \(\Gamma\). 如果对于任一点 \(q\in \Gamma\),在曲线族(15)中都有一条曲线 \(K(C^*)\)通过 \(q\)点且在该点与 \(\Gamma\)相切,而且 \(K(C^*)\)\(q\)点的某一领域内不同于 \(\Gamma\). 则称曲线 \(\Gamma\)为曲线族(15)的一支包络.

注:这个定义与微分几何中的定义稍有不同.

定理4.3:设微分方程 \[ F(x, y, \frac{\mathrm{d}y}{\mathrm{d}x})=0 \tag{16} \] 有通积分为 \[ U(x, y, C)=0 \tag{17} \] 又设(积分)曲线族(17)有包络为 \[ \Gamma: \qquad y=\varphi(x) \quad (x\in J) \] 则包络 \(y=\varphi(x)\)是微分方程(16)的奇解.

:任取 \(\Gamma\)上一点,由包络的性质即可.

由奇解的定义可知,奇解是通解的包络. 因此,由上定理可知,求微分方程的奇解归结到求它的通积分的包络. 通积分的包络就是原方程的奇解.

定理4.4:设 \(\Gamma\)是曲线族(15)的一支包络. 则它满足如下的C-判别式\[ V(x, y, C)=0, \quad V'_C(x, y, C)=0 \tag{18} \] 或(消去 \(C\),所得到的关系式) \[ \Omega(x, y)=0 \tag{19} \]

:设包络的参数方程为 \[ \begin{cases} x=f(C) \\ y=g(C) \end{cases} \] 根据切向量共线得出一个式子,然后对 \(V(x, y,C)=0\)\(C\)求导即可.

满足C-判别式的曲线未必是相应曲线族的包络,如 \[ (y-1)^{2}\biggl(\frac{\mathrm{d}y}{\mathrm{d}x}\biggr)^{2}=\frac{4}{9}y \] 解它的通积分和C-判别式得 \(y=0\)\(y=3\),然而只有前一个是曲线的包络.

下面的定理给出了包络的一个充分条件. 定理4.5:设由曲线族(15)的C-判别式 \[ V(x, y,C)=0 \quad V'_C(x, y, C)=0 \] 确定一支连续可微且不含于族(15)的曲线 \[ \Lambda: \quad x=\varphi(C), \quad y=\psi(C) \quad (C\in J) \] 而且它满足非蜕化性条件(即切向量不退化) \[ (\varphi'(C),\psi'(C))\neq (0,0), \quad (V'_x,V'_y)\neq (0,0) \] 其中 \(V'_x=V'_x(\varphi(C),\psi(C),C)\)\(V'_y=V'_y(\varphi(C),\psi(C),C)\). 则 \(\Lambda\)是曲线族(15)的一支包络. :在 \(\Lambda\)上任取一点,在这一点附近用隐函数定理找到一条曲线族中的曲线,然后证明它在这一点与 \(\Lambda\)相切,即其切向量共线.

该定理基本上只适用于能求出通积分的情形,使用时通常配合4.3,这是因为若能判定C-曲线不是积分曲线,则它不是解,从而不是奇解,从而不是包络.

奇解的存在定理

这是定理4.2的证明,可以参考教材117-120页 或者如下链接: ---

常微分方程学习笔记(4) https://zhuanlan.zhihu.com/p/97414951


下一章就进入高阶微分方程了.


常微分方程(4)
http://example.com/2022/04/06/常微分方程(4)/
Author
John Doe
Posted on
April 6, 2022
Licensed under