数值分析.南京电大 36讲 szfx32.docVIP

同样用泰勒展开的方法，将f展开到， 将展开到,四个,三个，六个共十 三个未知待定参数,共有十一个独立的方程的二 个自由度。 ， ，， ， ， ， ， 十三个未知量，十一个方程有无穷多组解。 这样就得到格式称为古典龙格－库塔格式，书上 称为标准龙格－库塔方法，其格式为： 得到的格式称为库塔法，格式为： 得到的方法称为吉尔方法。 总之,这三种方法的共同点是截断误差为 阶数为四阶精度,均称为显式四阶龙格－库塔方法， 在从来计算出均要计算四个f的值 一般一阶常微分方程初值问题均用四阶龙格－库 塔方法来计算，其精度均满足了实际问题的精度 要求。 数值例子： 考虑微分方程初值问题：。 解：其精确解为： ，用古典龙格－库塔方法,。 (1)、求，此时 (2)、求，此时 以下计算用表格列出： 误差 1.1 1.20027 1.20027 3(10-7 2.7(10-4 1.2 1.40182 1.40182 5(10-7 1.06(10-3 1.3 1.60524 1.60524 6(10-7 2.08(10-3 1.4 1.81079 1.71079 6(10-7 3.22(10-3 1.5 2.01856 2.01856 7(10-7 4.42(10-3 1.6 2.22855 2.22855 7(10-7 5.66(10-3 1.7 2.44070 2.44070 8(10-7 6.90(10-3 1.8 2.65497 2.65497 8(10-7 8.16(10-3 1.9 2.87127 2.87127 8(10-7 9.40(10-3 2.0 3.08953 3.08953 8(10-7 10.65(10-3 例：初值问题用四阶古典 Runge－Kutta方法，。 (1)、求， (2)、求， ， 0.2 1.1832293 1.1832160 1.3(10-5 0.4 1.3416803 1.3416408 4.0(10-5 0.6 1.4832838 1.4832397 4.4(10-5 0.8 1.6125172 1.6124515 6.6(10-5 1.0 1.7321463 1.7320508 9.6(10-5 §14.3 亚当斯方法(Adams) 单步方法和多步方法： 前面讲的方法：欧拉方法、改进欧拉方法、 龙格－库塔方法均是单步方法，即在每一步要计 算时,只要前面一个值已知的条件下秒可以 计算出了。 特点： (1)、可以自成系统进行直接计算，因为初始条 件只有一个已知，由可以计算，， ，不必借助于其它方法,这种称为单步 方法是自开始的。 (2)、如果格式简单如欧拉方法，则只有一阶精 度，如果提高精度，则计算很复杂. 如Runge－Kutta方法。 (3)、公式的构造推导也很复杂。 多步方法： 利用前面已知计算出来的，这 前面计算好的个值来计算，这样自由度的增 16 加来提高格式的精度，这样的方法称为多步方法， 利用k个值计算，称为k步方法。 多步方法的特点： (1)、 因初始条件只有一个，运用多步方法设法 开始，要借助高阶的单步方法来开始.例如,已知 用单步的四阶Runge－Kutta方法计算，再 计算，再由计算，用单步方法有 后运用四阶的四步方法，由计算；由 计算；由计算；一直下去 均匀可以用多步方法了，而且始终达到四阶精度. (2)、多步方法比较简单，只要在这四个点的函数 值的线性组合,而且每步中后三个函数值下一步还 可使用。 显式Adams方法： 1. 方法构造的思想： Adams方法是多步方法中的某一类，而不就是多 步方法，考虑微分方程初值问题， 将微分方程在上积分， 下面我们来推导四步显式Adams方法，即若已知 来计算，简记 用的拉格朗日插值多项式代替f