常微分方程初值问题的数值解法_第1页
常微分方程初值问题的数值解法_第2页
常微分方程初值问题的数值解法_第3页
常微分方程初值问题的数值解法_第4页
常微分方程初值问题的数值解法_第5页
已阅读5页,还剩5页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、常微分方程初值问题的数值解法,一阶常微分方程初值问题指的数值解法,1 欧拉法,1)欧拉法 根据导数的定义离散化之后,有,代入常微分方程,并构造迭代系列得,2)隐式欧拉法,代入常微分方程,并构造迭代系列得,隐式欧拉法的求解: 利用迭代的思路进行,中点欧拉法、改进的欧拉法,3)中点欧拉法,代入常微分方程,并构造迭代系列得,4)改进的欧拉法,欧拉法例题,在区间0,1上取步长h=0.2求解如下的常微分方程,解:按照欧拉法的计算公式,分别得到如下的函数值,建立高精度的单步递推格式,单步递推法的基本思想是从 ( xi , yi ) 点出发,以某一斜率沿直线达到 ( xi+1 , yi+1 ) 点。欧拉法及

2、其各种变形所能达到的最高精度为2阶,斜率 一定取K1 K2 的平均值吗,步长一定是一个h 吗,2 龙格 - 库塔法 /* Runge-Kutta Method *,2 Runge-Kutta Method,首先希望能确定系数 1、2、p,使得到的算法格式有2阶精度,即在 的前提假设下,使得,Step 1: 将 K2 在 ( xi , yi ) 点作 Taylor 展开,Step 2: 将 K2 代入第1式,得到,2 Runge-Kutta Method,Step 3: 将 yi+1 与 y( xi+1 ) 在 xi 点的泰勒展开作比较,要求 ,则必须有,这里有 个未知数, 个方程,3,2,存在无穷多个解。所有满足上式的格式统称为2阶龙格 - 库塔格式,注意到, 就是改进的欧拉法,Q: 为获得更高的精度,应该如何进一步推广,其中i ( i = 1, , m ),i ( i = 2, , m ) 和 ij ( i = 2, , m; j = 1, , i1 ) 均为待定系数,确定这些系数的步骤与前面相似,2 Runge-Kutta Method,最常用为四级4阶经典龙格-库塔法 /* Classical Runge-Kutta Method */,R-K例题,用R-K法求解方程的解

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论