Калькулятор метода Эйлера

О Калькулятор метода Эйлера

Этот Калькулятор Метода Эйлера численно решает любую задачу Коши для ОДУ первого порядка вида \( y' = f(x, y), \; y(x_0) = y_0 \), используя классический (прямой) метод Эйлера. Он выводит полную итерационную таблицу, строит полигон Эйлера поверх поля направлений, сравнивает решения при трех различных размерах шага для визуального контроля сходимости и — если указано точное решение — проводит анализ ошибок на каждом шаге.

Что такое метод Эйлера?

Метод Эйлера — это простейший алгоритм для аппроксимации решения задачи Коши. Начиная с известной точки \( (x_0, y_0) \) на кривой решения, он последовательно продвигается на малый шаг размером h вдоль локального наклона \( f(x, y) \):

Геометрически каждый шаг представляет собой короткий прямолинейный отрезок, наклон которого равен значению дифференциального уравнения в текущей точке. Результирующая ломаная линия — полигон Эйлера — является приближением к истинному (обычно криволинейному) решению.

Насколько он точен?

Метод Эйлера является методом первого порядка. Локальная ошибка аппроксимации на каждом шаге составляет \( O(h^2) \), а глобальная ошибка после интегрирования на фиксированном интервале — \( O(h) \). На практике:

• Уменьшение шага вдвое примерно вдвое уменьшает глобальную ошибку.
• Ошибка растет линейно с увеличением интервала интегрирования.
• Наибольшая ошибка возникает там, где решение имеет высокую кривизну.

Встроенное сравнение размеров шага (h, h/2, h/4) позволяет увидеть эту линейную сходимость напрямую: включите эту опцию и убедитесь, что три конечных значения приближаются к общему пределу, причем каждое последующее значение находится примерно в два раза ближе к пределу, чем предыдущее.

Как читать график

Визуализация накладывает четыре типа данных на одну координатную плоскость:

• Серое поле направлений — короткие отрезки, наклон которых равен \( f(x, y) \) в данной точке. Это "поток", диктуемый ОДУ. Любая кривая решения должна быть касательной к полю в каждой точке.
• Индиго полигон Эйлера — пошаговое численное решение. Каждый сегмент начинается в предыдущей точке сетки и направлен вдоль \( f(x_n, y_n) \) на расстояние h.
• Пунктирная зеленая точная кривая — отображается только при вводе точного решения. Вертикальные оранжевые пунктирные линии показывают локальные ошибки \( y_n - y_{\text{точно}}(x_n) \).
• Оранжевая и зеленая кривые сравнения — решение той же задачи при шагах h/2 и h/4 (при включенном сравнении шагов).

Как пользоваться этим калькулятором

1. Введите правую часть ОДУ в поле y' =. Используйте x и y в качестве переменных. Поддерживаются операторы + − × ÷ ^ и функции sin, cos, tan, asin, acos, atan, exp, ln, log, log10, log2, sqrt, abs.
2. Укажите начальные условия: начальное значение x₀, значение y₀ в этой точке, размер шага h (положительный для интегрирования вперед, отрицательный — назад) и количество шагов n.
3. (Опционально) Введите точное решение y(x), если оно вам известно. Калькулятор вычислит \( |y_n - y(x_n)| \) на каждом шаге и выведет максимальную и конечную ошибки.
4. Настройте визуализацию: поле направлений включено по умолчанию; сравнение шагов накладывает две дополнительные кривые при h/2 и h/4.
5. Нажмите "Запустить". В результатах появится статистика, график, панель сравнения сходимости и полная итерационная таблица. При наведении на строку подсвечивается соответствующая точка на графике (и наоборот).

Пример решения

Рассмотрим \( y' = x + y, \; y(0) = 1 \) при h = 0.1 и 10 шагах. Точное решение: \( y(x) = -x - 1 + 2e^x \). Применяя формулу Эйлера:

Конечная ошибка составляет около 0.249. Уменьшение h до 0.05 снижает ошибку примерно до 0.13, а при h = 0.025 она падает до 0.067 — чистая линейная сходимость, в полном соответствии с теорией.

Метод Эйлера vs другие численные методы

Когда метод Эйлера ошибается

Прямой метод Эйлера может работать некорректно в трех случаях:

• Слишком большой шаг — ломаная осциллирует или расходится. Решение — уменьшить h; сравнение h, h/2, h/4 делает это наглядным.
• Жесткие ОДУ — уравнения с одновременно быстро и медленно затухающими составляющими заставляют выбирать ничтожно малый h для устойчивости. Переходите на неявные методы или BDF.
• Особенности в f(x, y) — деление на ноль, sqrt от отрицательного числа или ln от неположительного остановят расчет. Калькулятор четко укажет на проблемный шаг.

Типовые применения

• Физика — второй закон Ньютона как система первого порядка, радиоактивный распад \( \dot{N} = -\lambda N \), закон охлаждения Ньютона.
• Биология и эпидемиология — логистический рост \( \dot{y} = r\,y(1 - y/K) \), компартментные SIR-модели.
• Экономика — непрерывное начисление процентов, простые модели роста Солоу.
• Химия — кинетика реакций первого порядка \( \dot{c} = -k c \).
• Обучение — введение в концепцию численного интегрирования перед изучением RK4 или адаптивных решателей.

Часто задаваемые вопросы

Что такое метод Эйлера?

Метод Эйлера — это простейшая численная процедура решения задачи Коши y' = f(x, y), y(x0) = y0. На каждом шаге решение продвигается по формуле y_{n+1} = y_n + h · f(x_n, y_n), следуя наклону в текущей точке. Он имеет первый порядок точности, то есть глобальная ошибка пропорциональна h.

Насколько точен метод Эйлера?

Он имеет локальную ошибку O(h²) и глобальную O(h). Уменьшение шага вдвое примерно вдвое снижает ошибку. Это наглядно показано в панели сравнения сходимости нашего калькулятора.

Когда метод Эйлера не работает?

Он становится неустойчивым для жестких задач или при слишком больших шагах. Решение может начать колебаться или уйти в бесконечность. В таких случаях нужно уменьшать шаг или использовать неявные методы.

Как выбрать размер шага?

Начните с такого h, чтобы на интервале получилось 10–50 шагов. Если ломаная сильно отклоняется от поля направлений, уменьшите h. Используйте сравнение h/2 и h/4, чтобы убедиться в сходимости кривых.

Чем метод Эйлера отличается от Рунге-Кутты (RK4)?

RK4 оценивает наклон в четырех точках на шаг, обеспечивая точность O(h⁴). Это на порядки точнее Эйлера при том же числе шагов, но метод Эйлера проще для понимания основ численного анализа.

Можно ли решать системы ОДУ?

Этот калькулятор предназначен для одного скалярного уравнения. Системы или уравнения высших порядков требуют специализированных решателей или покомпонентного сведения к уравнениям первого порядка.

Можно ли интегрировать в обратном направлении?

Да — введите отрицательное значение шага h. Калькулятор выполнит расчет от x₀ в сторону уменьшения x, что полезно для реконструкции прошлых состояний системы.

Дополнительные материалы

• Метод Эйлера — Википедия
• Методы Рунге — Кутты — Википедия
• Поле направлений — Википедия
• Жесткие ОДУ — Википедия