常微分方程求解(2)

一、全微分方程定义与性质

我们将一阶方程改写为对称的形式:

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

如果上式的左边恰好是某一个二元函数 u(x,y) 的全微分,即:

M(x,y)dx+N(x,y)dy=du(x,y)

则称之为全微分方程,u(x,y) 称为它的一个原函数。

y=y(x) 是全微分方程的一个解,且 y(x0)=y0。则:

M(x,y(x))dx+N(x,y(x))dy(x)0

即:

du(x,y(x))0

于是:

u(x,y(x))c0=u(x0,y0)

y=y(x)u(x,y)=u(x0,y0) 满足 y(x0)=y0 所确定的隐函数。反之,设y=y(x)u(x,y)=u(x0,y0)满足y(x0)=y0所确定的隐函数,则有:

u(x,y(x))u(x0,y0)

求微分,得:

du(x,y(x))0

于是:

M(x,y(x))dx+N(x,y(x))dy(x)0

y=y(x)是全微分方程的的解。

由以上论证可知,如果已知方程是一个全微分方程,u(x,y)是它的一个原函数,则:

u(x,y)=u(x0,y0)

是方程满足条件:当x=x0y=y0的积分,其中(x0,y0)G.并且

u(x,y)=C

是方程的通积分,C是任意常数。
现在存在两个问题:

  1. 如何根据M(x,y)N(x,y)去判别方程为全微分方程。
  2. 当方程为全微分方程时,如何去求原函数u(x,y).

二、全微分方程的判别

为了判别一个方程是否为全微分方程,给出下面的定理:

设函数M(x,y)N(x,y)在单连通区域G内连续且有连续的一阶偏导数,则方程为全微分方程的充要条件是:

MyNx,(x,y)G

一个区域D中,任一闭曲线所围区域,都属于区域D,称为“单连通区域”。这就意味着,区域D是实心的,不能有空洞。

三、全微分方程求解

1、方法一

当满足充要条件,确定方程为全微分方程时,可以由下述与路径无关的曲线积分求得一个原函数:

u(x,y)=(x0,y0)(x,y)M(x,y)dx+N(x,y)dy

其中点(x0,y0) 可以取G 内任意一个定点。
如果G是边平行于坐标轴的矩形区域,则从点 (x0,y0) 到点(x,y) 可取边平行于坐标轴的折线ADDBACCB

其计算公式为:

u(x,y)=x0xM(ξ,y0)dξ+y0yN(x,η)dη

或:

u(x,y)=x0xM(ξ,y)dξ+y0yN(x0,η)dη

2、方法二

我们也可以不使用曲线积分,采用下述不定积分的方法求原函数u(x,y)。因为:

du(x,y)=M(x,y)dx+N(x,y)dy

所以:

ux=M(x,y),uy=N(x,y)

由前一式,有

u(x,y)=M(x,y)dx+φ(y)

这里积分中把y看作常量,φ(y)y的任意可微函数。再由第二式,有:

uy=yM(x,y)dx+φ(y)=N(x,y)

由此求得:

φ(y)=N(x,y)yM(x,y)dx

在方程为全微分方程的条件下,可证上式右边与x无关,从而可求出φ(y),代入即得u(x,y)

3、方法三

此方法的优点是不必事先验证该方程是否满足充要条件来判断是否为全微分方程,但难点是找原函数要采用“凑”的办法,有一点技巧,要熟练掌握求二元函数全微分的逆向思维。例如下面一些等式,从右到左是计算二元函数的全微分,而从左到右就是“凑”出原函数。

ydx+xdy=d(xy)

ydxxdyy2=d(xy)

ydx+xdyx2=d(yx)

ydx+xdyx2+y2=d(arctanyx)

ydxxdyx2y2=d(12ln|xyx+y|)

三、使用积分因子法配凑全微分方程

我们考虑一阶微分方程:

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

若该方程不是全微分方程,即:

MyNx

我们希望找到一个积分因子 μ,使得乘以它后变为全微分方程:

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

我们分别考虑 μ 仅为 xy 的函数的情况。

1、若积分因子为 μ(x),即积分因子只与 x 有关:

μ=μ(x),则方程变为:

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

为使其成为全微分方程,需满足:

y[μ(x)M(x,y)]=x[μ(x)N(x,y)]

即:

μ(x)My=μ(x)N(x,y)+μ(x)Nx

两边除以 μ(x)(假设 μ(x)0)得:

μ(x)μ(x)=MyNxN

这给出了一个微分方程:

dμdx=μ(x)MyNxN

若右边仅依赖于 x(即分子与 x 有关,N 不含 y),则此方程可解,进而可求得积分因子 μ(x)

在解决问题时,首先计算:

MyNxN

如果上式结果仅仅与 x 有关,那么就记上式为 φ(x),去寻找一个函数 μ(x)(解微分方程),满足:

dμdx1μ=φ(x)


方法二:积分因子为 μ(y)

μ=μ(y),则方程变为:

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

为使其成为全微分方程,需满足:

y[μ(y)M(x,y)]=x[μ(y)N(x,y)]

即:

μ(y)M(x,y)+μ(y)My=μ(y)Nx

两边除以 μ(y) 得:

μ(y)μ(y)=NxMyM

得到微分方程:

dμdy=μ(y)NxMyM

若右边仅依赖于 y,则可求得 μ(y)

在解决问题时,首先计算:

MyNxM

如果上式结果仅仅与 y 有关,那么就记上式为 φ(y),去寻找一个函数 μ(y)(解微分方程),满足:

dμdy1μ=φ(y)

打赏
评论区
头像
  • 高兴
  • 小怒
  • 脸红
  • 内伤
  • 装大款
  • 赞一个
  • 害羞
  • 汗
  • 吐血倒地
  • 深思
  • 不高兴
  • 无语
  • 亲亲
  • 口水
  • 尴尬
  • 中指
  • 想一想
  • 哭泣
  • 便便
  • 献花
  • 皱眉
  • 傻笑
  • 狂汗
  • 吐
  • 喷水
  • 看不见
  • 鼓掌
  • 阴暗
  • 长草
  • 献黄瓜
  • 邪恶
  • 期待
  • 得意
  • 吐舌
  • 喷血
  • 无所谓
  • 观察
  • 暗地观察
  • 肿包
  • 中枪
  • 大囧
  • 呲牙
  • 抠鼻
  • 不说话
  • 咽气
  • 欢呼
  • 锁眉
  • 蜡烛
  • 坐等
  • 击掌
  • 惊喜
  • 喜极而泣
  • 抽烟
  • 不出所料
  • 愤怒
  • 无奈
  • 黑线
  • 投降
  • 看热闹
  • 扇耳光
  • 小眼睛
  • 中刀
  • 呵呵
  • 哈哈
  • 吐舌
  • 太开心
  • 笑眼
  • 花心
  • 小乖
  • 乖
  • 捂嘴笑
  • 滑稽
  • 你懂的
  • 不高兴
  • 怒
  • 汗
  • 黑线
  • 泪
  • 真棒
  • 喷
  • 惊哭
  • 阴险
  • 鄙视
  • 酷
  • 啊
  • 狂汗
  • what
  • 疑问
  • 酸爽
  • 呀咩爹
  • 委屈
  • 惊讶
  • 睡觉
  • 笑尿
  • 挖鼻
  • 吐
  • 犀利
  • 小红脸
  • 懒得理
  • 勉强
  • 爱心
  • 心碎
  • 玫瑰
  • 礼物
  • 彩虹
  • 太阳
  • 星星月亮
  • 钱币
  • 茶杯
  • 蛋糕
  • 大拇指
  • 胜利
  • haha
  • OK
  • 沙发
  • 手纸
  • 香蕉
  • 便便
  • 药丸
  • 红领巾
  • 蜡烛
  • 音乐
  • 灯泡
  • 开心
  • 钱
  • 咦
  • 呼
  • 冷
  • 生气
  • 弱
  • 阿鲁
  • 泡泡
    头像
    路边的一头
      

    可以

      头像
      laoguanTX
        
文章目录