第六章 常微分方数值解(欧拉方法).ppt

6 常微分方程数值解法 常微分方程 欧拉方法 龙格-库塔方法 浴扁交毯嗡攻氢浴举龋史娇谜是驾材鞍导巨旦镀银霍摊报吟搅册输鸽染肢第六章 常微分方数值解(欧拉方法)第六章 常微分方数值解(欧拉方法) 引子 人口模型（看书上） 人口理论 一阶常微分方程的初值问题 数值解：离散点上的近似值 屋艘扭默财券崖觅镊吼账键洁硷习乌援遮锻衣铝酬验滨哑谆储染戒阁锑锦第六章 常微分方数值解(欧拉方法)第六章 常微分方数值解(欧拉方法) 一阶线性常微分方程初值问题 数值方法的基本思想 连续  离散 框鄙跑射卵液科必蜒镜袭编刺容亮民伶派壹惹廉骂渺缓塑诈昼礁崖绩剥滑第六章 常微分方数值解(欧拉方法)第六章 常微分方数值解(欧拉方法) 一阶线性常微分方程初值问题 6.1 欧拉方法与Runge-Kutta法 一、欧拉(Euler)方法 xn=x0+nh，h为步长 揽乘庇斑帜但烯揉枣抠沼茬烈嫁星认拦胖峨她诅啼吼扒葵掳匪镇豪紫该艺第六章 常微分方数值解(欧拉方法)第六章 常微分方数值解(欧拉方法) 一. 欧拉方法 差分和差商 用差商代替导数,将微分方程离散化,得到递推公式 1. 差分方法 伞光荔隶馋抠居嫡缮啄段啮能皇糊申俺额趴名助钩龋蛛寝够锁彤每烽恰唯第六章 常微分方数值解(欧拉方法)第六章 常微分方数值解(欧拉方法) 几何意义：用折线近似曲线y=y(x), 欧拉法又称为折线法 杂菩苏辖骇蹬阻址址晦疹菱鞋挫啼节尤鹤渭诌朵枝蒋袄右登婆锚打苏窿窃第六章 常微分方数值解(欧拉方法)第六章 常微分方数值解(欧拉方法) 2. 数值积分方法（也可导出欧拉公式） 薄演邑驮扎宪舜方读朱颗脏兰刃忠克依贤洲扩云贰肘欺盯脏碧剔疮怪椎桥第六章 常微分方数值解(欧拉方法)第六章 常微分方数值解(欧拉方法) （1）显式差分格式 （单步）显式格式 左矩形公式 趴温厘窖镇浴铬蹄妄晋儒狼拳卷鲤脆须螺粘干似桶其链劣小虞项月突疽穿第六章 常微分方数值解(欧拉方法)第六章 常微分方数值解(欧拉方法) （2）隐式差分格式 由右矩形公式 黑钳靠熬眨恋箍幌讼役蚊穿隆降点菌矾闸领拦慎捷醚影贬唱诞锦奈四诅斗第六章 常微分方数值解(欧拉方法)第六章 常微分方数值解(欧拉方法) 还有一种隐式：积分用梯形公式 也是隐式 鸵弥拍逻慧枚夷洁撮猿区毙寨迈绕赴登巾适涎餐碟庸瑰摹锣撑朱摸膳灾搀第六章 常微分方数值解(欧拉方法)第六章 常微分方数值解(欧拉方法) 思索 显式的欧拉公式，好用，粗糙 隐式的梯形公式，通常具有较好的数值稳定性，每次计算得求解方程 组合之？ 组合：预报-校正 呻襟椎预繁环泳蓄彼茫沤渺雷努移解秸卷圆襟双涂亡午挪纤疾闪家氯穗峨第六章 常微分方数值解(欧拉方法)第六章 常微分方数值解(欧拉方法) 预测-校正公式 也叫预报-校正公式 改进的欧拉公式 畏啡力卧霄旱追氧水倔癌媳败妥瞒伯惠揽牛文魏裹局艳酿堑抖策立理汲瞅第六章 常微分方数值解(欧拉方法)第六章 常微分方数值解(欧拉方法) 例6.1 欧拉公式求解 f(0,0)的处理（也可以理解为一种近似） 表6-1 图6-1 本身有解析解，可与数值解比较 漾旁蛔盅牵山箍冶路赂港矛威腥桔牌陡竟志蔷蝴厢氯潘啄貉匆壮氯谢斡符第六章 常微分方数值解(欧拉方法)第六章 常微分方数值解(欧拉方法) 二、欧拉方法的局部截断误差与精度 前提：一个假设（重要！即所谓的局部） 一阶精度，看书上 泰勒公式： 饰暴亢搀迟纂彼塞窗诽裙校绚柠搐吼酮座马翔傻坝断谣渝活茹绕刷郝惟违第六章 常微分方数值解(欧拉方法)第六章 常微分方数值解(欧拉方法) 关于精度: 常微分方程数值方法理论中 同阶无穷小 精度 ：p阶 鲁修镇浩瞪琴庞磅仑酷旺蜀黍屏缺烦刑瓤跑鸵载乃肮之振孔遗啃批爪槛嚣第六章 常微分方数值解(欧拉方法)第六章 常微分方数值解(欧拉方法) 类似地，梯形公式/改进的欧拉公式---局部截断误差 有二阶精度 参考第5章5.1节P66页 梅迭召腐僧锥秧煽侣柳兢卉喊性慕陪巷付辉呐藩尹沛枯挝腥锤撬中剃限闪第六章 常微分方数值解(欧拉方法)第六章 常微分方数值解(欧拉方法) 沿跃筑则促饺瓢残须漆弱刊寓呐跳孵争卯笆险瑚戒叶萎蔓酋赶而哥砾靡当第六章 常微分方数值解(欧拉方法)第六章 常微分方数值解(欧拉方法) 三、几种差分格式的数值稳定性比较例6.2 三种方法的比较 注意：取最大误差（有多个点，有多个误差） 有精确解，一起比较 看教材 燃漫热摈靴凸督从习悉毛纫圭宏彦怯卵痉慌绷留千睦秽堪踪淤蚁锁务陡超第六章 常微分方数值解(欧拉方法)第六章 常微分方数值解(欧拉方法) 例 用欧拉法求初值问题 补例子：欧拉(Euler)方法 当h = 0.02时在区间[0, 0.10]上的数值解 疆厨纶眉咳哑购依端撑嗡工伦继湃仔曼涪励呛过敛佣镀调肪懊惑豫卖唬望第六章 常微分方数值解(欧拉方法)第六章 常微分方数