各类微分方程的解法

各类微分方程的解法

1.可分离变量的微分方程解法

一般形式:g(y)dy=f(x)dx

直接解得∫g(y)dy=∫f(x)dx

设g(y)及f(x)的原函数依次为G(y)及F(x),则G(y)=F(x)+C 为微分方程的隐式通解

2.齐次方程解法

一般形式:dy/dx=φ(y/x)

令u=y/x 则y=xu,dy/dx=u+xdu/dx,所以u+xdu/dx=φ(u),即du/［φ(u)-u ］=dx/x 两端积分,得∫du/［φ(u)-u ］=∫dx/x

最后用y/x 代替u,便得所给齐次方程的通解

3.一阶线性微分方程解法

一般形式:dy/dx+P(x)y=Q(x)

先令Q(x)=0则dy/dx+P(x)y=0解得y=Ce

－∫P(x)dx ,再令y=u e －∫P(x)dx 代入原方程 解得u=∫Q(x) e

∫P(x)dx dx+C,所以y=e －∫P(x)dx ［∫Q(x)e ∫P(x)dx dx+C ］ 即y=Ce

－∫P(x)dx +e －∫P(x)dx ∫Q(x)e ∫P(x)dx dx 为一阶线性微分方程的通解 4.可降阶的高阶微分方程解法 ①y

(n)=f(x)型的微分方程 y

(n)=f(x) y

(n-1)= ∫f(x)dx+C 1 y (n-2)= ∫［∫f(x)dx+C 1］dx+C 2

依次类推,接连积分n 次,便得方程y

(n)=f(x)的含有n 个任意常数的通解

②y ”=f(x,y ’) 型的微分方程 令y ’=p 则y ”=p ’,所以p ’=f(x,p),再求解得p=φ(x,C 1)

即dy/dx=φ(x,C 1),所以y=∫φ(x,C 1)dx+C 2

③y ”=f(y,y ’) 型的微分方程

令y ’=p 则y ”=pdp/dy,所以pdp/dy=f(y,p),再求解得p=φ(y,C 1)

即dy/dx=φ(y,C 1),即dy/φ(y,C 1)=dx,所以∫dy/φ(y,C 1)=x+C 2

5.二阶常系数齐次线性微分方程解法

一般形式:y ”+py ’+qy=0,特征方程r 2

+pr+q=0

6.二阶常系数非齐次线性微分方程解法

一般形式: y”+py’+qy=f(x)

先求y”+py’+qy=0的通解y0(x),再求y”+py’+qy=f(x)的一个特解y*(x) 则y(x)=y0(x)+y*(x)即为微分方程y”+py’+qy=f(x)的通解

求y”+py’+qy=f(x)特解的方法:

①f(x)=P m(x)eλx

型

令y*=x k Q m(x)eλx［k按λ不是特征方程的根,是特征方程的单根或特征方程的重根依次取0,1或2］再代入原方程,确定Q m(x)的m+1个系数

②f(x)=eλx［Pｌ(x)cosωx+P n(x)sinωx］型

令y*=x k eλx

［Q m(x)cosωx+R m(x)sinωx］［m=max﹛ｌ,n﹜,k按λ+iω不是特征

方程的根或是特征方程的单根依次取0或1］再代入原方程,分别确定Q m(x)和R m(x)的m+1个系数

附

微分方程在物理学中的应用:

⑴找准合适的研究对象

⑵确定正确的数学模型

⑶联列合理的微分方程

⑷解出最佳的方程结果

执笔：缪张华

一阶常微分方程解法总结

第 一 章 一阶微分方程的解法的小结 ⑴、可分离变量的方程： ①、形如 )()(y g x f dx dy = 当0)(≠y g 时，得到 dx x f y g dy )() (=，两边积分即可得到结果； 当0)(0=ηg 时，则0)(η=x y 也是方程的解。 例1.1、 xy dx dy = 解：当0≠y 时，有 xdx y dy =，两边积分得到)(2ln 2为常数C C x y += 所以)(112 12 C x e C C e C y ±==为非零常数且 0=y 显然是原方程的解； 综上所述，原方程的解为)(12 12 为常数C e C y x = ②、形如0)()()()(=+dy y Q x P dx y N x M 当0)()(≠y N x P 时，可有 dy y N y Q dx x P x M ) () ()()(=，两边积分可得结果； 当0)(0=y N 时，0y y =为原方程的解，当0(0=） x P 时，0x x =为原方程的解。 例1.2、0)1()1(2 2 =-+-dy x y dx y x 解：当0)1)(1(2 2 ≠--y x 时，有 dx x x dy y y 1 12 2-=-两边积分得到 )0(ln 1ln 1ln 22≠=-+-C C y x ，所以有)0()1)(1(22≠=--C C y x ； 当0)1)(1(2 2 =--y x 时，也是原方程的解； 综上所述，原方程的解为)()1)(1(2 2 为常数C C y x =--。 ⑵可化为变量可分离方程的方程：

①、形如)(x y g dx dy = 解法：令x y u =，则udx xdu dy +=，代入得到)(u g u dx du x =+为变量可分离方程，得到 )(0),,(为常数C C x u f =再把u 代入得到)(0),,(为常数C C x x y f =。 ②、形如)0(),(≠+=ab by ax G dx dy 解法：令by ax u +=，则b du adx dy +=，代入得到)(1u G b a dx du b =+为变量可分离方程， 得到)(0 ),,(为常数C C x u f =再把u 代入得到)(0),,(为常数C C x by ax f =+。 ③、形如 )(2 221 11c y b x a c y b x a f dx dy ++++= 解法：01、 02 2 11=b a b a ，转化为 )(by ax G dx dy +=，下同①； 0 2、 022 1 1≠b a b a ，?? ?=++=++00 222111c y b x a c y b x a 的解为),(00y x ，令???-=-=00y y v x x u 得到，)()( )(221 12211u v g u v b a u v b a f v b u a v b u a f du dv =++=++=，下同②； 还有几类：xy u dy xy xg dx xy yf ==+,0)()( xy v xy f dx dy x ==),(2 22),(x y w x y xf dx dy == θθsin ,cos ,0))(,())(,(r y r x ydx xdy y x N ydy xdx y x M ===-++ 以上都可以化为变量可分离方程。 例2.1、 2 5 --+-=y x y x dx dy 解：令2--=y x u ，则du dx dy -=，代入得到u u dx du 71+=- ，有dx udu 7-= 所以)(72 2 为常数C C x u +-=，把u 代入得到)(72 22 为常数） （C C x y x =+--。 例2.2、 1 212+-+-=y x y x dx dy

全微分方程及积分因子

全微分方程及积分因子

全微分方程及积分因子 内容：凑微分法，全微分方程的判别式，全微分方程的公式解，积分因子的微分方程，只含一个变量的积分因子和其他特殊形式的积分因子。由于有数学分析多元微积分的基础，本节的定理1可以简化处理。对课本中第三块知识即全微分方程的物理背景可以留到后面处理，对第四块知识增解和失解的情况要分散在本章各小节，每次都要重视这个问题。关于初等积分法的局限性可归到学习近似解法时一起讲解。 重点：全微分方程的公式解和积分因子的计算，难点为凑微分法和积分因子的计算。 习题1(1,3,5)，2，3 思考题：讨论其他特殊形式的积分因子。 方程：0),(),(=+dy y x N dx y x M 判定：全微分?x N y M ??≡?? 解法：C dy y x N dx y x M y y x x =+??00),(),(0 初值问题0=C 积分因子：x N y M y M x N ??-??=? ???????-??μμμ1

)(x μ： N x N y M dx d ?? -??=μμ1 )(y μ： M x N y M dy d ??- ??-=μμ1 1．解下列方程： 1）0)(222=-+dy y x xydx 解：x N y M ?? ≡??=x 2 ??=-+x y C dy y xydx 002 )0(2既 C y y x =-3/32 2）0)2(=+---dy xe y dx e y y 解：x N y M ??≡??=y e -- ??=-+-y x y C dy y dx e 00)2(既C y xe y =--2 3）0)1(222=---+dy y x dx y x x 解：x N y M ??≡??=y x --221 ??=---+x y C dy y dx y x x 002)1(2 C y y y x x =-+---+23 232322)(32 )(32 )(32 既C y x x =-+23 2 2)(32 4）0)ln (3 =++dy x y dx x y

第二节 几类简单微分方程及其解法

第二节 几类简单微分方程及其解法 本节将介绍可分离变量的微分方程、齐次方程以及一阶线性微分方程等一阶微分方程的解法. 一阶微分方程是微分方程中最基本的、最常见的一类方程.它的一般形式可表示为： 0)',,(=y y x F 或),('y x F y =， 其中)',,(y y x F 为,,'x y y 的已知函数，),(y x F 为,x y 的已知函数. 一、可分离变量的微分方程 如果一阶微分方程),('y x F y =的等式右端能分解为： )()(),(y g x f y x F =， 即)()('y g x f y = (7.2.1) 则称方程(7.2.1)为可分离变量的微分方程. 设)(y g ≠0，则方程(6.2.1)改写为： dx x f dy y g )() (1=， 上式两边积分，可得 ??=dx x f dy y g )()(1. 上述将微分方程化成分离变量形式求解的方法，称为分离变量法. 注：在分离变量时，未知函数y 的函数和微分要写在等式的左边. 例1 求微分方程)3(2'+=y x y 的通解. 解1： 原方程可改写为)3(2+=y x dx dy . 分离变量，两边积分,得,23 1??=+xdx dy y ,3ln 12c x y +=+即.312-±=+c x e y 记1c e c ±=，则微分方程的通解为 32 -=x ce y （c 为任意常数）. 解2：

原方程可改写为)3(2+=y x dx dy . 分离变量，两边积分,得,23 1??=+xdx dy y ,ln )3ln(2c x y +=+即,3ln 2x c y =+23x ce y =+ 则微分方程的通解为 32 -=x ce y （c 为任意常数）. 注：为了简化运算，规定： （1） 微分方程中出现形为 ?u du 的积分时，可不按不定积分基本积分公式表写成 ln du u c u =+?，而是写成ln du u u =?； （2） 不定积分等式中至少有一个形为?u du 的积分时，任意常数不写成c ，而写成c ln 并放在等式右侧. 例2 求微分方程y xy ='的通解. 解: 分离变量，两边积分, 得 ,dy dx y x =?? c x y ln ln ln += cx ln = 则微分方程的通解为cx y = （c 为任意常数）. 例3 求微分方程dx e x dy x e y y )1(2)1(2+=+的通解. 解: 分离变量，两边积分, 得 dx x x dy e e y y ??+=+2121， c x e y ln )1ln()1ln(2++=+ )1(ln 2x c +=， ).1(12x c e y +=+ 则微分方程的通解为 ]1)1(ln[2-+=x c y （c 为任意常数）. 例4 求微分方程)'('2 y y a xy y +=-的通解.

微分方程数值解法

《微分方程数值解法》 【摘要】自然界与工程技术中的很多现象，可以归结为微分方程定解问题。其中，常微分方程求解是微分方程的重要基础内容。但是，对于许多的微分方程，往往很难得到甚至不存在精确的解析表达式，这时候，数值解提供了一个很好的解决思路。，针对于此，本文对常微分方程数值解法进行了简单研究，主要讨论了一些常用的数值解法，如欧拉法、改进的欧拉法、Runge —Kutta 方法、Adams 预估校正法以及勒让德谱方法等，通过具体的算例，结合MA TLAB 求解画图，初步给出了一般常微分方程数值解法的求解过程。同时，通过对各种方法的误差分析，让大家对各种方法的特点和适用范围有一个直观的感受。 【关键词】 常微分方程 数值解法 MA TLAB 误差分析 引言 在我国高校，《微分方程数值解法》作为对数学基础知识要求较高且应用非常广泛的一门课程，不仅 在数学专业，其他的理工科专业的本科及研究生教育中开设这门课程．近四十年来，《微分方程数值解法》不论在理论上还是在方法上都获得了很大的发展．同时，由于微分方程是描述物理、化学和生物现象的数学模型基础，且它的一些最新应用已经扩展到经济、金融预测、图像处理及其他领域 在实际应用中，通过相应的微分方程模型解决具体问题，采用数值方法求得方程的近似解，使具体问题迎刃而解。 2 欧拉法和改进的欧拉法 2.1 欧拉法 2.1.1 欧拉法介绍 首先，我们考虑如下的一阶常微分方程初值问题 ???==0 0)() ,('y x y y x f y （2--1） 事实上，对于更复杂的常微分方程组或者高阶常微分方程，只需要将x 看做向量，（2--1）就成了一个一阶常微分方程组，而高阶常微分方程也可以通过降阶化成一个一阶常微分方程组。 欧拉方法是解常微分方程初值问题最简单最古老的一种数值方法，其基本思路就是把（2--1）中的导数项'y 用差商逼近，从而将一个微分方程转化为一个代数方程，以便求解。 设在[]b a ,中取等距节点h ，因为在节点n x 点上，由（2--1）可得：

一阶常微分方程解法总结

页脚内容1 第 一 章 一阶微分方程的解法的小结 ⑴、可分离变量的方程： ①、形如 )()(y g x f dx dy = 当0)(≠y g 时，得到 dx x f y g dy )()(=，两边积分即可得到结果； 当0)(0=ηg 时，则0)(η=x y 也是方程的解。 例1.1、xy dx dy = 解：当0≠y 时，有xdx y dy =，两边积分得到)(2ln 2为常数C C x y += 所以)(11212 C x e C C e C y ±==为非零常数且 0=y 显然是原方程的解； 综上所述，原方程的解为)(1212 为常数C e C y x = ②、形如0)()()()(=+dy y Q x P dx y N x M 当0)()(≠y N x P 时，可有dy y N y Q dx x P x M ) ()()()(=，两边积分可得结果； 当0)(0=y N 时，0y y =为原方程的解，当0(0=） x P 时，0x x =为原方程的解。 例1.2、0)1()1(22=-+-dy x y dx y x

页脚内容2 解：当0)1)(1(22≠--y x 时，有dx x x dy y y 1 122-=-两边积分得到 )0(ln 1ln 1ln 22≠=-+-C C y x ，所以有)0()1)(1(22≠=--C C y x ； 当0)1)(1(22=--y x 时，也是原方程的解； 综上所述，原方程的解为)()1)(1(22为常数C C y x =--。 ⑵可化为变量可分离方程的方程： ①、形如)(x y g dx dy = 解法：令x y u = ，则udx xdu dy +=，代入得到)(u g u dx du x =+为变量可分离方程，得到)(0),,(为常数C C x u f =再把u 代入得到)(0),,(为常数C C x x y f =。 ②、形如)0(),(≠+=ab by ax G dx dy 解法：令by ax u +=，则b du adx dy +=，代入得到)(1u G b a dx du b =+为变量可分离方程，得到)(0),,(为常数C C x u f =再把u 代入得到)(0),,(为常数C C x by ax f =+。 ③、形如 )(222111c y b x a c y b x a f dx dy ++++= 解法：01、02211 =b a b a ，转化为)(by ax G dx dy +=，下同①； 02、0221 1 ≠b a b a ，???=++=++00222111c y b x a c y b x a 的解为),(00y x ，令???-=-=00y y v x x u

各种类型的微分方程及其相应解法教程文件

各种类型的微分方程及其相应解法 专业班级：交土01班 姓名：高云 学号：1201110102 微分方程的类型有很多种，解题时先判断微分方程是哪种类型，可以帮助我们更快解题，所以我们有必要归纳整理一下各类型（主要是一阶和二阶）的微分方程及其相应解法。 一、一阶微分方程的解法 1.可分离变量的方程 dx x f dy y g )()(=,或)()(y g x f dx dy = 其特点是可以把变量x 和y 只分别在等式的两边，解法关键是把变量分离后两边积分。 例1.求微分方程ydy dx y xydy dx +=+2的通解. 解 先合并dx 及dy 的各项,得dx y dy x y )1()1(2-=- 设,01,012≠-≠-x y 分离变量得 dx x dy y y 1112-=- 两端积分??-=-dx x dy y y 1112得 ||ln |1|ln |1|ln 2 112C x y +-=- 于是 2212)1(1-±=-x C y 记,21C C ±=则得到题设方程的通解 .)1(122-=-x C y 2.齐次方程 （1）)(x y f dx dy = （2） )(c by ax f dx dy ++=（a ，b 均不等于0） 例2求解微分方程.2222xy y dy y xy x dx -=+- 解 原方程变形为=+--=2222y xy x xy y dx dy ,1222?? ? ??+--??? ??x y x y x y x y 令,x y u =则,dx du x u dx dy +=方程化为,1222u u u u dx du x u +--=+ 分离变量得?? ????-+--??? ??--112212121u u u u ,x dx du = 两边积分得 ,ln ln ln 2 1)2ln(23)1ln(C x u u u +=----

偏微分方程数值解法答案

1. 课本2p 有证明 2. 课本812,p p 有说明 3. 课本1520,p p 有说明 4. Rit2法，设n u 是u 的n 维子空间，12,...n ???是n u 的一组基底，n u 中的任一元素n u 可 表为1n n i i i u c ?==∑ ，则,11 11()(,)(,)(,)(,)22j n n n n n n i j i j j i j j J u a u u f u a c c c f ???=== -=-∑∑是12,...n c c c 的二次函数，(,)(,)i j j i a a ????=，令 () 0n j J u c ?=?，从而得到12,...n c c c 满足1 (,)(,),1,2...n i j i j i a c f j n ???===∑，通过解线性方程组，求的i c ，代入1 n n i i i u c ?==∑， 从而得到近似解n u 的过程称为Rit2法 简而言之，Rit2法：为得到偏微分方程的有穷维解，构造了一个近似解，1 n n i i i u c ?== ∑， 利用,11 11()(,)(,)(,)(,)22j n n n n n n i j i j j i j j J u a u u f u a c c c f ???===-=-∑∑确定i c ，求得近似解n u 的过程 Galerkin 法：为求得1 n n i i i u c ? == ∑形式的近似解，在系数i c 使n u 关于n V u ∈，满足(,)(,) n a u V f V =，对任 意 n V u ∈或（取 ,1j V j n ?=≤≤） 1 (,)(,),1,2...n i j i j i a c f j n ???===∑的情况下确定i c ，从而得到近似解1 n n i i i u c ?==∑的过程称 Galerkin 法为 Rit2-Galerkin 法方程： 1 (,)(,)n i j i j i a c f ???==∑ 5. 有限元法：将偏微分方程转化为变分形式，选定单元的形状，对求解域作剖分，进而构 造基函数或单元形状函数，形成有限元空间，将偏微分方程转化成了有限元方程，利用 有效的有限元方程的解法，给出偏微分方程近似解的过程称为有限元法。 6. 解：对求解区间进行网格剖分，节点01......i n a x x x x b =<<<<=得到相邻节点1,i i x x -

各种类型的微分方程及其相应解法

各种类型的微分方程及其相应解法 专业班级：交土01班 姓名：高云 学号：1201110102 微分方程的类型有很多种，解题时先判断微分方程是哪种类型，可以帮助我们更快解题，所以我们有必要归纳整理一下各类型（主要是一阶和二阶）的微分方程及其相应解法。 一、一阶微分方程的解法 1.可分离变量的方程 dx x f dy y g )()(=,或)()(y g x f dx dy = 其特点是可以把变量x 和y 只分别在等式的两边，解法关键是把变量分离后两边积分。 例1.求微分方程ydy dx y xydy dx +=+2的通解. 解 先合并dx 及dy 的各项,得dx y dy x y )1()1(2-=- 设,01,012≠-≠-x y 分离变量得 dx x dy y y 1112-=- 两端积分??-=-dx x dy y y 1112得 ||ln |1|ln |1|ln 2 112C x y +-=- 于是 2212)1(1-±=-x C y 记,21C C ±=则得到题设方程的通解 .)1(122-=-x C y 2.齐次方程 （1）)(x y f dx dy = （2） )(c by ax f dx dy ++=（a ，b 均不等于0） 例2求解微分方程.2222xy y dy y xy x dx -=+- 解 原方程变形为=+--=2222y xy x xy y dx dy ,1222?? ? ??+--??? ??x y x y x y x y 令,x y u =则,dx du x u dx dy +=方程化为,1222u u u u dx du x u +--=+ 分离变量得?? ????-+--??? ??--112212121u u u u ,x dx du = 两边积分得 ,ln ln ln 2 1)2ln(23)1ln(C x u u u +=----

各类微分方程的解法大全

创作编号：BG7531400019813488897SX 创作者：别如克* 各类微分方程的解法 1.可分离变量的微分方程解法 一般形式:g(y)dy=f(x)dx 直接解得∫g(y)dy=∫f(x)dx 设g(y)及f(x)的原函数依次为G(y)及F(x),则G(y)=F(x)+C为微分方程的隐 式通解 2.齐次方程解法 一般形式:dy/dx=φ(y/x) 令u=y/x则y=xu,dy/dx=u+xdu/dx,所以u+xdu/dx=φ(u),即du/［φ(u)-u］ =dx/x两端积分,得∫du/［φ(u)-u］=∫dx/x 最后用y/x代替u,便得所给齐次方程的通解 3.一阶线性微分方程解法 一般形式:dy/dx+P(x)y=Q(x) 先令Q(x)=0则dy/dx+P(x)y=0解得y=Ce－ ∫P(x)dx,再令y=u e－∫P(x)dx代入原方程 解得u=∫Q(x) e∫P(x)dx dx+C,所以y=e－∫P(x)dx［∫Q(x)e∫P(x)dx dx+C］ 即y=Ce－∫P(x)dx +e－ ∫P(x)dx∫Q(x)e∫P(x)dx dx为一阶线性微分方程的通解 4.可降阶的高阶微分方程解法 ①y(n)=f(x)型的微分方程 y(n)=f(x) y(n-1)= ∫f(x)dx+C1

y(n-2)= ∫［∫f(x)dx+C1］dx+C2 依次类推,接连积分n次,便得方程y(n)=f(x)的含有n个任意常数的通解②y”=f(x,y’) 型的微分方程 令y’=p则y”=p’,所以p’=f(x,p),再求解得p=φ(x,C1) 即dy/dx=φ(x,C1),所以y=∫φ(x,C1)dx+C2 ③y”=f(y,y’) 型的微分方程 令y’=p则y”=pdp/dy,所以pdp/dy=f(y,p),再求解得p=φ(y,C1) 即dy/dx=φ(y,C1),即dy/φ(y,C1)=dx,所以∫dy/φ(y,C1)=x+C2 5.二阶常系数齐次线性微分方程解法 一般形式:y”+py’+qy=0,特征方程r2+pr+q=0 6.二阶常系数非齐次线性微分方程解法 一般形式: y”+py’+qy=f(x) 先求y”+py’+qy=0的通解y0(x),再求y”+py’+qy=f(x)的一个特解y*(x) 则y(x)=y0(x)+y*(x)即为微分方程y”+py’+qy=f(x)的通解 求y”+py’+qy=f(x)特解的方法: ①f(x)=P m(x)eλx型 令y*=x k Q m(x)eλx［k按λ不是特征方程的根,是特征方程的单根或特征方程的重根依次取0,1或2］再代入原方程,确定Q m(x)的m+1个系数 ②f(x)=eλx［Pｌ(x)cosωx+P n(x)sinωx］型

微分方程的分类及其数值解法

微分方程的分类及其数值解法 微分方程的分类： 含有未知函数的导数，如dy/dx=2x 、ds/dt=0.4都是微分方程。 一般的凡是表示未知函数、未知函数的导数与自变量之间的关系的方程，叫做微分方程。未知函数是一元函数的，叫常微分方程；未知函数是多元函数的叫做偏微分方程。微分方程有时也简称方程。 一、常微分方程的数值解法： 1、Euler 法： 00d (,), (1.1)d (), (1.2) y f x y x y x y ?=???=? 001 (),(,),0,1,,1n n n n y y x y y hf x y n N +=??=+=-? (1.4) 其中0,n b a x x nh h N -=+=. 用(1.4)求解(1.1)的方法称为Euler 方法。 后退Euler 公式???+==+++),,(),(111 00n n n n y x hf y y x y y 梯形方法公式 )].,(),([2 111+++++=n n n n n n y x f y x f h y y 改进的Euler 方法11(,),(,),1().2p n n n c n n p n p c y y hf x y y y hf x y y y y ++?=+??=+???=+??? 2、Runge-Kutta 方法： p 阶方法 ： 1()O h -=?总体截断误差局部截断误差 二阶Runge-Kutta 方法 ??? ????++==++=+),,(),,(,2212 1211hk y h x f k y x f k k h k h y y n n n n n n

一阶常微分方程解法总结

第 一 章 一阶微分方程的解法的小结 ⑴、可分离变量的方程： ①、形如 )()(y g x f dx dy = 当0)(≠y g 时，得到 dx x f y g dy )() (=，两边积分即可得到结果； 当0)(0=ηg 时，则0)(η=x y 也是方程的解。 例1.1、 xy dx dy = 解：当0≠y 时，有xdx y dy =，两边积分得到)(2ln 2为常数C C x y += 所以)(112 12 C x e C C e C y ±==为非零常数且 0=y 显然是原方程的解； 综上所述，原方程的解为)(12 12 为常数C e C y x = ②、形如0)()()()(=+dy y Q x P dx y N x M 当0)()(≠y N x P 时，可有 dy y N y Q dx x P x M ) () ()()(=，两边积分可得结果； 当0)(0=y N 时，0y y =为原方程的解，当0(0=） x P 时，0x x =为原方程的解。 例1.2、0)1()1(22=-+-dy x y dx y x 解：当0)1)(1(22≠--y x 时，有 dx x x dy y y 1 122-=-两边积分得到 )0(ln 1ln 1ln 22≠=-+-C C y x ，所以有)0()1)(1(22≠=--C C y x ； 当0)1)(1(22=--y x 时，也是原方程的解； 综上所述，原方程的解为)()1)(1(22为常数C C y x =--。 ⑵可化为变量可分离方程的方程：

①、形如 )(x y g dx dy = 解法：令x y u =，则udx xdu dy +=，代入得到)(u g u dx du x =+为变量可分离方程，得 到)(0),,(为常数C C x u f =再把u 代入得到)(0),,(为常数C C x x y f =。 ②、形如)0(),(≠+=ab by ax G dx dy 解法：令by ax u +=，则b du adx dy +=，代入得到)(1u G b a dx du b =+为变量可分离方程， 得到)(0),,(为常数C C x u f =再把u 代入得到)(0),,(为常数C C x by ax f =+。 ③、形如 )(2 221 11c y b x a c y b x a f dx dy ++++= 解法：0 1、 02 2 11=b a b a ，转化为 )(by ax G dx dy +=，下同①； 02、 022 1 1≠b a b a ，???=++=++00 222111 c y b x a c y b x a 的解为),(00y x ，令???-=-=0 0y y v x x u 得到，)()( )(221 12211u v g u v b a u v b a f v b u a v b u a f du dv =++=++=，下同②； 还有几类：xy u dy xy xg dx xy yf ==+,0)()( 以上都可以化为变量可分离方程。 例2.1、 2 5--+-=y x y x dx dy 解：令2--=y x u ，则du dx dy -=，代入得到u u dx du 7 1+= - ，有dx udu 7-= 所以)(72 2 为常数C C x u +-=，把u 代入得到)(72 22 为常数） （C C x y x =+--。 例2.2、 1 212+-+-=y x y x dx dy 解：由???=+-=+-012012y x y x 得到?????=-=3131y x ，令?? ???-=+=3131y v x u ，有???==du dx dv dy ，代入得到

常微分方程数值解法

i.常微分方程初值问题数值解法 常微分方程初值问题的真解可以看成是从给定初始点出发的一条连续曲线。差分法是常微分方程初值问题的主要数值解法，其目的是得到若干个离散点来逼近这条解曲线。有两个基本途径。一个是用离散点上的差商近似替代微商。另一个是先对微分方程积分得到积分方程，再利用离散点作数值积分。 i.1 常微分方程差分法 考虑常微分方程初值问题：求函数()u t 满足 (,), 0du f t u t T dt =<≤ （i.1a ） 0(0)u u = (i.1b) 其中(,)f t u 是定义在区域G : 0t T ≤≤, u <∞上的连续函数，0u 和T 是给定的常数。我们假设(,)f t u 对u 满足Lipschitz 条件，即存在常数L 使得 121212(,)(,), [0,]; ,(,)f t u f t u L u u t T u u -≤-?∈∈-∞∞ （i.2） 这一条件保证了（i.1）的解是适定的，即存在，唯一，而且连续依赖于初值0u 。 通常情况下，（i.1）的精确解不可能用简单的解析表达式给出，只能求近似解。本章讨论常微分方程最常用的近似数值解法-差分方法。先来讨论最简单的Euler 法。为此，首先将求解区域[0,]T 离散化为若干个离散点： 0110N N t t t t T -=<< <<= （i.3） 其中n t hn =，0h >称为步长。 在微积分课程中我们熟知，微商（即导数）是差商的极限。反过来，差商就是微商的近似。在0t t =处，在（i.1a ）中用向前差商 10()()u t u t h -代替微商du dt ，便得 10000()()(,())u t u t hf t u t ε=++ 如果忽略误差项0ε，再换个记号，用i u 代替()i u t 便得到 1000(,)u u hf t u -= 一般地，我们有 1Euler (,), 0,1, ,1n n n n u u hf t u n N +=+=-方法： （i.4） 从(i.1b) 给出的初始值0u 出发，由上式可以依次算出1,,N t t 上的差分解1,,N u u 。

[整理]一阶微分方程解的存在定理.

第三章 一阶微分方程解的存在定理 [教学目标] 1. 理解解的存在唯一性定理的条件、结论及证明思路，掌握逐次逼近法，熟练近似解的误差估计式。 2. 了解解的延拓定理及延拓条件。 3. 理解解对初值的连续性、可微性定理的条件和结论。 [教学重难点] 解的存在唯一性定理的证明，解对初值的连续性、可微性定理的证明。 [教学方法] 讲授，实践。 [教学时间] 12学时 [教学内容] 解的存在唯一性定理的条件、结论及证明思路，解的延拓概念及延拓条件，解对初值的连续性、可微性定理及其证明。 [考核目标] 1.理解解的存在唯一性定理的条件、结论，能用逐次逼近法解简单的问题。 2.熟练近似解的误差估计式，解对初值的连续性及可微性公式。 3.利用解的存在唯一性定理、解的延拓定理及延拓条件能证明有关方程的某些性质。 §1 解的存在性唯一性定理和逐步逼近法 微分方程来源于生产实践际，研究微分方程的目的就在于掌握它所反映的客观规律，能动解释所出现的各种现象并预测未来的可能情况。在第二章介绍了一阶微分方程初等解法的几种类型，但是，大量的一阶方程一般是不能用初等解法求出其通解。而实际问题中所需要的往往是要求满足某种初始条件的解。因此初值问题的研究就显得十分重要，从前面我们也了解到初值问题的解不一定是唯一的。他必须满足一定的条件才能保证初值问题解的存在性与唯一性，而讨论初值问题解的存在性与唯一性在常微分方程占有很重要的地位，是近代常微分方程定性理论，稳定性理论以及其他理论的基础。 例如方程 dy dx =过点(0,0)的解就是不唯一，易知0y =是方程过(0,0)的解，此外，容易验证，2 y x =或更一般地，函数 2 0 0() c<1x c y x c x ≤≤?=?-≤? 都是方程过点(0,0)而且定义在区间01x ≤≤上的解，其中c 是满足01c <<的任一数。 解的存在唯一性定理能够很好地解释上述问题，它明确地肯定了方程的解在一定条件下的存在性 和唯一性。另外，由于能得到精确解的微分方程为数不多，微分方程的近似解法具有重要的意义，而解的存在唯一性是进行近似计算的前提，如果解本身不存在，而近似求解就失去意义；如果存在不唯一，不能确定所求的是哪个解。而解的存在唯一性定理保证了所求解的存在性和唯一性。 1．存在性与唯一性定理： （1）显式一阶微分方程 ),(y x f dx dy = （3.1） 这里),(y x f 是在矩形域：00:||,||R x x a y y b -≤-≤ （3.2）

一阶常微分方程的解法

一阶常微分方程的解法 摘要：常微分方程是微积分学的重要组成部分，广泛用于具体问题的研究中，在整个数学中占有重要的地位。本文对一阶常微分方程的解法作了简要的总结，并举例加以分析了变量可分离方程，线性微分方程，积分因子，恰当微分方程，主要归纳了一阶微分方程的初等解法，并以典型例题加以说明。 关键词：变量分离;积分因子;非齐次微分方程;常数变易法 Solution of first-order differential equation Abstract: Differential equations, important parts of calculus, are widely used in the research of practical problems, which also play important role in mathematics. The solution of a differential equation is summarized briefly, and illustrates the analysis of variable separable equation, linear differential equation, integral factor, exact differential equation, mainly summarizes the elementary solution of first order differential equations, and the typical examples to illustrate. Keywords: variable separation; integral factor; non-homogeneous differential equation; constant variation method 1. 引言 一阶常微分方程初等解法,就是把常微分方程的求解问题转化为积分问题, 能用这种方法求解的微分方程称为可积方程. 本文通过对一阶微分方程的初等解法的归纳与总结，以及对变量分离，积分因子，微分方程等各类初等解法的简要分析，同时结合例题把常微分方程的求解问题化为积分问题，进行求解. 2. 一般变量分离 2.1 变量可分离方程 形如 ()()dy f x g y dx = （1.1） 或 1122()()()()M x N y dx M x N y dy = （1.2） 的方程，称为变量可分离方程。分别称（1.1）、（1.2）为显式变量可分离方程和 微分形式变量可分离方程[1] . (1) 显式变量可分离方程的解法 在方程（1.1）中， 若()0g y ≠，（1.1）变形为 ()() dy f x dx g y =

微分方程数值解法答案

包括基本概念，差分格式的构造、截断误差和稳定性，这些内容是贯穿整个教材的主线。解答问题关键在过程，能够显示出你已经掌握了书上的内容，知道了解题方法。这次考试题目的类型：20分的选择题，主要是基本概念的理解，后面有五个大题，包括差分格式的构造、截断误差和稳定性。 习题一 1． 略 2． y y x f -=),(，梯形公式：n n n n n n y h h y y y h y y )121(),(2111+-+=+- =+++，所以0122)1(01])121[()121()121(y h h y h h y h h y h h n h h n n n +--+--+-+=+-+==+-+= ，当0→h 时， x n e y -→。 同理可以证明预报-校正法收敛到微分方程的解. 3． 局部截断误差的推导同欧拉公式; 整体截断误差: ? ++++++-++≤1 ),())(,(11111n n x x n n n n n n n dx y x f x y x f R εε 11)(++-++≤n n n y x y Lh R ε，这里R R n ≤ 而111)(+++-=n n n y x y ε，所以 R Lh n n += -+εε1)1(，不妨设1

微分方程数值解――

微分方程数值解―― 第二章 习题 1． 设)('x f 为)(x f 的一阶广义导数，试用类似办法定义)(x f 的k 阶广义导数) () (x f k （ ,2,1=k ）。 解：对一维情形，函数的广义导数是通过分部积分来定义的。 我们知，)(x f 的一阶广义导数位)(x g ，如果满足 dx x x f dx x x g b a b a )()()()('?? -=?? 类似的，)(x f 的k 阶广义导数为)()() (x f x g k =，如果有 dx x x f dx x x g b a k k b a )()()1()()()(?? -=?? 2． 试建立与边值问题 ?????====<<=+=) 2.1(0)()(,0)()() 1.1(,''44b u b u a u a u b x a f u dx u d Lu 等价的变分问题。 证明： 设}0)()(,0)()(),(|{' '2====∈=b v a v b v a v I H v v V 对方程)1.1(两边同乘以v ，再关于x 在),(b a 上积分)(V v ∈，得 ??=+b a b a fvdx vdx u dx u d )(44 其中 dx dx dv dx u d dx dx dv dx u d dx u d v dx u d d v vdx dx u d b a b a b a b a b a ???? -=-==33 33333344|)( dx dx v d dx u d dx dv dx u d dx u d d dx dv b a b a b a ??+-=-=22222222|)( dx dx v d dx u d b a ? = 2 222 (*) 记dx uv dx v d dx u d v u a b a ?+=)(),(2 222，?=b a fvdx v f ),(。于是我们得到以下等价变分问题的提法：

全微分方程的不定积分解法及其证明

全微分方程的不定积分解法及其证明 一个一阶微分方程写成 P (x,y ) dx + Q (x,y ) dy = 0 ⑴ 形式后，如果它的左端恰好是某一个函数u= u (x,y ) 的全微分： du (x,y ) = P (x,y ) dx + Q (x,y ) dy 那么方程⑴就叫做全微分方程。这里 5u 5x = P (x,y ）， 5u 5y = Q (x,y ) 方程⑴就是du (x,y ) = 0,其通解为： u (x,y ) = C(C 为常数） 可见，解全微分方程的关键在于求原函数u (x,y ）。因此，本文将提供一种求原函数u (x,y ) 的简捷 方法，并给出证明。 1引入记号 为了表述方便，先引入记号如下： 设M (x,y ) 为一个含有变量x,y 项的二元函数，定义： ⑴“M (x q ,y ) ”表示M (x,y ) 减去它里面含有变量x 的项; ⑵“M (x,y q ）”表示M (x,y ) 减去它里面含有变量y 的项； 注意：常数项看作既不含变量x 也不含变量y 的项。 现举一例如下： 设：M (x,y ) = xy + x ey+ x 1- x + sinx+ co sx co sy + y 2+ 1 按记号定义有： M (x q ,y ) = M (x,y ) - (x y + x ey + x 1 - x + sinx + co sx co sy ) = y 2 + 1 M (x,y q ）= M (x,y ) - (x y + x ey + co sx co sy + y 2) = x 1 - x

+ sinx + 1 2u (x,y ) 的简捷求法 引理设开区域G 是一个单连通域，函数P (x,y ),Q (x,y ) 在G 内具有一阶连续偏导数，则 P (x,y ) dx + Q (x,y ) dy 在G 内为某一函数u (x,y ) 的全微分的充分必要条件是等式 5P 5y = 5Q 5x

各类微分方程的解法大全

各类微分方程的解法 1.可分离变量的微分方程解法 一般形式:g(y)dy=f(x)dx 直接解得∫g(y)dy=∫f(x)dx 设g(y)及f(x)的原函数依次为G(y)及F(x),则G(y)=F(x)+C为微分方程的隐式通解 2.齐次方程解法 一般形式:dy/dx=φ(y/x) 令u=y/x则y=xu,dy/dx=u+xdu/dx,所以u+xdu/dx=φ(u),即du/［φ(u)-u］=dx/x 两端积分,得∫du/［φ(u)-u］=∫dx/x 最后用y/x代替u,便得所给齐次方程的通解 3.一阶线性微分方程解法 一般形式:dy/dx+P(x)y=Q(x) 先令Q(x)=0则dy/dx+P(x)y=0解得y=Ce－ ∫P(x)dx,再令y=u e－∫P(x)dx代入原方程解得u=∫Q(x) e∫P(x)dx dx+C,所以y=e－∫P(x)dx［∫Q(x)e∫P(x)dx dx+C］ 即y=Ce－∫P(x)dx +e－ ∫P(x)dx∫Q(x)e∫P(x)dx dx为一阶线性微分方程的通解 4.可降阶的高阶微分方程解法 ①y(n)=f(x)型的微分方程 y(n)=f(x) y(n-1)= ∫f(x)dx+C1 y(n-2)= ∫［∫f(x)dx+C1］dx+C2 依次类推,接连积分n次,便得方程y(n)=f(x)的含有n个任意常数的通解②y”=f(x,y’) 型的微分方程 令y’=p则y”=p’,所以p’=f(x,p),再求解得p=φ(x,C1) 即dy/dx=φ(x,C1),所以y=∫φ(x,C1)dx+C2 ③y”=f(y,y’) 型的微分方程

令y ’=p 则y ”=pdp/dy,所以pdp/dy=f(y,p),再求解得p=φ(y,C 1) 即dy/dx=φ(y,C 1),即dy/φ(y,C 1)=dx,所以∫dy/φ(y,C 1)=x+C 2 5.二阶常系数齐次线性微分方程解法 一般形式:y ”+py ’+qy=0,特征方程r 2+pr+q=0 6.二阶常系数非齐次线性微分方程解法 一般形式: y ”+py ’+qy=f(x) 先求y ”+py ’+qy=0的通解y 0(x),再求y ”+py ’+qy=f(x)的一个特解y*(x) 则y(x)=y 0(x)+y*(x)即为微分方程y ”+py ’+qy=f(x)的通解 求y ”+py ’+qy=f(x)特解的方法: ① f(x)=P m (x)e λx 型 令y*=x k Q m (x)e λx ［k 按λ不是特征方程的根,是特征方程的单根或特征方程的重根依次取0,1或2］再代入原方程,确定Q m (x)的m+1个系数 ② f(x)=e λx ［P ｌ(x)cos ωx+P n (x)sin ωx ］型 令y*=x k e λx ［Q m (x)cos ωx+R m (x)sin ωx ］［m=max ﹛ｌ,n ﹜,k 按λ+i ω不是特征方程的根或是特征方程的单根依次取0或1］再代入原方程,分别确定Q m (x)和R m (x)的m+1个系数