微分方程积分因子的研究

合集下载

从一道微分方程解题想到的积分因子变换

从一道微分方程解题想到的积分因子变换

由公式得:
dx x
∫ =y
2
e

1dx x
2
x3e
−2

1dx x
dx
+
c
( ) 通解 y = x2 x2 + C 。
解法三:观察法
对于 x2

=
2x ,
于是所求通解为:
y x2
=
x2
+C
( ) y = x2 x2 + C
观察法实际,在于构建全微分,其本质是积分因子法。 解法四:积分因子法 方程可化为:
1 问题背景
一阶微分方程的解法是多种多样的,我们在解题的过程中根据方
程的特点,通过变换可以将其上升为解题方法。
例 1 解微分方程
[1]。
解法一:常系数变易法
由线性齐次方程:
∂M − ∂N

令 µ= (x)
e∫
∂y ∂x dx
N=
e= ln x−2
1 ,有: x2
(2
ydx − x2
xdy)−
2
xdx
= 0
d (−
y x2
+x2 )
= 0
∴ 通解为:
y x2
=
x2
+C
得解
。设 y = C(x)x2 为原方程的解,则:
C′(x) = 2x , C(x) = x2 + C
即通解:
( ) y = x2 x2 + C
解法二:公式法
对于一阶线性齐次方程 dy − 2 =2x3
238
理论研究
从一道微分方程解题想到的积分因子变换
杨 冯 , 赵临龙 * (安康学院 数学与统计学院 , 陕西 安康 725000)

全微分方程与积分因子法

全微分方程与积分因子法

已构成全微分的项分出再把剩下的项凑成全微分.但这种方法
要求熟记一些简单二元函数的全微分,如
ydx+xdy=d(x,y)
ydx-xdy y2
=d(
x y

-ydx+xdy x2
=d(
x y

ydx-xdy =d(ιn| x |)
xy
y
ydx-xdy x2+y2
=d(arctg
x y

| | ydx-xdy x2-y2
的通解为
μ(x,y)=∫x0xP(x,y)dx+∫y0xQ(x,y)dy=C
(7)
其中点(x0,y0)可在与路径无关的单连通区域 G 内 任 意 取
得.很 多 情 况 下 都 选 (0,0)为 (x0,y0),只 有 当 点 (0,0)不 在 上 述
单连通区域 G 内,才考虑其他点作为曲线积分的始点.
坠p - 坠Q 坠y 坠x
-P
这里 φ 仅为 y 的函数.从而求得方程 (1)的一个积分因子 μ=
e 。 ∫φ(y)dy
例 4 试用公式法解线性微分方程(8)
解 : 将 (8)式 改 写 成 [Q(x)-P(X)Y]DX-DY=0
(10)
这时由公式,μ(x)=e∫p(x)dx.以 μ(x)=e∫p(x)dx 乘上(10)式得到
或 y=e-∫p(x)dx[∫Q(x)e∫p(x)dxdx+C]
2.公 式 法
由同一个方程
ydx-xdy=0
可以有不同的积分因子 1 y2

1 x2

1和 1 xy x2±y2
.可以证明,只要方程有解,则必有积分因子存在,
并且不是唯一的.因此,在具体解题过程中,由于求出的积分因

浅谈积分因子及首次积分

浅谈积分因子及首次积分

浅谈积分因子与首次积分摘要:本文先给出了微分方程中的积分因子、首次积分以及特征方程的相关定义并加深理解,后引出全微分方程积分因子存在的充要条件以及与之相关的两类重要命题,灵活的将用积分因子解微分方程的方法与偏微分方程首次积分联系起来,为求特殊积分因子提供了方便,最后应用性的求出了常见的几类微分方程的积分因子.关键词:微分方程;积分因子;首次积分;特征方程;偏微分:合分比Introduction to integral factor and the points for the first timeChen Xueyun(School of Mathematics and Statistics,Tianshui Normal University 741000) Abstract This paper firstly presents the definition of the integral factors ,first integral in differential equation and the characteristic equation and leads to the necessary and sufficient condition for the existence of all the integrating factor of differential equation as well as in connection with the two important types of proposition, Then it provides conveniences for special integral factor by combining the method of integral factor to solve differential equations with partial differential equation flexibly,Finally it finds out the integral factor of some types of differential equations via application.Keywords Differential equations,Integrating factor,For the first time points,Characteristic equation, Partial differential,points than目录0 引言 (1)1 相关概念 (1)1.1 积分因子 (1)1.2 首次积分 (2)1.3 特征方程(组) (2)2 预备知识 (2)3 问题的引入 (3)4 有关积分因子存在的充要条件的两类命题 (5)4.1 命题一 (5)4.2 命题二 (7)5 几类常见微分方程的积分因子 (8)5.1 可分离变量的微分方程 (8)5.2 齐次微分方程 (9)5.3 线性微分方程 (10)5.4 Bornoulli微分方程 (10)6 小结 (12)参考文献 (13)致谢 (14)浅谈积分因子与首次积分0 引言在微分方程的学习中,微分方程的求解是至关重要的,并且对于不同类型的微分方程可以给出不同的解法.Euler 曾在他的论文中指出:凡是可用分离变量法的地方都可以用积分因子法,这就启示我们要更多的了解和学习积分因子法.而这其中恰当微分方程就可以通过积分的方法求出,但并非所有的微分方程均为恰当微分方程,于是如何将一个非恰当微分方程转化为恰当微分方程,让求其通解变得简单?将是一个很显眼的问题,对此本文就特殊积分因子的求法引进了首次积分法,并加以论证运用.1 相关概念1.1 积分因子定义1 对于下面的一阶方程 (,)dy f x y dx= 现把它写成微分的形式(,)0f x y dx dy -=或者把,x y 平等的处理,写成具有对称形式的一阶微分方程,如下(,)(,)0,M x y dx N x y dy += (1.1.1)其中(,)M x y 和(,)N x y 是关于,x y 的连续函数且具有连续的一阶偏导数,(,)x y D ∈,D 为单连通区域.这样如果存在连续可微的函数(,)0x y μμ=≠,使得(,)(,)(,)(,)0x y M x y dx x y N x y dy μμ+=,为一恰当微分方程,即存在函数(,)u x y ,使 u u Mdx Ndy du dx dy x yμμ∂∂+≡=+∂∂, 则称(,)x y μ为方程(1.1.1)的积分因子.1.2 首次积分定义2 对于一般的常微分方程组 1112221212(,,,,),(,,,,),(,,,,).n n n n n dy f x y y y dx dy f x y y y dx dy f x y y y dx⎧=⎪⎪⎪=⎪⎨⎪⎪⎪=⎪⎩ (1.2.1)其中,右端函数12,,,n f f f 都在某个域G 内连续.设函数12(,,,,)n x y y y Φ=Φ在域G 内连续可微,且12(,,,,)n x y y y const Φ=Φ≠,如果以方程组的(1.2.1)的任一解()(1,2,)i y x i n =代入函数Φ时,使得函数12(,(),(),,())n x y x y x y x const Φ=Φ=(其中const 为任意常数),且此常数与x 无关,同时此常数值随不同解而异,则称表达式12=(,,,,)n x y y y const ΦΦ=为方程(1.2.1)的一个首次积分,有时也称12(,,,,)n x y y y Φ为首次积分.1.3 特征方程(组)定义3 对于一阶齐次线性偏微分方程 121(,,,)0n i n i iX x x x x μ=∂=∂∑, (1.3.1) 其中12,,,n x x x 是自变量,μ是12,,,n x x x 的未知函数,函数()(1,2,,)0i iX x i n x μ∂==∂为域'G 内相应的已知函数,则常微分方程组 1212n ndx dx dx X X X ===, (1.3.2) 称为一阶齐次线性偏微分方程(1.3.1)的特征方程组,也称特征方程.2 预备知识引理[1] 函数(,)x y μ是方程(,)(,)0M x y dx N x y dy +=的积分因子的充要条件是()()M N y xμμ∂∂=∂∂. (2.1) 现在对(2.1)式做如下变形 M N M N y y x xμμμμ∂∂∂∂+=+∂∂∂∂, M N N M x y y x μμμ⎛⎫∂∂∂∂-=- ⎪∂∂∂∂⎝⎭, 则可得到如下以μ为未知函数的一阶线性偏微分方程 ln ln M N N M x y y xμμ∂∂∂∂-=-∂∂∂∂. (2.2) 3 问题的引入通过对常微分方程的学习,我们了解了一阶微分方程的各种解法,其中对于恰当微分方程(,)(,)0M x y dx N x y dy +=我们可以通过积分求出它的通解,但是在学习研究中所遇到的微分方程大多并非全微分方程,于是我们引进了积分因子,以便更好的将一个非恰当微分微分方程转化为恰当微分方程,对于一般的积分因子我们可以通过观察法、分组法、公式法等来求解,那么对于特殊的又当如何呢?能否关联到有关偏微分方程的理论呢?下面的定理告诉我们这是肯定的.定理[2] 一阶微分方程(,)(,)0M x y dx N x y dy +=有可求的积分因子⇔方程(2.2)的特征方程有可求的首次积分.证 )⇒ 由引理我们知道如果存在函数(,)x y μ,使(,)(,)(,)(,)x y M x y d x x y N x y μμ+= 为全微分方程,则有方程(2.1)成立,而方程(2.1)又与偏微分方程(2.2)等价,所以想要求解微分方程(1.1.1)关键是要求出积分因子(,)x y μ,而要求出积分因子(,)x y μ的关键是求解(2.2)式即偏微分方程 ln ln M N NM x y y x μμ∂∂∂∂-=-∂∂∂∂,现在令 ln t μ=, M N T y x∂∂=-∂∂, (3.1) 则方程(2.2)可变形为如下偏微分方程 t t N M T x y∂∂-=∂∂, (3.2) 由定义3现在可以写出(3.2)的特征方程 dx dy dt N M T==-, (3.3) 所以求解积分因子(,)x y μ的关键就在于求出(3.3)的首次积分,假设现在求出特征方程(3.3)式的两个首次积分11ϕ(x,y,t )=c ,22ϕ(x,y,t )=c ,并且让它们相互独立,即雅可比行列式 1112220(,)x y D D x y x yϕϕϕϕϕϕ∂∂∂∂=≠∂∂∂∂(,), (3.4) 那么,若120ϕϕΦ=(,),并能从120ϕϕΦ=(,)中确定函数ln (,)t x y μϕ==, 则120ϕϕΦ=(,)为方程(2.2)的通解.)⇐ 假设特征方程(3.3)存在首次积分(,,)x y t c ϕ=()c 为常数,由 1ln (,)(,)t x y x y μϕ==可得 1(,)(,)x y x y e ϕμ=.把1(,)(,)x y x y e ϕμ=代入ln ln M N N M x y y xμμ∂∂∂∂-=-∂∂∂∂,经验证成立. 故(,)(,)(,)(,)0x y M x y dx x y N x y μμ+=是一个全微分方程,并且(,)x y μ为方程(,)(,)0M x y d x N x y d y +=的一个积分因子.例1 用首次积分求解方程组22,().()dx y dt y x dy x dt y x ⎧=⎪-⎪⎨⎪=⎪-⎩ 解 首先将两式相比可得 dx y dy x=, 即 0xdx ydy -=对上式积分我们就可以得到方程组的一个首次积分2211x y c ψ=-=.再次将两式作差得到 2()()()d x y x y dt x y ---=-, 即 ()()0dt x y d x y +--=.对上式关于(x-y)积分可得 2221()2t x y c ψ=+-=. 下面根据定理以及(3.4)式验证首次积分1ψ与2ψ的相互独立性,因为 1121222(,)2()0(,)x y D x y D x y x yψψψψψψ∂∂∂∂==--≠∂∂∂∂, 故首次积分1ψ与2ψ是相互独立的,因而原方程组的通解为 22122,1().2x y c t x y c ⎧-=⎪⎨+-=⎪⎩ 4 有关积分因子存在的充要条件的两类命题4.1 命题一对一阶微分方程(,)(,)0,M x y dx N x y dy += (,)x y D ∈(4.1.1) D 为单连通区域,其中(,)M x y 和(,)N x y 是关于,x y 的连续函数且具有连续的一阶偏导数,若存在只与x 有关的函数()P x ,使得M N N y x ⎛⎫∂∂- ⎪∂∂⎝⎭P(x)=成立, 则方程(4.1.1)存在只与x 有关的积分因子()()P x dx x e μμ⎰==,证 (方法一)对于方程(4.1.1),由题设条件知,若存在只与x 有关的积分因子()x μμ=,则 0yμ∂=∂. 这时方程 M N N M x y y x μμμ⎛⎫∂∂∂∂-=- ⎪∂∂∂∂⎝⎭, 变为 M N Nx y x μμ⎛⎫∂∂∂=- ⎪∂∂∂⎝⎭, 即 ()M N d y x dx P x dx Nμμ∂∂-∂∂==, (4.1.2) 进而得到 ()d P x dx μμ=, (4.1.3)这里()P x 仅为x 的函数,现在对方程(4.1.3)两边同时积分,可以求得方程(4.1.1)的一个积分因子 ()()P x dx x e μμ⎰==.(方法二) 由定理以及证明过程可知,设函数(,)x y μ为方程(4.1.1)的积分因子,那么它必满足偏微分方程(3.2)和特征方程 ln dx dy d M N N M y xμ==∂∂--∂∂, (4.1.4) 现在将条件M N N y x ⎛⎫∂∂- ⎪∂∂⎝⎭P(x)=变形为()*(,)M N P x N x y y x ∂∂-=∂∂并代入(4.1.4), 得 ln 1()dx d P x μ=, 对上式积分可得首次积分为ln ()p x dx μ=⎰.从而可得方程(4.1.1)的一个积分因子()()P x dx x e μμ⎰==.4.2 命题二对一阶微分方程(,)(,)0,M x y dx N x y dy += (,)x y D ∈(4.2.1) D 为单连通区域,其中(,)M x y 和(,)N x y 是关于,x y 的连续函数且具有连续的一阶偏导数,若存在只与y 有关的函数()Q y ,使得M N M y x ⎛⎫∂∂-- ⎪∂∂⎝⎭P(x)=成立,则方程(4.2.1)存在只与y 有关的积分因子()()Q y dy y e μμ⎰==.证 (方法一)对于方程(4.2.1),由题设条件知若存在只与y 有关的积分因子()y μμ=,则0x μ∂=∂, 这时方程 M N N M x y y x μμμ⎛⎫∂∂∂∂-=- ⎪∂∂∂∂⎝⎭变为 M N My y x μμ⎛⎫∂∂∂-=- ⎪∂∂∂⎝⎭, 即 ()M N d y x dy Q y dy Mμμ∂∂-∂∂==-, (4.2.2) 即 ()d Q y dy μμ=, (4.2.3)这里()Q y 仅为y 的函数,现在对方程(4.2.3)两边同时积分可以求得方程(4.2.1)的一个积分因子()()Q y dy y e μμ⎰==(方法二)由定理以及证明过程,现在设函数(,)x y μ为方程(4.2.1)的积分因子,那么它必满足偏微分方程(3.2)和特征方程 ln dx dy d M N N M y xμ==∂∂--∂∂, (4.2.4)现在将条件 ()M N y x Q y M∂∂-∂∂=-变形为 []()*(,)M N Q y M x y y x∂∂-=-∂∂并代入方程(4.2.4), 得 ln 1()dy d Q y μ=, 对上式积分可得首次积分为ln ()Q y dy μ=⎰,从而可得方程(4.2.1)的一个积分因子()()Q y dy y e μμ⎰==.例2 求解(2)0y y e dx x xy e dy -+=的通解.解 容易看出 (,)y M x y e =,(,)(2)y N x y x xy e =-+以及 y M e y ∂=∂, 4y N xy e x∂=--∂. 由于 2M N y x N x∂∂-∂∂=-, 即此方程存在只与x 有关的积分因子()x μ. 由命题一可知 2()21()dx x x e xμ-⎰==, 再用积分因子21()x x μ=乘原方程两端可得 220y ye e dx ydy dy x x--=, 即 20ye d dy x ⎛⎫--= ⎪⎝⎭, 于是,原方程的通解是 2ye y c x+=. 5 几类常见微分方程的积分因子5.1 可分离变量的微分方程设可分离变量方程为()()dy f x y dx ϕ=,其中()f x ,()y ϕ分别是,x y 的连续函数. 现在对上述方程做出变形得到更一般的形式1212()()()()0M x M y dx N x N y dy +=. (5.1.1)现在取211(,)()()x y M y N x μ=,同时乘到(5.1.1)式的两边,得到1212()()0()()M x N y dx dy N x M y +=. (5.1.2) 由于 1212()()0()()M x N y y N x x M y ⎛⎫⎛⎫∂∂== ⎪ ⎪∂∂⎝⎭⎝⎭,所以式(5.1.2)为恰当微分方程即全微分方程, 故211(,)()()x y M y N x μ=就是方程(5.1.1)的一个积分因子.5.2 齐次微分方程这里只考虑第一种形式的奇次微分方程.设齐次微分方程 ()dy yf dx x=, (5.2.1) 其中()f u 是u 的连续函数.作变量变换 yu x=,则 y ux =, (5.2.2)两边关于x 求导可得 dy du x u dx dx =+. (5.2.3) 现在将(5.2.1)式和(5.2.2)式代入(5.2.3)式可得到如下可分离变量的微分方程()duxu f u dx+= 或 []()0xdu u f u dx +-=, (5.2.4) 所以由前文对分离变量微分方程的讨论可知,微分方程(5.2.4)的积分因子是[]1(,)()x y x u f u μ=-. 现将[]1(,)()x y x u f u μ=-同时乘到(5.2.4)的两边得到 []110()du dx u f u x +=-, (5.2.5)由于 []110()u x x u f u ⎛⎫∂∂⎛⎫== ⎪ ⎪ ⎪∂∂-⎝⎭⎝⎭,所以(5.2.5)是全微分方程,也即 1(,)()x y y y xf xμ=- 是方程(5.2.1)的积分因子.5.3 线性微分方程设一阶线性微分方程()()dyP x y Q x dx=+, (5.3.1) 其中()P x ,()Q x 是关于x 的连续函数. 对上式进行简单变形可得如下对称形式 []()()0P x y Q x dx dy ++=, 其中 (,)()()M x y P x y Q x =+,(,)1N x y =.又由于 ()M N N P x y x ⎛⎫∂∂-= ⎪∂∂⎝⎭ (这里()P x 只是关于x 的函数),所以由前文的命题一结论可知此一阶线性微分方程有只与x 有关的积分因子()x μ,即 0y μ∂=∂,d x dxμμ∂=∂. (5.3.2)将(5.3.2)式代入(2.2)式,得到 ()M Nd y xdx Nμμ∂∂-∂∂= 即()d P x dx μμ=,对上式两边同时积分可得()()P x dx x e μ⎰=,即一阶线性微分方程(5.3.1)的积分因子是()()P x dx x e μ⎰=.5.4 Bornoulli 微分方程设Bornoulli 微分方程()()n dyP x y Q x y dx=+ (0,1)n ≠, (5.4.1) 其中()P x ,()Q x 为x 的连续函数.现在把(5.4.1)式变形为 ()()0n P x y Q x y dx dy ⎡⎤+-=⎣⎦,其中(,)()()n M x y P x y Q x y =+,(,)1N x y =-,并且 1()()n M P x ny Q x y -∂=+∂,0Nx∂=∂,于是可以得到它的特征方程为11(()())()()n n dx dy dtP x y Q x y P x ny Q x -==--++. (5.4.2) 对(5.4.2)式做一个变形可得1()()()()n n dx ndy ydtnP x y nQ x y P x y ny Q x -==++. (5.4.5) 利用分式的性质可得1(1)()dx ndy ydt n P x y+=-, 即 (1)()n P x ydx ndy ydt -=+, 1(1)()dt n P x dx ny dy -=--, 对上式积分可得首次积分为 (1)()ln t n P x dx n y =--⎰, 再由前文定理可知 (1)()ln ,(,)n P x dx n t x y y e μμ--⎰==.故Bornoulli 微分方程的积分因子是(1)()(,)n P x dx n x y y e μ--⎰=.例3 把方程32()()0x xy x y dx x y dy +++--=化为全微分方程. 解 由题目很容易知道32(,)M x y x xy x y =+++, (,)()N x y x y =--, 以及21M xy y ∂=+∂, 1Nx∂=-∂, 从而可以写出此方程的特征方程为32()()22dx dy dtx y x xy x y xy ==---++++,对上式进行简单变形可得233222(1)xdx ydy dtx xy x y xy xy y xy ==-+++-+,利用分式性质可得2222()(1)(1)xdx ydy dtx y xy xy +=++-+,即 2222()1d x y dtx y +=+-,对上式两边同时积分可得首次积分为 22ln()t x y =-+. 再由定理可知 ln t μ=,故221(,)x y x y μ=+,故所要化简方程的一个积分因子为221(,)x y x y μ=+,于是原方程可以化为一个全微分方程,即3222220x xy x y x ydx dy x y x y+++--=++, 其中 3222222222()x xy x y x y x y xyy x x y x y ⎛⎫⎛⎫∂+++∂---=-= ⎪ ⎪∂∂++⎝⎭⎝⎭. 6 小结通过本文使我们对微分方程的积分因子与偏微分方程有了一个初步的了解,尤其对偏微分方程中的特征方程,首次积分等概念有了一个简单的认识.其次就偏微分方程中的特征方程,首次积分在求解常微分方程组的积分因子中的应用有了进一步的认识,最后应用本文中的理论求出了几类常见的微分方程的积分因子,为今后能更好的研究与积分因子相关的理论提供帮助.参考文献[1] 王高雄,周之铭,朱思铭.常微分方程[M].第三版.北京:高等教育出版社,2006:50-60.[2] 李刚升.浅谈积分因子与偏微分方程[J].商丘科技职业技术学院,2004,(02).[3] 张奕河,郭文川.关于一阶常微分方程的积分因子的求解问题[J].四川理工学院学报,2009.22(6):11-13.[4] 徐安农,段复建.全微分方程与积分因子法[J].桂林电子工业学院学报,2002,22(2):11-12.[5] 李耀红,张海燕.几类微分方程积分因子存在定理[J].巢湖学院学报,2002,8(3):8-10.[6] 李伟鹏.求首次积分的几种方法[J].陇东学院学报.2007,17(1)22-24.致谢本论文的完成是在何万生老师的悉心指导下进行的,在何老师的指导下,我的各方面能力都有所提高,尤其是何老师渊博的知识,敏锐的学术思维,精益求精的工作态度以及诲人不倦的师者风范是我终生学习的楷模.何老师严谨求实和一丝不苟的治学态度和勤勉的工作态度也将永远鼓励我,至此,谨向何老师致以衷心的感谢和崇高的敬意.同时感谢所有教育过我的专业老师,你们传授的专业知识是我不断成长的源泉,也是我完成论文的基础.在此也感谢我同一组的组员和班里的同学,是你们在我遇到难题时帮我找大量相关资料,解决难题.再次真诚感谢所有帮过我的老师同学,另外,由于自己知识积累、能力有限和经验匮乏,论文中难免有许多考虑不周全的地方,希望各位老师多加指导.最后,我要向百忙之中抽出时间对本文进行审阅,评议和参与本人论文答辩的各位老师表示感谢.。

常微分方程的积分因子法

常微分方程的积分因子法

常微分方程的积分因子法在数学中,常微分方程是一种描述动态系统的重要工具。

在实际应用中,常微分方程模型广泛应用于物理、化学、生物学等领域,用于研究自然界中各种现象的演化规律。

常微分方程除了数值解法外,还有一种有力的解法——积分因子法。

积分因子法是通过引入一个特殊的乘数,将常微分方程转化为可积分的形式,从而求出它的通解。

1. 常微分方程与积分因子首先,我们需要了解什么是常微分方程。

简单来说,常微分方程是描述一个未知函数与其导数之间关系的方程。

比如,一阶常微分方程可以写成:$$\frac{dy}{dx}=f(x,y)$$其中,$y=y(x)$ 是未知函数,$f(x,y)$ 是已知函数。

解此方程的一般方法是使用分离变量法或者变量代换法,但是有些方程并不方便通过这些方法求解。

这时候,就需要借助积分因子法。

积分因子法是常微分方程中的一种特殊解法,通过引入一个特殊的函数,将原方程乘上这个函数,使它变为可积分的形式。

其必要条件是,乘上这个函数后,原方程满足以下形式:$$\mu(x,y,z)\frac{\partial f(x,y,z)}{\partialx}+\mu(x,y,z)\frac{\partial g(x,y,z)}{\partialy}+\mu(x,y,z)\frac{\partial h(x,y,z)}{\partialz}+\mu(x,y,z)P(x,y,z)=0$$其中,$\mu(x,y,z)$ 是引入的积分因子。

这时,我们可以通过将这个新方程改写成完全微分形式来求解,从而得到原方程的通解。

2. 积分因子法的应用举例下面,我们来看一个实际的例子,说明积分因子法的应用。

考虑以下常微分方程:$$\frac{dy}{dx}+2y=xe^{-x}$$这是一个一阶线性非齐次微分方程,我们可以使用常见的解法——待定系数法或变量分离法,但这里我们要演示积分因子法的应用。

首先,我们需要找到这个方程的积分因子。

一阶微分方程积分因子探讨

一阶微分方程积分因子探讨

一阶微分方程积分因子的求法探讨数学与信息科学学院 数学与应用数学专业 指导教师:郑丽丽 职称:教授摘 要:针对满足某些条件的微分方程,本文将给出几种直接、有效地求积分因子的方法.关键词:一阶微分方程;积分因子The Solution of Integral Factor for the First Order OrdinaryDifferential EquationAbstract :This paper has made a special effort to study how to quadrate integral factors directly and efficiently .When the differential equations meet some conditions , therefore , the common method we can get from it .Key Words :the first order ordinary differential equation ;integral factor0前 言一阶微分方程的求解是整个微分方程求解的基础,一般的有两种处理方式:一是以变量可分离的方程为基础,通过适当的变量代换把一阶微分方程化为可积型方程;另外就是以全微分方程为基础,采取积分因子法把一个一阶微分方程化为全微分方程求.这里我们讨论了积分因子存在的充要条件,给出了确定若干特殊类型的积分因子的求法.1 积分因子的定义若对于一阶微分方程()(),,0M x y dx N x y dy += (1)其中(),M x y ,(),N x y 在矩形域内是,x y 的连续函数,且有连续的一阶偏导数.若存在连续可微的函数(),0x y μ≠,使得()()()(),,,,0x y M x y dx x y N x y dy μμ+≡,为一恰当方程,即存在函数V ,使M dx N dy dV μμ+=.则称(),x y μ为方程(1)的积分因子. 通过计算可得,函数(),x y μ为0M dx N dy +=积分因子的充要条件为:()()M N xyμμ∂∂=∂∂,即M N N Mxy yx μμμ⎛⎫∂∂∂∂-=- ⎪∂∂∂∂⎝⎭(2) 这是个以μ为未知数的一阶线性偏微分方程,要想通过解方程()2来求积分因子通常很困难,但在若干特殊情况下,求积分因子还是容易的,下面总结了几种可以方便求出特殊类型的积分因子的方法.2 积分因子存在的充要条件定理1[5] 方程()(),,0M x y dx N x y dy +=具有形如(),x y μμφ=⎡⎤⎣⎦的积分因子的充要条件为:()1,M N M N f x y y x x y φφφ-⎛⎫⎛⎫∂∂∂∂--=⎡⎤ ⎪ ⎪⎣⎦∂∂∂∂⎝⎭⎝⎭.证明 因为()(),,0M x y dx N x y dy +=有积分因子的充要条件为M N N Mxy y x μμμ⎛⎫∂∂∂∂-=- ⎪∂∂∂∂⎝⎭.令(),x y μμφ=⎡⎤⎣⎦,则有()d d M N N Md xd y y x μφμφμφφφ⎛⎫∂∂∂∂-=- ⎪∂∂∂∂⎝⎭, 即()1,d M N N M f x y yx x y μφφφμ-⎛⎫⎛⎫∂∂∂∂=--=⎡⎤ ⎪ ⎪⎣⎦∂∂∂∂⎝⎭⎝⎭.并由此得出其积分因子为()(),f d x y e φφμ⎰=.根据这个定理可以得出以下特殊类型的积分因子的充要条件. 2.1 具有()x μμ=形式的积分因子[1]方程0M dx N dy +=具有特殊因子()x μμ=的充要条件为()M N yxx N ϕ∂∂-∂∂=,这里()x ϕ仅为x 的函数.于是积分因子为()x d xe ϕμ⎰=.例1[2] 求()20y x dx xdy --=的积分因子.解 因为2M y x =-,N x =-,且1M y∂=∂,1N x∂=-∂,则2MN yxN x∂∂-∂∂=-,于是积分因子为22dxx exμ--⎰==.2.2 具有()y μμ=形式的积分因子[1]方程0M dx N dy +=具有特殊因子()y μμ=的充要条件为()M N yxy Mψ∂∂-∂∂=-,这里()f y 仅为y 的函数.于是积分因子为()y dye ψμ⎰=.例2[5]求()()cos sin sin cos 0y x x x dx y x x x dy -++=的积分因子.解 因为cos sin M y x x x =-,sin cos N y x x x =+,且1MN yxM∂∂-∂∂=-,于是积分因子为(),dyy x y e e μ⎰==.2.3 具有()x y μμ=±形式的积分因子[8]方程0M dx N dy +=具有特殊因子()x y μμ=±的充要条件为()()1M N M N fx y y x -⎛⎫∂∂-±=± ⎪∂∂⎝⎭.例3[3] 求方程()()3223322223230x x y y y dx y xy x x ++-+++-=的积分因子. 解 因为322323M x x y y y=++-, 322223N y xy x x =++-,且()12M N N M yx x y-⎛⎫∂∂--=-⎪∂∂+⎝⎭,只与x y +有关,于是有积分因子()()22,d x y x y x y ex yμ-++⎰==+.2.4 具有()22x y μμ=±形式的积分因子[8]方程0M dx N dy +=具有特殊因子()22x y μμ=±的充要条件为()()122M N N x M y f xyy x -⎛⎫∂∂-±=± ⎪∂∂⎝⎭.例4[3] 求方程()220x y y dx xdy ++-=的积分因子. 解 因为22M x y y =++, N x =-,且()1221M N N x M y yx x y-⎛⎫∂∂--=-⎪∂∂+⎝⎭,于是积分因子为()22221221d xyx y ex yμ-++⎰==+.推广[7] 方程0M dx N dy +=具有特殊因子()x y αβμμ=+的充要条件是:()()111M N x N y M fxyyx αβαβαβ---⎛⎫∂∂--=+ ⎪∂∂⎝⎭.2.5 具有()x y αβμμ=形式的积分因子方程0M dx N dy +=具有特殊因子()x y αβμμ=的充要条件为()11M N N M fxyx yy x x y αβαβαβ-⎛⎫⎛⎫∂∂--= ⎪ ⎪∂∂⎝⎭⎝⎭.由此又可分为二种类型:()1 方程0M dx N dy +=具有特殊因子()xy μμ=的充要条件为()11M N N x M y yx xy -⎛⎫∂∂--=- ⎪∂∂⎝⎭;()2 方程0M dx N dy+=具有特殊因子x y μμ⎛⎫= ⎪⎝⎭的充要条件为12M N M N y y f yx x x x -⎛⎫∂∂⎛⎫⎛⎫--+= ⎪ ⎪⎪∂∂⎝⎭⎝⎭⎝⎭.例5[4] 求方程()()22324430y x y dx xy x dy +++=的积分因子. 解 设积分因子为p q x y ,于是有()()2232443p q p q x y y x y x y xy x y x ∂∂⎡⎤⎡⎤+=+⎣⎦⎣⎦∂∂,或写成()()()()121222414133pqpq p qpqq x yq xy p x yp xy+++++++=+++.上式对任意x 和y 都满足时,必须有()()2241q p +=+,()()4133q p +=+,解之得1p =,2q =.于是有积分因子2xy μ=.注 此种类型中α,β的确定可用待定系数法.以上所讨论的是微分方程具有特殊因子的求法.而有些方程具有特殊结构,我们可根据其特殊结构求出其积分因子.3 特殊结构方程的积分因子[6]定理2 方程()()()()0M x N y dx P x Q y dy +=有积分因子:()()1Ny P x μ=.定理3 如果0xM yN +≠,而M 和N 皆为m 次齐次函数,则方程0M dx N dy +=有积分因子:1xM yNμ=+.4 分组求积分因子法[9]对于一些复杂的方程,往往不容易直接求出它们的积分因子,这是可以把它的左边分组,分别求出各组的积分因子,然后再求总的式子的积分因子.例如分成两组:()()11220M dx N dy M dx N dy +++=(3) 可分别求出各组的积分因子1μ和2μ,也就是如果有1u ,2u 使:11111M dx N dy du μμ+=,22222M dx N dy du μμ+=.于是借助1μ,2μ常可求得0M dx N dy +=得积分因子.定理4[4] 如果μ是0M dx N dy +=的一个积分因子,且M dx N dy du μμ+=,则()u μϕ也是0M d x N d y +=的积分因子.此处()u ϕ是u 的任一连续函数.而()()()()()()u Mdx u Ndy u Mdx Ndy u du d u μϕμϕϕμμϕφ+=+==,其中φ是ϕ的一个原函数.据此知,对任意的函数()u ϕ,()u ψ,()11u μϕ及()22u μψ都分别是()3的第一组和第二组的积分因子.函数ϕ、ψ有着广泛选择的可能性,若能选择ϕ、ψ使:()()1122u u μμϕμψ==,则μ就既是()3的第一组也是第二组的积分因子.因而也就是0M dx N dy +=的积分因子.例6[9] 求方程()32420x y y dx x dy -+=的积分因子. 解 原方程改写为()34220x ydx x dy y dx +-=,显然131xμ=,1u xy =,221yμ=,2u x=.为使()()123211g xy g x xy=,只需取()()121g xy xy =,()251g x x=.于是求的原方程的一个积分因子:521x yμ=.综上所述,该文介绍一些特殊类型的积分因子的求法及部分特殊结构的微分方程的积分因子的求法,只要掌握这几种方法,就能很容易的解出一些方程的积分因子,将大大提高解微分方程的效率和可操作性.参考文献:[1] 焦宝聪,王在洪,时红廷.常微分方程[M].北京:清华大学出版社,2008.[2] 孙清华,李金兰.常微分方程内容、方法和技巧[M].武汉:华中科技大学出版社,2006.[3] 钱伟长.常微分方程的理论及其解法[M].北京:国防工业出版社.1992.[4] 丁同仁,李承浩.常微分方程教程[M].北京:高等教育出版社.1991.[5] 王高雄,周之铭,王寿松.常微分方程[M].北京:高等教育出版社.2006.[6] 潘鹤鸣.几种特殊类型积分因子的求法及在解微分方程中的应用[J].巢湖学院学报,2003(3):18-22.[7] 李德新.两类特殊微分方程的积分因子解法[J].福建农林大学学报,2004,33(2):269-271.[8] 李君士.积分因子的求法[J].九江师专学报:自然科学版,1989,8(2):64-68[9] 吴淼生.关于非恰当方程0M dx N dy+=积分因子的求法[J].宜春师专学报,1994,2(2):15-23.。

2.2-线性微分方程(积分因子法)

2.2-线性微分方程(积分因子法)

s
x Q(s)ex P(t)dt ds
x0
(3)
二 伯努利(Bernoulli )方程
形如 dy p(x) y Q(x) yn n 0,1是常数 dx
的方程,称为伯努利方程. 这里P(x), Q(x)为x的连续函数 。
解法: 10 引入变量变换 z y1n ,方程变为
dz (1 n)P(x)z (1 n)Q(x) dx 20 求以上线性方程的通解
§2.2 线性微分方程与积分因子法
一阶线性微分方程的一般形式为
a(x) dy b(x) y c(x) dx
在a(x) 0的区间上可写成 dy P(x) y Q(x) (1)L L 标准形式 dx
这里假设P(x), Q(x)在考虑的区间上是 x的连续函数 若Q(x) 0,则(1)变为
dy P(x) y 0 (2)L L 齐次线性方程 dx
若Q(x) 0,则(1)称为非齐次线性方程。
一 一阶线性微分方程的解法-----积分因子法
dy P(x) y Q(x) (1) dx
求解思想:方程两边乘一个函数,使得左边变成 一个函数的导数
y e p(x)dx ( Q(x)e p(x)dxdx c)
(3)
注: (i) 求(1)的通解可直接用公式(3) (ii) 课本用的是常数变易法,方程整理的形式不同
解: 原方程不是未知函数 y的线性方程 ,但将它改写为
dx 2x y2
dy y

dx 2 x ,
故其通解为 x e p( y)dy ( Q( y)e p( y)dydy c)
e
2 y
dy
(
(
y)e
2 y
dy
dy
c)

一阶微分方程积分因子的探讨

一阶微分方程积分因子的探讨

一阶微分方程积分因子的探讨一阶微分方程是常见的数学问题,其中积分因子是一种常用的解法之一。

积分因子是指一个函数,可以乘到微分方程的两边,使其变成可积的形式。

本文将讨论一阶微分方程中积分因子的概念、求法以及应用。

一、积分因子的概念对于一阶微分方程dy/dx + P(x)y = Q(x),其中P(x)和Q(x)是已知函数,如果存在一个函数μ(x),使得μ(x)dy/dx + μ(x)P(x)y = μ(x)Q(x)成立,则称μ(x)为该微分方程的积分因子。

二、积分因子的求法求解积分因子的方法有许多种,以下介绍两种常用的方法。

1.利用P(x)的特殊形式求积分因子当P(x)具有一定的特殊形式时,可以通过变形求得积分因子。

例如,当P(x)是一个常数时,我们可以将μ(x)设为e^(kx),其中k 为待定常数,代入μ(x)dy/dx + μ(x)P(x)y = μ(x)Q(x)中化简,得到:e^(kx)dy/dx + ke^(kx)y = e^(kx)Q(x)将该式两边同时乘以e^(kx),得到:e^(2kx)dy/dx + ke^(kx)dy/dx = e^(kx)Q(x)e^(kx) 对于左边的式子,我们可以发现它是一个乘积的求导形式,因此可以利用乘积法则进行求导,得到:d/dx(e^(kx)y) = e^(kx)Q(x)e^(kx)将该式两边同时积分,得到:e^(kx)y = ∫e^(kx)Q(x)e^(kx)dx + C因此,积分因子μ(x) = e^(kx) = e^(∫P(x)dx)。

2.利用常数变易法求积分因子常数变易法是一种常用的求解积分因子的方法。

具体步骤如下:(1)将微分方程写成dy/dx + P(x)y = Q(x)的形式。

(2)设积分因子为μ(x),则将μ(x)dy/dx + μ(x)P(x)y = μ(x)Q(x)两边同时积分,得到:μ(x)y = ∫μ(x)Q(x)dx + C其中C为常数。

一阶常系数微分方程的积分因子研究

一阶常系数微分方程的积分因子研究

文章编号:1673G2103(2018)05G0006G05一阶常系数微分方程的积分因子研究∗亓洪胜(蚌埠学院理学院,安徽蚌埠233030)摘㊀要:在现代电力和工程计算实例应用中,很多数学模型均是以一阶常系数微分方程的形式构建起来的.高于一阶常系数微分方程的初值或是边值问题均可等效转换为一阶常微分方程组问题,研究一阶常系数微分方程组的数值求解具有特殊的基础性意义.引入积分因子的概念,提出一阶常系数微分方程求解的基本方法.推导出四种常见的一阶常系数微分方程的积分因子的一般形式,并给出不同形式积分因子存在的充要条件,使求解过程更简单,更清晰.该方法可用于求解当前常规方法无法求解的微分方程.关键词:一阶常系数;微分方程;积分因子;充要条件中图分类号:O 241.81㊀文献标识码:A引言一阶常系数微分方程是运动学的重要组成部分,主要反映现实世界运动过程中量与量之间存在的各种数学模型,这些模型可满足大量一阶常系数微分法方程关系式的求解,通过对数学模型的求解可获取函数性质[1].一阶常系数微分方程的应用十分广泛,可以解决各种与导数相关的问题,如物理中的动力学和运动学等问题,都可以使用微分方程求解[2~3],积分因子法是解一阶常系数微分方程的有效的方法,它从一阶常系数微分方程的原理入手对方程进行求解[4~5].文献[6]利用M n c h 不动点定理对微分方程进行求解.文献[7]利用S c h a u d e r 不动点定理方法对微分方程进行求解,文献[8]采用再生核理论以及配置方法对微分方程进行求解,上述方法得到的都是一阶常系数微分方程的通解,但是很多特殊类型的微分方程往往无法给出解的具体表达式.1㊀一阶常系数微分方程的积分因子研究1.1㊀微分方程的积分因子的充要条件通过对积分方程进行求解,得到其该方程的积分充要条件.将一阶方程表示为:d xd y=f (x ,y )(1)用微分的形式表示为:f (x ,y )d x -d y =0㊀(2)假设将x ,y 看作是平等的,则对称的微分方程为:M (x ,y )d x +N (x ,y )d y =0(3)假设方程的左端是任意一个二次函数u (x ,y )的全微分,则有:6第40卷第5期V o l .40㊀N o .5㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀菏泽学院学报J o u r n a l o fH e z eU n i v e r s i t y㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀2018年10月O c t .㊀2018∗收稿日期:2018G02G10基金项目:安徽省自然科学基金项目(1808085MA 02),省级教学质量工程项目(2015j y x m 386,2014z y 141),蚌埠学院教研项目(2017J Y X M L 16,2017jx t d 2);蚌埠学院科研项目(2011Z R 09,2017Z R 10z d )作者简介:亓洪胜(1978-),男,山东莱芜人,硕士研究生,讲师,研究方向:微分方程数值解.M (x ,y )d x +N (x ,y )d y =d u (x ,y )=∂u ∂x +∂u∂x(4)即方程(4)为恰当方程.为了方便进行验证,将方程的通解表示为:u (x ,y )+c (5)上式中,c 为任意常数.假设方程是恰当方程,则有:∂u∂x=M (6)∂u∂y=N (7)利用式(6)与式(7)分别对x ㊁y 进行偏导数求解,得:∂2u ∂y ∂x =∂M ∂y ,∂2u ∂x ∂y =∂N∂x (8)上式中,恰当方程的必要条件为∂M ∂y=∂N ∂x .假设存在连续可微的函数μ=μ(x ,y )ʂ0,即:μ(x ,y )M (x ,y )d x +μ(x ,y )d y =0(9)假设式(9)为恰当方程,若存在函数v ,则:μM d x +μM d y =d v (10)若满足上式条件,则μ(x ,y )为该方程的积分因子[9~11].则该方程的积分因子的充要条件为:∂(μM )∂y=∂(μN )∂x (11)积分因子的形式多种多样[12~13],以下给出积分因子存在的充要条件:1)假设方程只与x 有关,则积分因子的充要条件为:∂M ∂y -∂N ∂x N=φ(x )(12)求得该方程的一个积分因子为:μ(x )=e ʏφ(x )d x (13)2)假设方程与y 有关,则积分因子的充要条件为:∂M ∂y -∂N∂x -M=φ(y )(14)求得该方程的一个积分因子为:μ(y )=e ʏφ(y )d y (15)3)如果方程为M (x ,y )d x +N (x ,y )d y =0,则u (x +y )的积分因子充要条件为:1N -M d M d y -d N d x æèçöø÷=f (x +y )(16)求得该方程的一个积分因子为:u =e x p(ʏ(x +y )d (x +y ))(17)4)如果方程为M (x ,y )d x +N (x ,y )d y =0,则u (x -y )的积分因子充要条件为:1N -M d M d y-d N d x æèçöø÷=f (x -y )(18)72018年㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀亓洪胜:一阶常系数微分方程的积分因子研究㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀第5期求得该方程的一个积分因子为:u =e x p(ʏ(x -y )d (x -y ))(19)1.2㊀一阶常系数微分方程的积分因子求解法在一阶常系数微分方程中,利用下式给出总微分方程的对称形式[14~15]为:P (x ,y )d x +Q (x ,y )d y =0(20)将一阶常系数微分方程的补充条件设为:∂p ∂y =∂Q ∂x 时,可将公式(20)进行转换为:μ(x ,y )P (x ,y )d x +μ(x ,y )Q (x ,y )d y =0(21)以公式(21)为基础,根据判别准则,将μ(x ,y )作为公式(22)的补充条件,则:p∂y =Q ∂x (22)1)第一种一阶常系数微分方程积分因子求解法假设,一阶常系数线性方程为:d y d x+p (x )y =f (x )(23)以公式(23)为基础,将一阶常系数线性方程改写成方程对称形式:p (x )y -f (x )[]d x +d y =0(24)利用公式(22)积分因子的补充条件,可将公式(24)转换为具有和x 有关的特殊性质积分因子μ(x ),μ(x )为:μ(x )=e f (x )d x ∂x (25)将公式(25)代入公式(22)中,可得一阶常系数线性方程两边积分:d μμ=p (x )d x(26)以公式(26)为基础,给出一阶常系数线性方程的积分因子:μ(x )=e ʏp (x )d x (27)例1㊀对d y d x +y x =s i n x x的积分因子求解解㊀一阶常系数微分方程积分因子μ(x )=e ʏp (x )d x =e l n x =x ,利用x 乘一阶常系数微分方程两端可得:x d y +yd x -s i n x d x =0(28)在一阶常系数微分方程两边积分,则解为:x y +c o s x =C (29)2)第二种一阶常系数微分方程积分因子求解法根据一阶常系数微分方程积分因子可以获取出贝努利的积分因子,利用下式给出方程:d y d x+p (x )y =f (x )yn(30)将公式(30)两边积分同时乘y -n ,当z =y -n 得:d z d x+(1-n )p (x )z =(1-n )f (x )(31)以公式(31)为基础,求出贝努利积分因子为:μ(x )=e (1-n )ʏp (x )d x(32)例2㊀对d y d x+y =(c o s x -s i n x )y 2的积分因子进行求解82018年㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀菏㊀泽㊀学㊀院㊀学㊀报㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀第5期解㊀贝努利积分因子μ(x ,y )=y -n e (1-n )ʏp (x )d x =y -2e -ʏd x ,微分方程两边同时乘以y -2e -x ,可得:y 2e -x d y +y -1e -x dx =(c o s x -s i n x )e -x d x (33)则贝努利方程两边积分解为:e -x s i n x +e -x y -1=C(34)3)第三种一阶常系数微分方程积分因子求解法利用观察法对可分离变量方程进行分析,得如下方程:p 1(x )q 1(y )d x +p 2(x )q 2(y )d y =0(35)将公式(35)两边同时乘以1p 2(x )q 1(y )得:p 1(x )p 2(x )d x +q 2(y )q 1(y )d y =0(36)公式(36)为全微分方程,利用下式给出,可分离变量公式(35)的积分因子为:μ(x ,y )=1p 2(x )q 1(y )(37)例3㊀对(x y -x )d x +(x y +x -y -1)d y =0的积分因子进行求解解㊀将微分方程进行变换为x (y -1)d x +(x -1)(y +1)d y =0,积分因子为:μ(x ,y )=1p 2(x )q 1(y )=1(x -1)(y -1),将变换后的微分方程两边同时乘以上述积分因子,可得:x x -1d x +y +1y -1d y =0(38)以公式(38)为基础,微分方程两边积分可解为:x +y +l n (x -1)(y -1)2=C (39)4)第四种一阶常系数微分方程积分因子求解法根据上述方程的积分因子导出齐次方程的积分因子,利用下式给出齐次方程:p (x ,y )d x +Q (x ,y )d y =0(40)将方程两边同时乘以1Q (x ,y),利用下式给出微分方程的对称式为:u -f (u )[]d x +x d u =0(41)将u =y x 乘以1Q (x ,y)代入到公式(41)中,求出积分因子为:μ(x ,y )=1x p (x ,y )+y q (x ,y )(42)例4㊀对(1+2e x y )d x +2e xy (1-x y )d y =0的积分因子进行求解解㊀将μ(x ,y )=1x p (x ,y )+y q (x ,y)这个微分方程两边同时乘以积分因子,可得:1+2e xy x +2y e x yd x +2e xy (1-x y )x +2y e xyd y =0(43)则方程两边积分解为:x +2y e xy =0(44)对于一阶常系数微分方程,可以利用上述四种求解方法进行积分因子求解.2㊀结束语本文引入积分因子的概念,阐述了一阶常系数微分方程的积分因子方法,提出了一阶常系数微分方程的92018年㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀亓洪胜:一阶常系数微分方程的积分因子研究㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀第5期2018年㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀菏㊀泽㊀学㊀院㊀学㊀报㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀第5期求解基本方法.(1)给出一阶常系数微分方程不同形式积分因子存在的充要条件.(2)推导出四种常见的一阶常系数微分方程的积分因子求解方法.该方法可用于求解当前常规方法无法求解的微分方程,使求解过程更简单,更清晰.参考文献:[1]何文正,陈晨,文竞舟.预应力圆弧曲梁振动微分方程推导及求解[J].人民长江,2016,47(17):88G92.[2]盛钰婷.三阶常微分方程的某些非线性特征值问题的正解[J].科学中国人,2016,23(5):513G519.[3]林文贤.一类具阻尼项的偶阶中立型微分方程的振动性分析[J].吉林大学学报:理学版,2017,55(5):1073G1076.[4]张滑,刘春凤.分数阶RGL微分方程初值问题解的存在唯一性[J].数学的实践与认识,2016,46(3):267G272.[5]蔡钢.B a n a c h空间上有限时滞退化微分方程的适定性[J].数学物理学报,2018,38(2):264G275.[6]王信峰,陈玉花,张莉.B a n a c h空间非线性二阶混合型脉冲积分G微分方程的解[J].数学的实践与认识,2017,47(4):266G272.[7]姚慧丽,刘婷.一类随机泛函积分微分方程的pG期望伪概自守温和解[J].哈尔滨理工大学学报,2017,22(4):135G142.[8]刘宏影,吕学琴.再生核结合配置法求解一类带有积分边值条件的四阶非线性微分方程[J].数学物理学报,2016,36(5):928G936.[9]王玥.复高阶微分方程的解[J].数学学报(中文版),2017,60(4):651G660.[10]曹玉林,马建萍.基于微分方程的移动自组网病毒传播模型研究[J].计算机工程,2017,43(1):172G177.[11]姚慧丽,刘婷.一类随机泛函积分微分方程的pG期望伪概自守温和解[J].哈尔滨理工大学学报,2017,22(4):135G142.[12]魏会贤,董文雷,郭彦平.一类时标上三阶微分方程的渐进性态[J].数学的实践与认识,2017,47(11):307G312.[13]徐鸿鹏,尹社会,张勇.一个超混沌类L o r e n z系统的非线性动力学行为及计算机仿真[J].电子设计工程,2016,24(10):38G41.[14]王平,范磊,梁铃,等.基于改进微分进化算法的发电机调速系统参数辨识[J].自动化与仪器仪表,2018(1):47G50.[15]朱鹏程,程浪,蔡栩.基于伪微分算法的电液稳定平台研究[J].自动化与仪器仪表,2017(10):29G32.O n t h e I n t e g r a l F a c t o r s o f F i r s tGo r d e rD i f f e r e n t i a l E q u a t i o n sw i t hC o n s t a n t C o e f f i c i e n t sQ IH o n gGs h e n g(C o l l e g e o f S c i e n c e,B e n g b uU n i v e r s i t y,B e n g b uA n h u i233030,C h i n a)A b s t r a c t:I nt h ea p p l i c a t i o no f m o d e r ne l e c t r i c p o w e ra n de n g i n e e r i n g c a l c u l a t i o ne x a m p l e s,m a n y m a t h e m a t i c a lm o d e l s a r e c o n s t r u c t e d i nt h e f o r mo f f i r s tGo r d e r c o n s t a n t c o e f f i c i e n td i f f e r e n t i a l e q u a t i o n s.T h e i n i t i a l v a l u e o r b o u n d a r y v a l u e p r o b l e m s o f h i g h e r t h a n t h e f i r s tGo r d e r d i f f e r e n t i a l e q u a t i o n sw i t h c o nGs t a n t c o e f f i c i e n t s c a nb e c o n v e r t e d i n t o t h e f i r s tGo r d e r o r d i n a r y d i f f e r e n t i a l e q u a t i o n s.I t i s o f s p e c i a l f u n d aGm e n t a l s i g n i f i c a n c e t o s t u d y t h e n u m e r i c a l s o l u t i o n o f f i r s t o r d e r d i f f e r e n t i a l e q u a t i o n sw i t h c o n s t a n t c o e f f iGc i e n t s.B y i n t r o d u c i n g t h e c o n c e p t o f i n t e g r a l f a c t o r,t h i s p a p e r p r o p o s e s ab a s i cm e t h o d f o r s o l v i n g f i r s tGo r d e r d i f f e r e n t i a l e q u a t i o n sw i t hc o n s t a n t c o e f f i c i e n t s.T h e g e n e r a l f o r m s o f i n t e g r a l f a c t o r s f o r f o u r c o mGm o n f i r s tGo r d e r d i f f e r e n t i a l e q u a t i o n sw i t hc o n s t a n t c o e f f i c i e n t s a r ed e d u c e d,a n d t h en e c e s s a r y a n ds u f f iGc i e n t c o n d i t i o n s f o r t h e e x i s t e n c e o f d i f f e r e n t f o r m s o f i n t e g r a l f a c t o r s a r e g i v e n.T h i sm e t h o d c a nb eu s e d t o s o l v e t h e d i f f e r e n t i a l e q u a t i o n sw h i c hc a n n o t b e d o n eb y c o n v e n t i o n a lm e t h o d s.K e y w o r d s:f i r s tGo r d e r c o n s t a n t c o e f f i c i e n t;d i f f e r e n t i a l e q u a t i o n s;i n t e g r a l f a c t o r;n e c e s s a r y a n d s u f f iGc i e n t c o n d i t i o n01。

微分方程的积分因子求解法

微分方程的积分因子求解法

常微分方程的积分因子求解法内容摘要:本文给出了几类特殊形式的积分因子的求解方法,并推广到较一般的形式。

关键词: 全微分方程,积分因子。

一、 基本知识定义1.1 对于形如0),(),(=+dy y x N dx y x M (1.1)的微分方程,如果方程的左端恰是x ,y 的一个可微函数),(y x U 的全微分,即d ),(y x U = dy y x N dx y x M ),(),(+,则称(1。

1)为全微分方程.易知,上述全微分方程的通解为 ),(y x U =C , (C 为任意常数).定理1。

1 (全微分方程的判别法)设),(y x M ,),(y x N 在x ,y 平面上的单连通区域G 内具有连续的一阶偏导数,则(1.1)是全微分方程的充要条件为xy x N y y x M ∂∂=∂∂),(),( (1。

2) 证明见参考文献[1].定义1。

2 对于微分方程(1。

1),如果存在可微函数),(y x μ,使得方程),(y x μ0),(),(),(=+dy y x N y x dx y x M μ (1.3)是全微分方程,则称),(y x μ为微分方程(1。

1)的积分因子。

定理1。

2 可微函数),(y x μ为微分方程(1.1)的积分因子的充要条件为x y x y x N ∂∂),(ln ),(μ-y y x y x M ∂∂),(ln ),(μ=xy x N y y x M ∂∂-∂∂),(),( (1。

4) 证明:由定理1。

1得,),(y x μ为微分方程(1。

1)的积分因子的充要条件为xy x N y x y y x M y x ∂∂=∂∂)),(),(()),(),((μμ, 展开即得: x y x y x N ∂∂),(),(μ—y y x y x M ∂∂),(),(μ=),(),(),(y x x y x N y y x M μ⎪⎪⎭⎫⎝⎛∂∂-∂∂。

上式整理即得(1。

积分因子法在求解常微分方程中的应用

积分因子法在求解常微分方程中的应用
基本原理


( 4 ) 齐 次 方 程 老= g ( ÷ ) .
( 1 . 1 . 1 )
1 . 如果微分形式 的一 阶方程 M( x , y ) d x+Ⅳ( , y ) d y=0
的左端恰好是一个二元 函数 U ( x , Y )的全微分 , 即
d U ( x , Y )=M( x , y ) d x+Ⅳ( , ) d y , ( 1 . 1 . 2 ) 得
[ M1 ( ) ⅣI ( y ) +M2 ( ) N 2 ( y ) d y ]=0 ,
+ 。 ,
则称式 ( 1 . 1 . 1 )是全微分方 程 , 而函数 u ( , Y )称为微 分式 ( 1 . 1 . 2 )的原函数 , 方程 ( 1 . 1 . 1 )的通解 为 ( , Y )= C ( c为任意常数 ) . 判别式 ( 1 . 1 . 1 )为全微分方 程 的充 要条
件 是

f I 一 1 J
( 【 1 J . 一 2 . 2 )
只与 有关 , 且此 时有
u =g ( ) 代入原方程 , 有 u+ d u= g ( u ) 即 d



u.

( ):。 J ( 等 一

( 1 . 2 . 3 )
分离变量得= 7

掣 … ,

( 一 ) :。 一

件为 d : d ; 求式 ( 1 . 1 . 1 )为 全 微 分 方 程 时 的 原 函数
U ( x , , , )有 如 下 公 式 :
U ( x , Y ) =J M ( x , y ) d x + J N ( x o ~ Y d y .

微分方程积分因子的研究

微分方程积分因子的研究

20 年 9月 O7
文 章编 号 :0 7—28 (0 7 0 —00 0 10 9 5 2 0 ) 5 0 6— 5
微 分 方 程 积 分 因 子 的 研 究
刘 俊 , 文 娟 汤
65 1) 50 1
( 曲靖师范学 院数学 系 , 云南 曲靖 摘
要: 对微分 方程 的积分 因子进行 了研 究, 找到 了几类微分方程 的积 分因子 文献标识码 : A
( Y M( y d ( Y N( y d =0 , ) , ) x+ , ) , ) y
为一 恰 当方程 , 称 ( Y 为方程 ()的积 分 因子 . 则 , ) 1
引 … 数 (, 为 1式 积 字 充 条 是 理1 函 y ( 的 分因 的 要 件 ) )
Ⅳ 一 M = ( d d
这就 证 明了 = ( y ] 方程 的积 分 因子 的充 要条件 为 , ) 为
( 一 )Ⅳ ( 一 ] .
并 由此得 出积 分 因子 为 ( y , ):

根据 定理 1文献 [ —1]中对 恰 当方 程 () 出 的仅 有关 于 和 Y的积 分 因子 的充要条 件 , , 1 5 1给 就变 成 了
( 一 ) 千 ( = y) 2,
积分因子为 :e ) 『 d ( (

(M O
d ’ 一 ,
推 论 3 方 程 () 1 具有积 分 因子 = (y x )的充要 条 件为
O N ) y N (

mx 一 )1=f x ) (y ,
积分因子为 :e d 『) (( .
微分方程积分因子作深入的研究 . 笔者在此基础上 , 对微分方程 的积分 因子进行深入研究 , 找到 了几类微

常微分方程初等积分法解法研究常微分方程及积分因子

常微分方程初等积分法解法研究常微分方程及积分因子

常微分方程初等积分法解法研究常微分方程及积分因子初等积分法解常微分方程的关键在于求解不定积分。

不定积分是解微分方程的主要手段,通过找到合适的积分因子,可以将一个一阶微分方程转化为一个可积的方程。

在本文中,将对常微分方程及积分因子进行研究。

dy/dx = f(x, y)其中,f(x,y)是已知函数。

解这个方程的方法之一就是通过积分来找到y。

我们需要将这个方程转化成一个可积的形式。

考虑一个形式为 dy/dx + P(x)y = Q(x) 的一阶常微分方程。

要将这个方程转化为可积的形式,需要找到一个因子M(x),使得通过乘以M(x)可以使得原方程的左侧变为一个可积的形式。

这个因子M(x)被称为积分因子。

要找到积分因子,通常通过求解方程 M(x) = 1/M dM/dx = P(x) 来确定。

最常见的积分因子是指数函数,即M(x) = e^(∫P(x)dx)。

通过乘以这个积分因子,原方程可以变为积分形式:d/dx (M(x)y) = M(x)Q(x)通过对上式两边进行不定积分,可以求解出y。

举个例子来说明。

考虑一阶常微分方程 dy/dx + xy = x^2、我们需要找到一个积分因子。

通过解方程 M(x) = 1/M dM/dx = x,可以得到M(x) = e^(1/2 x^2)。

d/dx (e^(1/2 x^2) y) = x e^(1/2 x^2)对上式两边不定积分,得到:e^(1/2 x^2) y = ∫x e^(1/2 x^2) dx通过不定积分求解上式,可以得到y。

通过求解积分因子,我们可以将一阶常微分方程转化为可积的形式。

这种方法适用于一阶线性常微分方程。

对于高阶常微分方程,可以通过转化为一组一阶微分方程来求解。

总结起来,常微分方程及积分因子的研究是通过寻找积分因子来将一阶常微分方程转化为可积的形式。

通过解不定积分,可以求解出未知函数。

初等积分法解常微分方程是一种常用的方法,对于一阶线性常微分方程特别适用。

积分因子法在常微分方程中的应用 开题报告

积分因子法在常微分方程中的应用  开题报告

积分因子法在常微分方程中的应用开题报告开题报告积分因子法在常微分方程中的应用一、选题的背景、意义在许多科学领域中,常常需要研究常微分方程的理论和其解是否存在.常微分方程的理论包括解的存在性和唯一性、奇解、定性理论等等.其中解的讨论也尤为重要,求解方法有很多种,例如,常数变易法、叠加法、积分因子法.求得常微分方程的解能使常微分方程在其他的科学领域有更好的应用.常微分方程在微积分概念出现后即已出现,对常微分方程的研究可分为以下几个阶段.发展初期是针对具体的常微分方程,希望能用初等函数或超越函数表示其解,属于“求通解”的时代.刘维尔在1841年证明了里卡蒂方程不存在一般的初等解,同时柯西又提出了初值问题.因此,早期的常微分方程的求解热潮中断了,而常微分方程从“求通解”时代转向“求定解”时代.19世纪末,常微分方程的研究从“求定解”时代转向“求所有解”的新时代.那是由天体力学中的太阳系稳定性问题需要研究常微分方程解的大范围性态引起的.20世纪末六七十年代以后,常微分方程在计算机技术发展的促进下,从“求所有解”时代转入“求特殊解”时代.求常微分方程的通解在历史上曾作为微分方程的主要目标,一旦求出通解的表达式,就能容易地求出问题所需要的特解;根据通解的表达式可以了解其对某些参数的依赖情况,便于参数取值,使它对应的解具有所需要的性能,也有助于解的其他研究.虽然通过求通解的方法可以求出方程的解,但是有些时候会比较复杂.因此,我们要寻找更为简便的求解方法.对常微分方程的求解.积分因子法是一种很好的求解方法,它能将复杂的计算简单化.二、研究的基本内容与拟解决的主要问题本课题主要对积分因子法进行归纳总结,旨在应用积分因子法来求解常微分方程.本课题的主要目的是通过查阅各种相关文献,寻找各种相关信息,来得到并了解用积分因子法求解常微分方程的一些计算技巧,达到化难为易的目的.先从定义出发,介绍相关的一些基本概念,如微分方程、常微分方程、全微分方程、解、积分因子等以及一些相关的定理和充要条件.接着归纳总结积分因子法:积分因子的求法在求积分因子之前,要对常用的一些简单函数的全微分形式比较熟悉,这样能更快地求出积分因子.(1)观察法求积分因子对于一些形式比较简单的微分方程,可以直接观察出方程的积分因子.如:方程,根据,可以直接观察出它的积分因子为.(2)分组凑微分法对于一些相对复杂的微分方程,可以对其进行分组,然后根据一些简单函数的全微分形式对其进行凑微分,得到其积分因子.(3)重新组合法对于一些相对复杂,不易观察出其积分因子的微分方程,可以将其各项重新组合,再根据一些简单函数的全微分形式通过观察来求得其积分因子.(4)指数待定法求积分因子如果微分方程中是的多项式,则可以找到形式的积分因子.(5)公式法求积分因子对一些非全微分方程可以用上面提到的四种方法求得它们的积分因子,但还有一些非全微分方程用上述四种方法不太容易得到它们的积分因子,这时就可以用一些公式来求解.不同的公式都有其相对应的条件需要满足.积分因子巧解常微分方程(1)观察法对于简单形式的微分方程,可以根据一些简单函数的全微分形式直接观察出方程的积分因子,再将积分因子乘到原方程的两边形成全微分方程进行求解.(2)分组凑微分法将微分方程重新分组,化成易求得积分因子的形式,求得其积分因子,再将积分因子乘到原方程的两边形成全微分方程进行求解.(3)重新组合法将微分方程进行重新组合,化成易求得积分因子的形式,求得其积分因子,再将积分因子乘到原方程的两边形成全微分方程进行求解.(4)指数待定法对符合特定条件的微分方程,用指数待定的方法求得其积分因子,再将积分因子乘到原方程的两边形成全微分方程进行求解.(5)公式法针对不同的微分方程,运用相对应的公式求得其积分因子,再将积分因子乘到原方程的两边形成全微分方程进行求解.积分因子法在一阶常微分方程中的应用(1)在可分离变量微分方程中的应用如果一阶微分方程可变化为的形式,则称这个方程为可分离变量方程.运用积分因子法求得这类方程的积分因子,将方程转化为全微分方程进行求解.(2)在齐次微分方程中的应用方程是齐次方程.运用积分因子法求得这类方程的积分因子,将方程转化为全微分方程进行求解.(3)在一阶线性微分方程中的应用设一阶线性微分方程为将其成对称的形式若方程有一个仅依赖于的积分因子,则,其中;反之,若仅依赖于,则是方程的一个积分因子.(4)在贝努力方程中的应用将贝努力方程令,可以将方程化为一阶线性微分方程然后用积分因子求解此方程积分因子法在二阶常微分方程中的应用二阶线性微分方程,当时,此方程为齐次方程;而当时,此方程为非齐次方程.运用积分因子法对二阶线性微分方程进行求解.积分因子法的其他应用证明一些初等公式或一些命题.三、研究的方法与技术路线、研究难点,预期达到的目标本课题归纳总结的主要内容是积分因子法在常微分方程中的应用.利用积分因子法来解决常微分方程的一些复杂的计算问题,使计算过程更加简单易理解.并且积分因子是不唯一,有简单也有复杂.不管它如何,它在常微分方程的计算中都有着不简单的力量.通过参考一些文献资料,以及自己对文献资料的理解和自己掌握的知识,并经过自己的努力,在最后可以用积分因子法解决一些常微分方程的计算.常微分方程的解本来就是一个难点,又由于对积分因子的了解不是很深,在之前学习的只是最基础的.因此,对于它的应用还是有一定的难度的.尽管这个课题有一定的难度,但是我相信不管困难是什么,总能找出方法来解决的.应用积分因子法可以使很多常微分方程的计算得到简化,能够达到化难为易的目的.常微分方程的研究与其他学科领域的结合,使得各种新的研究分支出现.相信常微分方程会在更多的科学领域有更好的应用,并会有更好的发展,做出更大的贡献.四、论文详细工作进度和安排2011-02-21至2011-03-20完成初稿;2011-03-21至2011-04-20在导师的指导下完成第一次修改;2011-04-21至2011-05-20在导师的指导下完成第二次修改并定稿;2011-05-21至2011-05-23准备论文答辩.五、主要参考文献:[1]时宝,黄朝炎.微分方程基础及其应用[M].北京:科学出版社.2007:2-3.[2]丁同仁,李承治.常微分方程教程[M].北京:高等教育出版社.2004,19:32-33,46-47.[3]试析一阶微分方程的积分因子[J].许昌师专学报.1993,312:9,35-39.[4]杨雨民.积分因子咋一阶线性微分方程中的应用[J]. 辽宁省交通高等专科学校学报.1997,51:30-33.[5]James?Stewart?Calculus:Early?Transcendentals5thed?[M].北京:高等教育出版社.2004:598-601,641-643.[6]张奕河,郭文川.关于一阶常微分方程的积分因子求解问题[J].四川理工学院学报(自然科学版).2009,226:11-13.[7]Ma?Yuan-jing.Runge-Kutta?type?iterative?method?for?nonlinear?e quations[J]JOUJOURNAL?OF?NATURAL?SCIENCE?OF?HEI?LONG?JIANG?UNIVERSITY.2009,264:431-435[8]潘鹤鸣.几种特殊类型积分因子的求法及在解微分方程中的应用[J].巢湖学院学报.2003,53:18-22.[9]徐安农,段复建.全微分方程与积分因子法[J].桂林电子工业学院学报.2002,222:10-12.[10]温启军,张丽静.关于积分因子的讨论[J]长春大学学报.2006,165:17-20.[11]龚雅玲.求解微分方程的积分因子法[J].黔南昌教育学院学报.2007,221:32-35.[12]张凤然,马金江.二阶变系数线性微分方程的积分因子解法[J].2008,628:13-15.。

c 微分方程积分因子的研究

c   微分方程积分因子的研究

第28卷第5期吉首大学学报(自然科学版)Vol.28No.5 2007年9月J ournal of J ishou University(Natural Science Edi ti on)Sept.2007文章编号:1007-2985(2007)05-0006-05微分方程积分因子的研究X刘俊,汤文娟(曲靖师范学院数学系,云南曲靖655011)摘要:对微分方程的积分因子进行了研究,找到了几类微分方程的积分因子.关键词:微分方程;积分因子;通解中图分类号:O175.14文献标识码:A恰当方程可以通过公式法求出它的通解,因此如何通过积分因子将一个非恰当方程化为恰当方程就显得尤其重要.文献[1-15]对恰当方程M(x,y)d x+N(x,y)d y=0仅给出了关于x和y的积分因子的充要条件,介绍了用熟记简单的二元函数全微分公式、用分项组合的方法求解微分方程的解通,但尚未对微分方程积分因子作深入的研究.笔者在此基础上,对微分方程的积分因子进行深入研究,找到了几类微分方程的积分因子的求法,以弥补这一工作.1预备知识定义1[1]若方程M(x,y)d x+N(x,y)d y=0(1)的左端恰好是某个二元函数u(x,y)的全微分,则称(1)式为恰当方程.定义2[1]若存在连续可微的函数L=L(x,y)X0,使得L(x,y)M(x,y)d x+L(x,y)N(x,y)d y=0为一恰当方程,则称L(x,y)为方程(1)的积分因子.引理1[1]函数L(x,y)为(1)式的积分因子的充要条件是9(L M)9y=9(L N)9x,即N 9L9x-M9L9y=(9M9y-9N9x)L.2主要结论定理1方程M(x,y)d x+N(x,y)d y=0具有形如L=L[<(x,y)]的积分因子的充要条件为(9M9y-9N9x)(N9<9x-M9<9y)-1=f[<(x,y)].证明方程M(x,y)d x+N(x,y)d y=0有积分因子的充要条件为X收稿日期:2007-05-14基金项目:国家自然科学基金资助项目(10361007);云南省自然科学基金资助项目(2006A0082M);云南省教育厅科学研究项目(5Z0030A)作者简介:刘俊(1963-),男,云南昆明人,曲靖师范学院数学系教授,主要从事微分方程、动力系统研究.N 9L9x-M9L9y=(9M9y-9N9x)L.令L=L[<(x,y)],则有N d Ld<#9<9x-M d Ld<#9<9y=(9M9y-9N9x)L(<),d L L=(9M9y-9N9x)(N9<9x-M9<9y)-1d<.这就证明了L=L[<(x,y)]为方程的积分因子的充要条件为(9M9y-9N9x)(N9<9x-M9<9y)-1=f[<(x,y)].并由此得出积分因子为L(x,y)=e Q f(<)d<.根据定理1,文献[1-15]中对恰当方程(1)给出的仅有关于x和y的积分因子的充要条件,就变成了定理1的特殊情形.而且可得到如下结论:推论1方程(1)具有积分因子L=L(x?y)的充要条件为(9,9y-9N9x)(NºM)-1=f(x?y),积分因子为L=eQ f(x?y)d(x?y).推论2方程(1)具有积分因子L=L(x2?y2)的充要条件为(9M9y-9N9x)(NxºMy)-1=f(x2?y2),积分因子为L=e Q(x2?y2)d(x2?y2).推论3方程(1)具有积分因子L=L(xy)的充要条件为(9,9y-9N9x)(Ny-Mx)-1=f(xy),积分因子为L=eQ(xy)d(xy).推论4方程(1)具有积分因子L=L(yx)的充要条件为-(9M9y-9N9x)(Mx+Mx2)-1=f(yx),积分因子为L=e Q(y x)d(y x).推论5方程(1)具有特殊积分因子L=L(x A y B)的充要条件为1 x A y B (9M9y-9N9x)(A Nx-B My)-1=f(x A y B),积分因子为L=eQ f(x A y B)d(x A y B).注推论5中A,B可用待定系数法确定.定理2假设(1)式中M(x,y)和N(x,y)存在以下关系:9M 9y-9N9x=Nf(x)-Mg(y),其中f(x),g(y)分别为x,y的连续函数.则方程(1)的积分因子为L(x,y)=e Q f(x)d x+Q g(y)d y.证明由M(x,y)和N(x,y)存在的关系得到9M9y+Mg(y)=9N9x+Nf(x).两边同乘以e Q f(x)d x+Q g(y)d y,得e Q f(x)d x+Q g(y)d y[9M9y+Mg(y)]=eQ f(x)d x+Q g(y)d y[9N9x+Nf(x)],即7第5期刘俊,等:微分方程积分因子的研究9(e Q f(x)d x+Q g(y)d y#M)9y=9(e Q f(x)d x+Q g(y)d y#N)9x.则(eQ f(x)d x+Q g(y)d y#M)d x=(e Q f(x)d x+Q g(y)d y#N)d y=0为恰当方程.所以L(x,y)=e Q f(x)d x+Q g(y)d y为方程M(x,y)d x+N(x,y)d y=0的积分因子.定理3若方程(1)为齐次方程,则积分因子为L=1xM+yN,其中M(x,y),N(x,y)是m次齐次函数.证明作变换z=yx,有d y=x d z+z d x.由于M(x,y)=M(x,xz)=x m#M(1,z),N(x,y)=N(x,xz)=x m#N(1,z),将其代入方程(1),得到变量分离方程x m[M(1,z)+N(1,z)]d x+x m+1#N(1,z)d z=0,其积分因子为L(x,z)=1x m+1[M(1,z)+N(1,z)#z].代回原变量,得L(x,y)=1xM(x,y)+yN(x,y),即为齐次方程的积分因子,其中xM+yN X0.定理4变量分离方程M(x)N(y)d x+P(x)Q(y)d y=0有积分因子L=1N(y)P(x).证明用L=1N(y)P(x)乘以变量分离方程两端,得M(x)P(x)d x+Q(y)N(y)d y=0.显然,这个方程是恰当方程.因此变量分离方程的积分因子为L=1N(y)P(x).定理5设函数f(xy)+g(xy)连续可微且f(xy)X g(xy),则方程yf(xy)d x+xg(xy)d y=0(2)有积分因子L=1xy[f(xy)-g(xy)].证明令xy=u,则原方程可化为u[f(u)-g(u)]xd x+g(u)d u=0.两边乘以1u[f(u)-g(u)],得d xx-g(u)u[f(u)-g(u)]d u=0.显然,这是一个恰当方程.故(2)式有积分因子1u[f(u)-g(u)],因而原方程有积分因子(xy[f(xy)-g(xy)])-1.定理6若r,s,m,n,r,s,k,l是任意常数,且ml-kn X0,则方程x r y s(my d x+nx d y)+x Q y R(ky d x+lx d y)=0有积分因子L=x A y B(A,B是常数),且通解为mA+r+1x A+r+1y B+s+1+kA+Q+1xA+Q+1y B+R+1=c,8吉首大学学报(自然科学版)第28卷或者n B +s +1x A +r +1y B +s +1+l B +R +1x A +Q +1y B +R +1=c .证明 (1)先证:若L (x ,y )是方程M (x ,y )d x +N (x ,y )d y =0的积分因子,则L 2(x ,y )也是该方程的积分因子.由定理2知L (x ,y )=e Q f (x )d x +Q g (y )d y 是方程M d x +N d y =0的积分因子,而L 2(x ,y )=(eQ f (x )d x +Qg (y )d x)2=e 2(Q f (x )d x +Q g (y )d y ).因为L 2M d x +L 2N d y =0,所以9(L 2M )9y -9(L 2N )9x =(2L M 9L 9y +L 29M 9y )-(2L N 9L 9x +L 29N9x)=2L (M9L 9y -N 9L 9x )+L 2(9M 9y -9N9x).因为f (x ),g (y )分别为x ,y 的连续函数,而9M 9y -9N 9x =2Nf (x )-2Mg (y ),9L 9y =L g (y ),9L9x =L f (x ),9(L 2M )9y -9(L 2N )9x=2L 2[Mg (y )-Nf (x )]+2L 2[Nf (x )-Mg (y )]=0]9(L 2M )9y =9(L 2N )9x, 所以L 2M d x +L 2N d y =0为恰当方程.(2)再证:形如x E y D(my d x +nx d y )=0的一阶微分方程的积分因子为L (x ,y )=x km -E -1ykn -D -1.令M (x ,y )=L (x ,y )mx E yD +1,N (x ,y )=L (x ,y )nxE +1y D.因为9M 9y =mx E y D [9L 9y y +L (D +1)],9N 9x =nx E y D [9L 9xx +L (E +1)],由定理2可知L (x ,y )=e Q f (x )d x +Q g (y )d y ,于是9L 9x =L f (x ),9L 9y=L g (y ),所以9M 9y =m L x E y D [g (y )y +D +1],9N 9x=n L x E y D [f (x )x +E +1].又因为Mg (y )=m L x E yD +1g (y ),Nf (x )=n L xE +1y Df (x ),所以Nf (x )-Mg (y )=L x E y D[nf (x )x -mg (y )y ],9M 9y-9N9x =L x E y D {m [g (y )y +D +1]-n [f (x )x +E +1]},9M 9y -9N9x=Nf (x )-Mg (y ).因此2f (x )x +E +1m =2g (y )y +D +1n.令2f (x )x +a +1m =2g (y )y +b +1n=k ,则f (x )=km -E -12x ,g (y )=kn -D -12y,Qf (x )d x =km -E -12ln x +c 1,Qg (y )d y =kn -D -12ln x +c 1.所以 L = c xkm -E -12ykn -D -12.令L = L 2,则L (x ,y )=xkm -E -1ykn -D -1(k I c )是原方程的积分因子.故方程x r y s(my d x +nx d y )=0有积分因子L =x ma -r -1yn a -s -1,方程x Q y R(ky d x +lx d y )=0有积分因子L =9第5期 刘 俊,等:微分方程积分因子的研究x kb-Q-1y lb-R-1.选择A,B使这2个积分因子重合,即ma-r-1=kb-Q-1=A,na-s-1=lb-R-1=B.则A,B需满足条件:当nk X ml时,n(A+r+1)=m(B+s+1),l(A+Q+1)=k(B+R+1)..用x A,y B乘原方程两端,得(mx A+r y B+s+1d x+nx A+r+1y B+s d y)+(kx A+Q y B+R+1d x+lx A+Q+1y B+R d y)=0.积分后,有mA+r+1x A+r+1y B+s+1+kA+Q+1xA+Q+1y B+R+1=c,或者nB+s+1x A+r+1y B+s+1+lB+R+1xA+Q+1y B+R+1=c.参考文献:[1]王高雄,周之铭,朱思铭,等.常微分方程(第3版)[M].北京:高等教育出版社,2006.[2][德]W.瓦尔特(著).邹应(译).常微分方程导论[M].北京:高等教育出版社,1988.[3][美]R.布朗森(著).王祖城,沈功(译).现代微分方程的理论和习题[M].北京:中国铁道出版社,1984.[4]王兴涛.常微分方程[M].哈尔滨:哈尔滨工业大学出版社,2003.[5]丁同仁.常微分方程基础[M].上海:上海科学技术出版社,2003.[6]东北师范大学数学系微分方程教研室.常微分方程[M].北京:高等教育出版社,2004.[7]周义仓,靳祯,秦军林.常微分方程及其应用)))方法、理论、建模、计算机[M].北京:科学出版社,2003.[8]都长清,焦宝聪,焦炳照.常微分方程[M].北京:北京师范学院出版社,1993.[9]王奎实,杨德全,陈锡华.常微分方程[M].长春:吉林教育出版社,1987.[10]温耀华.常微分方程[M].重庆:西南师范大学出版社,1989.[11]刘志汉.常微分方程[M].西安:陕西师范大学出版社,1987.[12]蔡燧林.常微分方程[M].杭州:浙江大学出版社,2002.[13]高素志,马遵路,曾昭.常微分方程[M].北京:北京师范大学出版社,2000.[14]庄万,黄启宇,丛树凡.常微分方程[M].济南:山东科学技术出版社,1998.[15]王光发,吴克乾,邓宗琦.常微分方程[M].长沙:湖南教育出版社,1988.Research of the Integral Factor of Differential EquationsLI U Jun,TANG Wen-juan(Mathematics Department,Qujing Normal College,Qujing655011,Yunan China)Abstract:The inte gral factors of the differential equations are studied,and the integral factors of categories of differen-tial equations are obtained.Key words:differential equation;integral factor;general solution(责任编辑向阳洁) 10吉首大学学报(自然科学版)第28卷。

关于对常微分方程中积分因子的研究毕业论文

关于对常微分方程中积分因子的研究毕业论文

毕业论文论文题目:关于对常微分方程中积分因子的研究姓名:学科专业:数学教育指导教师:完成时间:2011年 5 月20 日摘 要研究了四种一阶微分方程的积分因子存在的充要条件,主要通过一些特殊变形的方法来得到这四种类型的常微分方程积分因子的通解,具体可以分为:()(),,y x G y x +=μ()(),,by ax G y x +=μ()()b a y x G y x -=,μ及()()xy G y x =,μ, 四种类型各有特点但又互有联系。

关键词:常微分方程;积分因子;充要条件目录第一章引言 (1)第二章四类常微分方程积分因子的存在充要条件 (2)()()yμ (2),=y x+Gx()()byμ (3)=x+,axGy()()ba yμ (4)x-,=yxG()()xyμ (5),Gx=y参考文献 (7)第一章 引言一阶微分方程()()0,,=+dy y x N dx y x M ()的积分因子的形式,其一,其求解方法是根据类型确定解法,其中一类是全微分方程,所谓全微分方程就是方程()的左端恰为某个方程的全微分。

我们知道方程()是全微分的充要条件是,xNy M ∂∂=∂∂当此不满足的时候,方程()就不是全微分方程,此时若有一个恰当的函数()0,≠y x μ使方程(1)两端乘以后()0,≠y x μ所得的方程()()()()0,,,,=+y x N y x y x M y x μμ ()为全微分方程,则称函数()y x ,μ为方程()的积分因子。

积分因子存在的充要条件:如何求方程()的积分因子?一下就是关于()y x ,μ为积分因子的充要条件,微分方程()()()()0,,,,=+y x N y x y x M y x μμ ()为全微分方程的充要条件是: ()()()()()(),,,,,xy x N y x y y x M y x ∂∂=∂∂μμ即()()()()()()()()()()()xy x y x N y y x N y x x y x y x M y y x M y x ∂∂+∂∂=∂∂+∂∂,,,,,,,,μμμμ () (),,μμ=y x (),,M y x M = (),,N y x N =上式可以整理到μμμ⎪⎪⎭⎫⎝⎛∂∂-∂∂=∂∂-∂∂x N y M y M x N() 所以()y x ,μ为方程()的积分因子的充要条件是()y x ,μ为方程的解。

微分方程的积分因子法

微分方程的积分因子法

微分方程是数学中重要的研究对象,它通过描述变量之间的关系,可以用来解释许多自然现象和物理规律。

微分方程的求解是数学分析的重要方法之一,其中积分因子法是一种常用且有效的求解微分方程的方法。

首先,我们来了解什么是微分方程。

微分方程是包含未知函数及其导数的方程,一般形式为dy/dx = f(x,y),其中y是未知函数,f(x,y)是已知的函数。

微分方程可以分为常微分方程和偏微分方程两类,常微分方程中只包含一个自变量,而偏微分方程中包含多个自变量。

解微分方程要找出满足方程的函数形式,而积分因子法是一种特殊的方法用来解决一类形式为M(x,y)dx + N(x,y)dy = 0的一阶常微分方程。

积分因子法的思想是通过引入一个适当的积分因子来改变微分方程的形式,从而使其变得可积。

具体步骤如下:1.将方程化为其标准形式:M(x,y)dx + N(x,y)dy = 0,其中M(x,y)和N(x,y)为已知函数。

2.判断方程是否是恰当微分方程。

若满足∂M/∂y = ∂N/∂x,则该方程为恰当微分方程,直接求解即可;若不满足,则进行下一步。

3.求取积分因子。

积分因子可以通过通解公式I(x) = e^(∫P(x)dx),其中P(x)为方程的系数。

4.将积分因子乘到方程上,得到恰当微分方程:I(x)M(x,y)dx +I(x)N(x,y)dy = 0。

5.求解恰当微分方程。

由于恰当微分方程是可积的,可以直接求出其解。

通过这样的步骤,利用积分因子法可以将一些常见的非恰当微分方程转化为恰当微分方程,从而能够更方便地求解微分方程。

需要指出的是,积分因子法并不适用于所有的微分方程,只适用于一些具有特定形式的微分方程。

对于其他形式的微分方程,可能需要使用其他的求解方法。

总结来说,积分因子法是一种求解常微分方程的有效方法,它通过引入适当的积分因子,将非恰当微分方程转化为恰当微分方程,从而更容易求解。

使用积分因子法需要熟悉方程的形式及其特点,才能正确选择和应用积分因子。

积分因子法在求解常微分方程中的应用

积分因子法在求解常微分方程中的应用

积分因子法在求解常微分方程中的应用常微分方程作为现代数学的重要分支,其应用范围广泛,涉及到物理、计算机科学等领域。

求解常微分方程是常微分方程理论的核心,而积分因子法作为其中的重要方法之一,常常被应用于常微分方程的求解中。

1. 什么是积分因子法?积分因子法是利用一个与方程解相关的因子来将常微分方程转化为可积的形式的一种方法。

在求解常微分方程时,为了保证方程解的双曲性或椭圆性,我们可能需要乘上一个符合要求的函数因子使其可以进行精确积分,这个函数因子就被称为积分因子。

2. 如何应用积分因子法?应用积分因子法的关键是需要找到符合要求的积分因子。

一般来说,积分因子需要满足以下条件:(1)积分因子最好能够求得,即它可以具体的表达式表示出来。

(2)积分因子必须非零。

(3)积分因子的乘积与微分方程的系数的组合必须是可积的。

(4)积分因子在微分方程所考虑的区域上必须是连续的。

(5)积分因子应该是一种容易求得的函数形式。

找到符合要求的积分因子后,我们就可以将常微分方程乘上这个因子,从而将其转化为一个可积的形式。

通过对等式两边的乘积进行积分,最终获得方程的解析解。

3. 积分因子法在求解实际问题时的应用积分因子法在求解实际问题时的应用有很多。

例如在物理学中,通过应用积分因子法可以求解出多个物理系统的行为规律。

在这种情况下,微分方程主要描述物理量的变化,而积分因子则为了提高求解的准确度和精度。

在计算机科学领域,积分因子法的应用同样非常广泛。

在进行数值计算时,我们经常需要通过微分方程来描述系统的行为规律。

但由于数值方法的固有误差,我们得出的解往往不够精确。

而在这种情况下,我们可以通过引入一个积分因子来提高求解的精度。

总的来说,积分因子法在求解常微分方程中起着重要的作用。

它可以帮助我们获得更加精确和准确的解析解,而这些解析解在现代数学和其它学科领域中有着广泛的应用。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
吉首大学学报(自然科学版) JOURNAL OF JISHOU UNIVERSITY(NATURAL SCIENCE EDITION) 2007,28(5) 0次
定理3 若方程(1)为齐次方程,则积分因子为p=五—丽,其中肘(石,y),Ⅳ(并,,,)是m次齐次
函数.
证明 作变换z=孝,有dy=菇dz+zd菇.由于
M(菇,Y)=M(茗,愆)=戈“·M(1,z),N(z,Y)=N(戈,舷)=髫“·N(1,:), 将其代入方程(1),得到变量分离方程
菇“[M(1,名)+N(1,z)]dx+戈“+1·N(1,z)dz=0,
由定理主可知∥(z,y)=小枷#吖a(Y)dY,于是罄=朋菇),丝ay=昭(y),所以
又因为 所以
警=mttx y‘[g(Y)Y+¨1]aaN石=n,ux"yS Ef(石)菇+e+1].
Ug(y)=m1.tx"y8+19(y),彤(石)=nlzx"”/y(x), Nf(菇)一坳(),)=tzx",[nf(菇)菇一mg(y)y],
1 预备知识 , 定义1【l】 若方程
M(茗,y)dx+N(z,y)dy=0
的左端恰好是某个二元函数u(石,Y)的全微分,则称(1)式为恰当方程. 定义2… 若存在连续可微的函数户=产(髫,Y)≠O,使得
卢(菇,Y)M(髫,y)dx+∥(菇,Y)N(x,y)dy=0
为一恰当方程,则称岸(z,Y)为方程(1)的积分因子.
z’,(mydx+nxdy):0有积分因子口:髫一“1广。~,方程扩,(bdx+lxdy):0有积分因子肛:
万方数据
10
吉首大学学报(自然科学版)
xkb-纠Yn”1.选择a,卢使这2个积分因子重合,即
{m阳。a一-。r一-。1:=。k6b一-口p一-。1』爱’
则a,p需满足条件:当nk≠ml时,
一(筹一掣)(警+扩2=八》 、a,,
a茄7、搿 ’戈7
一J、髫7’
积分因子为弘:e『(詈)d(詈). 推论5 方程(1)具有特殊积分因子P=P(z4,)的充要条件为
南矿(,碧、一筹a,)(,警一一趔a菇)~,:、以z扩广一),Y 7 一.,、“了7’
积分因子为卢:小,,)d(,^.
注 推论5中a,p可用待定系数法确定. 定理2 假设(1)式中肘(菇,Y)和Ⅳ(省,Y)存在以下关系:
筹一筹=∥广{耐g(,,)y+艿+1]一凡If(茹)茹+e+1]},
晋一掣=瞅名)一us(y).
因此
2丛苎2篓±兰±:!一29(z2z±查±!
m’’

/'t
令丝巫坠m旦盟陬:枷纽龙∽:型掣以』l,l省n丛戈:)+=|c|掣l},,,f则gg(,(,,),d,)=,掣:,掣ln戈+c1.
所以五:施盘2甲,,盘学.令∥:∥--2,则产(菇,y):石h—z广小t(||}∈c)是原方程的积分因子.故方程
摘要:对微分方程的积分因子进行了研究,找到了几类微分方程的积分因子
关键词:微分另程;积分因子;通解
中图分类号:0175.14
文献标识码:A
V01.28 No.5 Sept.2007
恰当方程可以通过公式法求出它的通解,因此如何通过积分因子将一个非恰当方程化为恰当方程就 显得尤其重要.文献[1—15]对恰当方程M(x,y)dx+N(x,y)dy=0仅给出了关于菇和y的积分因子的 充要条件,介绍了用熟记简单的二元函数全微分公式、用分项组合的方法求解微分方程的解通,但尚未对 微分方程积分因子作深入的研究.笔者在此基础上,对微分方程的积分因子进行深入研究,找到了几类微 分方程的积分因子的求法,以弥补这一工作.
第28卷
参考文献: [1]王高雄,周之铭,朱思铭,等.常微分方程(第3版)[M].北京:高等教育出版社,2006. [2][德]w.瓦尔特(著).邹应(译).常微分方程导论[M].北京:高等教育出版社,1988. [3][美]R.布朗森(著).王祖城,沈功(译).现代微分方程的理论和习题[M].北京:中国铁道出版社,1984. [4] 王兴涛.常微分方程[M].哈尔滨:哈尔滨工业大学出版社,2003. [5] 丁同仁.常微分方程基础[M].上海:上海科学技术出版社,2003. [6] 东北师范大学数学系微分方程教研室.常微分方程[M].北京:高等教育出版社,2004.

万方数据

吉首大学学报(自然科学版)
第28卷
鱼(d!:!竺:』!!竺:丝)一鱼— (d!:!!::』!!:!竺:盟2
ay
a戈


(办枷*+胁m.M)d菇:(出枷*+如m.N)dy:0
为恰当方程.所以卢(石,),):小枷*+』s(,)d,为方程M(龙,y)d菇+Ⅳ(搿,,,)d,,:0的积分因子.
鱼(拦:丝)一鱼(笆!盟2
a,,

a算

所以,u2Mdx+p2Ndy=0为恰当方程.
(2)再证:形如矿广(mydx+nxdy)=0的一阶微分方程的积分因子为产(算,Y)=少。。1广。一.令 M(x,Y)=产(髫,,,)眦‘广“,Ⅳ(并,Y)=产(戈,,,)/1St;z“广.因为
警=臌。y.E警y),+产(艿+1)],筹=∥广[静+产(e+1)],
推论2 方程(1)具有积分因子p=产(戈2±Y2)的充要条件为
(警一掣)(眠干嘶)~=,(戈2。y2),
积分因子为口:e『(x2ty2)d(z2t^.
推论3 方程(1)具有积分因子∥=产(形)的充要条件为
(票一m。N菇)(M一胝)~=f(xy),
积分因子为卢:J(枷(圳.
推论4 方程(1)具有积分因子P=∥(詈)的充要条件为
g(形)])~.
定理6 若r,s,m,n,r,s,k,l是任意常数,且ml—kn≠0,则方程
xrys(mydx+nxdy)+扩厂(.砂dx+Ixdy)=0
:南石。+r1广+;+1南a+p1 有积分因子卢=矿广(口,p是常数),且通解为 。 口+r+
口+p+ 1广+口。 +1=c,’
万方数据
第5期
Key words:differential equation;integral factor;general solution
(责任编辑向阳洁)
万方数据
微分方程积分因子的研究
作者: 作者单位: 刊名:
英文刊名: 年,卷(期): 被引用次数:
刘俊, 汤文娟, LIU Jun, TANG Wen-juan 曲靖师范学院数学系,云南,曲靖,655011
第28卷第5期 2007年9月
吉首大学学报(自然科学版) Joumal of Jishou University(Natural Science Edition)
文章编号:1007—2985(2007)05一O006—05
微分方程积分因子的研究
刘俊,汤文娟
(曲靖师范学院数学系,云南曲靖655011)
[7] 周义仓,靳祯,秦军林.常微分方程及其应用一方法、理论、建模、计算机[M].北京:科学出版社,2003
[8]都长清,焦宝聪,焦炳照.常微分方程[M].北京:北京师范学院出版社,1993. [9]王奎实,杨德全,陈锡华.常微分方程[M].长春:吉林教育出版社,1987. [10] 温耀华.常微分方程[M].重庆:西南师范大学出版社,1989. [11] 刘志汉.常微分方程[M].西安:陕西师范大学出版社,1987. [12]蔡燧林.常微分方程[M].杭州:浙江大学出版社,2002. [13] 高素志,马遵路,曾昭.常微分方程[M].北京:北京师范大学出版社,2000. [14]庄万,黄启宇,丛树凡.常微分方程[M].济南:山东科学技术出版社,1998. [15]王光发,吴克乾,邓宗琦.常微分方程[M].长沙:湖南教育出版社,1988.
引理1… 函数p(龙,,,)为(1)式的积分因享的充要条件是丛g盟:丛筻堕,即


aV
a并
Ⅳo警Z 一M警dV=(筹一d’,警)d产髫’ .
2 主要结论
定理1 方程M(x,y)dx十N(x,y)dy=0具有形如户=P[≯(菇,Y)]的积分因子的充要条件为
(筹一掣)(Ⅳ墓一肘器)-1刊№∽].
证明 方程M(x,y)dx+N(x,y)dy=0有积分因子的充要条件为
警一瓦DN=Nf(菇)一ug(y),
其中以菇),g(y)分别为菇,y的连续函数.则方程(1)的积分因子为p(茹,y):小枷;小(,)打.
证明 由M(菇,Y)和N(菇,Y)存在的关系得到
碧+魄(y)=碧+Nf(菇).
两边同乘以出枷*+fs(,)“,得
小枷小Ⅲ,[警+坛(y)]=办枷小棚7[掣+Nf(菇)],
刘俊,等:微分方程积分因子的研究

Ⅳd誊并 一M警d'=,(警Yo 一筹)a产茗。.
令p=卢[5&(菇,Y)],贝0有
Ⅳ嚣·譬一M嚣·琴=(筹一DaN髫)P(≯),警=(DaMy一掣)(Ⅳ髦一M dDye).1 dj5.
这就证明了卢=∥[声(菇,y)]为方程的积分因子的充要条件为
(DaM,,一謦)(Ⅳ蓑一M磬)-1=∥㈠y)].
定理5 设函数八xy)+g(xy)连续可微且以xy)≠g(xy),则方程
2万D巧习之硼· 有积分因子P
xf(xy)dx+ 令xy=“,则原方程可!生化[为.£生些立_掣d菇+g(M)dⅡ:o.
两边乘以可蕊百产硼,得
显然,这是一个恰当方程·故(2)警髫 式有一一积面M欧分毋因u)兰子一可g(翮M灭)di]饷M叫=,o2叭 .因而原方程有积分因子(拶[八夥)一
『n(口+r+1)=m(p+s+1),
【Z(a+ID+1)=k(卢+盯+1).‘
用菇。,,乘原方程两端,得
(m,x”r/”“dx+rtx””1广”dy)+(础+9广”“d茹+舻+P+1广”dy)=0.
积分后,有
a+r+1广+J+1+南戈 ———m———1算
4+9+1,+4+1=c,
相关文档
最新文档