第9章常微分方程初值问题数值解法试卷.ppt

* * 消去 ，上式简化为： （9.39） （9.40） 上述方法同样可以用来处理三阶或更高阶的微分方程（或方程组）的初值问题 * * 例9.7 求解下列二阶微分方程的初值问题 取步长h=0.1 解:先作变换：令 ，代入上式，得一阶方程组 用四阶龙格-库塔方法求解,按式(9.37)及(9.38)进行计算： 取步长 , , , 时 * * * * 然后计算 时的 y2和z2；依此类推，直到i=9时的y10和z10，即可得到其数值解。 * * 本章小结 本章介绍了常微分方程初值问题的基本数值解法。包括单步法和多步法。单步法主要有欧拉法、改进欧拉法和龙格—库塔方法。多步法是亚当姆斯法。 它们都是基于把一个连续的定解问题离散化为一个差分方程来求解，是一种步进式的方法。用多步法求常微分方程的数值解可获得较高的精度。 实际应用时，选择合适的算法有一定的难度，既要考虑算法的简易性和计算量，又要考虑截断误差和收敛性、稳定性。 * * 龙格-库塔法较为常用，适用于多步方法中作初值计算和函数f(x,y)较为简单的场合。四阶标准龙格—库塔法精度高，程序简单，易于改变步长，比较稳定，也是一个常用的方法，但计算量较大。当函数f(x,y)较为复杂，可用显式亚当姆斯方法或亚当姆斯预测—校正方法，不仅计算量较小，稳定性也比较好，但不易改变步长。 一般采用龙格—库塔法提供初值y1, y2, y3，然后用亚当姆斯外推公式求得预测值 ，再由亚当姆斯内插值求得校正值yi+1，如此求得的值近似程度好且节省计算量，是一种较好的方法。 * * 改进的欧拉公式 * * 四阶经典 R—K公式 （9.20） * * 作业 习题九 P316 1、2、5(1) * * Thank you very much! 人有了知识，就会具备各种分析能力， 明辨是非的能力。 所以我们要勤恳读书，广泛阅读， 古人说“书中自有黄金屋。 ”通过阅读科技书籍，我们能丰富知识， 培养逻辑思维能力； 通过阅读文学作品，我们能提高文学鉴赏水平， 培养文学情趣； 通过阅读报刊，我们能增长见识，扩大自己的知识面。 有许多书籍还能培养我们的道德情操， 给我们巨大的精神力量， 鼓舞我们前进。 This is a placeholder for the demo. It reminds you when to switch over to the demo and it tells the audience why you are going to show them what you are showing both before the demo and when you switch back to the slides. * * * 例9.4 取步长 h=0.2，用经典 R-C 格式求解初值问题 解: 由四阶龙格-库塔公式可得 * * 可同样进行其余 yi 的计算。本例方程的解为 ，从表中看到所求的数值解具有4位有效数字。 龙格—库塔法的推导基于Taylor展开方法，因而它要求所求的解具有较好的光滑性。如果解的光滑性差，那么，使用龙格—库塔方法求得的数值解，其精度可能反而不如改进的欧拉方法。在实际计算时，应当针对问题的具体特点选择合适的算法。 * * 以经典的四阶龙格-库塔法(9.20)为例。从节点xi出发，先以h为步长求出一个近似值，记为 ，由于局部截断误差为 ，故有 当h值不大时，式中的系数c可近似地看作为常数。 9.3.6 变步长的龙格-库塔法 在微分方程的数值解中，选择适当的步长是非常重要的。单从每一步看，步长越小，截断误差就越小；但随着步长的缩小，在一定的求解区间内所要完成的步数就增加了。这样会引起计算量的增大，并且会引起舍入误差的大量积累与传播。因此微分方程数值解法也有选择步长的问题。 * * 然后将步长折半,即以为 步长,从节点xi出发,跨两步到节点xi+1,再求得一个近似值 ,每跨一步的截断误差是 ,因此有 这样 由此可得 这表明以 作为 的近似值，其误差可用步长折半前后两次计算结果的偏差 来判断所选步长是否适当 * * 当要求的数值精度为ε时： （1）如果Δε，反复将步长折半进行计算，直至Δε为止,并取其最后一次步长的计算结果作为 （2）如果Δε，反复将步长加倍，直到Δε为止，并以上一次步长的