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

常微分方程(二)

2022-10-09 20:21 作者:啊啊啊每当想起你  | 我要投稿

前言:

本节介绍第二章:一阶常微分的初等解法.

有一说一其实笔者之前学过高数,所以对一阶常微分方程的解法中的分离变量和常数变异法这两种方法多少有点印象,本章除了介绍这两种方法外还介绍了一种积分因子法(针对全微分方程的).


2.1变量分离方程与变量变换

2.1.1变量分离方程

    变量分离方程:形如

    %5Cfrac%7Bdy%7D%7Bdx%7D%3Df(x)%CF%86(y)                                                                                                    (2.1)

称为变量分离方程,这里f(x)%EF%BC%8C%CF%86(y)分别是关于x,y的连续函数.

如果%CF%86(y)%5Cneq0,可将方程改写(即分离变量)为:

    %5Cfrac%7Bdy%7D%7B%CF%86(y)%7D%3Df(x)dx                                                                                                    

两边积分得:

    %5Cint%5Cfrac%7Bdy%7D%7B%CF%86(y)%7D%3D%5Cint%7Bf(x)dx%7D%2Bc                                                                                  (2.2)

此处常数c的取值必须要使(2.3)有意义.

【P.S. 需要注意,如果存在y_0使得%CF%86(y_0)%3D0,则直接验证y%3Dy_0也是(2.1)的解,因此还必须寻求%CF%86(y)%3D0的,当y%3Dy_0不包括在方程(2.2)中时,必须补上特解y%3Dy_0.】

(之前看到丁版ODE时也看到了这点,不过是以例题的方式出现的;而大一学同济版高数的时候没有这点,包括很多课后习题也没有提到)

2.1.2可化为变量分离方程的类型(变量代换)

2.1.2.1齐次微分方程

    齐次微分方程:形如:

    %5Cfrac%7Bdy%7D%7Bdx%7D%3Dg(%5Cfrac%7By%7D%7Bx%7D)                                                                                                            (2.3)

的方程称为齐次微分方程,这里g(u)为u的连续函数.

    作变量代换

    u%3D%5Cfrac%7By%7D%7Bx%7D                                                                                                                      (2.4)

y%3Dux,两边对x求导得:

    %5Cfrac%7Bdy%7D%7Bdx%7D%3Dx%5Cfrac%7Bdu%7D%7Bdx%7D%2Bu                                                                                                    (2.5)

(2.4)与(2.5)式代入(2.3)式整理得:

    %5Cfrac%7Bdu%7D%7Bdx%7D%3D%5Cfrac%7Bg(u)-u%7D%7Bx%7D                                                                                                    (2.6)

    此时(2.6)式为一个变量分离方程,可用(2.2)式进行求解,之后再将(2.4)式代入所得解中的u即可

2.1.2.2u%3Dax%2Bby%2Bc 型

    定义:形如:

    %5Cfrac%7Bdy%7D%7Bdx%7D%3Df(ax%2Bby%2Bc)                                                                                         (2.7)

的方程称为齐次微分方程,这里f(u)为u的连续函数.

    作变量代换

    u%3Dax%2Bby%2Bc                                                                                                   (2.8)

y%3Dux,两边对x求导得:

    %5Cfrac%7Bdu%7D%7Bdx%7D%3Da%2Bb%5Cfrac%7Bdy%7D%7Bdx%7D                                                                                                     (2.9)

(2.8)与(2.9)式代入(2.7)式整理得:

    %5Cfrac%7Bdu%7D%7Bdx%7D%3Da%2Bbf(u)                                                                                                 (2.10)

    此时(2.10)式为一个变量分离方程,可用(2.2)式进行求解,之后再将(2.8)式代入所得解中的u即可

2.1.2.3f(%5Cfrac%7Ba_1x%2Bb_1y%2Bc_1%7D%7Ba_2x%2Bb_2y%2Bc_2%7D) 型

    定义:形如:

    %5Cfrac%7Bdy%7D%7Bdx%7D%3Df(%5Cfrac%7Ba_1x%2Bb_1y%2Bc_1%7D%7Ba_2x%2Bb_2y%2Bc_2%7D)                                                                               (2.11)

的方程称为齐次微分方程,这里f(u)为u的连续函数.

    情形一:c_1%3Dc_2%3D0,则%5Cfrac%7Bdy%7D%7Bdx%7D%3Df(%5Cfrac%7Ba_1x%2Bb_1y%7D%7Ba_2x%2Bb_2y%7D)%3Df(%5Cfrac%7Ba_1%2Bb_1%5Cfrac%7By%7D%7Bx%7D%7D%7Ba_2%2Bb_2%5Cfrac%7By%7D%7Bx%7D%7D)%3Dg(%5Cfrac%7By%7D%7Bx%7D)  ,此时为齐次方程(2.3).

    情形二:c_1%2Cc_2不全为0,则讨论方程(2.11)右端分子分母的一次多项式:

    %5Cbegin%7Bequation%7D%0A%5Cleft%5C%7B%0A%5Cbegin%7Baligned%7D%0Aa_1x%2Bb_1y%2Bc_1%3D0%20%5C%5C%0Aa_2x%2Bb_2y%2Bc_2%3D0%5C%5C%0A%5Cend%7Baligned%7D%0A%5Cright.%0A%5Cend%7Bequation%7D%0A  即%5Cbegin%7Bequation%7D%0A%5Cleft%5C%7B%0A%5Cbegin%7Baligned%7D%0Aa_1x%2Bb_1y%3D-c_1%20%5C%5C%0Aa_2x%2Bb_2y%3D-c_2%5C%5C%0A%5Cend%7Baligned%7D%0A%5Cright.%0A%5Cend%7Bequation%7D%0A                                          (2.12)

对系数行列式D%3D%0A%5Cbegin%7Bvmatrix%7D%0Aa_1%20%26%20b_1%20%5C%5C%0Aa_2%20%26%20b_2%0A%5Cend%7Bvmatrix%7D%0A分两种情况进行讨论:

(此处设D_1%3D%0A%5Cbegin%7Bvmatrix%7D%0A-c_1%20%26%20b_1%20%5C%5C%0A-c_2%20%26%20b_2%0A%5Cend%7Bvmatrix%7D%0A%3D-%0A%5Cbegin%7Bvmatrix%7D%0Ac_1%20%26%20b_1%20%5C%5C%0Ac_2%20%26%20b_2%0A%5Cend%7Bvmatrix%7D%0AD_2%3D%0A%5Cbegin%7Bvmatrix%7D%0Aa_1%20%26%20-c_1%20%5C%5C%0Aa_2%20%26%20-c_2%0A%5Cend%7Bvmatrix%7D%0A%3D-%0A%5Cbegin%7Bvmatrix%7D%0Aa_1%20%26%20c_1%20%5C%5C%0Aa_2%20%26%20c_2%0A%5Cend%7Bvmatrix%7D

    (1)D%5Cneq0%0A,则由Cramer法则可解得%5Cbegin%7Bequation%7D%0A%5Cbegin%7Bcases%7D%0Ax%20%3D%20%5Calpha%20%5C%5C%0Ay%20%3D%20%5Cbeta%20%0A%5Cend%7Bcases%7D%0A%5Cend%7Bequation%7D%0A,其中%5Calpha%20%3D%20%5Cfrac%7BD_1%7D%7BD%7D%0A%5Cbeta%20%3D%20%5Cfrac%7BD_2%7D%7BD%7D%0A.

    若令 %5Cbegin%7Bequation%7D%0A%5Cleft%5C%7B%0A%5Cbegin%7Baligned%7D%0AX%3Dx-%5Calpha%20%5C%5C%0AY%3Dy-%5Cbeta%5C%5C%0A%5Cend%7Baligned%7D%0A%5Cright.%0A%5Cend%7Bequation%7D%0A,则(2.12)式可化为:%5Cbegin%7Bequation%7D%0A%5Cleft%5C%7B%0A%5Cbegin%7Baligned%7D%0Aa_1X%2Bb_1Y%3D0%20%5C%5C%0Aa_1X%2Bb_1Y%3D0%20%5C%5C%0A%5Cend%7Baligned%7D%0A%5Cright.%0A%5Cend%7Bequation%7D%0A

即(2.11)式可化为:%5Cfrac%7BdY%7D%7BdX%7D%3D%0A%5Cfrac%7Bdy%7D%7Bdx%7D%3D%0Af(%5Cfrac%7Ba_1%2Bb_1%5Cfrac%7BY%7D%7BX%7D%7D%7Ba_2%2Bb_2%5Cfrac%7BY%7D%7BX%7D%7D)%3Dg(%5Cfrac%7BY%7D%7BX%7D)                     (2.13)

(2.13)式为齐次微分方程,求解代入原变量即可得方程(2.11)的解

    (2)D%3D0%0A,分以下三种情况

    (i)a_1%3Db_1%3D0%0A时,(2.11)式变为%5Cfrac%7Bdy%7D%7Bdx%7D%3Df(%5Cfrac%7Bc_1%7D%7Ba_2x%2Bb_2y%2Bc_2%7D) ,而a_2%3Db_2%3D0%0A时,(2.11)式变为%5Cfrac%7Bdy%7D%7Bdx%7D%3Df(%5Cfrac%7Ba_1x%2Bb_1y%2Bc_1%7D%7Bc_2%7D),二者均用(2.7)计算.

    (ii)a_1%3Da_2%3D0%0A时,(2.11)式变为%5Cfrac%7Bdy%7D%7Bdx%7D%3Df(%5Cfrac%7Bb_1y%2Bc_1%7D%7Bb_2y%2Bc_2%7D) ,而b_1%3Db_2%3D0%0A时,(2.11)式变为%5Cfrac%7Bdy%7D%7Bdx%7D%3Df(%5Cfrac%7Ba_1x%2Bc_1%7D%7Ba_2x%2Bc_2%7D),二者均为变量可分离方程(2.1)

    (iii)%5Cfrac%7Ba_1%7D%7Ba_2%7D%3D%5Cfrac%7Bb_1%7D%7Bb_2%7D%3Dk%5Cneq0%0A时,令u%3Da_2x%2Bb_2y ,此时%5Cfrac%7Bdu%7D%7Bdx%7D%3Da_2%2Bb_2%5Cfrac%7Bdy%7D%7Bdx%7D,(2.11)式变为f(%5Cfrac%7Bk(a_2x%2Bb_2y)%2Bc_1%7D%7Ba_2x%2Bb_2y%2Bc_2%7D)%3Df(%5Cfrac%7Bku%2Bc_1%7D%7Bu%2Bc_2%7D)%3Dg(u) .方程(2.17)化为%5Cfrac%7Bdu%7D%7Bdx%7D%3Da_2%2Bb_2g(u),属于变量分离方程(2.1)


2.2 线性微分方程与常数变易法

2.2.1 一阶线性微分方程

    定义:形如:

    %5Cfrac%7Bdy%7D%7Bdx%7D%3DP(x)y%2BQ(x)                                                                                        (2.14)的方程称为一阶线性常微分方程,其中P(x)%2CQ(x)在考虑的区间是x的连续函数.

    若Q(x)%3D0,则(2.14)式变为:

    %5Cfrac%7Bdy%7D%7Bdx%7D%3DP(x)y                                                                                                        (2.15)

(2.15)称为一阶齐次线性微分方程.若Q(x)%5Cneq0,则(2.14)称为一阶非齐次线性微分方程方程.

2.2.2常数变易法

    定义:将齐次微分方程中的常数c变易成待定函数c(x).

    解法:对于(2.15)式可用变量分离法解得其通解为y%3Dce%5E%7B%5Cint%7BP(x)dx%7D%7D                    (2.16)

这里c是任意常数.现在讨论(2.14)式的通解,利用常数变易法,将任意常数c换成待定函数c(x),得到(2.14)的形式求解,即:y%3Dc(x)e%5E%7B%5Cint%7BP(x)dx%7D%7D                                      (2.17)

    对(2.17)式求导得:

    %5Cfrac%7Bdy%7D%7Bdx%7D%3D%5Cfrac%7Bdc(x)%7D%7Bdx%7De%5E%7B%5Cint%7BP(x)dx%7D%7D%2Bc(x)P(x)e%5E%7B%5Cint%7BP(x)dx%7D%7D%0A%3DQ(x)%2Bc(x)P(x)e%5E%7B%5Cint%7BP(x)dx%7D%7D 

    %5Cfrac%7Bdc(x)%7D%7Bdx%7D%0A%3DQ(x)e%5E%7B-%5Cint%7BP(x)dx%7D%7D 

积分后得:

    c(x)%3D%5Cint%7BQ(x)e%5E%7B-%5Cint%7BP(x)dx%7D%7D%7D%2Bc'                                                                       (2.18)

(2.18)式代入(2.17)得通解为: 

    y%3De%5E%7B%5Cint%7BP(x)dx%7D%7D(%5Cint%7BQ(x)e%5E%7B-%5Cint%7BP(x)dx%7D%7D%7D%2Bc')                                                         (2.19)

【P.S. 在同济版高数大家见过的都是如下形式:

    %5Cfrac%7Bdy%7D%7Bdx%7D%2BP(x)y%3DQ(x) 

它的通解为:y%3De%5E%7B-%5Cint%7BP(x)dx%7D%7D(%5Cint%7BQ(x)e%5E%7B%5Cint%7BP(x)dx%7D%7D%7D%2Bc')

其实与(2.19)式一样的,将等式左边的P(x)移至右侧,得到%5Cfrac%7Bdy%7D%7Bdx%7D%3D-P(x)y%2BQ(x) ,便可以代入至(2.19)式进行求解,得到上面的通解,只是要注意这个差别】

2.2.3 Bernoulli微分方程

    定义:形如:

    %5Cfrac%7Bdy%7D%7Bdx%7D%3DP(x)y%2BQ(x)y%5En                                                                                   (2.20)

    解法:将(2.20)式两边同乘以y%5E%7B-n%7D,得:

    y%5E%7B-n%7D%5Cfrac%7Bdy%7D%7Bdx%7D%3DP(x)y%5E%7B1-n%7D%2BQ(x)                                                                          (2.21)

设:

    z%3Dy%5E%7B1-n%7D                                                                                                                (2.22)

求导得:

    %5Cfrac%7Bdz%7D%7Bdx%7D%3D(1-n)y%5E%7B-n%7D%5Cfrac%7Bdy%7D%7Bdx%7D                                                                                       (2.23)

代入(2.20)式得:

    %5Cfrac%7Bdz%7D%7Bdx%7D%3D(1-n)P(x)z%2B(1-n)Q(x)                                                           (2.24)

(2.24)式即可按一阶线性微分方程求解,然后再代回原变量得到(2.20)式的通解.

(P.S.S需要注意的是,n>0时还需考虑y=0的解)


(由于专栏投稿只让放100张图片,公式也算图片的,所以剩余部分包括整章感想将放到下一节去)


常微分方程(二)的评论 (共 条)

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