计算方法—-常微分方程的数值解法3.16.pptVIP

主 要 内 容 常微分方程表示方法 常微分方程解法回顾 在微分方程中, 自变量的个数只有一个, 称为常微分方程. 自变量的个数为两个或两个以上的微分方程叫偏微分方程。 简单的数值方法与基本概念 构造高阶单步法的直接方法 由Taylor公式： 当h充分小时，略去Taylor公式余项，并以yi、yi+1分别代替y(xi)、y(xi+1)，得到差分方程： 其局部截断误差为： Runge?Kutta方法 1. 基本思想 因为 其中K?=f(?,y(?))称为y(x)在[xi,xi+1]上的平均斜率。 若取K1=f(xi,y(xi))——Euler公式 取K2=f(xi+1,y(xi+1))——向后Euler公式 一阶精度 取 —— 梯形公式 二阶精度 猜想：若能多预测几个点的斜率，再取其加权平均作为K?，可望得到较高精度的数值解，从而避免求f 的高阶导数。 2. R?K公式 其中Kj为y = y(x)在xi + ajh (0 ? aj ? 1)处的斜率预测值。 aj，bjs，cj为特定常数。 3. 常数的确定 确定的原则是使精度尽可能高。 以二阶为例：