Метод Рунге-Кутты 4 порядка: численное решение ОДУ

Метод Рунге-Кутты 4 порядка (часто пишут РК4 или RK4) - самый популярный одношаговый метод численного решения задачи Коши для обыкновенного дифференциального уравнения. На каждом шаге он вычисляет четыре пробных наклона $k_1, k_2, k_3, k_4$ и берёт их взвешенное среднее, добиваясь глобальной погрешности $O(h^4)$ при четырёх вычислениях правой части. Этот баланс точности и стоимости и сделал метод Рунге-Кутты 4 порядка стандартом «по умолчанию» в инженерных расчётах и учебных задачах.

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

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

Самая грубая оценка - метод Эйлера: $y_{n+1} = y_n + h\, f(x_n, y_n)$. Он берёт наклон только в начале отрезка и продлевает его на весь шаг, накапливая глобальную погрешность $O(h)$. Метод Рунге-Кутты 4 порядка устраняет эту грубость, опрашивая наклон в четырёх точках отрезка.

Tool: один шаг метода Рунге-Кутты 4 порядка

Считать $k_1, k_2, k_3, k_4$ вручную для $y' = x^2 + y^2$ при $h = 0{,}1$ - занятие, в котором легко ошибиться арифметически: на каждом шаге четыре подстановки в правую часть и аккуратное взвешивание. Ниже - мини-форма: выбираешь уравнение из списка (или вписываешь своё), задаёшь $x_0$, $y_0$, шаг $h$ и число шагов - получаешь пошаговый расчёт всех четырёх коэффициентов и таблицу значений $y_n$.

Формулы коэффициентов k1–k4

Суть метода Рунге-Кутты 4 порядка - оценить наклон решения в четырёх точках отрезка $[x_n, x_{n+1}]$ и усреднить их по правилу Симпсона. Коэффициенты вычисляются последовательно, каждый следующий использует предыдущий:

$k_1 = f(x_n,\, y_n),$ $k_2 = f\!\left(x_n + \tfrac{h}{2},\; y_n + \tfrac{h}{2} k_1\right),$ $k_3 = f\!\left(x_n + \tfrac{h}{2},\; y_n + \tfrac{h}{2} k_2\right),$ $k_4 = f\!\left(x_n + h,\; y_n + h\, k_3\right).$

Здесь $k_1$ - наклон в начале шага (как у Эйлера), $k_2$ и $k_3$ - два уточнённых наклона в середине отрезка, $k_4$ - наклон в конце, посчитанный уже с учётом продвижения по $k_3$. Итоговое приращение - взвешенное среднее:

$y_{n+1} = y_n + \frac{h}{6}\bigl(k_1 + 2 k_2 + 2 k_3 + k_4\bigr).$

Веса $1, 2, 2, 1$ при сумме $6$ в знаменателе - не случайны: это в точности квадратурные веса метода Симпсона на отрезке с серединой. Если правая часть не зависит от $y$, то есть $y' = f(x)$, метод Рунге-Кутты 4 порядка вырождается ровно в формулу Симпсона для интеграла $\int_{x_n}^{x_{n+1}} f(x)\,dx$. Об этой квадратуре подробнее - в статье про метод Симпсона.

Геометрический смысл усреднения

Метод Эйлера ошибается, потому что наклон в начале отрезка не отражает, как кривая решения изгибается дальше. Пробные точки $k_2$ и $k_3$ заглядывают в середину: $k_2$ оценивает наклон, шагнув на полшага по эйлеровому $k_1$, а $k_3$ - повторно, но уже по уточнённому $k_2$. Эти два средних наклона ловят кривизну решения. Финальный $k_4$ проверяет наклон в конце шага. Взвешенное среднее $\tfrac{1}{6}(k_1 + 2k_2 + 2k_3 + k_4)$ гасит систематическую ошибку: ряды Тейлора приближения и точного решения совпадают вплоть до члена $h^4$ включительно.

Именно это совпадение тейлоровских разложений до четвёртого члена и означает «4 порядок»: каждый дополнительный согласованный член ряда поднимает порядок на единицу. Эйлер согласует только линейный член - отсюда его первый порядок.

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

Различают две погрешности. Локальная погрешность усечения - ошибка за один шаг при условии, что предыдущее значение точное; у метода Рунге-Кутты 4 порядка она составляет $O(h^5)$. Глобальная погрешность - накопленное отклонение $|y_n - y(x_n)|$ на всём отрезке; шагов порядка $1/h$, поэтому глобальная погрешность на порядок ниже локальной и равна $O(h^4)$. Отсюда и название: «4 порядка» - это порядок именно глобальной точности.

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

Сравнение с методом Эйлера и порядок методов

Семейство явных методов Рунге-Кутты параметризуется числом стадий $s$ (вычислений $f$ за шаг). Метод Эйлера - это РК первого порядка с одной стадией. Метод средней точки и метод Хойна (улучшенный Эйлер) - второго порядка, две стадии. Классический РК4 - четыре стадии, четвёртый порядок. До $s = 4$ порядок равен числу стадий, что и делает РК4 «золотой серединой»: дальше за каждую дополнительную единицу порядка приходится платить более чем одной стадией (барьер Бутчера), и выигрыш окупается только при очень высоких требованиях к точности.

Системы уравнений и уравнения высших порядков

Метод Рунге-Кутты 4 порядка переносится на системы без изменений: $y$ становится вектором, $f$ - вектор-функцией, а $k_1, \ldots, k_4$ - векторами; формулы те же. Уравнение второго порядка $y'' = g(x, y, y')$ сводят к системе из двух уравнений первого порядка заменой $u = y'$: получают $y' = u$, $u' = g(x, y, u)$ и решают совместно. Так РК4 покрывает почти все учебные дифференциальные уравнения, не имеющие аналитического решения. Если же интегральное преобразование всё-таки даёт замкнутый ответ, его стоит предпочесть - см. преобразование Лапласа.

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

• Путают $k_2$ и $k_3$. Оба считаются в середине $x_n + h/2$, но $k_2$ использует сдвиг $y_n + \tfrac{h}{2}k_1$, а $k_3$ - сдвиг $y_n + \tfrac{h}{2}k_2$. Подстановка $k_1$ вместо $k_2$ в третий коэффициент - самая частая арифметическая ошибка.
• Забывают полушаг в аргументе $x$. В $k_2$ и $k_3$ сдвигается не только $y$, но и $x$ на $h/2$; в $k_4$ - на полный $h$. Если $f$ зависит от $x$, пропуск этого сдвига портит порядок.
• Делят на 4 вместо 6. Веса $1, 2, 2, 1$ дают сумму $6$, поэтому в знаменателе стоит именно $6$, а не число коэффициентов.
• Считают $O(h^4)$ локальной погрешностью. Локальная - $O(h^5)$, глобальная - $O(h^4)$. Их путаница искажает оценку при выборе шага.
• Берут слишком большой $h$ для жёсткой задачи. Явный РК4 условно устойчив: на жёстких уравнениях при большом $h$ решение «разлетается». Это решается уменьшением шага или неявными методами.

FAQ

Почему именно 4 порядок, а не выше? До четырёх стадий порядок равен числу вычислений $f$, дальше за каждый дополнительный порядок нужно более одной лишней стадии. РК4 - точка наилучшего соотношения точности и стоимости для большинства гладких задач.

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

Подходит ли РК4 для систем и уравнений второго порядка? Да. Система решается покомпонентно теми же формулами с векторными $k_i$, а уравнение высшего порядка предварительно сводят к системе первого порядка заменой переменных.

Коротко

Метод Рунге-Кутты 4 порядка решает задачу Коши $y' = f(x, y)$, $y(x_0) = y_0$, вычисляя на каждом шаге четыре наклона $k_1$–$k_4$ и переходя по правилу $y_{n+1} = y_n + \tfrac{h}{6}(k_1 + 2k_2 + 2k_3 + k_4)$. Веса $1,2,2,1$ повторяют квадратуру Симпсона, разложения Тейлора совпадают до члена $h^4$, отсюда глобальная погрешность $O(h^4)$ - на три порядка лучше Эйлера при четырёх вычислениях правой части за шаг. Метод без изменений переносится на системы и, через сведение к системе, на уравнения высших порядков, что и делает его рабочим стандартом численного интегрирования ОДУ.