常微分方程的数值解及实验.pptVIP

y(xi+1)≈y(xi)+hy′(xi) =y(xi)+hf(xi，y(xi)) =yi+hf(xi，yi) 若取前三项，可得到截断误差为O(h3)的公式 第三十页，共六十九页，2022年，8月28日 这里 y′(xi)=f(xi，y(xi)) y″(xi)=f′x(xi，y(xi))+f′y(xi，y(xi))y′(xi) =f′x(xi，y(xi))+f(xi，y(xi))f′y(xi，y(xi)) 类似地，若取前k项作为y(xi+1)的近似值，便得到截断误差为O(hk)的数值计算公式。这些公式的计算必须依赖于求y(xi)的k阶导数，除非f(x，y)足够简单，否则直接用泰勒展开法求解较为复杂。但是泰勒级数展开法的基本思想是许多数值方法的基础。 第三十一页，共六十九页，2022年，8月28日 3.2 龙格库塔法 前面已经知道，初值问题(4―1)等价于 龙格库塔法的基本思想是：用f(x，