Конечная разница - Википедия - Finite difference

А конечная разница математическое выражение вида ж (Икс + б) − ж (Икс + а). Если конечная разность делится на ба, каждый получает коэффициент разницы. Приближение производные конечными разностями играет центральную роль в методы конечных разностей для числовой решение дифференциальные уравнения, особенно краевые задачи.

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

Сегодня термин «конечная разница» часто используется как синоним конечно-разностные аппроксимации деривативов, особенно в контексте численные методы.[1][2][3] Конечно-разностные аппроксимации - это конечно-разностные отношения в терминологии, использованной выше.

Конечные различия были введены Брук Тейлор в 1715 г., а также изучались как абстрактные самостоятельные математические объекты в работах Джордж Буль (1860), Л. М. Милн-Томсон (1933), и Кароли Джордан (1939). Конечные различия уходят корнями в одну из Йост Бюрги алгоритмы (c. 1592) и работа других, в том числе Исаак Ньютон. Формальное исчисление конечных разностей можно рассматривать как альтернативу исчислению бесконечно малые.[4]

Основные типы

Три типа конечных разностей. Центральная разница относительно x дает наилучшее приближение производной функции в точке x.

Обычно рассматривают три основных типа: вперед, назад, и центральный конечные разности.[1][2][3]

А форвардная разница является выражением формы

В зависимости от области применения интервал час может быть переменным или постоянным. Если опущено, час принимается равным 1: Δ [ж ](Икс) = Δ1[ ж ](Икс).

А обратная разница использует значения функции в Икс и Иксчас, вместо значений в Икс + час иИкс:

Наконец, центральная разница дан кем-то

Связь с производными

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

В производная функции ж в какой-то момент Икс определяется предел.

Если час имеет фиксированное (ненулевое) значение вместо того, чтобы приближаться к нулю, то правая часть приведенного выше уравнения будет записана

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

Та же формула верна для обратной разницы:

Однако центральная (также называемая центрированной) разность дает более точное приближение. Если ж дважды дифференцируема,

Главная проблема[нужна цитата ] Однако в случае метода центральной разности осциллирующие функции могут давать нулевую производную. Если ж (нэ) = 1 за п странно, и ж (нэ) = 2 за п даже тогда ж ′(нэ) = 0 если рассчитывается по центральной разностной схеме. Это особенно неприятно, если домен ж дискретно. Смотрите также Симметричная производная

Авторы, для которых конечные разности означают аппроксимации конечных разностей, определяют прямые / обратные / центральные разности как коэффициенты, приведенные в этом разделе (вместо использования определений, данных в предыдущем разделе).[1][2][3]

Различия высшего порядка

Аналогичным образом можно получить конечно-разностные аппроксимации производных высших порядков и дифференциальных операторов. Например, используя приведенную выше формулу центральной разности для ж ′(Икс + час/2) и ж ′(Иксчас/2) и применяя формулу центральной разности для производной от ж ′ в Икс, получаем центрально-разностное приближение второй производной ж:

Центральный

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

Нападающий второго порядка
Второй порядок назад

В более общем плане ппорядок вперед, назад и по центру различия даны соответственно

Вперед

или для час = 1,

Назад
Центральная

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

Обратите внимание, что центральная разница будет для нечетных п, имеют час умноженное на нецелые числа. Это часто является проблемой, потому что это означает изменение интервала дискретизации. Проблема может быть решена в среднем δп[ ж ](Иксчас/2) и δп[ ж ](Икс + час/2).

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

Связь этих разностей более высокого порядка с соответствующими производными очевидна,

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

приблизительно ж ′(Икс) до срока заказа час2. Это можно доказать, расширив приведенное выше выражение до Серия Тейлор, или используя исчисление конечных разностей, объясненное ниже.

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

Ядра произвольного размера

Используя линейную алгебру, можно построить конечно-разностные аппроксимации, которые используют произвольное количество точек слева и (возможно, различное) количество точек справа от точки оценки для производной любого порядка. Это включает решение такой линейной системы, что Расширение Тейлора суммы этих точек вокруг точки оценки наилучшим образом аппроксимирует разложение Тейлора желаемой производной. Такие формулы могут быть представлены графически на гексагональной или ромбовидной сетке.[5]

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

Подробности изложены в этих Примечания.

В Калькулятор конечно-разностных коэффициентов строит конечно-разностные аппроксимации для нестандартных (и даже нецелочисленных) шаблонов по произвольному шаблону и желаемому порядку производной.

Характеристики

  • Для всех положительных k и п

В дифференциальных уравнениях

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

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

Серия Ньютона

В Серия Ньютон состоит из условий Уравнение прямой разности Ньютона, названный в честь Исаак Ньютон; по сути, это Формула интерполяции Ньютона, впервые опубликованный в его Principia Mathematica в 1687 г.,[6] а именно дискретный аналог непрерывного разложения Тейлора,

что справедливо для любого многочлен функция ж и для многих (но не для всех) аналитические функции (Не действует, когда ж является экспоненциальный тип . Это легко увидеть, поскольку функция синуса обращается в нуль при целых кратных ; соответствующий ряд Ньютона тождественно равен нулю, поскольку все конечные разности в этом случае равны нулю. Однако ясно, что синусоидальная функция не равна нулю.) Здесь выражение

это биномиальный коэффициент, и

это "падающий факториал "или" нижний факториал ", а пустой продукт (Икс)0 определяется равным 1. В данном конкретном случае предполагается, что изменение значений Икс, час = 1 обобщения ниже.

Обратите внимание на формальное соответствие этого результата Теорема Тейлора. Исторически это, как и Тождество Чу – Вандермонда,

(вытекающая из нее и соответствующая биномиальная теорема ), входят в наблюдения, созревшие до системы темный камень.

Чтобы проиллюстрировать, как можно использовать формулу Ньютона на практике, рассмотрим несколько первых членов удвоения Последовательность Фибоначчи ж = 2, 2, 4, ... Можно найти многочлен который воспроизводит эти значения, сначала вычисляя таблицу различий, а затем подставляя различия, соответствующие Икс0 (подчеркнут) в формулу следующим образом:

Для случая неравномерных шагов по значениям Икс, Ньютон вычисляет разделенные различия,

серия товаров,

и полученный многочлен является скалярное произведение,[7]

.

В анализе с п-адические числа, Теорема Малера утверждает, что предположение, что ж является полиномиальной функцией, можно ослабить вплоть до предположения, что ж просто непрерывно.

Теорема Карлсона обеспечивает необходимые и достаточные условия единственности ряда Ньютона, если он существует. Однако рядов Ньютона, как правило, не существует.

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

В сжатом и несколько более общем виде и равноудаленных узлах формула выглядит так:

Исчисление конечных разностей

Разницу вперед можно рассматривать как оператор, называется оператор разницы, который отображает функцию ж к Δчас[ ж ].[8][9] Этот оператор составляет

куда Тчас это оператор смены с шагом час, определяется Тчас[ ж ](Икс) = ж (Икс + час), и я это оператор идентификации.

Конечная разность высших порядков может быть определена рекурсивным образом как Δп
час
≡ Δчасп − 1
час
)
. Другое эквивалентное определение: Δп
час
= [Тчася]п
.

Оператор разницы Δчас это линейный оператор, как таковая удовлетворяет Δчас[αf + βg](Икс) = α Δчас[ ж ](Икс) + β Δчас[грамм](Икс).

Он также удовлетворяет особому Правило Лейбница указанное выше,Δчас(ж (Икс)грамм(Икс)) = (Δчасж (Икс)) грамм(Икс+час) + ж (Икс) (Δчасграмм(Икс)). Аналогичные утверждения справедливы для обратных и центральных различий.

Формально применяя Серия Тейлор относительно час, дает формулу

куда D обозначает оператор производной континуума, отображение ж к производной ж ′. Расширение действительно, когда обе стороны действуют на аналитические функции, для достаточно малых час. Таким образом, Тчас = еHD, и формально обращая экспоненциальные доходности

Эта формула верна в том смысле, что оба оператора дают один и тот же результат при применении к многочлену.

Даже для аналитических функций не гарантируется сходимость ряда справа; это может быть асимптотический ряд. Однако его можно использовать для получения более точных приближений для производной. Например, сохранение первых двух членов ряда дает приближение второго порядка к ж ′(Икс) упоминается в конце раздел Различия высшего порядка.

Аналогичные формулы для обратного и центрального разностных операторов имеют вид

Исчисление конечных разностей связано с темный камень комбинаторики. Это удивительно систематическое соответствие связано с идентичностью коммутаторы теневых величин к их континуальным аналогам (час → 0 пределы),

Большое количество формальных дифференциальных соотношений стандартного исчисления с функциями ж (Икс) таким образом систематически отображать в мрачные конечно-разностные аналоги с участием ж (xT−1
час
)
.

Например, темный аналог одночлена Иксп является обобщением указанного падающего факториала (Почхаммер k-символ ),

так что

отсюда вышеприведенная формула интерполяции Ньютона (путем согласования коэффициентов в разложении произвольной функции ж (Икс) в таких символах) и так далее.

Например, теневой синус

Как и в непрерывном пределе, собственная функция Δчас/час тоже бывает экспоненциальной,

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

Так, например, Дельта-функция Дирака соответствует своему мрачному корреспонденту, функция кардинального синуса,

и так далее.[11] Разностные уравнения часто можно решить с помощью техник, очень похожих на методы решения дифференциальные уравнения.

Обратный оператор прямого разностного оператора, следовательно, умбральный интеграл - это неопределенная сумма или оператор антиразличия.

Правила исчисления конечно-разностных операторов

Аналогично правила нахождения производной, у нас есть:

Все вышеперечисленные правила одинаково хорошо применимы к любому разностному оператору, включая относительно Δ.

или же

См. Ссылки.[12][13][14][15]

Обобщения

  • А обобщенная конечная разность обычно определяется как

куда μ = (μ0,… μN) - его вектор коэффициентов. An бесконечная разница является дальнейшим обобщением, в котором конечная сумма выше заменена на бесконечная серия. Другой способ обобщения - это сделать коэффициенты μk зависит от точки Икс: μk = μk(Икс), таким образом учитывая взвешенная конечная разность. Также можно сделать шаг час зависит от точки Икс: час = час(Икс). Такие обобщения полезны для построения различных модуль непрерывности.

  • Обобщенную разность можно рассматривать как кольца многочленов р[Тчас]. Это приводит к разностным алгебрам.
  • Оператор разности обобщается на Инверсия Мёбиуса через частично заказанный набор.
  • Как оператор свертки: с помощью формализма алгебры инцидентности, разностные операторы и другое обращение Мёбиуса могут быть представлены как свертка с функцией на poset, называемой Функция Мёбиуса μ; для разностного оператора, μ - последовательность (1, −1, 0, 0, 0, ...).

Многомерные конечные разности

Конечные различия могут рассматриваться более чем в одной переменной. Они аналогичны частные производные в нескольких переменных.

Некоторые приближения частных производных:

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

поскольку единственные значения для вычисления, которые еще не нужны для предыдущих четырех уравнений: ж (Икс + час, у + k) и ж (Иксчас, уk).

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

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

  1. ^ а б c Пол Уилмотт; Сэм Ховисон; Джефф Дьюинн (1995). Математика финансовых деривативов: введение для студентов. Издательство Кембриджского университета. п.137. ISBN  978-0-521-49789-3.
  2. ^ а б c Питер Олвер (2013). Введение в уравнения с частными производными. Springer Science & Business Media. п. 182. ISBN  978-3-319-02099-0.
  3. ^ а б c М Ханиф Чаудри (2007). Открытый канал потока. Springer. п. 369. ISBN  978-0-387-68648-6.
  4. ^ Йордан, соч. соч., стр. 1 и Милн-Томсон, стр. xxi. Милн-Томсон, Луи Мелвилл (2000): Исчисление конечных разностей (Chelsea Pub Co, 2000) ISBN  978-0821821077
  5. ^ Фрейзер, Дункан К. (1 января 1909 г.). "О графическом обозначении формулы интерполяцииæ". Журнал института актуариев. 43 (2): 235–241. Дои:10.1017 / S002026810002494X. Получено 17 апреля, 2017.
  6. ^ Ньютон, Исаак, (1687). Principia, Книга III, лемма V, случай 1
  7. ^ Рихтмейер, Д. и Мортон К.В. (1967). Методы различия для задач с начальным значением, 2-е изд., Wiley, New York.
  8. ^ Буль, Джордж, (1872). Трактат об исчислении конечных разностей, 2-е изд., Macmillan and Company. В сети. Также [Dover edition 1960]
  9. ^ Джордан, Чарльз, (1939/1965). «Исчисление конечных разностей», Chelsea Publishing. В сети: [1]
  10. ^ Захос, К. (2008). «Темные деформации в дискретном пространстве-времени». Международный журнал современной физики A. 23 (13): 2005–2014. arXiv:0710.2306. Bibcode:2008IJMPA..23.2005Z. Дои:10.1142 / S0217751X08040548.
  11. ^ Curtright, T. L .; Захос, К. К. (2013). "Umbral Vade Mecum". Границы физики. 1: 15. arXiv:1304.0429. Bibcode:2013FrP ..... 1 ... 15C. Дои:10.3389 / fphy.2013.00015.
  12. ^ Levy, H .; Лессман, Ф. (1992). Конечно-разностные уравнения. Дувр. ISBN  0-486-67260-3.
  13. ^ Эймс, У. Ф. (1977). Численные методы для уравнений с частными производными., Раздел 1.6. Academic Press, Нью-Йорк. ISBN  0-12-056760-1.
  14. ^ Хильдебранд, Ф. Б., (1968). Конечно-разностные уравнения и моделирование, Раздел 2.2, Прентис-Холл, Энглвуд Клиффс, Нью-Джерси.
  15. ^ Флажолет, Филипп; Седжвик, Роберт (1995). «Преобразования Меллина и асимптотика: конечные разности и интегралы Райса» (PDF). Теоретическая информатика. 144 (1–2): 101–124. Дои:10.1016 / 0304-3975 (94) 00281-М..
  • Ричардсон, К. Х. (1954): Введение в исчисление конечных разностей (Ван Ностранд (1954) онлайн-копия
  • Миккенс, Р. Э. (1991): Разностные уравнения: теория и приложения (Чепмен и Холл / CRC) ISBN  978-0442001360

внешняя ссылка