欢迎光临散文网 会员登陆 & 注册

常微分方程笔记(三)

2023-03-14 23:41 作者:啊啊啊每当想起你  | 我要投稿

前言:本节继续介绍一阶ODE的解法.本节介绍了恰当微分方程及积分因子法,以及一阶隐式微分方程常见形式的解法.


2.3恰当微分方程及积分因子

2.3.1恰当微分方程

    定义:形如:

    M(x%2Cy)dx%2BN(x%2Cy)dy%3D0                                                                          (2.25)

其中M(x%2Cy)%2CN(x%2Cy)在考虑矩形域内是x,y的连续函数,且具有一阶偏导数.

    如果方程(2.25)的左边刚好是某个二元函数u(x%2Cy)的全微分,即:

    M(x%2Cy)dx%2BN(x%2Cy)dy%3Ddu(x%2Cy)%3D%5Cfrac%7B%5Cpartial%7Bu%7D%7D%7B%5Cpartial%7Bx%7D%7Ddx%2B%5Cfrac%7B%5Cpartial%7Bu%7D%7D%7B%5Cpartial%7By%7D%7Ddy                         (2.26)

那么就称其为恰当微分方程(全微分方程).易知其通解为:

    u(x%2Cy)%3Dc                                                                                                            (2.27)

    解法:(1)用恰当方程的判定求解

    定理1:二元函数M(x%2Cy)%2CN(x%2Cy)在某单连通域内是x,y的连续函数,且具有连续的一阶偏导数,则方程M(x%2Cy)dx%2BN(x%2Cy)dy%3D0 为恰当微分方程的充要条件是  

    %5Cfrac%7B%5Cpartial%7BM%7D%7D%7B%5Cpartial%7By%7D%7D%3D%5Cfrac%7B%5Cpartial%7BN%7D%7D%7B%5Cpartial%7Bx%7D%7D

    该解法即定理1的充分性证明【必要性可利用(2.26)式证明】:

    把(2.25)式两边对x积分得

    %5Cint%7BM(x%2Cy)dx%7D%2B%CF%86(y)%3Du(x%2Cy)                                                                     (2.28)

这里%CF%86(y)式关于y的任意可导函数,与x无关.两边对y求导得:

    %5Cfrac%7B%5Cpartial%7Bu%7D%7D%7B%5Cpartial%7By%7D%7D%0A%3D%5Cfrac%7B%5Cpartial%7D%7B%5Cpartial%7By%7D%7D%5Cint%7BM(x%2Cy)dx%7D%2B%5Cfrac%7Bd%CF%86(y)%7D%7Bdy%7D%3DN(x%2Cy)

因此

    %5Cfrac%7Bd%CF%86(y)%7D%7Bdy%7D%3DN(x%2Cy)-%0A%5Cfrac%7B%5Cpartial%7D%7B%5Cpartial%7By%7D%7D%5Cint%7BM(x%2Cy)dx%7D                                                        (2.29)

再将两边对x求偏导,显然右端恒等于0.事实上有:

    %5Cfrac%7B%5Cpartial%7D%7B%5Cpartial%7Bx%7D%7D%0A%5Cbegin%7Bbmatrix%7D%0AN(x%2Cy)-%5Cfrac%7B%5Cpartial%7D%7B%5Cpartial%7By%7D%7D%5Cint%7BM(x%2Cy)dx%7D%0A%5Cend%7Bbmatrix%7D%0A%3D%5Cfrac%7B%5Cpartial%7BN%7D%7D%7B%5Cpartial%7Bx%7D%7D-%5Cfrac%7B%5Cpartial%7D%7B%5Cpartial%7Bx%7D%7D%0A%5Cbegin%7Bbmatrix%7D%0A%5Cfrac%7B%5Cpartial%7D%7B%5Cpartial%7By%7D%7D%5Cint%7BM(x%2Cy)dx%7D%0A%5Cend%7Bbmatrix%7D(2.30)

又由于二元函数M(x%2Cy)%2CN(x%2Cy)在某单连通域内是x,y的连续函数,且具有连续的一阶偏导数,因此(2.30)中对x,y的求导顺序可以交换.再利用%5Cfrac%7B%5Cpartial%7BM%7D%7D%7B%5Cpartial%7By%7D%7D%3D%5Cfrac%7B%5Cpartial%7BN%7D%7D%7B%5Cpartial%7Bx%7D%7D,可得:

%5Cfrac%7B%5Cpartial%7BN%7D%7D%7B%5Cpartial%7Bx%7D%7D-%5Cfrac%7B%5Cpartial%7D%7B%5Cpartial%7Bx%7D%7D%0A%5Cbegin%7Bbmatrix%7D%0A%5Cfrac%7B%5Cpartial%7D%7B%5Cpartial%7By%7D%7D%5Cint%7BM(x%2Cy)dx%7D%0A%5Cend%7Bbmatrix%7D%0A%3D%5Cfrac%7B%5Cpartial%7BN%7D%7D%7B%5Cpartial%7Bx%7D%7D-%5Cfrac%7B%5Cpartial%7D%7B%5Cpartial%7By%7D%7D%0A%5Cbegin%7Bbmatrix%7D%0A%5Cfrac%7B%5Cpartial%7D%7B%5Cpartial%7Bx%7D%7D%5Cint%7BM(x%2Cy)dx%7D%0A%5Cend%7Bbmatrix%7D%0A%3D%5Cfrac%7B%5Cpartial%7BN%7D%7D%7B%5Cpartial%7Bx%7D%7D-%5Cfrac%7B%5Cpartial%7BM%7D%7D%7B%5Cpartial%7By%7D%7D%3D0

    于是(2.29)式右端为y的一元函数,两边积分可得:

    %CF%86(y)%3D%5Cint%0A%5Cbegin%7Bbmatrix%7D%0AN(x%2Cy)-%5Cfrac%7B%5Cpartial%7D%7B%5Cpartial%7By%7D%7D%5Cint%7BM(x%2Cy)dx%7D%0A%5Cend%7Bbmatrix%7D%0Ady                                              (2.31)

将(2.31)代入至(2.28)式可得公式:

    u(x%2Cy)%3D%5Cint%7BM(x%2Cy)dx%7D%2B%5Cint%0A%5Cbegin%7Bbmatrix%7D%0AN(x%2Cy)-%5Cfrac%7B%5Cpartial%7D%7B%5Cpartial%7By%7D%7D%5Cint%7BM(x%2Cy)dx%7D%0A%5Cend%7Bbmatrix%7D%0Ady          (2.32)

(2)利用曲线积分求u(x%2Cy)

    命题1:设二元函数M(x%2Cy)%2CN(x%2Cy)在某单连通区域D内是x,y的连续函数,且具有连续的一阶偏导数,则对D内任一按段光滑曲线L,曲线积分在区域D内的积分%5Cint_L%7BM(x%2Cy)dx%2BN(x%2Cy)dy%7D与路径无关的充要条件是%5Cfrac%7B%5Cpartial%7BM%7D%7D%7B%5Cpartial%7By%7D%7D%3D%5Cfrac%7B%5Cpartial%7BN%7D%7D%7B%5Cpartial%7Bx%7D%7D.

    利用命题1从定点A(x_0%2Cy_0)到动点B(x%2Cy)选择简单路径进行积分:

    u(x%2Cy)%3D%5Cint%5E%7Bx%7D_%7Bx_0%7D%7BM(x%2Cy_0)dx%7D%2B%5Cint%5E%7By%7D_%7By_0%7D%7BN(x%2Cy)dy%7D                                             (2.33)

    u(x%2Cy)%3D%5Cint%5E%7By%7D_%7By_0%7D%7BN(x_0%2Cy)dy%7D%2B%5Cint%5E%7Bx%7D_%7Bx_0%7D%7BM(x%2Cy)dx%7D                                             (2.34)

(3)分项组合方法(凑“全微分”)

    需要熟记简单的二元函数全微分:如:

    ydx%2Bxdy%3Dd(xy)     %5Cfrac%7Bydx-xdy%7D%7By%5E2%7D%3Dd(%5Cfrac%7Bx%7D%7By%7D)     %5Cfrac%7Bxdy-ydx%7D%7Bx%5E2%7D%3Dd(%5Cfrac%7By%7D%7Bx%7D)

   %5Cfrac%7Bydx-xdy%7D%7Bxy%7D%3Dd%0A%5Cleft(%5C%0A%5Cbegin%7Bmatrix%7D%0Aln%0A%5Cbegin%7Bvmatrix%7D%0A%5Cfrac%7Bx%7D%7By%7D%0A%5Cend%7Bvmatrix%7D%0A%5Cend%7Bmatrix%7D%0A%5Cright)%5C%5Cfrac%7Bydx-xdy%7D%7Bx%5E2%2By%5E2%7D%3Dd(arccot%5Cfrac%7By%7D%7Bx%7D)%0A%3Dd(arctan%5Cfrac%7By%7D%7Bx%7D)

   %5Cfrac%7Bydx-xdy%7D%7Bx%5E2-y%5E2%7D%3D%5Cfrac%7B1%7D%7B2%7Dd(ln%0A%5Cbegin%7Bvmatrix%7D%0A%5Cfrac%7Bx-y%7D%7Bx%2By%7D%0A%5Cend%7Bvmatrix%7D%0A)

2.3.2积分因子

    定义:如果不是一个恰当微分方程,那么可以通过将M(x%2Cy)%2CN(x%2Cy)都乘以一个连续可微函数%CE%BC(x%2Cy)%5Cneq0,使得

    %CE%BC(x%2Cy)M(x%2Cy)dx%2B%CE%BC(x%2Cy)N(x%2Cy)dy%3D0                                               (2.35)

为一恰当微分方程,即存在光滑的函数u(x%2Cy),使得

    %CE%BC(x%2Cy)M(x%2Cy)dx%2B%CE%BC(x%2Cy)N(x%2Cy)dy%3Ddu(x%2Cy)                                  (2.36)

则称%CE%BC(x%2Cy)为方程(2.25)的积分因子.

    此时有:%5Cfrac%7B%5Cpartial%7B(%CE%BCP)%7D%7D%20%7B%5Cpartial%7By%7D%7D%3D%0A%5Cfrac%7B%5Cpartial%7B(%CE%BCQ)%7D%7D%20%7B%5Cpartial%7Bx%7D%7D%0A%5Cbigg(%20%3D%20%5Cfrac%7B%5Cpartial%5E2%7B%CE%A6%7D%7D%20%7B%5Cpartial%7Bx%7D%5Cpartial%7By%7D%7D%20%5Cbigg)                                                      (2.37)

    由(2.37)式得:寻求可微函数%CE%BC(x%2Cy)等价于求解一阶偏微分方程

    P(x%2Cy)%5Cfrac%7B%5Cpartial%7B%5Cmu%7D%7D%20%7B%5Cpartial%7By%7D%7D-%0AQ(x%2Cy)%5Cfrac%7B%5Cpartial%7B%5Cmu%7D%7D%20%7B%5Cpartial%7Bx%7D%7D%3D%0A%5Cbigg(%20%5Cfrac%7B%5Cpartial%7BQ%7D%7D%20%7B%5Cpartial%7Bx%7D%7D-%0A%5Cfrac%7B%5Cpartial%7BP%7D%7D%20%7B%5Cpartial%7By%7D%7D%20%5Cbigg)%5Cmu(x%2Cy)                                (2.38)

    另外,同一个方程可以有不同的积分因子.可以证明,只要方程有解存在,则必有积分因子存在,并且不唯一,从而使通解有不同的形式.

    这里一种特殊形式的积分因子:

    (1)偏微分方程(2.38)有一个只依赖x的解%CE%BC(x)的充要条件是,由下式定义的函数G只依赖于x:G%3D%20-%20%5Cfrac%7B%20%5Cfrac%7B%5Cpartial%7BQ%7D%7D%20%7B%5Cpartial%7Bx%7D%7D-%0A%5Cfrac%7B%5Cpartial%7BP%7D%7D%20%7B%5Cpartial%7By%7D%7D%20%7D%7BQ(x%2Cy)%7D                                                                                    (2.39)

此时有:%CE%BC(x)%3De%5E%7B%5Cint%5E%7Bx%7D_%7Bx0%7DG(t)dt%7D                                                                                      (2.40)

    (2)偏微分方程(2.38)有一个只依赖y的解%CE%BC(y)的充要条件是,由下式定义的函数H只依赖于x:H%3D%20%5Cfrac%7B%20%5Cfrac%7B%5Cpartial%7BQ%7D%7D%20%7B%5Cpartial%7Bx%7D%7D-%0A%5Cfrac%7B%5Cpartial%7BP%7D%7D%20%7B%5Cpartial%7By%7D%7D%20%7D%7BP(x%2Cy)%7D                                                                                       (2.41)

此时有:%CE%BC(y)%3De%5E%7B%5Cint%5E%7By%7D_%7By0%7DH(s)ds%7D                                                                                     (2.42)

    (3)方程(2.25)有形如%CE%BC%3D%CE%BC(%5Cphi%20(x%2Cy))的积分因子的充分必要条件是:

    %20%5Cfrac%7B%20%5Cfrac%7B%5Cpartial%7BP%7D%7D%20%7B%5Cpartial%7By%7D%7D-%0A%5Cfrac%7B%5Cpartial%7BQ%7D%7D%20%7B%5Cpartial%7Bx%7D%7D%20%7D%0A%7BP%5Cfrac%7B%5Cpartial%7B%5Cphi%20%7D%7D%20%7B%5Cpartial%7By%7D%7D-%0AQ%5Cfrac%7B%5Cpartial%7B%5Cphi%20%7D%7D%20%7B%5Cpartial%7Bx%7D%7D%7D%0A%3Df(%5Cphi%20(x%2Cy))                                                                             (2.43)

其中f为某个一元函数.


(P.S 求解恰当微分方程时,如果使用公式只需要把原函数积分出来即可,比如M(x%2Cy)对x积分时,不需要带上仅含y的函数%CF%86(y)


2.4一阶隐式微分方程与参数表示

    一阶隐式微分方程形式为F(x%2Cy%2Cy')%3D0

    本节主要讨论下面四种类型:

    y%3Df(x%2Cy')x%3Df(y%2Cy')F(x%2Cy')%3D0F(y%2Cy')%3D0

2.4.1可以解出y(或x)的方程

(1)首先讨论形如

    y%3Df(x%2Cy')                                                                                                          (2.37)

的方程,这里假设f(x%2Cy')具有连续偏导数.

    引入参数p,令y'%3Dp.则方程(2.37)变为

     y%3Df(x%2Cp)                                                                                                          (2.38)

    将两边对x求导,并将y'%3Dp代入,得

    p%3D%5Cfrac%7B%5Cpartial%7Bf%7D%7D%7B%5Cpartial%7Bx%7D%7D%2B%5Cfrac%7B%5Cpartial%7Bf%7D%7D%7B%5Cpartial%7Bp%7D%7D%5Cfrac%7Bdf%7D%7Bdx%7D                                                                                              (2.39)

    (2.39)式即可按照2.1-2.3节解之.

    若(2.39)式的通解为p%3D%CF%86(x%2Cc),代入(2.38)式得y%3Df(x%2C%CF%86(x%2Cc))即为(2.37)式的通解.

    若(2.39)式的通解为x%3D%CF%88(p%2Cc),则(2.37)式的通解为参数方程形式:

    %5Cbegin%7Bcases%7D%0Ax%3D%CF%88(p%2Cc)%20%5C%5C%0Ay%3Df(%CF%88(p%2Cc)%2Cp)%0A%5Cend%7Bcases%7D(其中p为参数,c为任意常数).

    若(2.39)式的通解为%CE%A6(x%2Cp%2Cc)%3D0,则(2.37)式的通解为参数方程形式:

    %5Cbegin%7Bcases%7D%0A%CE%A6(x%2Cp%2Cc)%3D0%20%5C%5C%0Ay%3Df(x%2Cp)%0A%5Cend%7Bcases%7D(其中p为参数,c为任意常数).

(2)讨论形如

    x%3Df(y%2Cy')                                                                                                          (2.40)

的方程,这里假设f(y%2Cy')具有连续偏导数.

    引入参数p,令y'%3Dp,将(2.40)式两边对y求导,并将y'%3Dp代入,得

    %5Cfrac%7B1%7D%7Bp%7D%3D%0A%5Cfrac%7B%5Cpartial%7Bf%7D%7D%7B%5Cpartial%7By%7D%7D%2B%5Cfrac%7B%5Cpartial%7Bf%7D%7D%7B%5Cpartial%7Bp%7D%7D%5Cfrac%7Bdp%7D%7Bdy%7D                                                                                             (2.41)

    (2.41)式是关于y,p的一阶微分方程,可按照2.1-2.3节解之.

    若(2.40)式的通解为%CE%A6(y%2Cp%2Cc)%3D0,则(2.37)式的通解为参数方程形式:

    %5Cbegin%7Bcases%7D%0A%CE%A6(y%2Cp%2Cc)%3D0%20%5C%5C%0Ax%3Df(y%2Cp)%0A%5Cend%7Bcases%7D(其中p为参数,c为任意常数).

2.4.2不显含y(或x)的方程

(3)讨论形如

    F(x%2Cy')%3D0                                                                                                         (2.42)

的方程.令y'%3Dp,方程化为F(x%2Cp)%3D0,代表Oxp平面上的一条曲线.如果把曲线表示成参数方程形式:

    %5Cbegin%7Bcases%7D%0Ax%3D%CF%86(t)%20%5C%5C%0Ap%3D%CF%88(t)%0A%5Cend%7Bcases%7D                                                                                                           (2.43)

t为参数.代入y'%3Dp,得dy%3D%CF%88(t)%CF%86'(t)dt.两边积分,得y%3D%5Cint%7B%CF%88(t)%CF%86'(t)dt%7D%2Bc.代入(2.42)即可得原方程参数形式的通解%5Cbegin%7Bcases%7D%0Ax%3D%CF%86(t)%20%5C%5C%0Ay%3D%5Cint%7B%CF%88(t)%CF%86'(t)dt%7D%2Bc%0A%5Cend%7Bcases%7D,c为任意常数

(4)讨论形如

    F(y%2Cy')%3D0                                                                                                         (2.42)

的方程.其解法与(2.42)的解法相似.

    令y'%3Dp,方程化为F(y%2Cp)%3D0,代表Oyp平面上的一条曲线.如果把曲线表示成参数方程形式:

    %5Cbegin%7Bcases%7D%0Ay%3D%CF%86(t)%20%5C%5C%0Ap%3D%CF%88(t)%0A%5Cend%7Bcases%7D                                                                                                           (2.43)

t为参数.代入dy%3Dpdx,得dx%3D%5Cfrac%7B%CF%86'(t)%7D%7B%CF%88(t)%7Ddt.两边积分,得x%3D%5Cint%7B%5Cfrac%7B%CF%86'(t)%7D%7B%CF%88(t)%7Ddt%7D%2Bc.代入(2.42)即可得原方程参数形式的通解%5Cbegin%7Bcases%7D%0Ax%3D%5Cint%7B%5Cfrac%7B%CF%86'(t)%7D%7B%CF%88(t)%7Ddt%7D%2Bc%20%5C%5C%0Ay%3D%CF%86(t)%0A%5Cend%7Bcases%7D,c为任意常数.

    当y'%3Dp%3D0时,F(y%2C0)%3D0有实根y%3Dk,即y%3Dk也是方程的解.



常微分方程笔记(三)的评论 (共 条)

分享到微博请遵守国家法律