Трёхстадийный метод Гаусса–Лежандра – это неявный метод Рунге–Кутты порядка
p=6
, построенный на квадратурной формуле Гаусса с тремя узлами.
Он относится к классу:
- симметричных,
- A-устойчивых,
- симплектических методов.
Метод особенно известен тем, что хорошо сохраняет динамику консервативных систем:
- колебания,
- гамильтоновы инварианты,
- фазовую структуру.
Задача
Рассматривается задача Коши:
y' = f(t,y), \:\:\: y(t_0) = y_0
(1.1)
Метод имеет:
- число стадий
s=3, - порядок точности
p=6.
Узлы Гаусса
Узлы на интервале
[0,1]
:
\begin{aligned} c_{1} & = \frac{1}{2} - \frac{\sqrt{15}}{10}, \:\:\: c_{2} & = \frac{1}{2}, \:\:\: c_{3} & = \frac{1}{2} + \frac{\sqrt{15}}{10}. \\
\end{aligned}
(1.2)
Таблица Батчера
\renewcommand{\arraystretch}{1.5}
\begin{array}{c|ccc}
\frac12-\frac{\sqrt{15}}{10}&\frac{5}{36}&\frac{2}{9}-\frac{\sqrt{15}}{15}&\frac{5}{36}-\frac{\sqrt{15}}{30}\\
\frac12&\frac{5}{36}+\frac{\sqrt{15}}{24}&\frac{2}{9}&\frac{5}{36}-\frac{\sqrt{15}}{24}\\
\frac12+\frac{\sqrt{15}}{10}&\frac{5}{36}+\frac{\sqrt{15}}{30}&\frac{2}{9}+\frac{\sqrt{15}}{15}&\frac{5}{36}\\[0.4em]
\hline
&\frac{5}{18}&\frac{4}{9}&\frac{5}{18}
\end{array}\:\:\:\:\:\rightarrow\:\:\:\:\:
\begin{matrix} \textbf{c} & \vline & \textbf{A} \\
\hline \, & \vline & \textbf{b}^T
\end{matrix}
(1.3)
Формулы метода
Стадии:
\begin{aligned} Y_{i} & = y_{n} + h \sum^{\substack{3}}_{\substack{j=1}}a_{i}f(t_{n}+c_{j}h, Y_{j}), \\
\end{aligned}
(1.4)
где:
-
Y_{i}– внутренние стадии, -
a_{ij}– коэффициенты таблицы Батчера.
Переход:
\begin{aligned} y_{n+1} & = y_{n} + h \sum^{\substack{3}}_{\substack{j=1}}b_{i}f(t_{n}+c_{i}h, Y_{i}). \\
\end{aligned}
(1.5)
Функция устойчивости
Для тестового уравнения:
\begin{aligned} y' & = \lambda y \\
\end{aligned}
(1.6)
получается рациональная функция устойчивости:
\begin{aligned} R(z) & = \frac{1+\frac{1}{2}z + \frac{1}{10}z^2 + \frac{1}{120}z^3}{1-\frac{1}{2}z+\frac{1}{10}z^2-\frac{1}{120}z^3}, \:\:\:\:\;\: z = h\lambda . \\
\end{aligned}
(1.6)
Это Pade-аппроксимация:
\begin{aligned} R(z) & = P_{[3/3]}(e^z). \\
\end{aligned}
(1.7)
Динамические характеристики
Как мы говорили ранее,
Динамической характеристикой
\textbf{D}
численного метода дискретизации ОДУ называется логарифм натуральный его функции устойчивости:
\textbf{D} = ln\:R(z) .
(1.7)
Так как, в общем случае,
z
– комплексная величина, то функцию
\textbf{D}
можно представить в виде комплексной динамической харатреристики и рассматривать вещественную:
\textbf{D}_{R} = Re\:\textbf{D},
(1.8)
и мнимую:
\textbf{D}_{I} = Im\:\textbf{D}
(1.9)
динамические характеристики.
Вещественная характеристика
На следующем рисунке представлена вещественная динамическая характеристика 3-стадийного метода Гаусса-Лежандра.
\textbf{D}_{R}
).
Если представить
z
в виде комплексного числа
z = \mu + i\nu
, то координатная ось, отрисованная синим цветом, будет вещественной, а ось, отрисованная зеленым – будет отражать
мнимую координату точки
z
.
Ось, представленная красным цветом, показывает величину
\textbf{D}_{R}
.
Клавишей <P> (на клавиатуре) можно переключать перспективный и ортогональный режимы камеры, а колесом мыши – отдалять или приближать изображение.
Идеальная вещественная динамическая характеристика представляет собой плоскость, которой принадлежит
мнимая ось
Im(z) = \nu
, и повернутая на 45 градусов вокруг оси
Im(z) = \nu
.
Фактически, отклонение вещественной динамической характеристики
\textbf{D}_{R}
от идеальной, представляет собой амплитудную погрешность численного метода
(см. Рис. 3.1).
Мнимая характеристика
Мнимая динамическая характеристика 3-стадийного метода Гаусса-Лежандра представлена на рис.2:
\textbf{D}_{I}
).
Идеальной мнимой динамической характеристикой является плоскость, которой принадлежит
вещественная ось
Re(z) = \mu
, и повернутая на 45 градусов вокруг оси
Re(z) = \mu
.
Фактически, отклонение мнимой динамической характеристики
\textbf{D}_{I}
от идеальной, представляет собой частотную погрешность численного метода
(см. Рис. 3.1).
Области динамической согласованности
На рис.3 представлена область (зеленый контур), внутри которой относительные погрешности метода не превыщают одного процента. Вычислить границы данной области легко по формулам (3.25) и (3.25), представленными выше Определения 2, динамической теории погрешностей численных методов.
Кроме того, на рис.3 представлены области (красные контуры), гарантированно содержащие нули функции устойчивости метода, а также области (синие контуры), гарантированно содержащие полюса функции устойчивости метода.