第二章、一阶微分方程的初等解法[教学目标]1. 理解变量分离方程以及可化为变量分离方程的类型(齐次方程),熟练掌握变量分离方程的解法。
2. 理解一阶线性微分方程的类型,熟练掌握常数变易法及伯努力方程的求解。
3. 理解恰当方程的类型,掌握恰当方程的解法及简单积分因子的求法。
4. 理解一阶隐式方程的可积类型,掌握隐式方程的参数解法。
[教学重难点] 重点是一阶微分方程的各类初等解法 ,难点是积分因子的求法以及隐式方程的解法。
[教学方法] 讲授,实践。
[教学时间] 14学时[教学内容] 变量分离方程,齐次方程以及可化为变量分离方程类型,一阶线性微分方程及其常数变易法,伯努利方程,恰当方程及其积分因子法,隐式方程。
[考核目标]1.一阶微分方程的初等解法:变量分离法、一阶线性微分方程的常数变易法、恰当方程与积分因子法、一阶隐方程的参数解法。
2.会建立一阶微分方程并能求解。
§1 变量分离方程与变量变换1、 变量分离方程1) 变量分离方程形如()()dy f x g y dx= (或1122()()()()0M x N y dx M x N y dy +=) (2.1)的方程,称为变量分离方程,其中函数()f x 和()g y 分别是,x y 的连续函数.2) 求解方法如果()0g y ≠,方程(2.1)可化为, ()()dy f x dx g y = 这样变量就分离开了,两边积分,得到()()dy f x dx c g y =+⎰⎰ (2.2) 把,()()dy f x dx g y ⎰⎰分别理解为1,()()f x y ϕ的某一个原函数. 容易验证由(2.2)所确定的隐函数(,)y x c ϕ=满足方程(2.1).因而(2.2)是(2.1)的通解.如果存在0y 使0()0g y =,可知0y y =也是(2.1)的解.可能它不包含在方程的通解(2.2)中,必须予以补上.3) 例题例1 求解方程dy x dx y=- 解 将变量分离,得到ydy xdx =-两边积分,即得 22222y x c =-+ 因而,通解为22x y c += 这里的c 是任意的正常数.或解出显式形式y =例2 解方程 2cos dy y x dx= 并求满足初始条件:当0x =时.1y =的特解.解 将变量分离,得到2cos dy xdx y = 两边积分,即得 1sin x c y-=+因而,通解为 1sin y x c=-+ 这里的c 是任意的常数.此外,方程还有解0y =.为确定所求的特解,以0x =.1y =代入通解中确定常数c ,得到 1c =-因而,所求的特解为 11sin y x =- 例3 求方程()dy P x y dx = (2.3)的通解,其中()P x 是x 的连续函数.解 将变量分离,得到()dy P x dx y = 两边积分,即得 ln ()y P x dx c =+⎰这里的c 是任意常数.由对数的定义,即有 ()P x dx c y e +⎰= 即()P x dx c y e e ⎰=±令c e c ±=,得到 ()P x dx y ce ⎰=(2.4) 此外,0y =也是(2.3)的解.如果在(2.4)中允许0c =,则0y =也就包括在(2.4)中,因而,(2.3)的通解为(2.4),其中c 是任意常数.注: 1.常数c 的选取保证(2.2)式有意义.2.方程的通解不一定是方程的全部解,有些通解包含了方程的所有解,有些通解不能包含方程的所有解.此时,还应求出不含在通解中的其它解, 即将遗漏的解要弥补上.3.微分方程的通解表示的是一族曲线,而特解表示的是满足特定条件00()y x y =的一个解,表示的是一条过点00(,)x y 的曲线.2、可化为变量分离方程的类型1).形如dy y g dx x ⎛⎫= ⎪⎝⎭ (2.5)的方程,称为齐次方程,这里的()g u 是u 的连续函数.另外,ⅰ)对于方程 (,)(,)dy M x y dx N x y = 其中函数(,)M x y 和(,)N x y 都是x 和y 的m 次齐次函数,即对0t >有(,)(,)m M tx ty t M x y ≡ (,)(,)m N tx ty t N x y ≡ 事实上,取1t x=,则方程可改写成形如(2.5)的方程. (1,)(1,)(1,)(1,)m m y y x M M dy x x y y dx x N N x x== ⅱ)对方程 (,)dy f x y dx= 其中右端函数(,)f x y 是x 和y 的零次齐次函数,即对0t >有(,)(,)f tx ty f x y =则方程也可改写成形如(2.5)的方程(1,)dy y f dx x= 对齐次方程(2.5)利用变量替换可化为变量分离方程再求解.令y u x=(2.6)即y ux =,于是dy du x u dx dx =+ (2.7)将(2.6)、(2.7)代入(2.5),则原方程变为 ()du xu g u dx += 整理后,得到()du g u u dx x -= (2.8)方程(2.8)是一个可分离变量方程,按照变量分离法求解,然后将所求的解代回原变量,所得的解便是原方程(2.5)的解.例4 求解方程dy y y tg dx x x=+ 解 这是齐次方程,以,y dy du u x u x dx dx ==+代入,则原方程变为 du xu u tgu dx+=+ 即du tgu dx x= (2.9)分离变量,即有 dx ctgudu x=两边积分,得到 ln sin ln u x c =+这里的c 是任意的常数,整理后,得到sin u cx = (2.10)此外,方程(2.9)还有解0tgu =,即sin 0u =. 如果(2.10)中允许0c =,则sin 0u =就包含在(2.10)中,这就是说,方程(2.9)的通解为(2.10).代回原来的变量,得到原方程的通解为sin y cx x=例5 求解方程(0).dy x y x dx+=< 解 将方程改写为(0)dy y x dx x =< 这是齐次方程,以,y dy du u x u x dx dx==+代入,则原方程变为du xdx=(2.11)分离变量,得到dx x = 两边积分,得到(2.11)的通解ln()x c =-+即2[ln()](ln()0)u x c x c =-+-+>(2.12) 这里的c 是任意常数.此外,(2.11)还有解0u =注意,此解不包括在通解(2.12)中.代回原来的变量,即得原方程的通解2[ln()](ln()0)y x x c x c =-+-+>及解0y =.原方程的通解还可表为2[ln()],ln()0,0,x x c x c y ⎧-+-+>=⎨⎩它定义于整个负半轴上.注:1.对于齐次方程dy y g dx x ⎛⎫= ⎪⎝⎭的求解方法关键的一步是令y u x =后,解出y ux =,再对两边求关于x 的导数得dy du u x dx dx=+,再将其代入齐次方程使方程变为关于,u x 的可分离方程. 2.齐次方程也可以通过变换x v y =而化为变量分离方程.这时x vy =,再对两边求关于y 的导数得dx dv v y dy dy =+,将其代入齐次方程dx x f dy y ⎛⎫= ⎪⎝⎭使方程变为,v y 的可分离方程 小结:这一讲我们主要讲解了一阶微分方程的可分离变量法和齐次方程的dy y g dx x ⎛⎫= ⎪⎝⎭形状的解法.而这一齐次方程通过变量替换任然可化为可分离方程,因而,一定要熟练掌握可分离方程的解法.2)形如111222a x b y c dy dx a x b y c ++=++ (2.13)的方程经变量变换化为变量分离方程,这里的121212,,,,,a a b b c c 均为常数.分三种情况来讨论(1)120c c ==情形.这时方程(2.13)属齐次方程,有 1122a x b y dy y g dx a x b y x +⎛⎫== ⎪+⎝⎭此时,令y u x=,即可化为变量可分离方程. (2)11220a b a b =,即1122a b a b =的情形. 设1122a b k a b ==,则方程可写成22122222()()()k a x b y c dy f a x b y dx a x b y c ++==+++ 令22a x b y u +=,则方程化为22()du a b f u dx =+ 这是一变量分离方程.(3)1112220,a b c c a b ≠及不全为零的情形. 这时方程(2.13)右端的分子、分母都是,x y 的一次式,因此11122200a x b y c a x b y c ++=⎧⎨++=⎩ (2.14)代表xy 平面上两条相交的直线,设交点为(,)αβ.显然,0α≠或0β≠,否则必有120c c ==,这正是情形(1)(只需进行坐标平移,将坐标原点(0,0)移至(,)αβ就行了,若令X x Y y αβ=-⎧⎨=-⎩ (2.15)则(2.14)化为 112200a X bY a X b y +=⎧⎨+=⎩ 从而(2.13)变为1122a X bY dY Y g dX a X b Y X +⎛⎫== ⎪+⎝⎭ (2.16)因此,得到这种情形求解的一般步骤如下:(1)解联立代数方程(2.14),设其解为,x y αβ==;(2)作变换(2.15)将方程化为齐次方程(2.16);(3)再经变换Y u X=将(2.16)化为变量分离方程; (4)求解上述变量分离方程,最后代回原变量可得原方程(2.13)的解.上述解题的方法和步骤也适用于比方程(2.13)更一般的方程类型 111222a x b y c dy f dx a x b y c ⎛⎫+== ⎪++⎝⎭()dy f ax by c dx++ ()()0y xy dx xg xy dy += 2()dy x f xy dx=2dy y xf dx x⎛⎫= ⎪⎝⎭ 以及 (,)()(,)()0M x y xdx ydy N x y xdy ydx ++-=(其中,M N 为,x y 的齐次函数,次数可以不相同)等一些方程类型,均可通过适当的变量变换化为变量分离方程.例6 求解方程 13dy x y dx x y -+=+- (2.17)解 解方程组 1030x y x y -+=⎧⎨+-=⎩ 得1, 2.x y ==令12x X y Y =+⎧⎨=+⎩ 代入方程(2.17),则有dY X Y dX X Y -=+ (2.18)再令 Y u X =即 Y uX = 则(2.18)化为2112dX u du X u u+=-- 两边积分,得 22ln ln 21X u u c =-+-+22(21)c X u u e +-=±记1,c e c ±=并代回原变量,就得2212Y XY X c +-=221(2)2(1)(2)(1)y x y x c -+----=此外,易验证2210u u +-=即2220Y XY X +-=也就是(2.18)的解.因此方程(2.17)的通解为22262y xy x y x c +---=其中c 为任意的常数.3、 应用举例例7 电容器的充电和放电如图(2.1)所示的R C -电路,开始时电容C 上没有电荷,电容两端的电压为零.把开关K 合上“1”后,电池E 就对电容C 充电,电容C 两端的电压C u 逐渐升高,经过相当时间后,电容充电完毕,再把开关K 合上“2”,这时电容就开始放电过程,现在要求找出充、放电过程中,电容C 两端的电压C u 随时间t 的变化规律.解 对于充电过程,由闭合回路的基尔霍夫第二定理,c u RI E += (2.19)对于电容C 充电时,电容上的电量Q 逐渐增多,根据C Q Cu =,得到 ()C C du dQ d I Cu C dt dt dt === (2.20)将(2.20)代入(2.19),得到c u 满足的微分方程 c c du RCu E dt+= (2.21)这里R 、C 、E 都是常数.方程(2.21)属于变量分离方程.将(2.21)分离变量,得到C C du dtu E RC=-- 两边积分,得到11ln C u E t c RC-=-+ 即1112t t c RCRCC u E e e c e---=±=这里12c c e =±为任意常数.将初始条件:0t =时,0C u =代入,得到2c E =-. 所以 1(1)t RC C u E e -=-(2.22)这就是R C -电路充电过程中电容C 两端的电压的变化规律.由(2.22)知道,电压C u 从零开始逐渐增大,且当t →+∞时,C u E →,在电工学中,通常称RC τ=为时间常数,当3t τ=时,0.95C u E =,就是说,经过3τ的时间后,电容C 上的电压已达到外加电压的95%.实用上,通常认为这时电容C 的充电过程已基本结束.易见充电结果C u E =.对于放电过程的讨论,可以类似地进行.例8 探照灯反射镜面的形状在制造探照灯的反射镜面时,总是要求将点光源射出的光线平行地射出去,以保证照灯有良好的方向性,试求反射镜面的几何形状.解 取光源所在处为坐标原点,而x 轴平行于光的反射方向,设所求曲面由曲线()y f x z =⎧⎨=⎩(2.23)绕x 轴旋转而成,则求反射镜面的问题归结为求xy 平面上的曲线()y f x =的问题,仅考虑0y >的部分,过曲线()y f x =上任一点(,)M x y 作切线NT ,则由光的反射定律:入射角等于反射角,容易推知12αα= 从而OM ON = 注意到2dy MP tg dx NPα==及,,OP x MP y OM ===就得到函数()y f x =所应满足的微分方程式dy dx =(2.24)这是齐次方程.由2.12知引入新变量xu y=可将它化为变量分离方程.再经直接积分即可求得方程的解.对于方齐次方程(2.24)也可以通过变换xv y=而化为变量分离方程也可由x yv =得dx dvv y dy dy=+代入(2.24)得到sgn dvv y v y dy+=+于是sgn dy y y =(2.25)积分(2.25)并代回原来变量,经化简整理,最后得2(2)y c c x =+(2.26)其中c 为任意常数.(2.26)就是所求的平面曲线,它是抛物线,因此,反射镜面的形状为旋转抛物面22(2)y z c c x +=+ (2.27)小结: 本节我们主要讨论了一阶可分离微分方程和齐次微分方程的求解问题.将各种类型的求解步骤记清楚的同时要注意对解的讨论.§2 线性方程与常数变易法1、一阶线性微分方程()()()0dya xb x yc x dx++= 在()0a x ≠的区间上可以写成()()dyP x y Q x dx=+ (2.28)对于()a x 有零点的情形分别在()0a x ≠的相应区间上讨论.这里假设(),()P x Q x 在考虑的区间上是x 的连续函数.若()0Q x ≡,(2.28)变为 ()dyP x y dx= (2.3)称为一阶齐线性方程.若()0Q x ≠,(2.28)称为一阶非齐线性方程.2、常数变易法(2.3)是变量分离方程,已在例3中求得它的通解为 ()P x dxy ce ⎰=(2.4)这里c 是任意的常数.下面讨论一阶非齐线性方程(2.28)的求解方法.方程(2.3)与方程(2.28)两者既有联系又有区别,设想它们的解也有一定的联系,在(2.4)中c 恒为常数时,它不可能是(2.28)的解,要使(2.28)具有形如(2.4)的解, c 不再是常数,将是x 的待定函数()c x ,为此令 ()()P x dxy c x e ⎰=(2.29) 两边微分,得到()()()()()P x dxP x dx dy dc x e c x P x e dx dx⎰⎰=+ (2.30)将(2.29)、(2.30)代入(2.28),得到()()()()()()()()()P x dxP x dx P x dx dc x e c x P x e P x c x e Q x dx⎰⎰⎰+=+ 即()()()P x dx dc x Q x e dx-⎰= 积分后得到()()()P x dxc x Q x e dx c -⎰=+⎰(2.31)这里c 是任意的常数..将(2.31)代入(2.29),得到()()()()()() =()P x dxP x dx P x dx P x dx P x dxy e Q x e dx c ce e Q x e dx--⎛⎫⎰⎰=+ ⎪⎝⎭⎰⎰⎰+⎰⎰(2.32)这就是方程(2.28)的通解.这种将常数变易为待定函数的方法,通常称为常数变易法.实际上常数变易法也是一种变量变换的方法.通过变换(2.29)可将方程(2.28)化为变量分离方程.注: 非齐线性方程的通解是它对应的齐线性方程的通解与它的某个特解之和. 例1 求方程1(1)(1)x n dyx ny e x dx++-=+的通解,这里的n 为常数. 解 将方程改写为 (1)1x n dy n y e x dx x -=++ (2.33)先求对应的齐次方程01dy n y dx x -=+ 的通解,得(1)n y c x =+令 ()(1)n y c x x =+ (2.34) 微分之,得到()(1)(1)()n dy dc x x n x c x dx dx=+++ (2.35)以(2.34)、(2.35)代入(2.33),再积分,得 ()x c x e c =+ 将其代入公式(2.34),即得原方程的通解 (1)()n x y x e c =++ 这里c 是任意的常数. 例2 求方程22dy ydx x y=-的通解. 解 原方程改写为2dx x y dy y=- (2.36)把x 看作未知函数,y 看作自变量,这样,对于x 及dxdy来说,方程(2.36)就是一个线性方程了.先求齐线性方程2dx x dy y= 的通解为2x cy = (2.37)令2()x c y y =,于是 2()2()dx dc y y c y y dy dy=+ 代入(2.36),得到()ln c y y c =-+ 从而,原方程的通解为2(ln )x y c y =-这里c 是任意的常数,另外0y =也是方程的解. 特别的,初值问题00()()()dyP x y Q x dxy x y ⎧=+⎪⎨⎪=⎩ 的解为00()()()=()xxsx x x P d P d P d xx y ceeQ s eds ττττττ-⎰⎰⎰+⎰例3 试证(1)一阶非齐线性方程(2.28)的任两解之差必为相应的齐线性方程(2.3)之解;(2)若()y y x =是(2.3)的非零解,而()y y x =是(2.28)的解,则(2.28)的通解可表为()()y cy x y x =+,其中c 为任意常数.(3)方程(2.3)任一解的常数倍或两解之和(或差)仍是方程(2.3)的解. 证 (1)设12,y y 是非齐线性方程的两个不同的解,则应满足方程使1122()(1)()(2)dy py Q x dxdy py Q x dx=+=+(1)—(2)有1212()()d y y p y y dx-=-说明非齐线性方程任意两个解的差12y y -是对应的齐次线性方程的解.(2)因为(()())()()(()()()()d cy x y x dy x d y x c p cy p y Q x p cy y Q x dx dx dx+=+=++=++故结论成立.(3)因为12121212()()()(),(),()d y y d y y d cy p cy p y y p y y dx dx dx+-==+=- 故结论成立.3、Bernoulli 方程形如()()n dyP x y Q x y dx=+ ( 0,1n ≠) (2.38)的方程,称为伯努利(Bernoulli )方程,这里(),()P x Q x 为x 连续函数.利用变量变换可将伯努利方程化为线性方程来求解.事实上,对于0y ≠,用n y -乘(2.38)两边,得到1()()n n dyy y P x Q x dx--=+ (2.39) 引入变量变换1n z y -= (2.40) 从而(1)ndz dyn y dx dx-=- (2.41)将(2.40)、2.41)代入(2.39),得到 (1)()(1)()dzn P x z n Q x dx=-+- (2.42)这是线性方程,用上面介绍的方法求得它的通解,然后再代回原来的变量,便得到(2.38)的通解.此外,当0n >时,方程还有解0y =. 例4 求方程26dy yxy dx x=-的通解 解 这是2n =时的伯努利方程,令 1z y -=,得 2dz dyy dx dx-=- 代入原方程得到6dz z x dx x=-+ 这是线性方程,求得它的通解为268c x z x =+代回原来的变量y ,得到2618c x y x =+或者688x x c y -= 这是原方程的通解. 此外,方程还有解0y =. 例5 求方程331dy dx xy x y=+的解 解 将方程改写为33dxyx y x dy=+ 这是一个自变量为y ,因变量为x 的伯努利方程.解法同上.例6 求方程23y dy e x dx x+=的通解 这个方程只要做一个变换,令,y ydu dyu e e dx dx==,原方程改写为 22231du x u u dx x x=+ 便是伯努利方程.小结;这次主要讨论了一阶线性微分方程的解法.其核心思想是常数变易法.即将非齐线性方程对应的齐线性方程解的常数变易为待定函数,使其变易后的解函数代入非齐次线性方程,求出待定函数()c x ,求出非齐次方程的解.我们还讨论了伯努利方程,求解过程为,先变换,将原方程化为非齐线性方程,再求解.§3 恰当方程与积分因子1、恰当方程的定义 将一阶微分方程(,)dyf x y dx= 写成微分的形式(,)0f x y dx dy -=把,x y 平等看待,对称形式的一阶微分方程的一般式为(,)(,)0M x y dx N x y dy += (2.43)假设(,),(,)M x y N x y 在某区域G 内是,x y 的连续函数,而且具有连续的一阶偏导数. 如果存在可微函数(,)u x y ,使得(,)(,)du M x y dx N x y dy =+ (2.44)即 (,), (,)u u M x y N x y x y∂∂==∂∂ (2.45)则称方程(2.43)为恰当方程,或称全微分方程.在上述情形,方程(2.43)可写成(,)0du x y ≡,于是 (,)u x y C ≡就是方程(2.43)的隐式通解,这里C 是任意常数(应使函数有意义). 2、 恰当方程的判定准则定理1设(,),(,)M x y N x y 在某区域G 内连续可微,则方程(2.43)是恰当方程的充要条件是, (,)M Nx y G y x∂∂=∈∂∂ (2.46)而且当(2.46)成立时,相应的原函数可取为 00(,)(,)(,)xyx y u x y M s y ds N x t dt =+⎰⎰(2.47) 或者也可取为0(,)(,)(,)yxy x u x y N x t dt M s y ds =+⎰⎰(2.48)其中00(,)x y G ∈是任意取定的一点.证明 先证必要性.因为(2.43)是恰当方程,则有可微函数(,)u x y 满足(2.45), 又知(,),(,)M x y N x y 是连续可微的,从而有22M u u Ny y x x y x∂∂∂∂===∂∂∂∂∂∂ 下面证明定理的充分性,即由条件(2.46),寻找函数(,)u x y ,使其适合方程(2.45).从(2.47)可知(,)uN x y y∂=∂ 000000(,)(,) =(,)(,) =(,)(,)(,)yy y x y yy y u M x y N x t dt x x M x y N x t dtM x y M x t dt M x y ∂∂=+∂∂++=⎰⎰⎰即(2.45)成立,同理也可从(2.48)推出(2.45). 例1. 解方程21()02x xydx dy y++=(2.49)解 这里21, =()2x M xy N y=+,则y x M x N ==,所以(2.49)是恰当方程.因为N 于0y =处无意义,所以应分别在0y >和0y <区域上应用定理2.3,可按任意一条途径去求相应的原函数(,)u x y .先选取00(,)(0,1)x y =,代入公式(2.47)有 22011()ln 22xyx x u xdx dy y y y =++=+⎰⎰再选取00(,)(0,1)x y =-,代入公式(2.47)有22011()()ln()22xyx x u x dx dy y y y -=-++=+-⎰⎰ 可见不论0y >和0y <,都有2ln ||2x u y y =+ 故方程的通解为2ln ||2x y y C +=. 3、恰当方程的解法上述定理已给出恰当方程的解法,下面给出恰当方程的另两种常用解法. 解法1. 已经验证方程为恰当方程,从(,)x u M x y =出发,有2(,)(,)()()2x u x y M x y dx y y y φφ≡+=+⎰ (2.50)其中()y φ为待定函数,再利用(,)y u N x y =,有221()22x x y y φ'+=+ 从而1()y yφ'= 于是有 ()ln ||y y φ=只需要求出一个(,)u x y ,因而省略了积分常数.把它代入(2.50)便得方程的通解为2ln ||2x u y y C =+= 解法2. 分项组合的方法 对(2.49)式重新组合变为21()02x xydx dy dy y ++= 于是 2()ln ||02x d y d y +=从而得到方程的通解为 2ln ||2x y y C += 4、积分因子的定义及判别对于微分形式的微分方程(,)(,)0M x y dx N x y dy +=(2.43)如果方程(2.43)不是恰当方程,而存在连续可微的函数(,)0x y μμ=≠,使得(,)(,)0M x y dx N x y dy μμ+= (2.51)为一恰当方程,即存在函数(,)v x y ,使(,)(,)M x y dx N x y dy dv μμ+≡则称(,)x y μ是方程(2.43)的积分因子.此时(,)v x y C =是(2.51)的通解,因而也就是(2.43)的通解.如果函数(,),(,)M x y N x y 和(,)x y μ都是连续可微的,则由恰当方程的判别准则知道, (,)x y μ为(2.43)积分因子的充要条件是M Ny xμμ∂∂=∂∂ 即 ()M N N M x y y xμμμ∂∂∂∂-=-∂∂∂∂ (2.52)5、积分因子的求法方程(2.52)的非零解总是存在的,但这是一个以μ为未知函数的一阶线性偏微分方程,求解很困难,我们只求某些特殊情形的积分因子. 定理2 设(,),(,)M M x y N N x y ==和(,)x y ϕϕ=在某区域内都是连续可微的,则方程(2.43) 有形如((,))x y μμϕ=的积分因子的充要条件是:函数(,)(,)(,)(,)(,)(,)y x x y M x y N x y N x y x y M x y x y ϕϕ--(2.53)仅是(,)x y φ的函数,此外,如果(2.53)仅是(,)x y φ的函数((,))f f x y ϕ=,而()()G u f u du =⎰,则函数((,))G x y e ϕμ= (2.54)就是方程(2.43)的积分因子.证明 因为如果方程(2.43)有积分因子()μμϕ=,则由(2.52)进一步知()()d M N N M d x y y xμϕϕμϕ∂∂∂∂-=-∂∂∂∂ 即y x x yM N d d N M μϕμϕϕ-=-由()μμϕ=可知左端是ϕ的函数,可见右端y x x yM N N M ϕϕ--也是ϕ的函数,即()y x x yM N f N M ϕϕϕ-=-,于是,有()d f d μϕϕμ=, 从而 ()()f d G e e ϕϕϕμ⎰==反之,如果(2.53)仅是ϕ的函数,即()y x x yM N f N M ϕϕϕ-=-,则函数(2.54)是方程(2.52)的解.事实上,因为()()()()G x y y x NM N M f e M N x yϕμμϕϕϕμ∂∂-=-=-∂∂ 因此函数(2.54)的确是方程(2.43)的积分因子. 为了方便应用这个定理,我们就若干特殊情形列简表如下: 例解2(3y -解这里2231,M y xy N xy x =-+=-,注意y x M N y x -=-所以方程不是恰当的,但是1y xM N Nx-=它仅是依赖与x ,因此有积分因子1dx xe x μ⎰≡=给方程两边乘以因子x μ=得到2223(3)()0xy x y x dx x y x dy -++-=从而可得到隐式通解22321122u x y x y x C ≡-+=例3. 解方程2()(1)0xy y dx xy y dy ++++=解 这里2,1M xy y N xy y =+=++方程不是恰当的.但是1y xM N My-=-- 它有仅依赖于y 的积分因子 11dyy eyμ-⎰≡=方程两边乘以积分因子1y μ=得到 1()(1)0x y dx x dy y++++= 从而可得到隐式通解21ln ||2u x xy y y C ≡+++= 另外,还有特解0y =.它是用积分因子乘方程时丢失的解. 例4. 解方程 223(2)()0y x y dx xy x dy +++=解 这里2232,M y x y N xy x =+=+,不是恰当方程.设想方程有积分因子()x y αβμμ=,其中α,β是待定实数.于是2112111()(2)y xM N y x x N y M x y y x x y x y αβαβαβαβαβαβ----⋅=⋅=--+-只须取3,2αβ==.由上述简表知原方程有积分因子32x y μ=从而容易求得其通解为:446313u x y x y C ≡+=六、积分因子的其他求法以例4为例,方程的积分因子也可以这样来求:把原方程改写为如下两组和的形式:223()(2)0y dx xydx x ydx x dy +++=前一组有积分因子11yμ=,并且 21()()y dx xydy d xy y+=后一组有积分因子21xμ=,并且 2321(2)()x ydx x dy d x y x+= 设想原方程有积分因子211()()xy x y y xαβμ== 其中α,β是待定实数.容易看出只须3,2αβ==,上述函数确实是积分因子,其实就是上面找到一个.例5. 解方程 1212()()()()0M x M y dx N x N y dy += 其中1M ,2M ,1N ,2N 均为连续函数.解 这里12()()M M x M y =,12()()N N x N y =.写成微商形式就形式上方程是变量可分离方程,若有0y 使得20()0M y =,则0y y =是此方程的解;若有0x 使得10()0N x =,则0x x =是此方程的解;若21()()0M y N x ≠,则有积分因子211()()M y N x μ=并且通解为1212()()()()M x N y u dx dy N x M y ≡+⎰⎰ 例6、试用积分因子法解线性方程(2.28).解 将(2.28)改写为微分方程[()()]0P x y Q x dx dy +-=(2.55)这里()(),1M P x y Q x N =+=-,而()M Ny xP x N∂∂-∂∂=- 则线性方程只有与x 有关的积分因子()P x dxe μ-⎰= 方程(2.55)两边乘以()P x dxe μ-⎰=,得()()()()()0P x dx P x dx P x dxxP x e ydx e dy Q x e dx ---⎰⎰⎰-+= (2.56)(2.56)为恰当方程,又分项分组法()()()()0P x dx P x dxd ye Q x e dx --⎰⎰-=因此方程的通解为()()()P x dx P x dxye Q x e dx c --⎰⎰-=⎰即()()[()]P x dxP x dxy e Q x e dx c -⎰⎰=+⎰与前面所求得的结果一样.注:积分因子一般不容易求得可以先从求特殊形状的积分因子开始,或者通过观察法进行“分项分组”法求得积分因子.§4 一阶隐方程与参数表示1、一阶隐方程一阶隐式微分方程的一般形式可表示为:(,,)0F x y y '=如果能解出(,)y f x y '=,则可化为显式形式,根据前面的知识求解. 例如方程2()()0y x y y xy ''-++=,可化为y x '=或y y '=但难以从方程中解出y ',或即使解出y ',而其形式比较复杂,则宜采用引进参数的方法求解.一般隐式方程分为以下四种类型:1) (,)y f x y '= 2) (,)x f y y '= 3) (,)0F x y '= 4)(,)0F y y '=2、求解方法Ⅰ)可以解出y (或)x 的方程1) 讨论形如(,)y f x y '= (2.57)的方程的解法,假设函数(,)f x y '有连续的偏导数,引进参数y p '=,则方程(2.57)变为 (,)y f x p = (2.58)将(2.58) 的两边对x 求导数,得到 f f dp p x y dx∂∂=+∂∂ (2.59)方程(2.59)是关于,x p 的一阶微分方程,而且属于显式形式.若求得(2.59)的通解形式为(,)p x c ϕ=,将其代入(2.58),于是得到(2.57)通解为 (,(,))y f x x c ϕ=若求得(2.59)的通解形式为(,)x p c ψ=,于是得到(2.57)的参数形式的通解为(,)((,),)x p c y f p c p ψψ=⎧⎨=⎩其中p 为参数, c 是任意常数.若求得(2.59)的通解形式为(,,)0x p c Φ=,于是得到(2.57)的参数形式的通解为(,,)0(,)x p c y f x p Φ=⎧⎨=⎩其中p 为参数, c 是任意常数. 例1 求方程3()20dy dyx y dx dx+-= 的解 解 令dyp dx=,于是有32y p xp =+ (2.60)两边对x 求导数,得到2322dp dp p p x p dx dx=++ 即 2320p dp xdp pdx ++= 当0p ≠时,上式有积分因子p μ=,从而32320p dp xpdp p dx ++= 由此可知4234p xp c += 得到42223344c pc x p p p -==- 将其代入(2.60),即得43342()c p y p p-=+故参数形式的通解为22334(0) 212c x p p p c y p p ⎧=-⎪⎪≠⎨⎪=-⎪⎩当0p =时,由(2.60)可知0y =也是方程的解.例2 求方程22()2dy dy x y xdx dx =-+的解. 解 令dy p dx=,得到222x y p xp =-+(2.61)两边对x 求导数,得到2dp dp p px p x dx dx =--+ 或 (2)(1)0dpp x dx--= 由10dpdx-=,解得p x c =+,于是得到方程的通解为222x y cx c =++ (2.62)由20p x -=,解得2xp =,于是得到方程的一个解为24x y =(2.63)特解(2.63)与通解(2.62)中的每一条积分曲线均相切,因此称为方程的奇解.2) 讨论形如(,)dyx f y dx= (2.64)的方程的求解方法,方程(2.64)与方程(2.57)的求解方法完全类似,假定函数(,)f y y ' 有连续偏导数. 引进参数dyp dx=,则(2.64) 变为 (,)x f y p = (2.65)将(2.65) 的两边对y 求导数,得到 1f f dp p y x dy∂∂=+∂∂ (2.66)方程(2.66))是关于,y p 的一阶微分方程,而且属于显式形式.设其通解为 (,,)0y p c Φ= 则(2.64)的通解为(,,)0(,)y p c x f y p Φ=⎧⎨=⎩Ⅱ)不显含y (或)x 的方程 3) 讨论形如(,)0F x y '=(2.67) 的方程的解法.记dyp y dx'==,此时(,)0F x p =表示的是xp 平面上的一条曲线,设曲线用参数形式表示为()x t ϕ=,()p t ψ= (2.68)由于dy pdx =,进而()()dy t t dt ψϕ'= 两边积分,得到()()y t t dt c ψϕ'=+⎰于是得到方程(2.67)参数形式的解为()()()x t y t t dt c ϕψϕ=⎧⎪⎨'=+⎪⎩⎰c 是任意常数.例3 求解方程3330x y xy ''+-=解 令y p tx '==,则由方程得331tx t=+, 2331t p t =+ 于是 23339(12)(1)t t dy dt t -=+ 积分得到23333329(12)314(1)2(1)t t t y dt c c t t -+=+=+++⎰故原方程参数形式的通解为:3332313142(1)t x t t y c t ⎧=⎪+⎪⎨+⎪=+⎪+⎩4) 讨论形如(,)0F y y '=(2.69)的方程,其解法与方程(2.67)的求解方法类似.记dyp y dx'==,此时(,)0F y p =表示的是yp 平面上的一条曲线,设曲线用参数形式表示为()y t ϕ=,()p t ψ= 由关系式dy pdx =可知 ()()t dt t dx ϕψ'=,于是0p ≠时,有 ()()t dx dt t ϕψ'=, ()()t x dt c t ϕψ'=+⎰ 故方程(2.69)的参数形式的通解()()()t x dt c t y t ϕψϕ'⎧=+⎪⎨⎪=⎩⎰c 是任意常数.此外,不难验证,若(,0)0F y =有实根y k =,则y k =也是方程的解.例4 求解方程 22(1)(2)y y y ''-=-.解 令2y yt '-=,则有222(1')y y y t -=由此可以得2'1y t =-,1y t t=+ 代入1dx dy p=,得到 222111(1)1dx dt dt t t t=-+=-- 积分,得到1x c t=+故原方程参数形式的通解为 11x c t y t t ⎧=+⎪⎪⎨⎪=+⎪⎩其中c 是任意常数.此外, 当0y '=时原方程变为24y =,于是2y =±也是方程的解. 例5求解方程y '=解 令y p '=,则有p =,取,(,)22p tgt t ππ=∈-,则sin sec tgt x t t === 由dy pdx =得到cos sin dy tgt tdt tdt ==所以cos y t c =-+故原方程参数形式的通解为sin cos x t y t c =⎧⎨=-+⎩ 其中c 是任意常数.(注:本资料素材和资料部分来自网络,仅供参考。