Метод Рунге-Кутты четвертого порядка с автоматическим выбором шага интегрирования решения задачи Коши (стр. 3 из 7)

Правило 3/8

Классический метод Рунге-Кутты

1.5 «Оптимальные» формулы

Предпринималось много исследования, чтобы выбрать возможно «лучшие» из множества различных формул Рунге-Кутты 4-го порядка.

Первой попыткой в этом направлении был очень популярный метод, который в 1951 году предложил Гилл. Он преследовал цель уменьшить на сколько возможно количество требуемой машинной памяти («регистров»). Этот метод широко использовался на первых компьютерах в пятидесятых годах и представляет поэтому исторический интерес. Гилл заметил, что больше всего машинной памяти нужно при вычислении

, когда «требуются регистры для хранения в какой-либо форме» величин

.

Ясно, что для третьей стадии будет достаточно трех регистров, если подлежащие хранению величины линейно зависимы, то есть если

.

Гилл заметил, что это условие удовлетворяется для методов типа 3), если

. Получающийся метод можно в таком случае переформулировать следующим образом:

В настоящее время производительность компьютеров очень сильно возросла по сравнению с машинами 50-х годов, что привело к прекращению использования данного метода на практике при расчетах. Справедливости ради стоит отметить, что этот метод все-таки рационально употреблять в случаях, когда требуется решать системы дифференциальных уравнений очень высокой размерности со сложными функциями.

1.6 Условия порядков для методов Рунге-Кутты

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

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

Метод

(где на свободных местах должны стоять нули) имеет порядок

, если удовлетворяется уравнение

(2.6.1)

для каждого дерева

с корнем и не более чем с разветвлениями[1].

При

эти условия, обеспечивающие порядок 4, и соответствующие деревья имеют следующий вид:

	(2.6.2)
	(2.6.3)
	(2.6.4)
	(2.6.5)
	(2.6.6)
	(2.6.7)
	(2.6.8)
	(2.6.9)

Заметим, что для меньших значений

мы берем соответствующее подмножество этих условий, а для меньших оставляем лишь некоторые из указанных членов.

Из (2.9) видим, что действительно необходимо 4 этапа, так как если бы их было меньше, то был бы опущен единственный член в левой части этого уравнения. Для явных методов в общем случае выполняется неравенство

. Фактически (для тех значений, для которых это известно) минимальное значение для данного указано в следующей таблице:

	1	2	3	4	5	6	7
	1	2	3	4	6	7	9

Общие классы методов с этими значениями

и легко найти в случае .

Для

:

0

1

Это известный метод Эйлера.

Для

:

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

.

Для

имеется три семейства, из которых первые два таковы:

Каждое из них имеет один параметр

. Третье семейство имеет в качестве параметров и , причем

.

Вывод методов с

более сложен, но его можно упростить, положив

(2.6.10)

(что влечет равенство

), так как это позволяет опустить уравнения (2.6.3), (2.6.5), (2.6.8) и (2.6.9). Интересно также, что (2.6.10) является следствием (2.6.2) – (2.6.9).