蚂蚁文库
换一换
首页 蚂蚁文库 > 资源分类 > PDF文档下载
 

上海交通大学计算方法课件宋宝瑞CH8

  • 资源ID:55609261       资源大小:535.11KB        全文页数:15页
  • 资源格式: PDF        下载权限:游客/注册会员    下载费用:15积分 【人民币15元】
快捷注册下载 游客一键下载
会员登录下载
三方登录下载: 微信快捷登录 QQ登录  
下载资源需要15积分 【人民币15元】
邮箱/手机:
温馨提示:
支付成功后,系统会自动生成账号(用户名和密码都是您填写的邮箱或者手机号),方便下次登录下载和查询订单;
支付方式: 微信支付    支付宝   
验证码:   换一换

 
友情提示
2、PDF文件下载后,可能会被浏览器默认打开,此种情况可以点击浏览器菜单,保存网页到桌面,既可以正常下载了。
3、本站不支持迅雷下载,请使用电脑自带的IE浏览器,或者360浏览器、谷歌浏览器下载即可。
4、本站资源下载后的文档和图纸-无水印,预览文档经过压缩,下载后原文更清晰   

上海交通大学计算方法课件宋宝瑞CH8

第八章第八章常微分方程初值问题的数值解法常微分方程初值问题的数值解法 dy   f x, yx[x 0 ,b] ((1 1))dx  yx0  y0 ((1 1)的解解析解)的解解析解函数函数 y  yx, dy  f x, yxx[x 0 ,x N ]yx xx  y 0 0 dx 常微分方程课程中讨论了(常微分方程课程中讨论了(1 1)的解的存在性,唯一性条件)的解的存在性,唯一性条件 例如例如f x, yC[x0,b;,],且,且f x, y满足对满足对y的的 LipchitzLipchitz 条件条件 f x, y 1 f x, y2  L y 1  y 2 x[x 0 ,x n ] 则(则(1 1)的解存在,唯一)的解存在,唯一 以后我们总设以后我们总设f x, yLipLy 解析解不易求得,或太复杂。实际问题中归结出的方程主要用数值解,解析解不易求得,或太复杂。实际问题中归结出的方程主要用数值解, 即求即求y  yx在一系列离散点上的近似值,这些点是在一系列离散点上的近似值,这些点是 x 0 , x 1  x 0 h 1, , xn  x n1 h n 诸诸hi可以不同,为方便计算,设可以不同,为方便计算,设h i  h,i 1,2, 方法方法 ⊙据常微分方程理论,已知⊙据常微分方程理论,已知yxk  yk,则(,则(1 1)在)在[xk,b]上的解满足上的解满足 y  f x, y  yx  y k  k 1 提示我们从提示我们从x 0 出发,一步一步向前跨,得到出发,一步一步向前跨,得到y i  yx i , i  0,1,,n 初值问题初值问题 TaylorTaylor展式法展式法 数值积分法数值积分法 EulerEuler 折线法折线法 分点分点 x k  x 0 khk  0,1,2,.,nh  b x 0 n 给定(给定(1 1)) ,在,在x  x k 处将处将yx展成展成 TaylorTaylor展式展式 yx k1  yxk  yx k h 2 1 yh2 2 一般一般h很小,略去很小,略去h项,得项,得 y 1  y 0 hf x 0 , y 0 y 2  y 1 hf x 1, y1 一般地,一般地,yk yk1hf xk1, yk1 k 1,2,, n yx x  y k x x k f x k , y k x[x k ,x k1] 分段线性函分段线性函 数(数(EulerEuler 折线法名称的由来)折线法名称的由来) 如果如果yk yxk(没有误差)(没有误差) y k1 用用 EulerEuler 折线法求得折线法求得 h2 y T k1 局部截断误差局部截断误差则则yxk1 yk1 2 h2 T k1 yx k oh2 2 主项 EulerEuler 折线法算法简单,自开始,但精度差折线法算法简单,自开始,但精度差P.281P.281,表,表 9-19-1,几乎不单独,几乎不单独 2 用。用。 向后的向后的 EulerEuler 公式公式 y k1  y k hf x k1, yk1 h2h2 TaylorTaylor展开可得,展开可得,T k1   2 y, , 主项主项 2 yx k 隐式,隐式,y k1 可迭代求解,精度也不高。可迭代求解,精度也不高。 梯形公式(梯形公式(向前、向后向前、向后 EulerEuler 法,取算术平均)法,取算术平均) yy h k1  k  2 [ f x k , y k  f x k1, yk1] ○○ 1 1 平均斜率平均斜率○○ 2 2 消去截断误差中的消去截断误差中的h2项。提高精度项。提高精度T 3 k1 Oh 隐式,迭代方法隐式,迭代方法 y0 k1  y k hf x k , y k yn1 h n k1  y k  2 [ f x k , y k  f x k1, yk1] 迭代有限步,或迭代至收敛(收敛吗下证)迭代有限步,或迭代至收敛(收敛吗下证) 2-32-3 y n1 h k1  y k1  2 [ f x n k1, yk1 f xk1, yk1] LipchitzLipchitz 条件条件 y n hL k1  y k1  2 y n1 k1  y k1 当当h充分小,即充分小,即 hL 2 1时,方法收敛,缺点迭代次数无法控制。 时,方法收敛,缺点迭代次数无法控制。 如果只迭代一次,得到改进的如果只迭代一次,得到改进的 EulerEuler 公式公式 2 3 3 y k1  y k hf x k , y k T k1 Oh2 y h 3 k1  y k  2 [ f x k , y k  f x k1, yk1] T k1 Oh 预估预估校正法校正法 y k  f x k , y k   y k1  y k1 2 预估  hf x k , y k   y T k1 Oh3 k1  f x k1, yk1 h  校正 yk1  y k  2 y k1  y 3 k T k1 Oh y k1  y k hy k  h2 2 y k Oh3 2 说明说明y   h Oh3 k1  y k hy k 2 y k y k1  y k1 2hf x k , y k Oh3 T k1 Oh3 优点优点 预估与校正精度相同;不需迭代,精度较高。预估与校正精度相同;不需迭代,精度较高。 问题问题 已知已知y0, y1才可起步,要用其它方法做“表头”才可起步,要用其它方法做“表头” EulerEuler 法的整体误差法的整体误差 e n  yx n  y n ,,受第受第 1 1,,2 2,,第第n步截断误差的影响步截断误差的影响 记记yn yxn1hf xn1, yxn1,则,则 4 e n  yx n  y n  n  y  n  y n  yx n  y  T n  yx n1hf xn1, yxn1 yn1 hf x n1, yn1  T n  yx n1 yn1 h f x n1, yxn1 f xn1, yn1  T n 1hL e n1 反复应用上式

注意事项

本文(上海交通大学计算方法课件宋宝瑞CH8)为本站会员(sunhongz117)主动上传,蚂蚁文库仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知蚂蚁文库(发送邮件至2303240369@qq.com或直接QQ联系客服),我们立即给予删除!

温馨提示:如果因为网速或其他原因下载失败请重新下载,重复下载不扣分。




关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们


网站客服QQ:2303240369

copyright@ 2017-2027 mayiwenku.com 

网站版权所有  智慧蚂蚁网络

经营许可证号:ICP备2024020385号



收起
展开