Описания методов решения:, Метод Эйлера - Методы решения обыкновенных дифференциальных уравнений

Метод Эйлера

Расчетную формулу метода Эйлера можно получить, используя разложение функции U(X) в ряд Тейлора в окрестности некоторой точки XI:

. (1)

Если приращение H мало (то есть H<<XI), то члены ряда, начиная со слагаемого, включающего H во второй степени, могут быть отброшены как малые величины. Тогда из (3)в первом приближении получим

. (2)

Воспользуемся формулой (4), применив ее к единственной известной из условия задачи точкеX0. Найдем в X0 производную Du(X0)/Dx, подставив (2) в (1):

.

Подставив последнее выражение в (4) и полагая XI= X0, получим

Или, сокращая обозначения, в окончательном виде

.

Таким образом, (4) при известном значении функции U0= U(X0) в начальной точке X0 позволяет найти приближенное значение U1= U(X1) при малом смещении H от X0. На рис. 1 графически показан начальный шаг решения методом Эйлера.

метод эйлера

Рис. 1. Метод Эйлера

Решение можно продолжить, используя найденное значение функции U1 для вычисления следующего значения - U2. Распространяя эти рассуждения на последующие точки, запишем расчетную формулу метода Эйлера в виде

. (5)

Из рис. 1 видно, что ошибка метода Эйлера на шаге связана с используемой линейной аппроксимацией U(X). Хотя тангенс угла наклона касательной к кривой точного решения в точке (X0,U0) известен и равен Du(X0)/Dx, он изменяется при смещении от X0 до X1. Следовательно, при сохранении начального наклона касательной на всем интервале H расчет U1Выполняется с погрешностью.

Ошибка метода Эйлера на каждом шаге имеет порядок H2, так как члены, содержащие H во второй и более высоких степенях, отбрасываются - см. (3) и (4). Уменьшая H можно снизить Локальную ошибку на шаге.

Похожие статьи




Описания методов решения:, Метод Эйлера - Методы решения обыкновенных дифференциальных уравнений

Предыдущая | Следующая