WikiDer > Явные и неявные методы

Explicit and implicit methods

Явные и неявные методы подходы, используемые в числовой анализ для получения численных приближений решений нестационарных обычный и уравнения в частных производных, как требуется в компьютерное моделирование из физические процессы. Явные методы вычислить состояние системы в более позднее время из состояния системы в текущий момент, в то время как неявные методы найти решение, решив уравнение, включающее как текущее, так и более позднее состояние системы. Математически, если текущее состояние системы и состояние в более позднее время ( - малый шаг по времени), то для явного метода

а для неявного метода решается уравнение

найти

Неявные методы требуют дополнительных вычислений (решения вышеуказанного уравнения), и их может быть намного сложнее реализовать. Неявные методы используются потому, что многие проблемы, возникающие на практике, являются жесткий, для которых использование явного метода требует непрактично малых шагов по времени чтобы ошибка в результате оставалась ограниченной (см. числовая стабильность). Для таких задач для достижения заданной точности требуется гораздо меньше вычислительного времени для использования неявного метода с большими временными шагами, даже с учетом того, что необходимо решать уравнение вида (1) на каждом временном шаге. Тем не менее, следует ли использовать явный или неявный метод, зависит от решаемой проблемы.

Поскольку неявный метод не может быть реализован для каждого типа дифференциального оператора, иногда рекомендуется использовать так называемый метод разделения операторов, что означает, что дифференциальный оператор переписывается как сумма двух дополнительных операторов

в то время как один обрабатывается явно, а другой - неявно. Для обычных приложений неявный член выбирается линейным, а явный член может быть нелинейным. Эта комбинация первого метода называется Неявно-явный метод (сокращенно IMEX [1], [2]).

Иллюстрация с использованием прямого и обратного методов Эйлера

Рассмотрим обыкновенное дифференциальное уравнение

с начальным условием Рассмотрим сетку для 0 ≤k ≤ п, то есть шаг по времени и обозначим для каждого . Дискретность это уравнение с использованием простейших явных и неявных методов, которые являются нападающий Эйлер и обратный Эйлер методы (см. числовые обыкновенные дифференциальные уравнения) и сравните полученные схемы.

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

Нападающий Метод Эйлера

дает

для каждого Это явная формула для .

Обратный метод Эйлера

С обратный метод Эйлера

находится неявное уравнение

за (сравните это с формулой (3), где был задан явно, а не как неизвестное в уравнении).

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

В подавляющем большинстве случаев уравнение, которое необходимо решить при использовании неявной схемы, намного сложнее, чем квадратное уравнение, и аналитического решения не существует. Затем используется алгоритмы поиска корней, Такие как Метод Ньютона, чтобы найти численное решение.

Кривошипный метод Николсона

С Метод Кранка-Николсона

находится неявное уравнение

за (сравните это с формулой (3), где был задан явно, а не как неизвестное в уравнении). Это можно решить численно, используя алгоритмы поиска корней, Такие как Метод Ньютона, чтобы получить .

Крэнк Николсон можно рассматривать как форму более общего IMEX (Ясложный-Бывшийсложные) схемы.

Вперед-назад метод Эйлера
Результат применения как метода Эйлера вперед-назад, так и метода Эйлера вперед-назад и .

Чтобы применить IMEX-схему, рассмотрим несколько иное дифференциальное уравнение:

Следует, что

и поэтому

для каждого

Смотрите также

Источники

  1. ^ У.М. Ашер, С.Дж. Руут, Р.Дж. Спитери: Неявно-явные методы Рунге-Кутты для нестационарных дифференциальных уравнений с частными производными, Appl Numer Math, т. 25 (2-3), 1997 г.
  2. ^ Л. Парески, Г. Руссо: Неявно-явные схемы Рунге-Кутты для жестких систем дифференциальных уравнений, Последние тенденции в численном анализе, Vol. 3, 269-289, 2000