科学与工程计算科学计算.ppt_第1页
科学与工程计算科学计算.ppt_第2页
科学与工程计算科学计算.ppt_第3页
科学与工程计算科学计算.ppt_第4页
科学与工程计算科学计算.ppt_第5页
已阅读5页,还剩24页未读 继续免费阅读

下载本文档

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

文档简介

,科学与工程计算方法,主要内容常微分数值计算方法简介偏微分方程数值计算方法(差分法、有限元法、有限体积法)统计方法简介选讲内容计算流体力学方法简介计算电磁学方法简介,主要参考书科学与工程计算方法北京理工大学出版社数值计算方法北京理工大学出版社等偏微分方程数值解法清华大学出版社、华中理工大学出版社等微分方程的数值方法(英文)springer出版社实用多元统计方法与SAS系统北京大学出版社计算电磁学的数值方法清华大学出版社,常微分方程的初值问题第一节几种经典方法,求解常微分方程初值问题,的数值解,就是寻求准确解,在一系列离散节点,上的近似值,只有初值问题的解存在且唯一时,使用数值解法才有意义,这一前提条件由下面定理保证。,1:欧拉法(欧拉折线法),向前、向后欧拉(Euler)公式。以上方法称为欧拉法或欧拉折线法。,欧拉公式的几何意义,欧拉法的其它解释:,2:改进的欧拉法(梯形法),欧拉法形式简单,计算方便,但曲线的曲率较大时,欧拉法的效果差。于是,对欧拉法进行改进,用梯形公式计算,得到改进的欧拉公式,称为改进的欧拉法,也称为梯形法。,先用欧拉公式,3:预估一校正法,预估校正公式也常写成下列形式:,公式的局部截断误差,欧拉法的局部截断误差,近似解与准确解比较,欧拉法的结果大致只有两位有效数字,而预估校正法的结果则有3位有效数字。,四显式龙格库塔法,一泰勒级数法,设有初值问题,由泰勒展开式,即公式(2)为k阶方法。,所以这种方法实际上不能用来解初值问题。,二龙格库塔方法(R-K方法),R-K方法不是通过求导数的方法构造近似公式,而是通过计算不同点上的函数值,并对这些函数值作线性组合,构造近似公式,再把近似公式与解的泰勒展开式进行比较,使前面的若干项相同,从而使近似公式达到一定的阶数。先分析欧拉法与预估校正法。,与二阶龙格库塔公式的讨论方法类似,只需8个参数满足,满足条件(6)的一族公式(5)统称为三阶龙格库塔公式。,一个比较简单的三阶龙格龙塔公式是,(6),(6)式中13个待定常数需满足下列11个方程的方程组,和吉尔公式,例用标准四阶龙格库塔方法求解初值问题,解计算过程和结果,五、常用的线性多步公式,四阶Adams显式公式,四阶Adams隐式公式,1.Adams公式,(二)Milne公式,(三)Hamming公式,第二节:常微分方程组与高阶方程的数值解法经典方法的推广,前面几节介绍了一阶微分方程初值问题的各种数值解法,这些解法同样适用于一阶微分方程组与高阶方程。首先从简单的开始,下面仅就两个未知函数的方程组和二阶方程为例说明各种方法的计算公式。,3四阶标准龙格库塔公式,4四阶亚当斯显式公式,高阶方程,我们以二阶方程为例来说明高阶方程的数

温馨提示

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

评论

0/150

提交评论