已阅读5页,还剩72页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
第8章常微分方程初值问题数值解法,引言常用数值解法的导出与使用数值解中误差的积累外推极限法微分方程组与高阶方程解法,引言,在工程和科学计算中,所建立的各种常微分方程的初值或边值问题,除很少几类的特殊方程能给出解析解,绝大多数的方程是很难甚至不可能给出解析解的,其主要原因在于积分工具的局限性。因此,人们转向用数值方法去解常微分方程,并获得相当大的成功,讨论和研究常微分方程的数值解法是有重要意义的。,8.1常用数值解法的导出与使用,8.1.1数值微分法局部截断误差,8.1.1数值微分法局部截断误差,精确解:y=(1+2x)-0.45,8.1.2数值积分法隐式公式的使用,例8.2用梯形公式求解例8.1的问题,8.1.2数值积分法隐式公式的使用,预测值,校正值,初值问题的Euler方法,初值问题的Euler方法,误差概述,误差概述,误差概述,6.1.3数值稳定性分析,数值稳定性分析,定义6.1.3若某数值算法的绝对稳定性区域包含h平面上的左半平面Re(h)0,则称该方法是A稳定的。隐式Euler法是A稳定的。,6.2Runge-Kutta方法,Runge-Kutta方法,Runge-Kutta方法,Runge-Kutta方法,6.2.2四阶Runge-Kutta方法,四阶Runge-Kutta方法,6.2.3R-K法的稳定性,R-K法的稳定性,R-K法的稳定性,6.2.5隐式R-K法,隐式R-K法,隐式R-K法,隐式R-K法,隐式R-K法,6.3线形多步法,单步法主要依据yn的信息去计算yn+1。线性多步法是想依据yn,yn-1,yn-r(r1)的信息去计算yn+1。考虑到线性组合较为方便,因此,线性多步法一般形式可设为,6.3.1基于数值积分的方法,基于数值积分的方法,基于数值积分的方法,基于数值积分的方法,基于数值积分的方法,基于数值积分的方法,Adams预估校正法预估校正并取,6.3.2基于Taylar展开式的方法,基于Taylar展开式的方法,基于Taylar展开式的方法,6.4一阶常微分方程组数值解法,在许多实际问题中,常常出现高阶微分方程和高阶微分方程组,通过引入新的变量,总可化为一阶微分方程组。由此可知,讨论一阶常微分方程组的数值解法是很有意义的。,6.4.1解一阶常微分方程组的R-K方法,一阶常微分方程组的R-K方法,一阶常微分方程组的R-K方法,一阶常微分方程组的R-K方法,一阶常微分方程组的R-K方法,一阶常微分方程组的R-K算法,一阶常微分方程组的R-K方法,6.4.2刚性方程组,刚性方程组,刚性方程组,6.5常微分方程边值问题的数值解法,设二阶线性常微分方程为常见边界条件有三类:,6.5.1差分方程的建立,差分方程的建立,
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 初中数学竞赛专题5:分式
- 交叉设计在生物等效性试验中的EMA要求解析
- 骨科大手术VTE预防指南专家讲座
- 渔业行业的困境与改进对策设想
- 2007年上半年市直事业单位招聘工作人员职位表【模板】
- 临床路径标准化在单病种成本核算中的应用
- 信息资源检索与利用总复习题答案(部分)①
- 临床路径虚拟仿真系统的迭代优化策略
- 复发性流产病因分级筛查临床实践中国专家共识2026
- 风湿病的概述课件
- 物业清运垃圾车管理制度
- 手足口病教学课件
- 香料调味培训课件模板
- 国家职业标准 -碳排放管理员
- 肿瘤患者并发症管理
- 儿童口腔诊疗行为管理学
- 北京市2025学年高二(上)第一次普通高中学业水平合格性考试物理试题(原卷版)
- 十五五原材料工业发展规划
- 建筑工程质量员课件
- 驾考宝典三力测试考试试题及答案
- 2025年中国汽车零部件市场研究报告-2025-04-零部件
评论
0/150
提交评论