2012年微分方程数值解考试题_A_

1 / 2

北京邮电大学2011——2012学年第2学期

《微分方程数值解》期末考试试题A 卷

1.（15分）写出求解常微分方程初值问题0

00(,),,

()du

f t u t T dt u u ?=<≤???=?的欧拉（Euler）

格式，证明欧拉（Euler）法的局部截断误差的阶为2()O h .

2.（15分）若单步法1(,,)n n n n u u h t u h ?+=+中(,,)t u h ? 0(,t t T ≤≤00,h h ≤≤ (,))u ∈?∞∞关于u 满足Lipschitz 条件，即：

1212

(,,)(,,)t u h t u h L u u ?≤????

其中L ?为与u，t 无关的常数。证明该算法稳定。

3.（15分）用待定系数法确定求解常微分方程初值问题0

00(,),,

()du

f t u t T dt u u ?=<≤???=?的

四步四阶显式格式1123555937924

()n n n n n n h

u u f f f f +???=+?+?.

4.（15分）对于两点边值问题2222422101310122(),(,)

(),().d u

x u x x dx

u u e ??++=+∈??

??==+??

1）推导上微分方程的中心差分格式。

2）假设求解区间上共有5个分点，即N＝4，写出求解该问题的中心差分格

式的矩阵形式。

5. (15分）推导初边值问题22000000(,),,(,)(),(,)(,),u u

a f x t t T t x u x x x l u t u l t t T φ???=+<≤?????

的向前差分格式和向

后差分格式。

2 / 2

6. (15分）用有限体积法构造逼近方程

()u u k u k k f x x y y ??

???????????=?+=????????????????

第一边值问题|(,)u x y αΓ=的五点差分格式，这里0min (,)k k x y k =≥>. 7.(10分) 学习完本课程后，你最大的收获是什么？请结合实际谈谈你对本课程的看法以及将来学习的展望。

微分方程数值解习题

习题2 1．略 2．略 3．略 4．差分格式写成矩阵形式为： n n M n M n n n M n M n n e u u u u r t r r r t r r r t r r r t u u u u +????? ? ?? ?????????? ? ?-?--?--?--?-=???????? ??--+-+-++12211221121212121M O O O M αβαααβαααβαααβ 矩阵的特征值为：)cos(221M j r r t j π ααβλ+-?-=，要使格式稳定，则特征值须满足 t c j ?+≤1λ，即2 1≤r α 5．利用泰勒展式可以得到古典隐式差分格式的截断误差为)(2 h t O +?。古典隐式差分格式写成矩阵形式为： n n M n M n n n M n M n n e u u u u u u u u t r r r t r r r t r r r t r +??????? ? ?????????? ? ?=? ? ??? ? ? ? ?????????? ???++--?++--?++--?++--+-+ -++122 112211111121212121M O M O O O βαααβαααβαααβα 特征值为： 1 ))cos( 221(--+?++=M j r r t j πααβλ，即： )(1))2( cos 41(1 2t o M j r t j ?+≤+?++=-παβλ，所以无条件稳定。 6．由Von-Neumann 方法，令mh i n l n m e u β?=，代入差分格式得到增长因子为： )2 ( sin 41),(2h r i t G βωβ-=?，所以1)]2 ( sin 4[1),(22≥+=?h r t G βωβ，恒不稳定。 7． n m n m u v =+1，则原三层格式等价于： ??=-+=+--+++-+++n m n m n m n m n m n m n m n m u v v u u u u r u 111111)21()2()1(θθθ，令mh i n l n l n m n m e v u βη???? ? ??=???? ??，

常微分方程期中考试题

常微分方程期中测试试卷(1) 一、填空 1 微分方程 ) (2 2= + - +x y dx dy dx dy n 的阶数是____________ 2 若 ) , (y x M和) , (y x N在矩形区域R内是) , (y x的连续函数,且有连续的一阶偏导数,则方程 ) , ( ) , (= +dy y x N dx y x M有只与y有关的积分因子的充要条件是 _________________________ 3 _________________________________________ 称为齐次方程. 4 如果 ) , (y x f___________________________________________ ,则 ) , (y x f dx dy = 存在唯一的解 ) (x y? =,定义于区间h x x≤ - 0上,连续且满足初始条件 ) ( x y? = ,其中 = h_______________________ . 5 对于任意的 ) , ( 1 y x,) , ( 2 y x R ∈ (R为某一矩形区域),若存在常数)0 (> N N使 ______________________ ,则称 ) , (y x f在R上关于y满足利普希兹条件. 6 方程 2 2y x dx dy + = 定义在矩形区域R:2 2 ,2 2≤ ≤ - ≤ ≤ -y x上 ,则经过点)0,0(的解的存在区间是 ___________________ 7 若 ) ,..... 2,1 )( (n i t x i = 是齐次线性方程的n个解,)(t w为其伏朗斯基行列式,则)(t w满足一阶线性方程 ___________________________________ 8若 ) ,..... 2,1 )( (n i t x i = 为齐次线性方程的一个基本解组, )(t x为非齐次线性方程的一个特解,则非齐次线性方程的所有解可表为 _________________________ 9若 ) (x ?为毕卡逼近序列{})(x n?的极限，则有≤ -) ( ) (x x n ? ? __________________ 10 _________________________________________ 称为黎卡提方程，若它有一个特解 ) (x y，则经过变换___________________ ，可化为伯努利方程．二求下列方程的解１ 3 y x y dx dy + = ２求方程 2 y x dx dy + = 经过 )0,0(的第三次近似解３讨论方程 2 y dx dy = ， 1 )1(= y的解的存在区间 4 求方程 1 ) (2 2= - +y dx dy 的奇解

常微分方程知识点总结

常微分方程知识点总结常微分方程知识点你学得怎么样呢？下面是的常微分方程知识点总结，欢迎大家阅读！微分方程的概念方程对于学过中学数学的人来说是比较熟悉的；在初等数学中就有各种各样的方程，比如线性方程、二次方程、高次方程、指数方程、对数方程、三角方程和方程组等等。这些方程都是要把研究的问题中的已知数和数之间的关系找出来，列出包含一个数或几个数的一个或者多个方程式，然后取求方程的解。但是在实际工作中，常常出现一些特点和以上方程完全不同的问题。比如：物质在一定条件下的运动变化，要寻求它的运动、变化的规律；某个物体在重力作用下自由下落，要寻求下落距离随时间变化的规律；火箭在发动机推动下在空间飞行，要寻求它飞行的轨道，等等。物质运动和它的变化规律在数学上是用函数关系来描述的，因此，这类问题就是要去寻求满足某些条件的一个或者几个函数。也就是说，凡是这类问题都不是简单地去求一个或者几个固定不变的数值，而是要求一个或者几个的函数。解这类问题的基本思想和初等数学解方程的基本思想很相似，也是要把研究的问题中已知函数和函数之间的关系找出来，从列出的包含函数的一个或几个方程中去求得函数的表达式。但是无论在方程

的形式、求解的具体方法、求出解的性质等方面，都和初等数学中的解方程有许多不同的地方。在数学上，解这类方程，要用到微分和导数的知识。因此，凡是表示函数的导数以及自变量之间的关系的方程，就叫做微分方程。微分方程差不多是和微积分同时先后产生的，苏格兰数学家耐普尔创立对数的时候，就讨论过微分方程的近似解。牛顿在建立微积分的同时，对简单的微分方程用级数来求解。后来瑞士数学家雅各布?贝努利、欧拉、法国数学家克雷洛、达朗贝尔、拉格朗日等人又不断地研究和丰富了微分方程的理论。常微分方程的形成与发展是和力学、天文学、物理学，以及其他科学技术的发展密切相关的。数学的其他分支的新发展，如复变函数、李群、组合拓扑学等，都对常微分方程的发展产生了深刻的影响，当前计算机的发展更是为常微分方程的应用及理论研究提供了非常有力的工具。牛顿研究天体力学和机械力学的时候，利用了微分方程这个工具，从理论上得到了行星运动规律。后来，法国天文学家勒维烈和英国天文学家亚当斯使用微分方程各自计算出那时尚未发现的海王星的位置。这些都使数学家更加深信微分方程在认识自然、改造自然方面的巨大力量。微分方程的理论逐步完善的时候，利用它就可以精确地表述事物变化所遵循的基本规律，只要列出相应的微分方程，有了解方程的方法。微分方程也就成了最有生命力的数学分支。

偏微分方程数值解复习题(2011硕士)

偏微分方程数值解期末复习（2011硕士）一、考题类型本次试卷共六道题目，题型及其所占比例分别为：填空题20%；计算题80% 二、按章节复习内容第一章知识点：Euler法、向前差商、向后差商、中心差商、局部截断误差、整体截断误差、相容性、收敛性、阶、稳定性、显格式、隐格式、线性多步法、第一特征多项式、第二特征多项式、稳定多项式、绝对稳定等；要求: 会辨认差分格式, 判断线性多步法的误差和阶; 第二章知识点：矩形网格、(正则,非正则)内点、边界点、偏向前(向后,中心)差商、五点差分格式、增设虚点法、积分插值法、线性椭圆型差分格式、极值原理、比较定理、五点差分格式的相容收敛和、稳定性等；要求: 建立椭圆型方程边值问题的差分格式, 极值原理; 第四章知识点：最简显格式、最简隐格式、CN格式、双层加权格式、Richardson 格式、网格比、传播因子法(分离变量法) 、传播因子、传播矩阵、谱半径、von Neumann条件、跳点格式、ADI格式、线性椭圆型差分格式、极值原理、比较定理、五点差分格式的相容收敛和稳定性等；要求: 建立抛物型方程边值问题的差分格式, 计算局部截断误差; 第五章知识点：左偏心格式、右偏心格式、中心格式、LF格式、LW格式、Wendroff 格式、跳蛙格式、特征线、CFL条件等；要求: 建立双曲型方程边值问题的差分格式, 计算局部截断误差; 第七章要求: 会用线性元(线性基)建立常微分方程边值问题的有限元格式

三练习题 1、已知显格式21131()22 n n n n u u h f f +++-=-，试证明格式是相容的，并求它的阶。 P39+P41 2、用Taylor 展开原理构造一元函数一阶导数和二阶导数的数值微分公式。提示：向前、向后和中心差商与一阶导数间关系，二阶中心差商与二阶导数之间的关系课件 3、用数值微分方法或数值积分方法建立椭圆型方程 2222(,),(,),u u f x y x y x y ??--=?∈Ω?? :01,01x y Ω≤≤≤≤ 内点差分格式。 P75+课件 4、构造椭圆型方程边值问题的差分格式. P101 (4)题 5、构建一维热传导方程220,(0)u u Lu a a t x ??=-=>??的数值差分格式(显隐格式等)。参考P132-135相关知识点 6、设有逼近热传导方程22(0)u u Lu a f a const t x ??≡-==>??的带权双层格式 ()()1111111122(1)2k k j j k k k k k k j j j j j j u u a u u u u u u h θθτ++++-+-+-??=-++--+?? 其中[0,1]θ∈，试求其截断误差。并证明当2 1212h a θτ=-时，截断误差的阶最高阶为24()O h τ+。 P135+P165+课件 7、传播因子法证明抛物型方程22(0)u u Lu a f a const t x ??≡-==>??的最简显隐和六点CN 格式稳定性。 P156+课件 8、对一阶常系数双曲型方程的初边值问题 0,0,0,0,(,0)(),0,(0,)(),0, u u a t T x a t x u x x x u t t t T φψ???+=<≤<<∞>?????=≤<∞??=≤≤?

《常微分方程》期末模拟试题

《常微分方程》模拟练习题及参考答案一、填空题（每个空格4分，共80分） 1、n 阶线性齐次微分方程基本解组中解的个数恰好是 n 个。 2、一阶微分方程 2=dy x dx 的通解为 2=+y x C （C 为任意常数），方程与通过点（2，3）的特解为 2 1=-y x ，与直线y=2x+3相切的解是 2 4=+y x ，满足条件3 3ydx =?的解为 22=-y x 。 3、李普希兹条件是保证一阶微分方程初值问题解惟一的必要条件。 4、对方程 2()dy x y dx =+作变换 =+u x y ，可将其化为变量可分离方程，其通解为 tan()=+-y x C x 。 5、方程过点共有无数个解。 6、方程 ''2 1=-y x 的通解为 42 12122=-++x x y C x C ，满足初始条件13|2,|5====x x y y 的特解为 4219 12264 =-++x x y x 。 7、方程无奇解。 8、微分方程2260--=d y dy y dx dx 可化为一阶线性微分方程组 6?=??? ?=+??dy z dx dz z y dx 。 9、方程的奇解是 y=0 。 10、35323+=d y dy x dx dx 是 3 阶常微分方程。 11、方程 22dy x y dx =+满足解得存在唯一性定理条件的区域是 xoy 平面。 12、微分方程22450d y dy y dx dx --=通解为 512-=+x x y C e C e ，该方程可化为一阶线性微分方程组 45?=??? ?=+??dy z dx dz z y dx 。 2 1d d y x y -=)1,2 (πx x y x y +-=d d y x y =d d

(完整版)常微分方程的大致知识点

= + ?x = + ?x = + ?x 常微分方程的大致知识点（一）初等积分法 1、线素场与等倾线 2、可分离变量方程 3、齐次方程（一般含有 x 或 y 的项） y x 4、一阶线性非齐次方程常数变易法，或 y = e ? a ( x )dx [? b (x )e -? a ( x )dx dx + C ] 5、伯努力方程令 z = y 1-n ，则 dz = (1 - n ) y -n dy ，可将伯努力方程化成一阶线性非齐次或一阶线性齐次 dx 6、全微分方程若?M ?y 若 ?M ?y dx = ?N ，则u (x , y ) = C ，（留意书上公式） ?x ≠ ?N ，则找积分因子，（留意书上公式） ?x f (x f ( y , （二）毕卡序列 x y 1 y 0 0 x f (x , y 0 )dx ， y 2 y 0 0 x f (x , y 1 )dx ， y 3 y 0 0 f (x , y 2 )dx ，其余类推（三）常系数方程 1、常系数齐次L (D ) y = 0 方法：特征方程 7、可降阶的二阶微分方程 d 2 y = , dy ) ，令 dy = d 2 y p ,则 = dy dx 2 d 2 y = dx dy ) ，令 dx dy = p ,则 dx 2 d 2 y dx = p dp dx 2 dx dx dx 2 dy 8、正交轨线族

? ? dy 单的实根, ， y = C e 1x + C e 2 x 1 2 1 2 单的复根1, 2 = ± i ， y = e x (C cos x + C 2 sin x ) 重的实根 = = ， y = (C + C x )e x 1 2 1 2 重的复根1, 2 = ± i ，3, 4 = ± i ， y = e x [(C + C 2 x ) c os x + (C 3 + C 4 x ) sin x ] 2、常系数非齐次L (D ) y = 方法：三部曲。 f (x ) 第一步求L (D ) y = 0 的通解Y 第二步求L (D ) y = f (x ) 的特解 y * 第三步求L (D ) y = f (x ) 的通解 y = Y + y * 如何求 y * ？当 f (x ) = P m (x )e x 时， y * = x k Q (x )e x 当 f (x ) = P m (x )e ux cos vx + Q (x )e ux sin vx 时， y * = x k e ux (R (x ) cos vx + S m (x ) sin vx ) 当 f (x ) 是一般形式时， y * = ? x W (x ,) f ()d ，其中 W(.)是郎斯基行列式 x 0 W () （四）常系数方程组方法：三部曲。第一步求 dX dt = A (t ) X 的通解， Φ(t )C 。利用特征方程 A - I = 0 ，并分情况讨论。第二步求 dX dt 第三步求 dX dt = A (t ) X + f (t ) 的特解， Φ(t )?Φ-1 (s ) f (s )ds ，（定积分与不定积分等价） = A (t ) X + f (t ) 的通解， Φ(t )C + Φ(t )?Φ-1 (s ) f (s )ds （五）奇点与极限环 ? dx = ax + b y dt ? ? = cx + dy 1、分析方程组? dt 的奇点的性质，用特征方程： A - I = 0 特征方程的根有 3 种情况：相异实根、相异复根、相同实根。第一种情况：相异实根，1 ≠ 2 1 1 m m m

常微分方程初值问题的数值解法

第七章常微分方程初值问题的数值解法 --------学习小结一、本章学习体会通过本章的学习，我了解了常微分方程初值问题的计算方法，对于解决那些很难求解出解析表达式的，甚至有解析表达式但是解不出具体的值的常微分方程非常有用。在这一章里求解常微分方程的基本思想是将初值问题进行离散化，然后进行迭代求解。在这里将初值问题离散化的方法有三种，分别是差商代替导数的方法、Taylor 级数法和数值积分法。常微分方程初值问题的数值解法的分类有显示方法和隐式方法，或者可以分为单步法和多步法。在这里单步法是指计算第n+1个y 的值时,只用到前一步的值，而多步法则是指计算第n+1个y 的值时，用到了前几步的值。通过对本章的学习，已经能熟练掌握如何用Taylor 级数法去求解单步法中各方法的公式和截断误差，但是对线性多步法的求解理解不怎么透切，特别是计算过程较复杂的推理。在本章的学习过程中还遇到不少问题，比如本章知识点多，公式多，在做题时容易混淆，其次对几种R-K 公式的理解不够透彻，处理一个实际问题时，不知道选取哪一种公式，通过课本里面几种方法的计算比较得知其误差并不一样，，这个还需要自己在往后的实际应用中多多实践留意并总结。二、本章知识梳理常微分方程初值问题的数值解法一般概念步长h ，取节点0,(0,1,...,)n t t nh n M =+=，且M t T ≤，则初值问题000 '(,),()y f t y t t T y t y =≤≤?? =?的数值解法的一般形式是 1(,,,...,,)0,(0,1,...,)n n n n k F t y y y h n M k ++==-

常微分方程期末考试练习题及答案

一，常微分方程的基本概念常微分方程：含一个自变量x，未知数y及若干阶导数的方程式。一般形式为：F（x，y，y，.....y(n)）=0 (n≠0). 1. 常微分方程中包含未知函数最高阶导数的阶数称为该方程的阶。如：f(x)（3）+3f(x)+x=f(x)为3阶方程。 2.若f（x）使常微分方程两端恒等，则f（x）称为常微分方程的解。 3.含有独立的任意个常数（个数等于方程的阶数）的方程的解称为常微分方程的通解。如常系数三阶微分方程F（t，x（3））=0的通解的形式为：x（t）=c1x（t）+c2x（t）+c3x（t）。 4.满足初值条件的解称为它的特解（特解不唯一，亦可能不存在）。 5.常微分方程之线性及非线性：对于F（x，y，y，......y(n)）=0而言，如果方程之左端是y，y，......y(n)的一次有理式，则次方程为n阶线性微分方程。（方程线性与否与自变量无关）。如：xy（2）-5y，+3xy=sinx 为2阶线性微分方程；y（2）+siny=0为非线性微分方程。注：a.这里主要介绍几个主要的，常用的常微分方程的基本概念。余者如常微分方程之显隐式解，初值条件，初值问题等概念这里予以略去。另外，有兴趣的同学不妨看一下教材23页的雅可比矩阵。 b.教材28页第八题不妨做做。二.可分离变量的方程 A.变量分离方程

1.定义：形如 dx dy =f （x)φ(y)的方程，称为分离变量方程。这里f （x ），φ（x ）分别是x ，y 的连续函数。 2.解法：分离变量法? ? +=c dx x f y dy )()(?. （*）说明： a 由于（*）是建立在φ（y ）≠0的基础上，故而可能漏解。需视情况补上φ（y ）=0的特解。（有时候特解也可以和通解统一于一式中） b.不需考虑因自变量引起的分母为零的情况。例1.0)4(2=-+dy x x ydx 解：由题意分离变量得：04 2=+-y dy x dx 即： 0)141(41=+--y dy dx x x 积分之，得：c y x x =+--ln )ln 4(ln 4 1 故原方程通解为：cx y x =-4)4( （c 为任意常数），特解y=0包含在通解中（即两者统一于一式中）。 *例2.若连续函数f （x ）满足 2 ln )2 ()(20 +=? dt t f x f x ，则f （x ）是？解：对给定的积分方程两边关于x 求导，得： )(2)('x f x f = （变上限求积分求导）分离变量，解之得：x Ce x f 2)(= 由原方程知： f （0）=ln2，代入上解析式得： C=ln2， B.可化为分离变量方程的类型。解决数学题目有一个显而易见的思想：即把遇到的新问题，结合已知

常微分方程的大致知识点

常微分方程的大致知识点Last revision on 21 December 2020

常微分方程的大致知识点（一）初等积分法 1、线素场与等倾线 2、可分离变量方程 3、齐次方程（一般含有x y y x 或的项） 4、一阶线性非齐次方程常数变易法，或])([)()(?+??=-C dx e x b e y dx x a dx x a 5、伯努力方程令n y z -=1，则dx dy y n dx dz n --=)1(，可将伯努力方程化成一阶线性非齐次或一阶线性齐次 6、全微分方程若x N y M ??=??，则C y x u =),(，（留意书上公式）若 x N y M ??≠??，则找积分因子，（留意书上公式） 7、可降阶的二阶微分方程 ),(22dx dy x f dx y d =，令dx dy dx y d p dx dy ==22,则 ),(22dx dy y f dx y d =，令dy dp p dx y d p dx dy ==22,则 8、正交轨线族（二）毕卡序列 ?+=x x dx y x f y y 0),(001，?+=x x dx y x f y y 0),(102，?+=x x dx y x f y y 0),(203，其余类推（三）常系数方程 1、常系数齐次0)(=y D L 方法：特征方程单的实根21,λλ，x x e C e C y 2121λλ+= 单的复根i βαλ±=2,1，)sin cos (21x C x C e y x ββα+= 重的实根λλλ==21，x e x C C y λ)(21+= 重的复根i βαλ±=2,1，i βαλ±=4,3，]sin )(cos )[(4321x x C C x x C C e y x ββα+++=

微分方程数值解习题(李立康)

常微分方程习题《李立康》习题 1.用Euler 方法求初值问题 ? ? ?=-='0)0(21u tu u 在1=t 时的近似解（取4 1= h ）。 2.初值问题 1 3 00 u u u()??'=? ?=? 有解32 23/u(t )t ?? = ? ?? 。但若用Euler 方法求解，对一切N T ,和H T h = ,都只能得到N t u t ,...,2,1,0==，试解释此现象产生的原因。 3.用Euler 方法计算 ?? ?=='1 )0(u u u 在1=t 处的值，取16 1 和41= h ，将计算结果与精确值e =）1（u 相比较。 4.设),(u t f 满足定理2.1的条件，对改进Euler 法（2.10）式证明：（1）其局部截断误差为)()(12 43 h O t u h -'''- ；（2）当1

?? ?=='1 )0(u u u 计算公式 m m h h u ??? ? ??-+=22 取4 1 = h 计算)1(u 的近似值，并与习题3的结果比较。 6.就初值问题 ?? ?=+='0 )0(u b at u 分别导出用Euler 方法和改进Euler 法求近似解的表达式，并与真解 bt t a u += 22 相比较。 7.证明改进Euler 法的绝对稳定区域是整个左半平面0)Re(

2018年电大第三版常微分方程答案知识点复习考点归纳总结参考

习题1.2 1．dx dy =2xy,并满足初始条件：x=0,y=1的特解。解：y dy =2xdx 两边积分有：ln|y|=x 2+c y=e 2x +e c =cex 2另外y=0也是原方程的解，c=0时，y=0 原方程的通解为y= cex 2,x=0 y=1时 c=1 特解为y= e 2x . 2. y 2dx+(x+1)dy=0 并求满足初始条件：x=0,y=1的特解。解：y 2dx=-(x+1)dy 2y dy dy=-11+x dx 两边积分: -y 1=-ln|x+1|+ln|c| y=|)1(|ln 1+x c 另外y=0,x=-1也是原方程的解 x=0,y=1时 c=e 特解：y=|)1(|ln 1 +x c 3．dx dy =y x xy y 321++ 解：原方程为：dx dy =y y 21+31x x + y y 21+dy=31x x +dx 两边积分：x(1+x 2)(1+y 2)=cx 2 4. (1+x)ydx+(1-y)xdy=0 解：原方程为： y y -1dy=-x x 1+dx 两边积分：ln|xy|+x-y=c 另外 x=0,y=0也是原方程的解。 5．（y+x ）dy+(x-y)dx=0 解：原方程为：

dx dy =- y x y x +- 令x y =u 则dx dy =u+x dx du 代入有： -1 12++u u du=x 1 dx ln(u 2+1)x 2=c-2arctgu 即 ln(y 2+x 2)=c-2arctg 2x y . 6. x dx dy -y+22y x -=0 解：原方程为： dx dy =x y +x x ||-2)(1x y - 则令x y =u dx dy =u+ x dx du 211u - du=sgnx x 1 dx arcsin x y =sgnx ln|x|+c 7. tgydx-ctgxdy=0 解:原方程为：tgy dy =ctgx dx 两边积分：ln|siny|=-ln|cosx|-ln|c| siny=x c cos 1=x c cos 另外y=0也是原方程的解，而c=0时，y=0. 所以原方程的通解为sinycosx=c. 8 dx dy +y e x y 32+=0 解：原方程为：dx dy =y e y 2e x 3 2 e x 3-3e 2y -=c. 9.x(lnx-lny)dy-ydx=0 解：原方程为： dx dy =x y ln x y 令 x y =u ,则dx dy =u+ x dx du

MATLAB求解常微分方程数值解

利用MATLAB求解常微分方程数值解

目录 1. 内容简介 (1) 2. Euler Method(欧拉法)求解 (1) 2.1. 显式Euler法和隐式Euler法 (2) 2.2. 梯形公式和改进Euler法 (3) 2.3. Euler法实用性 (4) 3. Runge-Kutta Method(龙格库塔法)求解 (5) 3.1. Runge-Kutta基本原理 (5) 3.2. MATLAB中使用Runge-Kutta法的函数 (7) 4. 使用MATLAB求解常微分方程 (7) 4.1. 使用ode45函数求解非刚性常微分方程 (8) 4.2. 刚性常微分方程 (9) 5. 总结 (9) 参考文献 (11) 附录 (12) 1. 显式Euler法数值求解 (12) 2. 改进Euler法数值求解 (12) 3. 四阶四级Runge-Kutta法数值求解 (13) 4.使用ode45求解 (14)

1.内容简介把《高等工程数学》看了一遍，增加对数学内容的了解，对其中数值解法比较感兴趣，这大概是因为在其它各方面的学习和研究中经常会遇到数值解法的问题。理解模型然后列出微分方程，却对着方程无从下手，无法得出精确结果实在是让人难受的一件事情。实际问题中更多遇到的是利用数值法求解偏微分方程问题，但考虑到先从常微分方程下手更为简单有效率，所以本文只研究常微分方程的数值解法。把一个工程实际问题弄出精确结果远比弄清楚各种细枝末节更有意思，因此文章中不追求非常严格地证明，而是偏向如何利用工具实际求解出常微分方程的数值解，力求将课程上所学的知识真正地运用到实际方程的求解中去，在以后遇到微分方程的时候能够熟练运用MATLAB得到能够在工程上运用的结果。文中求解过程中用到MATLAB进行数值求解，主要目的是弄清楚各个函数本质上是如何对常微分方程进行求解的，对各种方法进行MATLAB编程求解，并将求得的数值解与精确解对比，其中源程序在附录中。最后考察MATLAB中各个函数的适用范围，当遇到实际工程问题时能够正确地得到问题的数值解。 2.Euler Method(欧拉法)求解 Euler法求解常微分方程主要包括3种形式，即显式Euler法、隐式Euler法、梯形公式法，本节内容分别介绍这3种方法的具体内容，并在最后对3种方法精度进行对比，讨论Euler法的实用性。本节考虑实际初值问题使用解析法，对方程两边同乘以得到下式

常微分方程数值解

第四章常微分方程数值解 [课时安排]6学时 [教学课型]理论课 [教学目的和要求] 了解常微分方程初值问题数值解法的一些基本概念，如单步法和多步法，显式和隐式，方法的阶数，整体截断误差和局部截断误差的区别和关系等；掌握一阶常微分方程初值问题的一些常用的数值计算方法，例如欧拉（Euler）方法、改进的欧拉方法、龙贝－库塔（Runge-Kutta）方法、阿达姆斯（Adams）方法等，要注意各方法的特点及有关的理论分析；掌握构造常微分方程数值解的数值积分的构造方法和泰勒展开的构造方法的基本思想，并能具体应用它们导出一些常用的数值计算公式及评估截断误差；熟练掌握龙格－库塔（Ｒ－Ｋ）方法的基本思想，公式的推导，Ｒ－Ｋ公式中系数的确定，特别是能应用“标准四阶Ｒ－Ｋ公式”解题；掌握数值方法的收敛性和稳定性的概念，并能确定给定方法的绝对稳定性区域。[教学重点与难点] 重点：欧拉方法，改进的欧拉方法，龙贝－库塔方法。难点：R—K方法，预估-校正公式。 [教学内容与过程] 4.1 引言本章讨论常微分方程初值问题 (4.1.1) 的数值解法，这也是科学与工程计算经常遇到的问题，由于只有很特殊的方程能用解析方法求解，而用计算机求解常微分方程的初值问题都要采用数值方法.通常我们假定(4.1.1)中 f(x,y)对y满足Lipschitz条件，即存在常数L＞0，使对，有 (4.1.2) 则初值问题(4.1.1)的解存在唯一. 假定(4.1.1)的精确解为，求它的数值解就是要在区间上的一组离散点上求的近似.通常取，h称为步长，求(4.1.1)的数值解是按节点的顺序逐步推进求得.首先，要对方程做离散逼近，求出数值解的公式，再研究公式的局部截

常微分方程解题方法总结.doc

常微分方程解题方法总结来源：文都教育复习过半，课本上的知识点相信大部分考生已经学习过一遍 . 接下来，如何将零散的知识点有机地结合起来，而不容易遗忘是大多数考生面临的问题 . 为了加强记忆，使知识自成体系，建议将知识点进行分类系统总结 . 著名数学家华罗庚的读书方法值得借鉴，他强调读书要“由薄到厚、由厚到薄”，对同学们的复习尤为重要 . 以常微分方程为例，本部分内容涉及可分离变量、一阶齐次、一阶非齐次、全微分方程、高阶线性微分方程等内容，在看完这部分内容会发现要掌握的解题方法太多，遇到具体的题目不知该如何下手，这种情况往往是因为没有很好地总结和归纳解题方法 . 下面以表格的形式将常微分方程中的解题方法加以总结，一目了然，便于记忆和查询 . 常微分方程通解公式或解法 ( 名称、形式 ) 当 g( y) 0 时，得到 dy f (x)dx ， g( y) 可分离变量的方程 dy f ( x) g( y) 两边积分即可得到结果； dx 当 g( 0 ) 0 时，则 y( x) 0 也是方程的解 . 解法：令 u y xdu udx ，代入，则 dy 齐次微分方程 dy g( y ) x dx x u g (u) 化为可分离变量方程得到 x du dx 一阶线性微分方程 P ( x)dx P ( x) dx dy Q(x) y ( e Q( x)dx C )e P( x) y dx

伯努利方程解法：令 u y1 n，有 du (1 n) y n dy ， dy P( x) y Q( x) y n（n≠0,1）代入得到du (1 n) P(x)u (1 n)Q(x) dx dx 求解特征方程：2 pq 三种情况：二阶常系数齐次线性微分方程 y p x y q x y0 二阶常系数非齐次线性微分方程 y p x y q x y f ( x) (1)两个不等实根：1, 2 通解： y c1 e 1x c2 e 2x (2) 两个相等实根：1 2 通解： y c1 c2 x e x (3) 一对共轭复根：i , 通解： y e x c1 cos x c2 sin x 通解为 y p x y q x y 0 的通解与 y p x y q x y f ( x) 的特解之和. 常见的 f (x) 有两种情况： x （ 1）f ( x)e P m ( x) 若不是特征方程的根，令特解 y Q m ( x)e x；若是特征方程的单根，令特解 y xQ m ( x)e x；若是特征方程的重根，令特解 y*x2Q m (x)e x；（2）f (x) e x[ P m ( x) cos x p n ( x)sin x]

常微分方程的数值解

实验4 常微分方程的数值解【实验目的】 1．掌握用MATLAB软件求微分方程初值问题数值解的方法； 2．通过实例用微分方程模型解决简化的实际问题； 3．了解欧拉方法和龙格-库塔方法的基本思想和计算公式，及稳定性等概念。【实验内容】题3 小型火箭初始重量为1400kg，其中包括1080kg燃料。火箭竖直向上发射时燃料燃烧率为18kg/s，由此产生32000N的推力，火箭引擎在燃料用尽时关闭。设火箭上升时空气阻力正比于速度的平方，比例系数为m，求引擎关闭瞬间火箭的高度、速度、加速度，及火箭到达最高点的时的高度和加速度，并画出高度、速度、加速度随时间变化的图形。模型及其求解火箭在上升的过程可分为两个阶段，在全过程中假设重力加速度始终保持不变，g=s2。在第一个过程中，火箭通过燃烧燃料产生向上的推力，同时它还受到自身重力（包括自重和该时刻剩余燃料的重量）以及与速度平方成正比的空气阻力的作用，根据牛顿第二定律，三个力的合力产生加速度，方向竖直向上。因此有如下二式： a=dv/dt=/m=/(1400-18t) dh/dt=v 又知初始时刻t=0，v=0，h=0。记x(1)=h，x(2)=v，根据MATLAB 可以求出0到60秒内火箭的速度、高度、加速度随时间的变化情况。程序如下： function [ dx ] = rocket( t,x ) a=[*x(2)^2)/(1400-18*t)]; dx=[x(2);a]; end ts=0:1:60;

x0=[0,0]; [t,x]=ode45(@rocket,ts,x0); h=x(:,1); v=x(:,2); a=[*(v.^2))./(1400-18*t)]; [t,h,v,a]; 数据如下： t h v a 000

常微分方程数值解法

第七章常微分方程数值解法常微分方程中只有一些典型方程能求出初等解(用初等函数表示的解)，大部分的方程是求不出初等解的。另外，有些初值问题虽然有初等解，但由于形式太复杂不便于应用。因此，有必要探讨常微分方程初值问题的数值解法。本章主要介绍一阶常微分方程初值问题的欧拉法、龙格-库塔法、阿达姆斯方法，在此基础上推出一阶微分方程组与高阶方程初值问题的数值解法；此外，还将简要介绍求解二阶常微分方程值问题的差分方法、试射法。第一节欧拉法求解常微分方程初值问题 ?????==0 0)() ,(y x y y x f dx dy （1）的数值解，就是寻求准确解)(x y 在一系列离散节点 <<<<

微分方程数值解(学生复习题)

一．填空 1. Euler 法的一般递推公式为，整体误差为，局部截断误差为： .，改进Euler 的一般递推公式整体误差为，局部截断误差为：。 2. 线性多步法绝对稳定的充要条件是。 3.当，则单步法1(,,)0,1,2,,n n n n T u u h t u h n h ?+=+= ，稳定。 4. 一个相容，稳定的多步法若绝对稳定，则绝对稳定域在。 5. 若，则多步法是相容的。 6．所有内点，界点的差分方程组成一个封闭的线性代数方程组，其系数矩阵是。 7.刚性方程是： 8．Runge-Kutta 法的特征值为，相容的充要条件为： 8.二阶常微分方程边值问题：22,(), ()d u Lu qu f a x b dx u a u b αβ?=-+=<

4、一级Runge-Kutta 法的绝对稳定域（-2，0） 5、若差分方程满足相容条件，且按右端稳定，则差分解收敛至波动方程的解。 6、Euler 法非A 稳定。 7．对任意网比0r >，六点对称格式的解有收敛阶22()O h τ+ 8. 对任意网比12 r ≤，向前差分格式的解有收敛阶2()O h τ+。 9、相容，稳定的多步法一定绝对稳定。三．选择 1.抛物型方程的加权隐式差分格式的稳定性为（） A 绝对稳定 B 无条件稳定 C 条件稳定 D 非条件稳定 2.von Neumann 条件是差分格式稳定的（） A 充分条件 B 必要条件 C 充要条件 D 既非充分也非必要条件 3．实系数二次方程20b c λλ--=的根按模小于或者等于1的充要条件是（） A 12b c ≤-≤ B 1+2b c ≤≤ C 12c b ≤-≤ D 12c b ≤+≤ 4.若线性多步法A 稳定，则有（），其中1,2,,i i k λ= （）为()()0h ρλσλ-=的根。 A Re 01,1,2,,i h i k λ= B 1Re 0i h λ≥?≥ C Re 01,1,2,,i h i k λ≤?≤= D 1Re 0i h λ