Явный и неявный метод Эйлера: устойчивость и выбор шага

Явный и неявный метод Эйлера - две простейшие одношаговые схемы численного решения задачи Коши для обыкновенного дифференциального уравнения. Обе используют один и тот же наклон $f(x, y)$ и дают глобальную погрешность $O(h)$, но различаются тем, в какой точке шага этот наклон берётся: явный (прямой) метод Эйлера опирается на наклон в начале отрезка, неявный (обратный) - на наклон в его конце. Это, казалось бы, мелкое различие радикально меняет устойчивость схемы, и именно оно объясняет, почему на жёстких уравнениях явный метод «взрывается», а неявный спокойно считает с крупным шагом.

Постановка задачи Коши и сетка

Решаем задачу Коши: дано уравнение первого порядка $y' = f(x, y)$ и начальное условие $y(x_0) = y_0$. Когда аналитической первообразной нет, решение строят по точкам. Вводят равномерную сетку с шагом $h$: узлы $x_n = x_0 + n h$, и ищут приближения $y_n \approx y(x_n)$. Любая одношаговая схема даёт правило перехода от $y_n$ к $y_{n+1}$; вопрос лишь в том, какой наклон и в какой точке использовать для приращения за шаг.

Метод Эйлера - самый прямой ответ: заменим производную $y'$ в точке разностным отношением и продлим наклон на весь шаг. В зависимости от того, в каком узле вычисляется правая часть, получаются две разные схемы. Та же развилка «вперёд или назад» лежит в основе разностных схем для краевых задач.

Tool: один шаг явного и неявного метода Эйлера

Считать шаги вручную и для прямой, и для обратной схемы - особенно когда неявная формула $y_{n+1} = y_n + h f(x_{n+1}, y_{n+1})$ требует решать уравнение относительно $y_{n+1}$ - занятие кропотливое. Ниже мини-форма: выбираешь уравнение, начальную точку, шаг $h$, число шагов и схему (явная / неявная / обе для сравнения) - получаешь пошаговый расчёт и таблицу значений.

Явный (прямой) метод Эйлера

Явный метод Эйлера, он же прямой или forward Euler, берёт наклон в начале отрезка $[x_n, x_{n+1}]$ и продлевает его на весь шаг:

$y_{n+1} = y_n + h\, f(x_n, y_n).$

Формула явная в буквальном смысле: правая часть зависит только от уже известных $x_n$ и $y_n$, поэтому $y_{n+1}$ вычисляется прямой подстановкой, без решения уравнений. Геометрически мы движемся по касательной к решению, проведённой в текущей точке. Это самый дешёвый из всех методов - одно вычисление $f$ на шаг, - но и самый грубый: касательная отрывается от изгибающейся кривой решения, и ошибка накапливается.

Явный метод Эйлера естественно выводится из разложения Тейлора: $y(x_{n+1}) = y(x_n) + h\, y'(x_n) + \tfrac{h^2}{2} y''(\xi)$. Отбрасывая член с $h^2$ и заменяя $y'(x_n)$ на $f(x_n, y_n)$, получаем схему. Отброшенный член $\tfrac{h^2}{2} y''$ и есть источник локальной погрешности.

Неявный (обратный) метод Эйлера

Неявный метод Эйлера (обратный, backward Euler) берёт наклон в конце отрезка:

$y_{n+1} = y_n + h\, f(x_{n+1}, y_{n+1}).$

Здесь неизвестное $y_{n+1}$ входит и в левую, и в правую часть - это алгебраическое уравнение относительно $y_{n+1}$, которое в общем случае нелинейно. На каждом шаге его решают: для линейной по $y$ правой части - явной формулой, для нелинейной - итерациями Ньютона или простой итерацией. Цена шага выше, чем у явной схемы, но устойчивость несравнимо лучше.

Для линейного тестового уравнения $y' = \lambda y$ обратная схема даёт $y_{n+1} = y_n + h\lambda y_{n+1}$, откуда $y_{n+1} = y_n / (1 - h\lambda)$ - решается сразу, без итераций. Именно на этом уравнении удобно сравнивать устойчивость двух методов.

Локальная и глобальная погрешность O(h)

Обе схемы имеют одинаковый порядок точности. Локальная погрешность усечения - ошибка за один шаг при точном предыдущем значении - составляет $O(h^2)$: это и есть отброшенный член $\tfrac{h^2}{2} y''$. Шагов на отрезке порядка $1/h$, поэтому глобальная погрешность $|y_n - y(x_n)|$ на порядок ниже локальной и равна $O(h)$. Отсюда название «методы первого порядка».

Практический вывод из $O(h)$ суров: чтобы вдвое уменьшить глобальную ошибку, нужно вдвое уменьшить шаг - то есть вдвое больше вычислений. Для сравнения, у метода Рунге-Кутты 4 порядка ошибка падает как $O(h^4)$, и уменьшение шага вдвое снижает её в шестнадцать раз. Поэтому метод Эйлера в чистом виде применяют редко - но он остаётся фундаментом, на котором строят и понимают все более точные схемы.

Абсолютная устойчивость и тестовое уравнение

Главное различие между явным и неявным методом Эйлера - не точность, а устойчивость. Её исследуют на модельном уравнении $y' = \lambda y$ с $\mathrm{Re}\,\lambda < 0$, точное решение которого затухает: $y = y_0 e^{\lambda x} \to 0$. Хорошая схема должна тоже давать затухающую последовательность $y_n \to 0$.

Для явного метода $y_{n+1} = (1 + h\lambda) y_n$, то есть $y_n = (1 + h\lambda)^n y_0$. Чтобы $y_n \to 0$, нужен множитель $|1 + h\lambda| < 1$ - это условие задаёт ограниченную область абсолютной устойчивости (для вещественного $\lambda < 0$ это $h < 2/|\lambda|$). Превысишь шаг - и приближённое решение начнёт расти по модулю, хотя точное затухает: схема «взрывается».

Для неявного метода $y_{n+1} = y_n / (1 - h\lambda)$, то есть $y_n = (1 - h\lambda)^{-n} y_0$. При $\mathrm{Re}\,\lambda < 0$ знаменатель по модулю больше единицы при любом $h > 0$, поэтому $|y_n| \to 0$ всегда. Неявный метод Эйлера абсолютно устойчив (A-устойчив): его область устойчивости включает всю левую полуплоскость, и ограничения на шаг по соображениям устойчивости нет.

Жёсткие уравнения: где явный метод проигрывает

Жёсткое уравнение - то, где сосуществуют процессы с сильно разными масштабами времени (большие по модулю отрицательные $\lambda$), например быстро затухающие компоненты на фоне медленной динамики. Явный метод Эйлера вынужден держать шаг $h < 2/|\lambda_{\max}|$ из-за самой быстрой компоненты - даже после того, как она давно затухла и на решение уже не влияет. Это делает явную схему катастрофически медленной: тысячи мелких шагов там, где по точности хватило бы десятка крупных.

Неявный метод Эйлера от этого свободен: будучи A-устойчивым, он выбирает шаг исходя только из требуемой точности, а не из устойчивости. Дополнительная стоимость решения уравнения на каждом шаге окупается многократно меньшим числом шагов. Поэтому для жёстких задач - кинетика реакций, электрические цепи, теплопроводность - выбирают именно неявные схемы. Та же логика «явная схема дешева, но условно устойчива; неявная дороже на шаг, но устойчива безусловно» дословно повторяется в методе конечных разностей для уравнений в частных производных.

Частые ошибки

• Считают явный и неявный методы разными по точности. Оба первого порядка, глобальная погрешность у обоих $O(h)$. Различие - в устойчивости, а не в порядке.
• Решают неявную формулу как явную. В $y_{n+1} = y_n + h f(x_{n+1}, y_{n+1})$ нельзя просто подставить $y_n$ в правую часть - это превратит схему в явную. $y_{n+1}$ надо найти из уравнения.
• Путают локальную и глобальную погрешность. Локальная - $O(h^2)$, глобальная - $O(h)$. Имя «первый порядок» относится к глобальной.
• Берут крупный шаг для явной схемы на жёсткой задаче. При $h > 2/|\lambda|$ решение разлетается. Это не ошибка округления, а потеря устойчивости - лечится уменьшением шага или переходом к неявному методу.
• Думают, что A-устойчивость = высокая точность. Неявный Эйлер устойчив, но всё ещё первого порядка: при крупном шаге решение не «взрывается», но и точным не будет.

FAQ

В чём главное различие явного и неявного метода Эйлера? В точке вычисления наклона: явный берёт $f$ в начале шага (узел $n$), неявный - в конце (узел $n+1$). Из-за этого явный требует решения уравнения на $y_{n+1}$, зато неявный абсолютно устойчив.

Почему неявный метод предпочитают для жёстких уравнений? Он A-устойчив: шаг ограничивается только точностью, а не устойчивостью. Явный метод на жёсткой задаче вынужден держать крошечный шаг из-за быстро затухающей компоненты, что делает его непомерно медленным.

Какой порядок точности у метода Эйлера? Первый: глобальная погрешность $O(h)$, локальная $O(h^2)$ - одинаково у явной и неявной схемы. Для высокой точности переходят к методам Рунге-Кутты.

Коротко

Явный и неявный метод Эйлера решают задачу Коши $y' = f(x, y)$, $y(x_0) = y_0$, продлевая наклон на шаг $h$: явная схема $y_{n+1} = y_n + h f(x_n, y_n)$ берёт наклон в начале отрезка и считается прямой подстановкой, неявная $y_{n+1} = y_n + h f(x_{n+1}, y_{n+1})$ - в конце, требуя решения уравнения на $y_{n+1}$. Обе первого порядка с глобальной погрешностью $O(h)$, но различаются устойчивостью: явный метод условно устойчив ($h < 2/|\lambda|$ на тесте $y' = \lambda y$), неявный - абсолютно (A-устойчив). Поэтому для гладких нежёстких задач хватает дешёвой явной схемы, а на жёстких уравнениях выбирают неявную, окупающую стоимость шага свободой в выборе $h$.