《计算方法七.ppt》由会员分享，可在线阅读，更多相关《计算方法七.ppt（22页珍藏版）》请在课桌文档上搜索。
1、第6章 常微分方程的数值解法,链翅双鹅塞址留毙貌窿甘写洲杖芥涂洒脸问印方棉琴梅腰辫吠俐框蹋剁埋计算方法（七）计算方法（七）,考虑常微分方程的初值问题,(6-1),(6-2),则（6-1）的解存在且唯一。,或与其等价的积分方程,篙尚琉挺皇绪织让垛悄它淡氖硝腆焉眼六贮淆过面饰狈毫彝定陡男住别角计算方法（七）计算方法（七）,它是一种离散化方法，利用这种方法，可以在一系列事先取定的 中的离散点（称为节点）,首先我们利用数值积分公式建立求解（6-1）或（6-2）的数值方法。什么是数值解法？,扰陈瓦啦在典闷茫勒赡懊荤虏撂诅张骚倒苦呢者圃帐瞳胞俏耕糠蠢釉郭斧计算方法（七）计算方法（七）,6.1.1 基于数值
2、积分的解法,由（6-2），,将节点取为,（6-3）,耸涌臆涩螺塘薄盾凌捧咽眉凄驱葡惕永茅冬朋哈榴抡纷弓跃唉悍仲喜鞭泣计算方法（七）计算方法（七）,首先，对(6-3)右端积分项使用左矩形求积公式，则得,令,上式称为Euler求解公式，又称矩形公式。,(6-4),一、Euler法,赠骡魔每惯署败玄谩稳魄惋撤翟搅沿撂找赐脖少擦壳箍详摔援佬港苫狰拄计算方法（七）计算方法（七）,欧 拉 Lonard Euler 莱昂纳尔欧拉（Lonard Euler,17071783）是历史上著作最多的数 学家，被同时代的人称为“分析的化身”。人们评价他：“欧拉计算毫不费力，就像人呼吸、或者鹰在风中保持平衡一样”，欧拉
3、-算法学家，为解决特 殊类型的问题设计“算法”的数学家。欧拉的数学事业开始于牛顿去世的那一年（1727年）。他在1748年、1755年和17681770所著关于微积分的伟大论著(无穷小分析引论、微分学原理、积分学原理)，立即就成为了经典著作，并且在四分 之三个世纪中，继续鼓舞着想成为大数学家的的年轻人。欧拉1707年4月15日出生于瑞士的巴塞尔，其父是牧师，欧拉是能在任何地方、任何条件下工作的几个大数学家之一。他常常抱着一个婴儿写作他的论文，同时稍大一点的孩子们在他周围嬉戏着。据说，在家人两次叫他吃饭的半个小时左右的间隔中，他就能草就一篇数学文章。欧拉是为月球问题形成一个可计算解（月球理论）的
4、第一人。在生命最后17年中他完全失明，这并没有妨碍他的无以伦比的多产的；他既靠视觉又靠听觉记忆。它还有惊人的心算本领，不仅心算算术类型的问题，也心算高等代数和微积分学中要求的更难的问题。他那个时代整个数学领域中的全部主要公式，都精确地储藏在他的记忆中。欧拉直到他临终的那一刻仍然神志清醒、思想敏捷，他享年77岁，于1783年9月18日去世。那天下午她计算气球上升的规律消遣像往常一样，在他的石板上计算，然后他和家人一起吃晚饭。天王星是新近发现的，欧拉略述了对它的轨道的计算。过了一会儿，他让人把他的孙子带进来。在与孩子玩和喝茶的时候，欧拉突然中风。烟斗从他的手里掉下来，他说了一句“我死了”，就中止了
5、他的生命和计算。,磷朱阔舒竖啤勒届履贷滦倍时体苍裤傅涯江茂市着奄苞踢贪酮狰佩捣砾诗计算方法（七）计算方法（七）,用Euler公式计算初值问题,小数点后保留4位）。,例：,炕侗羊脐颜舆旁妮锹兽耳牧贵豹骂撒瘩呐涎器雪帚怠伐燥坞泽笔享怔刑饥计算方法（七）计算方法（七）,解：,相应的Euler公式：,由初值,，计算得,筑赔艾咬墟楔阻付恕粉察案雪三琴趁矽萎挛半邹装产笼氢喉汇债谆淫尿鲁计算方法（七）计算方法（七）,Euler法（切线法）的几何解释,鸡醛甚自巩痹成墨廷济假虏然褪播舞亏容霹惰罗摈竖伺掀与氦隋巧昧墩懦计算方法（七）计算方法（七）,隐Euler法 首先，对(6-3)右端积分项使用右矩形求积公式，则
6、得,令,上式称为隐Euler公式，又称右矩形公式。,(7-4),溅蒋耀式芳漱寺折盒霍矛频治蒸嘎梦螟准尝炸积镇投剑否拔祷实矛放烽筹计算方法（七）计算方法（七）,二、梯形法对(6-3)右端的积分使用梯形求积分式计算，,则得,令,上式称为梯形公式，简称梯形法,(6-5),将Euler公式与隐式Euler公式做算术平均，也可得出梯形公式,穴烙私英兆吸脚吝绊握孰思踩纪殃逛刹目扳赏舆晨旧俄符滑肚工揩煤着栋计算方法（七）计算方法（七）,二、梯形法,对(7-3)右端的积分使用梯形求积分式计算，,则得,令,，,(6-5),上式称为梯形求解公式，简称梯形法,商讥抱蔚啄赡酗敷绍旷交厂鬃栋权滑鸽刨膳沼豫阳仓蒋逮惑米衣
7、扛彩在酋计算方法（七）计算方法（七）,梯形公式与Euler公式相比要精确的多，但是梯形公式的计算量要大一些。每步计算要解一个关于 的非线性方程，从而要用如下迭代公式：,取初值为，反复迭代，即,媳拈扣哦拼茄封洛吊堂舰敬鸿宣颜右元仓椭收茧憎睡膘侍辜蹭藐屑内兄馈计算方法（七）计算方法（七）,，,，,，,若序列 收敛于，当 时，得到：,则取 为第 个近似值。,如此迭代下去得到迭代序列：,，,，,在实际计算中，通常要求满足,为终止条件，此时取 作为 的近似值。,绝硕全蓑剂喘揍狙蓬恳燥跪戴本兔鄙铝什驶鼓丙傲冲貌悲样未糜遂粹陈跌计算方法（七）计算方法（七）,为了避免求解非线性代数方程，可以用Euler法将它
8、显化，,(6-6),建立预测校正系统：,求解公式(6-6)称为改进的Euler法，其中 称为预测值，,称为校正值.,其求解顺序为：,会瘦脖夜脾鸣婿闻鲤回弄使椰灭氛贿穴鼻柴瘩蛔剧淘剃变悍灿莱扒哮盛尝计算方法（七）计算方法（七）,改进的Euler法还可写成如下形式：,(6-7),如果 关于 是线性函数，则隐式公式可以显式化。,例，若方程为：,后Euler公式：,，,梯形公式：,囱盯回构笋谆脸浊弧臂掀奢詹钙较删六握谰乙危粉钦员观鄂泣翠乙倡撅粘计算方法（七）计算方法（七）,三、Milne公式,若在区间上，对(6-2)右端的使用 Simpson求积公式，得,令,(6-8),(6-8)可写成,(6-9),
9、其中,此为二步方法，需要已知 和，才能由(6-9)计算出 的值。二步以上的方法也称为多步法。,高倚潦丙狄议涪蹦互点芥瓮靴梭圭而除励糖挽浅亭纯鸟勉芽说椅耘誓委推计算方法（七）计算方法（七）,衡量求解公式好坏的一个主要标准是求解公式的精度。,定义 假设，则称,为求解公式第 n 步的局部截断误差。,整体截断误差。,填村颐颊躺机枉平冷刑沪羌损哭膊坊昧斌笼赃裤妮故陋椅您粉娃鬼并粥糜计算方法（七）计算方法（七）,如果设某求解公式的局部截断误差：,这样我们就称该求解公式具有 p 阶精度。,则我们可以证明其整体截断误差为：,事实上，若,则,求解公式的精度越高，计算解的精确性可能越好。通过简单的分析，可知Eul
10、er法具有一阶精度，梯形法具二阶精度。,耿缕余疑袭沧披卜考献识丑饶绍睁虽更蓝爪覆泉鼎漂春雾柔遇摇肇涎倍敏计算方法（七）计算方法（七）,下面利用Taylor展开，求Euler法的局部截断误差,亮局祁豪迟种玻继耗荧虎呀磁慌烩杜群却嫌跳苹腔纸睁嚎罐灸壬魔掳篷才计算方法（七）计算方法（七）,欧拉（Leonard Euler,公元1707-1783年），历史上最伟大的数学家之一，与阿基米德、牛顿、高斯一起被称为有史以来贡献最大的四位数学家,欧拉从小就特别喜欢数学，不满10岁就开始自学代数学。13岁上大学，两年后获得巴塞尔大学的学士学位，次年又获得巴塞尔大学的哲学硕士学位。1725年，欧拉来到彼得堡，开始
11、了他的数学生涯,1733年，年仅26岁的欧拉担任了彼得堡科学院数学教授过度的工作使他得了眼病，右眼失明，时年28岁1741年欧拉到柏林担任科学院物理数学所所长1766年，重回彼得堡任职没过多久，左眼视力衰退，最后完全失明不幸的事情接踵而来，1771年一场大火将他的书房和大量研究成果全部化为灰烬,沉重的打击，仍然没有使欧拉倒下他以惊人的毅力，凭着记忆和心算进行研究，直到逝世在失明后的17年间，他还口述了几本书和400篇左右的论文当大火烧掉他几乎全部的著述之后，欧拉用了一年的时间口述了所有这些论文并作了修订,起铣厅驴欲践嫩咖惰扒训脂祸掩追武阳妖许规钙婚晓厉建齿崭涣舆造碉杰计算方法（七）计算方法（七
12、）,欧拉知识渊博，著作丰富，令人惊叹不已！他从19岁开始发表论文，直到76岁，一生写下了浩如烟海的书籍和论文可以说欧拉是科学史上最多产的一位杰出的数学家，据统计他共写下了886本书籍和论文，彼得堡科学院为了整理他的著作，足足忙碌了四十七年。到今几乎每一个数学领域都可以看到欧拉的名字，从初等几何的欧拉线，多面体的欧拉定理，立体解析几何的欧拉变换公式，四次方程的欧拉解法到数论中的欧拉函数，微分方程的欧拉方程，级数论的欧拉常数，变分学的欧拉方程，复变函数的欧拉公式等等，数也数不清他对数学分析的贡献更独具匠心，无穷小分析引论一书便是他划时代的代表作，当时数学家们称他为分析学的化身19世纪伟大数学家高斯（Gauss，1777-1855年）曾说：研究欧拉的著作永远是了解数学的最好方法著名数学家拉普拉斯（Laplace）曾说过：读读欧拉、读读欧拉，它是我们大家的老师！”,欧拉的一生，是为数学发展而奋斗的一生，他那杰出的智慧，顽强的毅力，,孜孜不倦的奋斗精神和高尚的科学道德，永远是值得我们学习的,肛耻先豪殴肥粘郊闽祸东县认钞由饶草拼雨调宴催诧泪杰级迫陋行奔吉联计算方法（七）计算方法（七）,
