六常微分方程及方程组的解法.pptx
《六常微分方程及方程组的解法.pptx》由会员分享,可在线阅读,更多相关《六常微分方程及方程组的解法.pptx(116页珍藏版)》请在咨信网上搜索。
1、浙江大学研究生学位课程16.1 常微分方程及其求解概述常微分方程及其求解概述6.2 初值问题解法初值问题解法 1.Euler方法方法 2.线性多步法线性多步法 3.Runge-Kutta法法 4.方程组及刚性问题的方程组及刚性问题的Gear方法方法6.3 边值问题解法边值问题解法 1.Shooting(试射法试射法)2.差分法差分法 浙江大学研究生学位课程26.1 常微分方程及求解概述常微分方程及求解概述(Ordinary Differential Equations,ODE)6.1.1 基本概念基本概念 描述自由落体的描述自由落体的ODE:浙江大学研究生学位课程3 只有一个自变量的微分方程为
2、只有一个自变量的微分方程为 ODE,否则称为偏微分方程否则称为偏微分方程 PDE。方程中未知函数导数的最高阶数称为方程中未知函数导数的最高阶数称为 方程的阶。方程的阶。(6-4)是二阶的)是二阶的 方程中关于未知函数及其各阶导数方程中关于未知函数及其各阶导数 均是一次的,则称为线性微分方程。均是一次的,则称为线性微分方程。和(和(6-1)都是线性二阶)都是线性二阶ODE。(6-2),(6-3)是是(6-1)的的初始条件初始条件。亦称定解。亦称定解 条件。条件。(6-1)(6-2)叫做叫做初值问题初值问题。6.1.1 浙江大学研究生学位课程46.1.1 (6-1),(6-3)叫做边值问题。在没有
3、给定解条件时。方程一般在没有给定解条件时。方程一般 有一族解曲线有一族解曲线 y(x,c)。如:。如:对任意的对任意的n阶阶ODE,如果能写成:,如果能写成:则称该方程为则称该方程为显式显式的。方程(的。方程(6-4)是)是 显式的。而下面方程是显式的。而下面方程是隐式隐式的。的。浙江大学研究生学位课程5 对于高阶显式方程。通过定义对于高阶显式方程。通过定义n-1个个 新变量,可以写成新变量,可以写成n维一阶方程组。维一阶方程组。即令:即令:6.1.1 浙江大学研究生学位课程6 在讨论初值问题时,我们从一阶方在讨论初值问题时,我们从一阶方 程开始:程开始:然后毫不费力地套用来解方程组。然后毫不
4、费力地套用来解方程组。当当 f(x,y)与与y无关时,无关时,f(x,y)=g(x)6.1.1 浙江大学研究生学位课程7 6.1.2 数值解及其重要性数值解及其重要性浙江大学研究生学位课程8 6.1.3 ODE数值解的基本思想和方法特点数值解的基本思想和方法特点 基本思想有两点基本思想有两点 1.离散化离散化 用Taylor级数,数值积分和差商 逼近导数等手段,把ODE转化为离散 的代数方程(称差分方程)。2.递推化递推化 在具有唯一解的条件下,通过 步进法逐步计算出解在一系列离散 点上的值。从而得到原ODE的数值 近似解。浙江大学研究生学位课程9 6.2 初值问题解法初值问题解法 我们讨论一
5、阶ODE,而高阶可 能化为一阶ODEs。一阶初值问题 可以一般地写成:6.2.1 欧拉(欧拉(Euler)方法方法 Euler方法是求解(6-8)最简单简单方法,但精度差精度差,故不实用不实用。然而对理论分 析很有用。浙江大学研究生学位课程10 6.2.1.1 方法原理及推导方法原理及推导 设初值问题设初值问题(6-8)满足:满足:6.2.1浙江大学研究生学位课程116.2.1 图图 6.1 常微分方程初值问题的数值解常微分方程初值问题的数值解浙江大学研究生学位课程126.2.1浙江大学研究生学位课程13欧拉方法的几何意义:欧拉方法的几何意义:h步长6.2.1 图图 6.2 Euler方法的几
6、何意义方法的几何意义浙江大学研究生学位课程146.2.1 浙江大学研究生学位课程156.2.1 浙江大学研究生学位课程166.2.1 多多步步法法单步法单步法自动起步自动起步显式显式隐式隐式半隐式半隐式图图 6.3 ODE求解方法的类型求解方法的类型浙江大学研究生学位课程17表表 6.1 自由落体运动方程的自由落体运动方程的Euler公式求解公式求解6.2.1浙江大学研究生学位课程18图图 6.4 运动轨迹运动轨迹6.2.1浙江大学研究生学位课程19图图 6.5 Euler公式的误差公式的误差 6.2.1.2 Euler方法的误差估计方法的误差估计 一般其它方法的误差估计也类似。这里误差是指截
7、断误差截断误差(算法理论误 差)而不是舍入误差舍入误差。后者由计算机字 长等决定,属于稳定性问题。i)几何分析几何分析6.2.1 浙江大学研究生学位课程206.2.1 局部截断局部截断误差,误差,浙江大学研究生学位课程216.2.1 整体截断误差:整体截断误差:设ym是Euler公式(6-9)精确解,而 y(x)是初值问题(6-8)的解。则 整体截断误差定义为 它是局部截断误差的积累。定理:定理:若f(x,y)关于y满足Lipschitz条件。则有估计式:浙江大学研究生学位课程226.2.1 浙江大学研究生学位课程236.2.1 浙江大学研究生学位课程246.2.1 浙江大学研究生学位课程25
8、6.2.1 注意注意:稳定性稳定性:浙江大学研究生学位课程266.2.1 浙江大学研究生学位课程276.2.1 浙江大学研究生学位课程286.2.1浙江大学研究生学位课程296.2.1 浙江大学研究生学位课程306.2.1浙江大学研究生学位课程316.2.1表表 6.2 予估校正求解结果对比予估校正求解结果对比浙江大学研究生学位课程326.2.1表表 6.3 Euler法与外推结果的比较法与外推结果的比较浙江大学研究生学位课程33 6.2.2 线性多步法线性多步法 浙江大学研究生学位课程346.2.2()()浙江大学研究生学位课程356.2.2 Adams 外插法(k=2)3阶3步 显式表表
9、6.4 外插系数外插系数bki值值图图 6.6 3阶阶3步外插法步外插法浙江大学研究生学位课程366.2.2 浙江大学研究生学位课程376.2.2 Adams 外插法(k=2)4阶3步 图图 6.7 4步步3阶阶Adams内插公式内插公式浙江大学研究生学位课程386.2.2 浙江大学研究生学位课程396.2.2 浙江大学研究生学位课程406.2.2 图图 6.8 一般化插值形式一般化插值形式浙江大学研究生学位课程416.2.2 浙江大学研究生学位课程426.2.2 浙江大学研究生学位课程436.2.2 浙江大学研究生学位课程446.2.2 浙江大学研究生学位课程456.2.2 浙江大学研究生学
10、位课程466.2.2 图图 6.9浙江大学研究生学位课程476.2.2 线性多步法的绝对稳定性:线性多步法的绝对稳定性:浙江大学研究生学位课程486.2.2 定义定义:绝对稳定绝对稳定。绝对稳定区域。绝对稳定区域。浙江大学研究生学位课程496.2.2 Milne浙江大学研究生学位课程506.2.2 表表 6.6 计算结果计算结果浙江大学研究生学位课程516.2.2 6.2.2.4 预估预估-校正方法校正方法 (Predictor-Corrector Method)浙江大学研究生学位课程526.2.2 浙江大学研究生学位课程536.2.2 注意:注意:一步校正的计算量 预估计算量。所以要适当选取
11、h才能发挥PC的优点。设绝对稳定区域:设绝对稳定区域:达到精度的校正次数为达到精度的校正次数为N,则,则h的选取,的选取,应满足:应满足:否则可用否则可用N步显式算法稳定达到目的步显式算法稳定达到目的。h浙江大学研究生学位课程546.2.2 Adams四步四阶四步四阶预估校正算法预估校正算法浙江大学研究生学位课程55 6.2.3 Runge-Kutta 方法方法浙江大学研究生学位课程566.2.3 浙江大学研究生学位课程576.2.3 浙江大学研究生学位课程586.2.3 浙江大学研究生学位课程596.2.3 六个未知数,二个自由六个未知数,二个自由,故可取故可取故:故:浙江大学研究生学位课程
12、606.2.3 N=4:四级四阶四级四阶R-K方法方法 -最常用的古典最常用的古典Runge-Kutta方法。方法。增加计算函数值的次数(级)与提高精增加计算函数值的次数(级)与提高精 度(阶)的关系见下表:度(阶)的关系见下表:表表 6.7 Runge-Kutta方法中方法中级级与与阶阶的关系的关系浙江大学研究生学位课程616.2.3 表表 6.8 各种解法在例题中的结果比较各种解法在例题中的结果比较浙江大学研究生学位课程626.2.3 (2).单步法,自动起步单步法,自动起步 (3).易改为变步长易改为变步长 (4).绝对稳定区域较同阶线性多步法大绝对稳定区域较同阶线性多步法大 (5).计
13、算工作量较大,有时大于隐式方法计算工作量较大,有时大于隐式方法 (6).估计误差不易估计误差不易 绝对稳定性讨论绝对稳定性讨论 浙江大学研究生学位课程636.2.3 表表 6.9 各级各级R-K方法的绝对稳定区域方法的绝对稳定区域浙江大学研究生学位课程646.2.3 表表 6.10 不同步长对精度的影响不同步长对精度的影响浙江大学研究生学位课程656.2.3 浙江大学研究生学位课程666.2.3 浙江大学研究生学位课程67 P阶图图 6.10 变步长Runge-Kutta方法框图6.2.3 浙江大学研究生学位课程686.2.3 浙江大学研究生学位课程696.2.4 出发值的计算出发值的计算浙江
14、大学研究生学位课程706.2.4 浙江大学研究生学位课程71质量控制质量控制 Runge-Kutta 步进程序步进程序SUBROUTINE RKQC(Y,DYDX,N,X,HTRY,EPS,YSCAL,HDID,HNEXT,DERIVS)PARAMETER(NMAX=10,PGROW=-0.20,PSHRINK=-0.25,FCOR=1./15.ONE=1.0,SAFETY=0.9,ERRCON=6.E-4EXTERNAL DERIVSDIMENSION Y(N),DYSX(N),YSCAL(N),YTEMP(NMAX),YSAV(NMAX),DYSAV(NMAX)XSAV=X 保留初值保留初
15、值DO 11 I=1,N YSAV(I)=Y(I)DYSAV(I)=DYDX(I)H=HTRYHH=0.5 HCALL RK4(YSAV,DYSAV,N,XSAV,HH,YTEMP,DERIVS)X=XSAV+HHCALL DERIVS(X,YTEMP,DYDX)CALL RK4(YTMP,DYDX,N,X,HH,Y,DERIVS)X=XSAV+HIF(X.EQ.XSAV)PAUS E 步长无意义步长无意义CALL RK4(YSAV,DYSAV,N,XSAV,HH,YTEMP,DERIVS)ERRMAX=0.DO 12 I=1,N YTEMP(I)=Y(I)-YTEMP(I)ERRMAX=MA
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 微分方程 方程组 解法
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,个别因单元格分列造成显示页码不一将协商解决,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【a199****6536】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【a199****6536】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。