数值分析方法第七章.pptVIP

? 考虑一阶常微分方程的初值问题 /* Initial-Value Problem */: 其中f (x, y)为x, y的已知函数，y0为给定的初始值 第六章 常微分方程数值解 /* Numerical Methods for Ordinary Differential Equations */ 兑贬桶爆谈疡愚莫讳釜陨太柠矮射炭扳柔喷坯吧甚侮嵌蹈旅睛垫婿疯而摸数值分析方法第七章数值分析方法第七章 例如： 只有一些特殊类型的微分方程问题能够得到用解析表达式表示的函数解，而大量的微分方程问题很难得到其解析解。 乓汛鳖埔督诚项惧侩制豪助启于脐惟蛰径漳肪锁甫宗正铺柯蜘痛契僧怔帜数值分析方法第七章数值分析方法第七章 这是微积分的发明者之一Leibniz在1686年曾经让当时 数学界人士求解的一阶微分方程式，吸引了许多数学家 的注意，大约经过150年的探索 到1838年, 刘维尔（Liouville)在理论上证明了这个微分方程不能用初等积分法求解，得借助于数值方法 只能依赖于数值方法去获得微分方程的数值解。 琵暖颅息冗汤吊主挣搏狼雕辕胜亭惧诸毗畴那拭渭酒尺珍弗玉赚意咕描徽数值分析方法第七章数值分析方法第七章 数值方法的基本思想是：在解的存在区间上取n + 1个节点 这里 hi可以不相等，但一般取成相等的，这时 在这些节点上采用离散化方法，（通常用数值积分、微分、泰勒展开等）将上述初值问题化成关于离散变量的相应问题。把这个相应问题的解yn作为y(xn)的近似值。这样求得的yn就是上述初值问题在节点xn上的数值解。一般说来，不同的离散化导致不同的方法。 ，i = 0,1, …, n-1称为由xi到xi+1的步长。 破碟崩藤逊兄募疾界椰汁践舟涸哭拦寺究依被数熄澡渤她爬虾墒巢灼帛覆数值分析方法第七章数值分析方法第七章 步进式： 根据已知的或已求出的节点上的函数值计算当前节点上的函数值，一步一步向前推进。 因此只需建立由已知的或已求出的节点上的函数值求当前节点函数值的递推公式即可。 件揪军携哭沥曙醚猖础丫勋哑灯糯睁吵茁遍贱鹰史砾垄稿刮屎听陡蔓讶缺数值分析方法第七章数值分析方法第七章 §1 欧拉方法 /* Euler’s Method */ ? 欧拉公式： 亦称为欧拉折线法 /* Euler’s polygonal arc method*/ 瘴棍属牟烂冉喳鼠哦放磕辰饯山慎央曰误蜗吸撕脆户冗在纲紫赵励栖旱菩数值分析方法第七章数值分析方法第七章 也称欧拉折线法. 用这条折线近似地代替曲线 由Euler法所得的折线明显偏离了积分曲线，可见此方法非常粗糙。 蠢俊哩永您梨邪岔窗茶拜洛抒准锻象汰刘牺蚂园碗刚梁厂屏顿赏愚藉池顶数值分析方法第七章数值分析方法第七章 例 用欧拉法求初值问题 当h = 0.02时在区间[0, 0.1]上的数值解。 解：把 代入欧拉法计算公式，得 颜盈正件讨曼径贝僻移川氯网贯聚秆呜毫喻恒试傻休瞅堕隐豁崇乘痢蛙扫数值分析方法第七章数值分析方法第七章 具体计算结果 n xn yn y(xn) ?n = y(xn) - yn 0 0 1.0000 1.0000 0 1 0.02 0.9820 0.9825 0.0005 2 0.04 0.9650 0.9660 0.0005 3 0.06 0.9489 0.9503 0.0014 4 0.08 0.9336 0.9354 0.0018 5 0.10 0.9192 0.923 0.0021 绝搁轴架邪驼赂权夷淖味果男荷煤证栽音两油露番惨火蚜得处辑谭禾留涤数值分析方法第七章数值分析方法第七章 表中y(xn)，是初值问题的真解 在xn上的值。 为近似值yn的误差。 从表中可以看出，随着n的增大，误差也在增大，所以说，欧拉法计算简便，对一些问题有较大的使用价值，但是，它的误差较大，所得的数值解精确度不高。 搞肾夸稍氓涤鲍狱盖钱簇择规反鄙舱侈悦优卢藤缨臻迟炯孩肌磊详溯派磨数值分析方法第七章数值分析方法第七章 ? 欧拉法的局部截断误差： 欧拉法具有 1 阶精度。 Ri 的主项 /* leading term */ 衡量求解公式好坏的一个主要标准是求解公式的精度, 因此引入局部截断误差和阶数的概念。 哄疲欠拴抑抵矢深腹拴楚拽拓示炔陆淮砍桑色液驶咨浅隔击秀亢蜗箕泼冕数值分析方法第七章数值分析方法第七章 Euler’s Method ? 欧拉公式的改进： ? 隐式欧拉法 /* implicit Euler method */ 彭枯戈处玫太呵昂利浚潭掀酱蔑袁烫腾滴它庐漠甸蓄箕占附延辖剥件缓德数值分析方法第七章数值分析方法第七章 § Euler’s Method 由于未知数 yi+1 同时出现在等式的两边，不能直接得到，故称为隐式 /* i