




已阅读5页,还剩21页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
解一阶常微分方程欧拉法局部截断误差与p阶精度Range-Kutta公式一阶常微分方程组和二阶方程线性多步法简介,数值分析23,17:20,常微分方程(ODE)的定解问题主要有初值问题(IVP)和边值问题(BVP)两大类。,急性传染病数学模型(SIR模型),常微分方程的初值问题是描述系统发展演变的重要工具。,洛伦兹吸引子lorenzgui,关于爱情的动力系统Loveraffairesanddifferentialequations,17:20,一阶常微分方程初值问题:,离散化(discretization)取定离散点x0x1x2xN,步长h,其中y=y(x)是未知函数,右端函数f(x,y)是已知函数,初值y0是已知数据。,求未知函数y(x)在离散点处的近似值y1,y2,y3,yN,考察点xn列出微分方程,Euler法与修正的Euler法,其中yn是y(xn)的近似。,足够多的点来近似连续对象,17:20,求解常微分方程初值问题的Euler方法,取定步长:h,记xn=x0+nh,(n=1,2,N)称计算格式:yn+1=yn+hf(xn,yn)为显式Euler公式。对应的求初值问题数值解的方法称为Euler方法。,例1.用Euler法求初值问题,的数值解。,解:记f(x,y)=yxy2,xn=nh(n=0,1,2,N)由Euler公式得:yn+1=yn+h(ynxnyn2)(n=0,1,N),17:20,取步长h=2/10,2/20,2/30,2/40,用Euler法求解的数值实验结果如下.,解析解:,o数值解-准确解,17:20,y=f(x,y),左矩形公式,用数值积分方法离散化常微分方程,显式Euler法,求解常微分方程初值问题的数值方法可以由很多,各种方法都可以看成使用不同的数值积分方法计算,17:20,梯形公式:,右矩阵公式:,隐式Euler法,17:20,预-校方法算法如下k1=f(xn,yn),k2=f(xn+1,yn+hk1),预报-校正方法(修正的Euler法):,预报,校正,17:20,预-校方法(h=0.2)误差最大值0.0123,欧拉方法(h=0.2)误差最大值0.1059,17:20,基本概念,单步法:即为了求得后一点xn+1上数值解yn+1,只要知道一点xn的数值解yn就可以了。,多步法:即计算yn+1时,要用到前面多点的数值解yn,yn-1,yn-r等信息。,显式格式:数值解yn+1可以用xn和yn解析表出。这类格式称为显式格式。反之称为隐式格式。对隐式格式而言,计算数值解yn+1形式上需要解方程(组)。,显式方法易于计算,但通常数值稳定性较差。而隐式方法通常稳定性较好,但计算上显然很不方便。将显式方法和隐式方法联合使用是重要的思想方法。,17:20,Tn+1=y(xn+1)-yn+1称为局部截断误差。设yn=y(xn),即,由Taylor公式,Euler公式:yn+1=yn+hf(xn,yn)的局部截断误差,Euler公式的局部截断误差记为O(h2),则称Euler公式具有1阶精度。,局部截断误差,仅考虑xn到xn+1的局部情况,并假定xn之前的计算没有误差,17:20,若局部截断误差为O(hp+1)则称显式单步法具有p阶精度。,例2.证明梯形方法具有2阶精度,从x0开始计算,如果考虑每一步产生的误差直到xn,则误差y(xn)-yn称为数值方法在节点xn处的整体误差。,粗略的讲,整体误差把许多步的局部截断误差加在一起,而步数正比于1/h(即步长倒数),因此整体误差O(hp)。,17:20,(Runge-Kutta)龙格-库塔法,改进的Euler方法这样理解:它用xn和xn+1两个点的斜率值K1和K2取算术平均作为平均斜率,而xn+1处的斜率值则利用已知信息xn来预报。,RK方法是一大类的方法,其基本思想是采用如下形式:,是否可以推广改进的Euler方法?,17:20,取点个数s称为龙格-库塔法的级数。一般的RK公式如下给出:,RK方法是单步法(为了求得后一点xn+1上数值解yn+1,只要知道一点xn的数值解yn)。,我记得我的朋友JohnvonNeumann曾经说过,用四个参数我可以拟合出一头大象,而用五个参数我可以让它的鼻子摆动-费米,下面以2级方法为例子具体介绍龙格-库塔法。,构造的基本思想是选择适当的系数使得方法的局部截断误差阶数尽可能高。,17:20,17:20,三阶Range-Kutta公式一般形式,yn+1=yn+hk1+4k2+k3/6k1=f(xn,yn),k2=f(xn+0.5h,yn+0.5hk1)k3=f(xn+h,ynhk1+2hk2),四阶Range-Kutta公式一般形式,yn+1=yn+hk1+2k2+2k3+k4/6k1=f(xn,yn),k2=f(xn+0.5h,yn+0.5hk1)k3=f(xn+0.5h,yn+0.5hk2),k4=f(xn+h,yn+hk3),17:20,例4,数值实验:几种不同求数值解公式的误差比较,17:20,常微分方程初值问题的MATLAB,数值求解命令:x,y=ode23(f,a,b,y0),f=inline(y-x.*y.2);x,y=ode23(f,0,2,1),符号求解命令:dsolve(eqn1,.),symsxydsolve(Dy=y-x*y2,y(0)=1,x),解析解:,参考:TheMatlabODEsuite,SIAMJScientificComputing,17:20,一阶常微分方程组,欧拉公式:,(n=0,1,2,N),一阶常微分方程组向量形式,记,17:20,修改的欧拉公式:,(n=0,1,N),经典龙格-库塔公式:,17:20,洛伦兹模型,取=8/3,=10,=28。,x(0)=0,y(0)=0,z(0)=0.01。,t0,80,,微分方程右端函数:,17:20,线性多步法,(n=0,1,),
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 虎门中学模拟考试题及答案
- 福建幼教考试题目及答案
- 行政自学考试试题及答案
- 2025年执业药师历届真题及答案
- 2025年抖音小店客服试题及答案
- 化学实验仪器使用方法判断试题
- 化学创新精神素养突破常规评价试题
- 古代大学考试题目及答案
- 公文考试题库及答案
- 工伤保险考试题及答案
- 2025年肺功能证考试题及答案
- 2026中国海洋石油集团有限公司秋季校园招聘备考考试题库附答案解析
- 2025年及未来5年中国羊奶粉行业市场调研分析及投资战略咨询报告
- 学校物业委托管理服务合同7篇
- 2025-2026学年人教版二年级上册数学第三单元测试卷(含答案)(三套)
- 《守望成长法治护航》法制教育主题班会
- 高中日语学习宣讲+课件
- 2023年中国通用技术(集团)控股有限责任公司招聘笔试题库含答案解析
- 生物技术与农业
- GB/T 5668-2017旋耕机
- 动物资源保护与利用
评论
0/150
提交评论