微分方程-高阶微分方程

高阶微分方程

对于高阶微分方程

F\left(x,y,\dfrac{\text{d}y}{\text{d}x},\cdots,\dfrac{\text{d}^ny}{\text{d}x^n}\right)=0

一般没有普遍的解法,求解告诫微分方程的基本思想是降阶,通过变量变换的方法将方程(2.79)化为阶数较低的方程来求解,从而将问题简化.


接下来看这几类特殊的方程

不显含函数 y 的方程

F\left(x,\dfrac{\text{d}^k}{\text{d}x^k},\cdots,\dfrac{\text{d}^ny}{\text{d}x^n}\right)=0\;(k\geqslant1)\quad(2.80)

p=\dfrac{\text{d}^ky}{\text{d}x^k},则方程(2.80)变为关于 pn-k 阶微分方程

F\left(x,p,\dfrac{\text{d}p}{\text{d}x},\cdots,\dfrac{\text{d}^{n-k}y}{\text{d}x^{n-k}}\right)=0\quad(2.81)

如果能够求出方程(2.81)的通解 p=\varphi(x,C_1,\cdots,C_{n-k}),则方程

\dfrac{\text{d}^ky}{\text{d}x^k}=\varphi(x,C_1,\cdots,C_{n-k})\quad(2.82)

经过 k 次积分后得到的通解

y=\psi(x,C_1,C_2,\cdots,C_n)

就是方程(2.80)的通解,这里的 C_j,\;j=1,2,\cdots,n 为任意常数.

解方程

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

Sol:

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

x\dfrac{\text{d}p}{\text{d}x}-2(p^2-p)=0\quad(2.83)

p\not=0p\not=1 时用分离变量法,可得

p=\dfrac{1}{1-ax^2}

其中 a 为任意非零常数.

a>0 时,令 a=C_1^2,则

\dfrac{\text{d}^2y}{\text{d}x^2}=\dfrac{1}{1-C_1^2x^2}

经过两次积分后解得

y=\dfrac{1}{2C_1^2}\ln|1+C_1x||1-C_1x|+\dfrac{x}{2C_1}\ln\left|\dfrac{1+C_1x}{1-C_1x}\right|+C_2x+C_3

a<0 时,令 a=-C_1^2,则

\dfrac{\text{d}^2y}{\text{d}x^2}=\dfrac{1}{1+C_1^2x^2}

经过两次积分后解得

y=\dfrac{x}{C_1}\arctan(C_1x)-\dfrac{1}{2C_1^2}\ln(1+C_1^2x^2)+C_2x+C_3

其中 C_1\not=0,\,C_2,\,C_3 为任意常数.

此外,由于常函数 p=0p=1 也是方程(2.83)的解,因此,函数 y=C_1x+C_2,\;y=\dfrac12x^2+C_1x+C_2 也是原方程的解,其中 C_1,\,C_2 为任意常数.


不显含自变量 x 的方程

F\left(y,\dfrac{\text{d}y}{\text{d}x},\cdots,\dfrac{\text{d}^ny}{\text{d}x^n}\right)=0\quad(2.84)

这种方程也被称为自治微分方程. 令 p=\dfrac{\text{d}y}{\text{d}x},则方程(2.84)可变为关于 pn-1 阶微分方程.

事实上,若 p=\dfrac{\text{d}y}{\text{d}x},则

\dfrac{\text{d}^2y}{\text{d}x^2}=\dfrac{\text{d}}{\text{d}x}\left(\dfrac{\text{d}y}{\text{d}x}\right)=\dfrac{\text{d}p}{\text{d}y}\cdot\dfrac{\text{d}y}{\text{d}x}=p\dfrac{\text{d}p}{\text{d}y}

用数学归纳法容易证明:对任意的 i<j\leqslant n,\dfrac{\text{d}^jy}{\text{d}x^j} 可以用

y,p,\dfrac{\text{d}p}{\text{d}y},\cdots,\dfrac{\text{d}^{j-1}p}{\text{d}y^{j-1}}

来表出. 把它们带入方程(2.84)就得到形如

G\left(y,p,\dfrac{\text{d}p}{\text{d}y},\cdots,\dfrac{\text{d}^{n-1}p}{\text{d}y^{n-1}}\right)=0\quad(2.85)

的关于 pn-1 阶微分方程,比方程(2.84)低了一阶.

解方程

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

Sol:

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

\dfrac{\text{d}^2y}{\text{d}x^2}=p\dfrac{\text{d}p}{\text{d}y}.

故原方程可化为

yp\dfrac{\text{d}p}{\text{d}y}-p^2-2p=0\quad(2.86)

易知 p=0p=-2 是方程(2.86)的解. 因此,y=C,\;y=-2x+C 是原方程的解,其中 C 为任意常数.

p\not=0p\not=-2 时,使用分离变量的方法,可得方程(2.86)的通解为

p=C_1y-2,

其中 C_1 为任意常数. 求解方程

\dfrac{\text{d}y}{\text{d}x}=C_1y-2

C_1\not=0 时,原方程的通解为

y=C_2e^{C_1x}+\dfrac{2}{C_1}

其中 C_2 为任意常数. 当 C_1=0 时,积为上面讨论过的 p=-2 的情况.

齐次方程

F\left(x,y,\dfrac{\text{d}y}{\text{d}x},\cdots,\dfrac{\text{d}^ny}{\text{d}x^n}\right)=0\quad(2.87)

其中左边是关于变量

y,\dfrac{\text{d}y}{\text{d}x},\cdots,\dfrac{\text{d}^ny}{\text{d}x^n}

的零元齐次函数,即

F\left(x,ty,t\dfrac{\text{d}y}{\text{d}x},\cdots,t\dfrac{\text{d}^ny}{\text{d}x^n}\right)=F\left(x,y,\dfrac{\text{d}y}{\text{d}x},\cdots,\dfrac{\text{d}^ny}{\text{d}x^n}\right),\;\forall t\not=0

显然,当 y\not=0,齐次方程(2.87)等价于

F\left(x,1,\dfrac{1}{y}\dfrac{\text{d}y}{\text{d}x},\cdots,\dfrac{1}{y}\dfrac{\text{d}^ny}{\text{d}x^n}\right)=0\quad(2.88)

若令

p=\dfrac{1}{y}\dfrac{\text{d}y}{\text{d}x}

并以它为新未知函数,则方程就可降低一阶. 事实上,在此所设的假定下,有

\dfrac{\text{d}y}{\text{d}x}=yp

\dfrac{\text{d}^2y}{\text{d}x^2}=\dfrac{\text{d}y}{\text{d}x}p+y\dfrac{\text{d}p}{\text{d}x}=y(p^2+\dfrac{\text{d}p}{\text{d}x})

用数学归纳法不难证明:对任意的 1<k\leqslant n,\;\dfrac{1}{y}\dfrac{\text{d}^ky}{\text{d}x^k} 可用

p,\dfrac{\text{d}p}{\text{d}x},\cdots,\dfrac{\text{d}^{k-1}p}{\text{d}x^{k-1}}

表出. 将这些表达式带入方程(2.88),可得形如

G\left(x,p,\dfrac{\text{d}p}{\text{d}x},\cdots,\dfrac{\text{d}^{n-1}p}{\text{d}x^{n-1}}\right)=0\quad(2.89)

的关于 pn-1 阶微分方程,比方程(2.87)低了一阶.

解方程

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

Sol:

p=\dfrac{1}{y}\dfrac{\text{d}y}{\text{d}x}

\dfrac{\text{d}y}{\text{d}x}=yp

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

故原方程可化为

x^2y^2\dfrac{\text{d}p}{\text{d}x}=y^2-2xy^2p\quad(2.90)

y\not=0 时,方程(2.90)等价于

x^2\dfrac{\text{d}p}{\text{d}x}=1-2xp

解得

p=\dfrac{1}{x}+\dfrac{C_1}{x^2}

其中 C_1 为任意常数. 因此原方程的通解为

y=C_2e^{\int p\text{d}x}=C_2xe^{-\frac{C_1}{x}}

其中 C_2 为任意常数.

此外,y=0 显然时原方程的一个特解,已经包含在上面的通解表达式之中(取 C_2=0 即可).

全微分方程

F\left(x,y,\dfrac{\text{d}y}{\text{d}x},\cdots,\dfrac{\text{d}^ny}{\text{d}x^n}\right)=0\quad(2.91)

其中左边是某个形如

\varPhi\left(x,y,\dfrac{\text{d}y}{\text{d}x},\cdots,\dfrac{\text{d}^{n-1}y}{\text{d}x^{n-1}}\right)

的表达式对 x 的全导数,即

F\left(x,y,\dfrac{\text{d}y}{\text{d}x},\cdots,\dfrac{\text{d}^ny}{\text{d}x^n}\right)=\dfrac{\text{d}}{\text{d}x}\varPhi\left(x,y,\dfrac{\text{d}y}{\text{d}x},\cdots,\dfrac{\text{d}^{n-1}y}{\text{d}x^{n-1}}\right)

这里 n+1 元函数 \varPhi\left(x_1,x_2,\cdots,x_{n+1}\right) 的对各变元的一阶偏导数都存在且连续,故方程(2.91)有形式:

F\left(x,y,\dfrac{\text{d}y}{\text{d}x},\cdots,\dfrac{\text{d}^ny}{\text{d}x^n}\right)=\dfrac{\partial \varPhi}{\partial x_1}+\dfrac{\partial \varPhi}{\partial x_2}\dfrac{\text{d}y}{\text{d}x}+\cdots+\dfrac{\partial \varPhi}{\partial x_{n+1}}\dfrac{\text{d}^ny}{\text{d}x^n}

其中函数 \varPhi\left(x_1,x_2,\cdots,x_{n+1}\right) 的对各变元的一阶偏导数在

(x_1,x_2,\cdots,x_{n+1})=\left(x,y,\dfrac{\text{d}y}{\text{d}x},\cdots,\dfrac{\text{d}^{n-1}y}{\text{d}x^{n-1}}\right)

处取值.

此时,方程(2.91)等价于

\varPhi\left(x,y,\dfrac{\text{d}y}{\text{d}x},\cdots,\dfrac{\text{d}^{n-1}y}{\text{d}x^{n-1}}\right)=C\quad(2.92)

其中 C 为任意常数. 而方程(2.92)是 n-1 阶的,这样就降低了方程(2.91)的阶数. 与一阶微分方程类似,有时方程(2.91)本身不是全微分方程,但有时乘以一个积分因子

\mu\left(x,y,\dfrac{\text{d}y}{\text{d}x},\cdots,\dfrac{\text{d}^{n-1}y}{\text{d}x^{n-1}}\right)

后,就变成了全微分方程.

解方程

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

Sol:

y'\not=0 时,由于

\dfrac{y''}{y'}-\dfrac{2yy'}{1+y^2}=(\ln|y'|-\ln(1+y^2))'

故原方程可化为

\ln|y'|-\ln(1+y^2)=C_0

其中 C_0 为任意常数,这等价于

y'=\tilde{C}(1+y^2)\quad(\tilde{C}=\pm e^{C_0}).

由此解得 y=\tan(\tilde{C}x+C_2),其中 C_2 为任意常数.

此外,当 y'=0 时,y=C 也是原方程的解,其中 C 为任意常数.

因此原方程的通解可统一表示为

y=\tan(C_1x+C_2)

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

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

推荐阅读更多精彩内容

  • 隐式方程 因为会遇到一些导数未接触的一阶微分方程. 这里讨论一阶隐式方程,其一般形式为 求解这类方程的基本思想是将...
    洛玖言阅读 6,397评论 0 1
  • 考试形式和试卷结构一、试卷满分及考试时间 试卷满分为150分,考试时间为180分钟 二、答题方式 答题方式为闭卷、...
    幻无名阅读 751评论 0 3
  • 本章涉及知识点1、微分方程的定义2、一阶线性微分方程的定义3、求齐次线性方程通解的算法4、求非齐次线性方程通解的算...
    PrivateEye_zzy阅读 39,957评论 1 16
  • 考试科目:高等数学、线性代数、概率论与数理统计 考试形式和试卷结构 一、试卷满分及考试时间 试卷满分为150分,考...
    Saudade_lh阅读 1,076评论 0 0
  • 2018.8.24.星期五,天气晴,(255)二年级一班刘紫涵 今天彻底结束了所有的辅导班了,早上睡了个懒觉...
    凤舞涵歆阅读 129评论 0 1