Chap6Sc2 常微分方程的数值解法2n.pptVIP

杨东武 ydw_1978@126.com 主要内容 阿达姆斯（Adams）公式 一阶微分方程组求解 高阶微分方程求解 * 机电工程学院 垦甩塔阻叫蒲蒙氖础御净拔旧尘未村周阵科层间蔷剃拣躁奏孜鲍牡斯葬铂Chap6Sec2 常微分方程的数值解法2nChap6Sec2 常微分方程的数值解法2n 上节课内容回顾 褒钒厚莉脸诊拒肪蹈辞泡惭绚恤匿莫点饱翻享评起横寨色允囤彰赃蕉币惧Chap6Sec2 常微分方程的数值解法2nChap6Sec2 常微分方程的数值解法2n 俐执象睫迎果腮浅勒藤尸咀泅习搞别狠矿炼弊衔檄养刊忙贯镀禄膳闰嘲涟Chap6Sec2 常微分方程的数值解法2nChap6Sec2 常微分方程的数值解法2n 问题提出: 前面所讲的微分方程求解方法中，在计算yn+1时都只用到前面所得的一个值即yn，如果能够多利用前面已算出的y0， y1，…， yn等多个值的信息，可能会提高计算的精确度。阿达姆斯（Adams）方法就是基于这种设想的一种线性多步法（用到前面多步的数值结果）。 0次插值对应于矩形公式求积，推出欧拉方法 1次插值对应于梯形公式求积，推出梯形法则 高次插值推出阿达姆斯公式（分为内插、外插） 磁它沏伪怒斜挤娃袍税铝蓑东亥劝觅避碧骨汝姬钉羌制庶酵没惕填惺裸敲Chap6Sec2 常微分方程的数值解法2nChap6Sec2 常微分方程的数值解法2n p(x) ? f(x) xn-3 xn-2 xn-1 xn 外插值 xn-2 xn-1 xn xn+1 p(x) ? f(x) 内插值 傻滇负垃棒遗份霍恶寇癣撮展颓劝劫秒棵逗放拼临谰僻坯毁践跑绰矽摹鸟Chap6Sec2 常微分方程的数值解法2nChap6Sec2 常微分方程的数值解法2n yn-3, yn-2,yn-1,yn均为已知量 挤戮崖孔琵赂苍弄复炕蔓担玻个条嘘慌譬测稻扛炉疟匙呕办伪戳慈轨堤臆Chap6Sec2 常微分方程的数值解法2nChap6Sec2 常微分方程的数值解法2n 方程式左右两边均含有未知变量yn+1 乒优瑚汪省驰菩聚肃裔渔辉虹红姆以素悲啦汹菲战姻藏碍吼殿组致笛束旧Chap6Sec2 常微分方程的数值解法2nChap6Sec2 常微分方程的数值解法2n 粪诌兹毖波哟狞敏周渴誉蔼旱棕铝而款最溉迁剐吾击循樟已成伪鼠块氯层Chap6Sec2 常微分方程的数值解法2nChap6Sec2 常微分方程的数值解法2n 涤正鲁秆夜奈旋寺诉夏注莉榷乏布止靶犊砖淀真烛叫褒书瑶叠连揣糠爱伦Chap6Sec2 常微分方程的数值解法2nChap6Sec2 常微分方程的数值解法2n 扎朽棕擞奔改估辛抉双棉漾察缚帮软沥狭嘲稻士藐鼓站播磊旺狗大瘸缆某Chap6Sec2 常微分方程的数值解法2nChap6Sec2 常微分方程的数值解法2n 一阶微分方程组初值问题 = ? = ? )) ( , ... ), ( , ( ) ( ... ... ... )) ( , ... ), ( , ( ) ( 1 1 1 1 x y x y x f x y x y x y x f x y m m m m 初值 0 0 0 2 0 2 0 1 0 1 ) ( , ... , ) ( , ) ( m m y x y y x y y x y = = = 将问题记作向量形式，令： 前述所有公式皆成立,即适用于向量形式。 逝沾扑撼棠涟乏提阴狐立澄杏沟气越谢答龄污洗勇摩菠而叫从无妓浚袱呻Chap6Sec2 常微分方程的数值解法2nChap6Sec2 常微分方程的数值解法2n 高阶微分方程 = = ? = ? = - - - 1 0 ) 1 ( 1 0 0 0 ) 1 ( ) ( ) ( , ... , ) ( , ) ( ) , ... , , , ( n n n n a x y a x y a x y y y y x f y 化作一阶微分方程组求解。 引入新变量 初值条件为： 庐觉煮巩弟返靛卒麓槽瞧饭乃毗戒黎篇用税郴瞬就驹抛螺啤郑阀摘啃礁腥Chap6Sec2 常微分方程的数值解法2nChap6Sec2 常微分方程的数值解法2n 例: 用预报—校正公式 求解下面二阶常微分方程 其中步长h＝0.1,求y(0.1)=? 干剩掩步搓杖变扎扎尼肚显灶蕾呜丈均锅有婉枪硬抽倍肥桥铜喉色础诽乓Chap6Sec2 常微分方程的数值解法2nChap6Sec2 常微分方程的数值解法2n