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

Самое большое распространение из всех численных методов решения дифференциальных уравнений с помощью ЭВМ получил метод Рунге-Кутта 4-го порядка.

В этом методе на каждом шаге интегрирования дифференциальных уравнений искомая функция y(x) аппроксимируется рядом Тейлора (4), содержащим члены ряда с h⁴:

В результате ошибка на каждом шаге имеет порядок h⁵.

Для сохранения членов ряда, содержащих h₂, h₃, h₄ необходимо определить вторую y ², третью y ² ¢ и четвертую y⁽⁴⁾ производные функции y (x). Эти производные аппроксимируем разделенными разностями второго, третьего и четвертого порядков соответственно.

В результате для получения значения функции y_i+1 по методу Рунге-Кутта выполняется следующая последовательность вычислительных операций:

Методы Рунге-Кутта относятся к так называемым одношаговым методам, поскольку для вычисления значения функции y (x) в точке x_{i+ 1} требуется знать только значение функции y (x) в одной предыдущей точке x_i.

Многошаговые методы построены путем интерполирования по нескольким соседним точкам; для их использования необходимо знать значение функции y (x)в нескольких предыдущих точках.

Эти методы численного интегрирования дифференциальных уравнений первого порядка были разработаны Адамсом в 1855 году.

Достоинство многошаговых методов состоит в том, что независимо от порядка метода для вычисления значения функции y (x) в одной точке требуется один раз вычислить функцию f (x, y).

Метод Адамса второго порядка записывается следующим образом

где .

Методы Адамса третьего и четвертого порядков имеют вид

Метод Рунге-Кутта четвертого порядка и метод Адамса четвертого порядка имеют одинаковую погрешность, но метод Адамса требует примерно вчетверо меньшего объема вычислений.

ПОРЯДОК ВЫПОЛНЕНИЯ РАБОТЫ

1. Ознакомиться с методами решения обыкновенных дифференциальных уравнений.

2. С помощью моделирующей программы решить указанные в задании дифференциальное уравнение по трем различным методам с одинаковым значением шага. Варианты заданий выбрать из табл. 12.

3. Пользуясь моделирующей программой решить дифференциальное уравнение одним из исследуемых методов при разных значениях шага.

Таблица 12

№№ вар.	Задание	№№ вар.	Задание
	,		,
	,		,
	,		,
	,		,
	,		,
	,		,

ОТЧЕТ О РАБОТЕ

Отчет должен содержать:

1. Исследуемое дифференциальное уравнение и краткое описание используемых методов решения.

2. Таблицу с вычисленными значениями решения дифференциального уравнения.

3. Таблицу и график зависимости решения дифференциального уравнения от величины шага.

КОНТРОЛЬНЫЕ ВОПРОСЫ

1. Что такое порядок метода?

2. Какие методы дают точное значение при решении дифференциального уравнения с линейной правой частью?

3. Что выгоднее – увеличивать порядок метода, или уменьшать величину шага?

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

СПИСОК РЕКОМЕНДУЕМОЙ ЛИТЕРАТУРЫ

1. Березин И. С., Жидков Н. П. Методы вычислений.- М.: Наука, 1966. -632 с.

2. Бодров В. И., Дворецкий С. И., Калинин В. Ф. Численные методы и программирование: Учебное пособие.- М.: МИХМ, 1986. - 92 с.

3. Демидович Б. П., Марон И. А. Основы вычислительной математики, - М.: Наука, 1970.- 664 с.

4. Копченова Н. В., Марон И. А. Вычислительная математика в примерах и задачах.- М.: Наука, 1972, - 337 с.

5. Крылов В. И., Бобков В. В., Монастырский П. И. Вычислительные методы.- М.: Наука, 1976.- 304 с.

6. Основные численные методы и их реализация на микрокалькуляторах/ И.. Сулима и др. – Киев: Виша школа, 1987. – 312 с.