Дифференциальных
Скачать 157.37 Kb.
|
(1.9)Таким образом, погрешность δi+1 отличается от погрешности δiна два слагаемых: O (h2) есть следствие погрешности аппроксимации (1.5), a hO (δi) есть следствие неточности значения yi. При нахождении y1 начальное значение y0 задается, как правило, точ- но: δ0 = 0. Отсюда Oh2 , Oh2 hOh2 Oh2 2 h Oh2 . 1 2 1 Отсюда видно, что последнее слагаемое в (1.9) можно отбросить. Уравнение примет вид: i1 Oh2 , i т. е. погрешность на каждом шаге увеличивается на величину O(h2). При нахождении решения в точке xn, отстоящей на конечном рас- стоянии L от точки x0, погрешность состоит из n слагаемых O (h2). Если учесть, что h=L/n, то для погрешности δnполучаем окончательное вы- ражение: nOh2 LOh2 Oh. (1.10)n h Отсюда следует, что метод Эйлера имеет первый порядок точности. Методы Рунге-КуттаРассмотренный метод Эйлера (1.5) является частным случаем мето- дов первого и второго порядков, относящихся к классу методов Рунге- Кутта. Эти методы применяют для вычисления значения yi+1, (i= 0, 1, …) через yiи f(x,y), определенных при некоторых специальным обра- зом выбираемых значениях x x, x и y(x). На их основе могут i i 1 быть построены разностные схемы разного порядка точности. Одним из наиболее часто используемых методов является метод Рунге-Кутта чет- вертого порядка. Алгоритм метода записывается в виде y y 1 k 2k 2k k, i 0,1,.., i1 i 6 0 1 2 3 k0 hfxi, yi, |