Методы Рунге-Кутта 2-го порядка

Для уменьшения погрешности методов, использующих тейлоровское разложение искомого решения, необходимо использовать большее количество членов ряда. Однако при этом возникает необходимость аппроксимации производных от правых частей ОДУ. Основная идея методов Рунге-Кутта заключается в том, что производные аппроксимируются через значения функции F(t_k, f(t_k)) в точках на интервале [t_k, t_k+h], которые выбираются из условия наибольшей близости алгоритма к ряду Тейлора.

В зависимости от старшей степени h, с которой учитываются члены ряда, построены схемы Рунге-Кутта разных порядков точности. Так, например, для 2-го порядка получено однопараметрическое семейство схем вида

f(t_k+h)=f(t_k+h[(1-L)F_k+LF(t_k+гh, t_k+гF_kh)]+O(h³), (*)

где 0<L≤1 – свободный параметр, F_k=F(t_k, f(t_k), г=1/2L, O(h³) – остаточная погрешность.

Для параметра L наиболее часто используют значения L=1/2 и L=1. В первом случае формула (*) приобретает вид

f(t_k+h)=f(t_k)+h[F_k+F(t_k+h, f(t_k))+hF_k)]/2

Вначале вычисляется приближенное решение ОДУ в точке t_k+h по формуле Эйлера

f_э=f_k+hF_k.

Затем определяется наклон интегральной кривой в найденной точке F(t_k+h, f_э) и после нахождения среднего наклона на шаге h находится уточненное значение f_RK=f(t_k+h).

Схемы подобного типа называют «прогноз – коррекция». С целью экономии памяти при программировании алгоритма, обобщенного на системы ОДУ, изменим его запись с учетом того, что f_k=f_э–hF_k:

f_i(t_k+h)=f_iэ+(1/2)h[F_i(t_k+h, f_iэ)],

где i – номер решения для системы ОДУ.

Теперь не придется держать в памяти массив начальных значений f_i0 – его можно забыть после вычисления значений эйлеровских приближений f_iэ, размещаемых на месте массива f_i0.

Во втором случае при L=1 от формулы (*) переходим к схеме

f(t_k+h)=f_k+hF(t_k+h/2, f_k+hF_k/2).

Здесь при прогнозе определяется методом Эйлера решение в точке (t_k+h/2):

f_1/2=f_k+hF_k/2,

а после вычисления наклона касательной к интегральной кривой в средней точке решение корректируется по этому наклону.

Таким образом, для методов Рунге-Кутта 2-го порядка функции j(t, f(t)) имеют вид:

j(t, f(t))=(1/2)[F(t, f(t))+F(t+h), f(t)+hF(t, f(t))]

или

j(t, f(t))=F(t_k+h/2, f_k+hF_k/2).