数值分析方法第七章.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