Главная страница

Дифференциальных


Скачать 157.37 Kb.
НазваниеДифференциальных
Дата08.06.2022
Размер157.37 Kb.
Формат файлаdocx
Имя файлаA.YU.-Krajnov-K.M.-Moiseeva-CHislennye-metody-resheniya-kraevyh-.docx
ТипДокументы
#579564
страница6 из 28
1   2   3   4   5   6   7   8   9   ...   28

(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) можно отбросить.

Уравнение примет вид:

i1

Oh2 ,


i
т. е. погрешность на каждом шаге увеличивается на величину O(h2).

При нахождении решения в точке xn, отстоящей на конечном рас- стоянии L от точки x0, погрешность состоит из n слагаемых O (h2). Если учесть, что h=L/n, то для погрешности δnполучаем окончательное вы-

ражение:




nOh2 LOh2 Oh.


(1.10)


n h

Отсюда следует, что метод Эйлера имеет первый порядок точности.

    1. Методы Рунге-Кутта



Рассмотренный метод Эйлера (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,..,



i1

i 6 0 1 2 3

k0 hfxi, yi,




1
khf xi

  • h, y

2 i

k



0 ,

2
1   2   3   4   5   6   7   8   9   ...   28


написать администратору сайта