Методы Рунге-Кутты

Методы Рунге-Кутты – это группа методов, широко применяемых на практике для решения ОДУ. В этих методах при вычислении значения искомой функции в очередной точке х_i+1 используется информация о предыдущей точке х_i, y_i. Методы различаются объемом вычислений и точностью результата.

Порядок метода Рунге-Кутты определяется кратностью вычисления значения производной искомой функции f(x,y ) на каждом шаге. В соответствии с этим метод Эйлера является методом Рунге-Кутты первого порядка, поскольку для получения очередного значения y_i+1 функция f(x) вычисляется один раз в предыдущей точке х_i, y_i. В методах Рунге-Кутты более высоких порядков для вычисления очередного значения искомой функции в точке х_i+1 значение правой части уравнения y’= f(x,y ) вычисляется несколько раз, количество которых и определяет порядок метода.

Метод Рунге-Кутты 2-го порядка (Усовершенствованный метод Эйлера ). Вычисление значения искомой функции в точке х_i+1 проводится в два этапа. Сначала вычисляют вспомогательную величину по методу Эйлера:

(6.5.3-1)

Затем значение производной искомой функции в точке (x_i+1,y_i+1) используется для вычисления окончательного значения функции:

(6.5.3-2)

Подставляя (6.5.3-1) в (6.5.3-2), окончательно получим расчетную формулу метода Рунге-Кутты 2 -го порядка:

(6.5.3-3)

Этот метод также называют методом прогноза и коррекций. Сначала находят грубое приближение по методу Эйлера (прогноз), а затем уточненное значение y_i+1 (коррекция).

В общем виде формулу (6.5.3-3) можно представить как

(6.5.3-4)

Метод Рунге-Кутты второго порядка имеет наглядную геометрическую интерпретацию (рис. 6.5.3-1). Построение проводится следующим образом: определяется пересечением перпендикуляра, восстановленного из точки x_i+1 c касательной L₁, проведенной к кривой y(x) в предыдущей точке (х_i,y_i). Затем в точке проводится прямая L₂ с тангенсом угла наклона, равным . Прямую проводят через точку под углом, тангенс которого находим усреднением значений тангенсов углов наклона L₁ и L₂. Прямая L проводится параллельно через точку (х_i,y_i). Ее пересечение с перпендикуляром, восстановленным из точки х_i+1, и дает уточненное значение y_i+1.

Рис. 6.5.3-1

Погрешность метода Рунге-Кутты второго порядка связана с величиной шага интегрирования отношением e₂ =C₂h³, где C₂– произвольная постоянная.

Пример 6.5.3-1. Решить методом Рунге-Кутты второго порядка ОДУ y¢= 2x/y с начальными условиями x₀ = 1 и y₀ = 1 на отрезке [1;1.4] и шагом h = 0.2.

Проводя дальнейшее обобщение формул Рунге-Кутты, для решения ОДУ первого порядка можно записать следующее:

где Ф – линейная функция аргументов x, y, h и f(x,y), которая может быть представлена как

(6.5.3-5)

Величина n в (6.5.3-4) определяется порядком метода, а коэффициентам a₂,a₃, …,a_n, Р₁, Р₂, …,P_n подбирают такие значения, которые обеспечивают минимальную погрешность. Так, для метода Рунге-Кутты четвертого порядка (n=4) получена расчетная формула при следующих коэффициентах: a₂= a₃=1/2, a₄=1, P₁ = P₄=1/6, P₂ = P₃ =2/6.

Подставив значения коэффициентов в (6.5.3-4), имеем

(6.5.3-6)

Геометрическая интерпретация этого метода очень сложна и потому не приводится.

Погрешность метода Рунге-Кутты четвертого порядка значительно меньше методов первого и второго порядков и пропорциональна величине h (e₄ =C₄h⁵).

Пример 6.5.3-2. Решить методом Рунге-Кутты четвертого порядка ОДУ y¢= 2x/y с начальными условиями x₀ = 1 и y₀ = 1 на отрезке [1;1.4] с шагом h = 0,2.

Сведем в таблицу результаты решения уравнения y¢=2x/y методами Рунге-Кутты, соответственно, первого (y¹_i), второго (y²_i) и четвертого (y⁴_i) порядков и сравним с результатами, полученными точным методом (y_i).

хi	y^1i	y^2i	y^4i	yi
1.2 1.4	1.4 1.74286	1.3714 1.7091	1.37115 1.7089	1.37113 1.7088

На практике для обеспечения требуемой точности (при использовании любого приближенного метода решения ОДУ) применяется автоматический выбор шага методом двойного просчета. При этом в каждой точке х_i по формуле, соответствующей выбранному методу, производится расчет y_i с шагом h (y_i^(h)) и с шагом h/2 (y_i^(h/2)). Цель двойного просчета состоит в том, чтобы для каждой точки численного решения эти значения отличались на величину, не превышающую заданной погрешности e. В этом случае общая формула для оценки погрешности решения ОДУ методами Рунге-Кутты имеет следующий вид:

где p – порядок метода Рунге-Кутты. Эта формула называется также правилом Рунге.

Если | y_i^(h)) - y_i^(h/2)|< e, то шаг для следующей точки выбирается равным h, иначе шаг уменьшается вдвое и продолжается уточнение y _i в точке х_i.

Дата добавления: 2014-12-27