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

Материал из Википедии — свободной энциклопедии

Методы Рунге — Кутты^[1] — важное семейство численных алгоритмов решения (систем) обыкновенных дифференциальных уравнений. Данные итеративные методы явного и неявного приближенного вычисления были разработаны около 1900 года немецкими математиками К. Рунге и М. В. Куттой.

[править] Классический метод Рунге — Кутты 4 порядка

Метод Рунге — Кутты 4 порядка столь широко распространен, что его часто называют просто метод Рунге — Кутты.

Рассмотрим задачу Коши $\textbf{y}'=\textbf{f}(x,\textbf{y}), \textbf{y}(x_0)=\textbf{y}_0$ . Тогда значение в следующей точке вычисляется по формуле:

$\textbf{y}_{n+1} = \textbf{y}_n + {h \over 6} (\textbf{k}_1 + 2\textbf{k}_2 + 2\textbf{k}_3 + \textbf{k}_4)$

где

$\textbf{k}_1 = \textbf{f} \left( x_n, \textbf{y}_n \right),$

$\textbf{k}_2 = \textbf{f} \left( x_n + {h \over 2}, \textbf{y}_n + {h \over 2} \textbf{k}_1 \right),$

$\textbf{k}_3 = \textbf{f} \left( x_n + {h \over 2}, \textbf{y}_n + {h \over 2} \textbf{k}_2 \right),$

$\textbf{k}_4 = \textbf{f} \left( x_n + h, \textbf{y}_n + h\textbf{k}_3 \right),$

$\textbf{ }h$ — величина шага сетки по $\textbf{ }x$ .

Этот метод имеет 4 порядок, т.е. ошибка на каждом шаге составляет $O (h 5)$ , а суммарная ошибка на конечном интервале интегрирования $O (h 4)$ .

[править] Прямые методы Рунге — Кутты

Семейство прямых методов Рунге — Кутты является обобщением метода Рунге — Кутты 4 порядка. Оно задается формулами

$\textbf{y}_{n+1} = \textbf{y}_n + h\sum_{i=1}^s b_i \textbf{k}_i,$

где

$\textbf{k}_1 = f(x_n, \textbf{y}_n), \,$

$\textbf{k}_2 = f(x_n+c_2h, \textbf{y}_n+a_{21}h\textbf{k}_1), \,$

$\textbf{k}_3 = f(x_n+c_3h, \textbf{y}_n+a_{31}h\textbf{k}_1+a_{32}h\textbf{k}_2), \,$

$\vdots$

$\textbf{k}_s = f(t_n+c_sh, \textbf{y}_n+a_{s1}h\textbf{k}_1+a_{s2}h\textbf{k}_2+\cdots+a_{s,s-1}h\textbf{k}_{s-1}).$

Конкретный метод определяется числом $s$ и коэффициентами $b i, a i j$ и $c i$ . Эти коэффициенты часто упорядочивают в таблицу

0
$c 2$	$a 21$
$c 3$	$a 31$	$a 32$
$\vdots$	$\vdots$		$\ddots$
$c s$	$a s 1$	$a s 2$	$\cdots$	$a s, s - 1$
	$b 1$	$b 2$	$\cdots$	$b s - 1$	$b s$

Для коэффициентов метода Рунге — Кутты должны быть выполнены условия $\sum_{j=1}^{i-1} a_{ij} = c_i$ для $i=2, \ldots, s$ . Если мы хотим, чтобы метод имел порядок $p$ , то следует так же обеспечить условие $\bar\textbf{y}(h+x_0)-\textbf{y}(h+x_0)=O(h^{p+1})$ , где $\bar\textbf{y}(h+x_0)$ — приближение полученное по методу Рунге — Кутты. После многократного дифференцирования это условие преобразуется в систему полиномиальных уравнений на коэффициенты метода.

[править] Примечания

↑ Правильно на русском языке название методов будет звучать именно «Методы Рунге — Кутты», т.к. фамилия последнего автора заканчивается на гласную «а» в им. падеже и соотвественно склоняется гласной «ы» в родительном.

Категория: Численные методы

See also ebooksgratis.com: no banners, no cookies, totally FREE.

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

Материал из Википедии — свободной энциклопедии

[править] Классический метод Рунге — Кутты 4 порядка

[править] Прямые методы Рунге — Кутты

[править] Примечания

Views

Навигация

Участие

Поиск

На других языках