Явный метод Эйлера

Явный метод Эйлера 1-го порядка (RK1) – самый простой явный метод. Он имеет:

  • порядок точности p=1
  • число стадий s=1
  • явную схему
  • область устойчивости конечного размера
  • хорошую устойчивость для не жёстких задач

Для задачи Коши:

y' = f(t,y), \:\:\: y(t_0) = y_0

(1.1)

при шаге h вычисляются один промежуточный наклон:

k_1

(1.2)

После чего берётся взвешенная комбинация.

Смысл метода:

  • k_1 – наклон в начале шага

Итог:

  • средняя комбинация аппроксимирует интеграл по шагу с ошибкой O(h^2)

Формулы метода

\begin{aligned} k_1 & = f(t_n,y_n). \\ \end{aligned}

(1.3)

Вычисление значения на следующем шаге:

\begin{aligned} y_{n+1} & = y_n + k_1. \\ \end{aligned}

(1.4)

Таблица Батчера

\left.\begin{matrix} 0 \\ \hline \, \end{matrix}\right| \begin{matrix} \: 0 \\ \hline 1 \\ \end{matrix}\:\:\:\:\:\rightarrow\:\:\:\:\: \begin{matrix} \textbf{c} & \vline & \textbf{A} \\ \hline \, & \vline & \textbf{b}^T \end{matrix}

(1.5)

Функция устойчивости

\begin{aligned} R(z) & = 1 + z. \\ \end{aligned}

(1.6)

Динамические характеристики

Как мы говорили ранее, Динамической характеристикой \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)

динамические характеристики.

Вещественная характеристика

На следующем рисунке представлена вещественная динамическая характеристика явного метода Эйлера.

Fig.1. Вещественная динамическая характеристика явного метода Эйлера ( \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).

Мнимая характеристика

Мнимая динамическая характеристика явного метода Эйлера представлена на рис.2:

Fig.2. Мнимая динамическая характеристика явного метода Эйлера ( \textbf{D}_{I} ).

Идеальной мнимой динамической характеристикой является плоскость, которой принадлежит вещественная ось Re(z) = \mu , и повернутая на 45 градусов вокруг оси Re(z) = \mu .

Фактически, отклонение мнимой динамической характеристики \textbf{D}_{I} от идеальной, представляет собой частотную погрешность численного метода (см. Рис. 3.1).

Области динамической согласованности

На рис.3 представлена область (зеленый контур), внутри которой относительные погрешности метода не превыщают одного процента. Вычислить границы данной области легко по формулам (3.25) и (3.25), представленными выше Определения 2, динамической теории погрешностей численных методов.

Fig.3. Область динамической согласованности метода Рунге-Кутты.

Кроме того, на рис.3 представлена область устойчивости явного метода Эйлера (контур цвета морской волны) и области, гарантированно содержащие нули функции устойчивости метода (красные контуры).