Метод Эйлера¶

Метод Эйлера — это численный метод первого порядка точности для решения задачи Коши для обыкновенных дифференциальных уравнений (ОДУ). Является простейшим представителем семейства одношаговых методов Рунге — Кутты. Метод позволяет приближённо находить значения искомой функции в дискретных точках, заменяя производную конечной разностью.

¶Основная идея

Пусть дано дифференциальное уравнение первого порядка:

\[ y' = f(x, y) \]

с начальным условием \( y(x_0) = y_0 \). Требуется найти функцию \( y(x) \) на отрезке \( [x_0, X] \). Метод Эйлера основан на разложении искомой функции в ряд Тейлора в окрестности точки \( x_n \):

\[ y(x_{n+1}) = y(x_n) + h y'(x_n) + \frac{h^2}{2} y''(\xi), \quad \xi \in [x_n, x_{n+1}] \]

Отбрасывая члены второго порядка и выше, получают приближённую формулу:

\[ y_{n+1} = y_n + h f(x_n, y_n) \]

где \( h = x_{n+1} - x_n \) — шаг интегрирования. Геометрически метод означает движение по касательной к интегральной кривой в каждой точке сетки.

¶Алгоритм

Вычислительный процесс для равномерной сетки с шагом \( h \) состоит из следующих шагов:

1. Задать начальные данные: \( x_0 \), \( y_0 \), шаг \( h \), количество шагов \( N \).
2. Для \( n = 0, 1, 2, \dots, N-1 \) выполнить:

• Вычислить \( k_1 = f(x_n, y_n) \).
• Вычислить \( y_{n+1} = y_n + h \cdot k_1 \).
• Вычислить \( x_{n+1} = x_n + h \).

1. Результатом является последовательность точек \( (x_n, y_n) \), аппроксимирующая истинное решение.

¶Погрешность и сходимость

Метод Эйлера имеет первый порядок точности. Локальная погрешность на одном шаге составляет \( O(h^2) \), а глобальная погрешность на всём отрезке — \( O(h) \). Это означает, что для уменьшения ошибки вдвое необходимо уменьшить шаг в два раза.

¶Факторы, влияющие на точность

• Шаг интегрирования: уменьшение \( h \) снижает погрешность, но увеличивает количество вычислений и накопление ошибок округления.
• Жёсткость системы: для жёстких дифференциальных уравнений метод Эйлера может требовать чрезвычайно малого шага для устойчивости.
• Гладкость правой части: если функция \( f(x, y) \) имеет большие производные, погрешность возрастает.

¶Устойчивость

Явный метод Эйлера является условно устойчивым. Для линейного тестового уравнения \( y' = \lambda y \) (с комплексным \( \lambda \)) область устойчивости — круг радиуса 1 с центром в точке \( (-1, 0) \) на комплексной плоскости \( h\lambda \). При нарушении этого условия численное решение начинает неограниченно расти, даже если истинное решение затухает.

¶Модификации метода

¶Неявный метод Эйлера

Формула неявного метода:

\[ y_{n+1} = y_n + h f(x_{n+1}, y_{n+1}) \]

Значение \( y_{n+1} \) находится из решения алгебраического (или трансцендентного) уравнения на каждом шаге. Неявный метод обладает абсолютной устойчивостью (A-устойчив) для линейных задач, что делает его пригодным для жёстких систем.

¶Усовершенствованный метод Эйлера (метод Хойна)

Для повышения точности используется двухэтапная схема:

\[ \begin{aligned} k_1 &= f(x_n, y_n) \\ k_2 &= f(x_n + h, y_n + h k_1) \\ y_{n+1} &= y_n + \frac{h}{2}(k_1 + k_2) \end{aligned} \]

Этот метод имеет второй порядок точности и является простейшим представителем методов Рунге — Кутты второго порядка.

¶Метод Эйлера с пересчётом

Вариант с коррекцией: сначала вычисляется предварительное значение \( \tilde{y}_{n+1} \) по явной схеме, затем уточняется с использованием среднего значения производной:

\[ y_{n+1} = y_n + \frac{h}{2} \left[ f(x_n, y_n) + f(x_{n+1}, \tilde{y}_{n+1}) \right] \]

¶Применение

Метод Эйлера используется в следующих областях:

• Образование: как введение в численные методы решения дифференциальных уравнений. Простота реализации позволяет студентам быстро понять принципы дискретизации.
• Предварительные расчёты: для грубой оценки поведения решения, когда высокая точность не требуется.
• Моделирование простых динамических систем: например, радиоактивный распад, охлаждение тел, простые модели популяционной динамики.
• Компьютерная графика и анимация: для численного интегрирования уравнений движения частиц в симуляциях физики (метод Эйлера часто используется в игровых движках из-за его скорости, хотя он менее точен, чем метод Верле или Рунге — Кутты).

¶Ограничения

• Низкая точность: для достижения приемлемой погрешности требуется малый шаг, что приводит к большому объёму вычислений.
• Неустойчивость для жёстких задач: явная схема может давать расходящиеся решения даже при малом шаге.
• Накопление ошибок: при большом числе шагов погрешность может стать неприемлемо большой.

¶Пример

Рассмотрим задачу Коши:

\[ y' = -2y, \quad y(0) = 1 \]

Точное решение: \( y(x) = e^{-2x} \). Применим метод Эйлера с шагом \( h = 0.1 \) на отрезке \( [0, 1] \).

Видно, что погрешность растёт с увеличением расстояния от начальной точки, но остаётся в пределах \( O(h) \). Уменьшение шага до \( h = 0.01 \) даёт погрешность около 0.003.

¶Историческая справка

Метод назван в честь швейцарского математика Леонарда Эйлера, который описал его в 1768 году в своей работе «Интегральное исчисление» (лат. «Institutiones calculi integralis»). Эйлер разработал метод для решения задач небесной механики, где аналитическое решение дифференциальных уравнений часто было невозможно. Несмотря на свою простоту, метод остаётся важным инструментом в вычислительной математике и служит основой для более сложных алгоритмов.

¶Связь с другими методами

Метод Эйлера является частным случаем:

• Методов Рунге — Кутты: явный метод Эйлера соответствует методу Рунге — Кутты первого порядка с одним этапом.
• Метода конечных разностей: применяется для дискретизации производной в дифференциальных уравнениях.
• Метода ломаных Эйлера: геометрическая интерпретация, где интегральная кривая заменяется ломаной линией.

¶Литература

• Эйлер Л. Интегральное исчисление. — М.: ГИТТЛ, 1956. — Т. 1.
• Бахвалов Н. С., Жидков Н. П., Кобельков Г. М. Численные методы. — М.: Бином, 2008.
• Самарский А. А., Гулин А. В. Численные методы. — М.: Наука, 1989.
• Хайрер Э., Нёрсетт С., Ваннер Г. Решение обыкновенных дифференциальных уравнений. Нежёсткие задачи. — М.: Мир, 1990.