微分方程-隐式通解

隐式方程

因为会遇到一些导数未接触的一阶微分方程. 这里讨论一阶隐式方程,其一般形式为

\displaystyle F\left(x,y,\dfrac{\text{d}y}{\text{d}x}\right)=0 (2.46)

求解这类方程的基本思想是将 p=\dfrac{\text{d}y}{\text{d}x} 看成独立的变量而考虑把由代数方程 F(x,y,p)=0 所定义的 \mathbb{R}^3 上的曲面的参数化,再通过变量替换的方法把方程(2.46)化为导数已解出的显式方程,然后用之前的方法求解.

一般求解的具体做法:

第一步

将曲面 F(x,y,p)=0 表示成参数形式

\begin{cases} x=\phi(s,t),\\ y=\psi(s,t),\\ p=\kappa(s,t).\quad(2.47) \end{cases}

第二步

对(2.47)求 x,y 的微分,用 p=\dfrac{\text{d}y}{\text{d}x} 给出 \text{d}y\text{d}x 的关系:

\begin{aligned} &\text{d}x=\dfrac{\partial \phi}{\partial s}\text{d}s+\dfrac{\partial \phi}{\partial t}\text{d}t,&(2.48)\\ &\text{d}y=\dfrac{\partial \psi}{\partial s}\text{d}s+\dfrac{\partial \psi}{\partial t}\text{d}t,&(2.49)\\ &\text{d}y=\dfrac{\text{d}y}{\text{d}x}\text{d}x=\kappa\text{d}x,&(2.50) \end{aligned}

第三步

将 (2.48)、(2.49)带入(2.50)得

\displaystyle\dfrac{\partial \psi}{\partial s}\text{d}s+\dfrac{\partial \psi}{\partial t}\text{d}t=\kappa\left(\dfrac{\partial \phi}{\partial s}\text{d}s+\dfrac{\partial \phi}{\partial t}\text{d}t\right)

合并得到

\displaystyle\left(\dfrac{\partial \psi}{\partial s}-\dfrac{\partial \phi}{\partial s}\kappa\right)\text{d}s+\left(\dfrac{\partial \psi}{\partial t}-\dfrac{\partial \phi}{\partial t}\kappa\right)\text{d}t=0\;\;(2.51)

从而化成了对成型是的微分方程.

第四步

如果用学过的方法求出了方程(2.51)的通解 s=w(t,C),则将他带入(2.47)就得到方程(2.46)的参数形式的解

\begin{cases} x=\phi(w(t,C),t),\\ y=\psi(w(t,C),t).\\ \end{cases}\quad(2.52)

其中 C 为任意常数. 如果方程(2.51)的通解是另一种形式 t=w(s,C),我们可以得到类似结果.


讨论集中特殊形式的方程


可以解出 y 的方程

y=f\left(x,\,\dfrac{\text{d}y}{\text{d}x}\right),\quad(2.53)

这里函数 f 有连续的一阶偏导数. 这时曲面 F(z,y,p)=0 的参数形式可为

\begin{cases} x=x,\\ y=f\left(x,p\right),\\ p=p. \end{cases}

其中 x,p=\dfrac{\text{d}y }{\text{d}x} 为参数. 对方程(2.53)两边关于 x 求导,得

\displaystyle p=\dfrac{\partial f}{\partial x}(x,p)+\dfrac{\partial f}{\partial p}(x,p)\dfrac{\text{d}p}{\text{d}x}

整理可得到对称形式的方程:

\displaystyle\left(p-\dfrac{\partial f}{\partial x}(x,p)\right)\text{d}x-\dfrac{\partial f}{\partial p}(x,p)\text{d}p=0\quad(2.54)


可以解出 x 的方程:

x=f\left(y,\dfrac{\text{d}y}{\text{d}x}\right),\quad(2.55)

这里函数 f 有连续的一阶偏导数. 类似地曲面 F(x,y,p)=0 的参数形式可为

\begin{cases} x=f(y,p),\\ y=y,\\ p=p. \end{cases}

其中 y,p=\dfrac{\text{d}y}{\text{d}x} 为参数. 对方程(2.55)两边关于 x 求导,得

\displaystyle1=p\dfrac{\partial f}{\partial y}(y,p)+p\dfrac{\partial f}{\partial p}(y,p)\dfrac{\text{d}p}{\text{d}y}.

由上式可解出 \dfrac{\text{d}p}{\text{d}y},从而得到如下规范形式的一阶微分方程:

\displaystyle\dfrac{\text{d}p}{\text{d}y}=\dfrac{\dfrac{1}{p}-\dfrac{\partial f}{\partial y}(y,p)}{\dfrac{\partial f}{\partial p}(y,p)}.\quad(2.57)


不显含 y 的隐式方程

F\left(x,\dfrac{\text{d}y}{\text{d}x}\right)=0\quad(2.58)

p=\dfrac{\text{d}y}{\text{d}x},这时代数方程 F(x,p)=0 代表 Oxp 平面上的一条曲线,设该曲线有参数表示

\begin{cases} x=\varphi(s),\\ p=\psi(s), \end{cases}\quad(2.59)

其中 s 为参数. 由微分关系得

\begin{cases} \text{d}y=p\text{d}x=\psi(s)\text{d}x,\\ \text{d}x=\varphi(s)\text{d}s. \end{cases}

因此

\text{d}y=\psi(s)\varphi'(s)\text{d}s,

这是一个变量分离的方程,其通解为

\displaystyle y(s)=\int\psi(s)\varphi'(s)\text{d}s+C

其中 C 为任意常数. 由此得方程(2.58)的参数形式的通解为

\begin{cases} x=\varphi(s),\\ y=\int\psi(s)\varphi'(s)\text{d}s+C. \end{cases}\quad(2.60)

不显含 x 的隐式方程

F\left(y,\dfrac{\text{d}y}{\text{d}x}\right)=0\quad(2.61)

p=\dfrac{\text{d}y}{\text{d}x},同样,代数方程 F(y,p)=0 代表 Oyp 平面上的一条曲线,设其参数表示为

\begin{cases} y=\varphi(s),\\ p=\psi(s). \end{cases}\quad(2.62)

由微分关系得

\begin{cases} \text{d}y=\varphi'(s)\text{d}s\\ \text{d}y=p\text{d}x=\psi(s)\text{d}x. \end{cases}

因此,

\text{d}x=\dfrac{\varphi'(s)}{\psi(s)}\text{d}s.

故方程(2.61)的参数形式的通解为

\begin{cases} x=\int\dfrac{\varphi'(s)}{\psi(s)}\text{d}s+C.\\ y=\varphi(s). \end{cases}\quad(2.63)

其中 C 为任意常数.


例子


1

\displaystyle\left(\dfrac{\text{d}y}{\text{d}x}\right)^3+2x\dfrac{\text{d}y}{\text{d}x}-y=0

Sol:

p=\dfrac{\text{d}y}{\text{d}x},则有

y=p^3+2xp

对方称两边关于 x 求导,得

p=3p^2\dfrac{\text{d}p}{\text{d}x}+2x\dfrac{\text{d}p}{\text{d}x}+2p

(3p^2+2x)\text{d}p+p\text{d}x=0

p\not=0 时,方程有积分因子 \mu=p,用 \mu 乘 上述方程的两端,得

3p^3\text{d}p+(2xp\text{d}p+p^2\text{d}x)=0.

由此求出放曾的隐式通解:

\dfrac{3p^4}{4}+xp^2=C_1

其中 C_1 为任意常数. 接触 x

x=\dfrac{C_1-\dfrac34p^4}{p^2}=\dfrac{C-3p^4}{4p^2},

其中 C=4C_1. 从而原微分方程得参数形式的解为

\begin{cases} x=\dfrac{C-3p^4}{4p^2},\\ y=p^3+2p\dfrac{C-3p^4}{4p^2}=\dfrac{C-p^4}{2p} \end{cases}\quad(p\not=0).

p=0 时,由(2.65)可直接推知 y=0 也是微分方程的解.


2

x^3+\left(\dfrac{\text{d}y}{\text{d}x}\right)^3-3x\dfrac{\text{d}y}{\text{d}x}=0

Sol:
\dfrac{\text{d}y}{\text{d}x}=p=tx,则由方程得

\begin{cases} x=\dfrac{3t}{1+t^3},\\ p=\dfrac{3t^2}{2+t^3}.\\ \end{cases}

于是

\text{d}y=\dfrac{9(1-2t^3)t^2}{(1+t^3)^3}\text{d}t.

两边积分得

\displaystyle y=\int\dfrac{9(1-2t^3)t^2}{(1+t^3)^3}\text{d}t=\dfrac{3(1+4t^3)}{2(1+t^3)^2}+C

因此,方程的参数形式得通解为

\begin{cases} x=\dfrac{3t}{1+t^3},\\ y=\dfrac{3(1+4t^3)}{2(1+t^3)^2}+C. \end{cases}

其中 C 为任意常数.


3

\left(\dfrac{\text{d}y}{\text{d}x}\right)^3-y^2\left(4-\dfrac{\text{d}y}{\text{d}x}\right)=0

Sol:
\dfrac{\text{d}y}{\text{d}x}=p,\,y=pt,\, 由方程可得

\begin{cases} p=\dfrac{4t^2}{1+t^2},\\ y=\dfrac{4t^3}{1+t^2}.\\ \end{cases}

p\not=0 时,由 \text{d}x=\dfrac{1}{p}\text{d}y,则

\displaystyle x=\int\dfrac{3+t^2}{1+t^2}\text{d}t=t+2\arctan t+C.

因此,该方程的参数形式的解为

\begin{cases} x=t+2\arctan t+C,\\ y=\dfrac{4t^3}{1+t^2}.\\ \end{cases}

其中 C 为任意尝试. 此外,当 p=0 时,易知 y=0也是方程的解.


在某些情况下我们可从隐式方程中接触 \dfrac{\text{d}y}{\text{d}x},因此可以将方程化为前两节讨论过的显式方程. 例如:若方程的形式为

\displaystyle\sum_{k=0}^na_k(x,y)\left(\dfrac{\text{d}y}{\text{d}x}\right)^k=0

若该方程关于 \dfrac{\text{d}y}{\text{d}x} 多项式有 s 个不同的实根 f_k(x,y),\,k=1,2,\cdots,s,s\leqslant n, 则对每个 k,该方程的求解问题都可归结为形式较简单的显式方程

\dfrac{\text{d}y}{\text{d}x}=f_k(x,y)

的求解问题. 例如

y'^2-\left(xy+\dfrac{1}{2y}\right)y'+\dfrac{x}{2}=0

可以写成

\left(y'-\dfrac{1}{2y}\right)\left(y'-xy\right)=0

由此得两个方程

y'=\dfrac{1}{2y},\,\;y'=xy.

对这两个方程分别用分离变量法求解,从而得到原方程的不同解为

x=y^2+C_1,

y=C_2e^{\frac{x^2}{2}}

其中 C_1,\,C_2 为任意常数.


若方程不显含 xy,即方程的形式为

F\left(\dfrac{\text{d}y}{\text{d}x}\right)=0\quad(2.68)

这时若方程(2.68)至少有一个实根 y'=p ,则有 y=px+C. 将 p=\dfrac{y-C}{x} 带入方程(2.68),即得方程(2.68)的隐式通解

F\left(\dfrac{y-C}{x}\right)=0

其中 C 为任意常数. 例如方程

\left(\dfrac{\text{d}y}{\text{d}x}\right)^7-2\left(\dfrac{\text{d}y}{\text{d}x}\right)^4+5\left(\dfrac{\text{d}y}{\text{d}x}\right)^3+5=0

由于方程的左边是一个关于 \dfrac{\text{d}y}{\text{d}x} 的 7 次多项式,因此该方程至少有一个实根,故有隐式通解

\left(\dfrac{y-C}{x}\right)^7-2\left(\dfrac{y-C}{x}\right)^4+5\left(\dfrac{y-C}{x}\right)^3+5=0

其中 C 为任意常数.

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 216,544评论 6 501
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 92,430评论 3 392
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 162,764评论 0 353
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,193评论 1 292
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,216评论 6 388
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,182评论 1 299
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,063评论 3 418
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,917评论 0 274
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,329评论 1 310
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,543评论 2 332
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,722评论 1 348
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,425评论 5 343
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,019评论 3 326
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,671评论 0 22
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,825评论 1 269
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,729评论 2 368
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,614评论 2 353

推荐阅读更多精彩内容

  • 考试科目:高等数学、线性代数、概率论与数理统计 考试形式和试卷结构 一、试卷满分及考试时间 试卷满分为150分,考...
    Saudade_lh阅读 1,077评论 0 0
  • 2017年考研数学一大纲原文 考试科目:高等数学、线性代数、概率论与数理统计 考试形式和试卷结构 一、试卷满分及考...
    SheBang_阅读 624评论 0 7
  • 考试形式和试卷结构一、试卷满分及考试时间 试卷满分为150分,考试时间为180分钟 二、答题方式 答题方式为闭卷、...
    幻无名阅读 753评论 0 3
  • 今天是什么日子? 4月10日,108自律行动营结束倒计时第12天,也是2019年的第100天 天气:雨心情:还不错...
    搬砖的小姐姐阅读 326评论 0 0
  • 很像鸡汤,但又很有道理, 我67岁,这些是我学到的最好的东西: 你头20年都在担心别人对你的看法。你在接下来的20...
    前端混合开发阅读 278评论 0 3