Формула обратной дифференциации - Backward differentiation formula - Wikipedia

В формула обратного дифференцирования (BDF) - это семейство неявных методов для численное интегрирование обыкновенных дифференциальных уравнений. Они есть линейные многоступенчатые методы чтобы для данной функции и времени аппроксимировать производную этой функции, используя информацию из уже вычисленных моментов времени, тем самым увеличивая точность приближения. Эти методы особенно используются для решения жесткие дифференциальные уравнения. Впервые методы были представлены Чарльз Ф. Кертисс и Джозеф О. Хиршфельдер в 1952 г.[1]

Общая формула

BDF используется для решения проблема начального значения

Общая формула для BDF может быть записана как [2]

куда обозначает размер шага и . С оценивается для неизвестного , Методы BDF скрытый и, возможно, потребуют решения нелинейных уравнений на каждом шаге. Коэффициенты и выбираются так, чтобы метод достигал порядка , что является максимально возможным.

Вывод коэффициентов

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

Конкретные формулы

В s-шаговые BDFs с s <7 являются:[3]

  • BDF1: (это обратный метод Эйлера )
  • BDF2:
  • BDF3:
  • BDF4:
  • BDF5:
  • BDF6:

Методы с s > 6 не нулевой стабильный поэтому их нельзя использовать.[4]

Стабильность

Устойчивость численных методов решения жесткие уравнения обозначается их областью абсолютной стабильности. Для методов BDF эти области показаны на графиках ниже.

В идеале область содержит левую половину комплексной плоскости, и в этом случае метод называется A-устойчивым. Тем не мение, линейные многоступенчатые методы с порядком больше 2 не может быть А-стабильный. Область устойчивости методов BDF более высокого порядка содержит большую часть левой полуплоскости и, в частности, всю отрицательную действительную ось. Методы BDF являются наиболее эффективными линейными многоступенчатыми методами такого рода.[4]

Рекомендации

Цитаты

  1. ^ Кертисс, К. Ф., и Хиршфельдер, Дж. О. (1952). Интегрирование жестких уравнений. Слушания Национальной академии наук, 38 (3), 235-243.
  2. ^ Ашер и Петцольд 1998, §5.1.2, с. 129
  3. ^ Изерлес 1996, п. 27 (для s = 1, 2, 3); Сюли и Майерс 2003, п. 349 (для всех s)
  4. ^ а б Сюли и Майерс 2003, п. 349

Упомянутые работы

  • Ascher, U.M .; Петцольд, Л. (1998), Компьютерные методы решения обыкновенных дифференциальных и дифференциально-алгебраических уравнений, СИАМ, Филадельфия, ISBN  0-89871-412-5.
  • Изерлес, Арье (1996), Первый курс численного анализа дифференциальных уравнений, Издательство Кембриджского университета, ISBN  978-0-521-55655-2.
  • Сули, Эндре; Майерс, Дэвид (2003), Введение в численный анализ, Издательство Кембриджского университета, ISBN  0-521-00794-1.

дальнейшее чтение

  • Методы BDF в вики SUNDIALS (SUNDIALS - это библиотека, реализующая методы BDF и подобные алгоритмы).