Метод Эйлера — разные подходы к построению

Учитывая ключевую позицию, которую занимает метод Эйлера в теории численных методов ОДУ, рассмотрим несколь­ко способов его вывода. При этом будем считать, что вычисления проводятся с расчетным шагом расчетными точками (узлами) служат точки хi = х0 + ih (i=0,1,..., п) про­межутка [x0,b]ицелью является построение таблицы

приближенных значений yi решения у = у(х) задачи (1)-(2) в расчетных точках хi.

Геометрический способ.Пользуясь тем, что в точке х0 из­вестно и значение решения у(х0) = у0 (согласно (2)), и значе­ние его производной у'(xо)=f(x0, y0) (согласно (1)), можно записать уравнение касательной к графику искомой функции y= у(х) в точке (х0;y0):

(6)

При достаточно малом шаге h ордината

(7)

этой касательной, полученная подстановкой в правую часть (6) значения по непрерывности должна мало от­учаться от ординаты у(х1) решения у(х) задачи (1)-(2). Следовательно, точка (х1, у1) пересечения касательной (6) с прямой x=x1 может быть приближенно принята за новую на­чальную точку. Через эту точку снова проведем прямую

которая уже приближенно отражает поведение касательной к у=у(х) в точке (х1,у(х1)). Подставляя сюда х=х2(=х1+h), иначе, пересекая эту «касательную» прямой х = х2, получим приближение значения у(х2) значением

и т.д. В итоге этого процесса, определяемого формулой

i = 0,1,..., п (8)

и называемого методом Эйлера, график решения у = у(х) дан­ной задачи Коши (1)-(2) приближенно представляется ло­маной, составленной из отрезков приближенных касатель­ных (рис. 1), откуда происходит другое название метода (8) — метод ломаных.

Рис. 1. Геометрическая интерпретация метода Эйлера

 

Применение формулы Тейлора.Описываемый здесь спо­соб вывода метода Эйлера тесно связан с предыдущим. Линеаризуя решение вокрестности начальной точки по формуле Тейлора, имеем

Отсюда при х = х1получаем

(9)

Точное равенство (9), переписанное в виде

говорит о том, что здесь мы имеем одновременно как саму фор­мулу Эйлера для вычисления значения (сравните с формулой (7)), так и ее остаточный член

(10)

где ξi — некоторая точка интервала (x0, x1).

Остаточный член (10) характеризует локальную (или, иначе, шаговую) ошибку метода Эйлера, т.е. ошибку, совершае­мую на одном шаге. Очевидно, что от шага к шагу, т.е. при мно­гократном применении формулы (8), возможно наложение ошибок. За п шагов, т.е. в точке b, образуется глобальная ошибка. Порядок глобальной ошибки (относительно шага h) на единицу ниже, чем порядок локальной ошибки, а по­рядком глобальной ошибки и определяется порядок соответст­вующего численного процесса решения задачи Коши. Таким обра­зом, локальная ошибка метода Эйлера, согласно (10), есть глобальная — 0(h), т.е. метод Эйлера относится к ме­тодам первого порядка.

Разностный способ.Рассматривая уравнение (1) в точке x0 с учетом (2) имеем равенство

Применяя к его левой части аппроксимацию производной пра­вым разностным отношением первого порядка точности

получаем

что идентично равенству (9), поставляющему формулу для вычисления у1 вида (7) и локальный остаточный член (10). Ясно, что для получения общей расчетной формулы (8) мож­но было сразу применить аппроксимацию по формуле (6.16) в равенстве

(11)

заменив неизвестное точное значение y(xi) известным прибли­женным значением yi.

Заметим, что порядок получающегося таким способом ме­тода численного интегрирования дифференциальной задачи (1)-(2) совпадает с порядком аппроксимации производ­ной в левой части уравнения (1).

Квадратурный способ.Как было показано на­чальную задачу для ОДУ (1)-(2) можно заменить эквива­лентным интегральным уравнением (3). При х=х1 из него получится равенство

(12)

 

Применение к интегралу в правой части равенства (12) про­стейшей (одноточечной) формулы левых прямоугольников дает приближенную формулу

правая часть которой, очевидно, совпадает с выражением (7) для подсчета значения у1. В общем случае расчетная форму­ла (8) метода Эйлера получается численным интегрированием посредством простейшей формулы левых прямоугольников в ра­венстве

(13)

в предположении, что на каждом i-м шаге в роли начальной точки (x0,y0) выступает точка (xi, yi). Зная точность исполь­зуемой в (13) квадратурной формулы, легко прийти к тому же выражению локальной ошибки метода Эйлера, что и при других способах его построения.

Существуют и другие подходы к выводу метода Эйлера. В частности, он будет возникать далее как частный случай некото­рых семейств численных методов решения задачи (1)-(2).