next up previous
Next: Методы Рунге-Кутта с интерполяцией Up: О применении i-гладкого анализа Previous: Методы Рунге-Кутта второго порядка

6. Методы Рунге-Кутта с кусочно-линейной интерполяцией предыстории и экстраполяцией предыстории модели по методу Эйлера

Согласно теореме 3.1 порядок сходимости методов Рунге-Кутта (2.1) - (2.3) при кусочно-постоянной интерполяции (2.4) и экстраполяции (2.5) определяется не только порядком аппроксимации метода, но и порядком модуля непрерывности точного решения. Даже если точное решение достаточно гладкое, например, выполняется условие 4.1, то модуль непрерывности имеет первый порядок и, поэтому, это обеспечивает лишь первый порядок сходимости при любом порядке аппроксимации метода (2.1) - (2.3). Ниже приводится модификация метода Рунге-Кутта, обеспечивающая второй порядок сходимости (при условии, что порядок аппроксимации метода (2.1) - (2.3) не ниже второго).

Введем функцию tex2html_wrap_inline1735, tex2html_wrap_inline1965, кусочно-линейную при tex2html_wrap_inline1737:
equation466
Таким образом tex2html_wrap_inline1969, если tex2html_wrap_inline1971, tex2html_wrap_inline1849, а между точками tex2html_wrap_inline1975 проведена линейная интерполяция при tex2html_wrap_inline1849.

Функция tex2html_wrap_inline1739 при tex2html_wrap_inline1741 продолжена вправо (проведена экстраполяция за точку tex2html_wrap_inline1815) по методу Эйлера:
displaymath1985

equation486
Заметим, что в отличие от (2.4), (2.5) функция tex2html_wrap_inline1735 в (6.1) - (6.2) является непрерывной.

Теорема 6.1. Пусть решение x(t) дважды непрерывно дифференцируемо на tex2html_wrap_inline1785, тогда если метод (2.1) - (2.3) , (6.1), (6.2) имеет порядок аппроксимации l, то он сходится, причем порядок сходимости не меньше минимума из l и 2.

Доказательство теоремы в основном повторяет доказательство теоремы 3.1, поэтому укажем лишь отличия.

В силу линейной интерполяции (6.1) имеем при tex2html_wrap_inline1999, k=1,...,n,
displaymath2003
где tex2html_wrap_inline2005.

Тогда вместо оценки (3.12) имеем оценку
displaymath2007

Для получения аналога оценки (3.14) проведем оценку
displaymath2009

displaymath2011

displaymath2013

displaymath2015

displaymath2017

displaymath2019

displaymath2021

displaymath2023
Поэтому, вместо (3.13) имеем
displaymath2025
и оценка (3.14) справедлива при так определенном tex2html_wrap_inline2027.

Вместо (3.18) имеем
displaymath2029
а вместо (3.19) получаем
displaymath2031
где
displaymath2033

displaymath2035

displaymath2037
Далее теорема 6.1 доказывается как и теорема 3.1, только вместо tex2html_wrap_inline2039 из (3.21) нужно взять
displaymath2041

Замечание 6.1. Теорема 6.1 справедлива также, если точное решение x(t) дважды непрерывно дифференцируемо не на всем отрезке tex2html_wrap_inline1785, а лишь на полуинтервалах tex2html_wrap_inline2047, tex2html_wrap_inline2049 при n=0,...,N-1.

Замечание 6.2. Интерполяцию предыстории модели (6.1) и ее экстраполяцию (6.2) имеет смысл применять, когда используется метод Рунге-Кутта второго порядка аппроксимации (см. раздел 5). В этом случае общий порядок сходимости равен 2. Экстраполяция предыстории модели по методу Эйлера (6.2) играет роль прогноза, а роль коррекции играет сам метод Рунге-Кутта второго порядка.


next up previous
Next: Методы Рунге-Кутта с интерполяцией Up: О применении i-гладкого анализа Previous: Методы Рунге-Кутта второго порядка