logo search
флеров

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

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

Пусть дана функция , определенная на множестве действительных (возможно целых) чисел и принимающая действительные значения. Определим новую функцию (x), называемую первой разностью , формулой

.

Оператор  называется разностным оператором первого порядка; кратко и очень упрощенно можно определить исчисление конечных разностей как исследование оператора  . Можно применить оператор  k раз и получить k-ый разностный оператор

.

Число k(x) называется k-ой разностью  в точке x (k(0) называется k-ой разностью  в 0).

Определим другой оператор E, называемый оператором сдвига, формулой

.

Таким образом, =E  I, где I означает единичный оператор:

.

Тогда первая разность функции может быть записана в виде:

.

Разности более высоких порядков определяются рекуррентным соотношением

Откуда получаем выражение для n -ой разности:

В частности,

, (1.6)

что дает явную формулу для k-ой разности в терминах значений (0), (1), ... , (k). Нетрудно обратить формулу ( 1 .6) и выразить (n) через i (0). Именно,

. (1.7)

Напишем теперь в строку значения

... (-2) (-1) (0) (1) (2) (3)...

Если внизу написать между каждой парой последовательных членов (i), (i +1) их разность (i +1)  (i) = (i), то получим последовательность

...(-2) (-1) (0) (1) (2) ... .

Повторение этой процедуры приводит к таблице разностей функции , k-ая строка которой состоит из значений k(n). Диагональ, начинающаяся в (0) и идущая направо вниз, состоит из разностей k (0) в 0. Например, пусть (n) = n4 . Таблица разностей (начинающаяся с (0) ) выглядит так:

0

1

16

81

256

625

...

1

15

65

175

369

14

50

110

194

36

60

84

24

24

0


Из формулы ( 1 .7) следует, что

(1.8)

В этом случае, так как n4 - многочлен четвертой степени и при фиксированном k есть многочлен степени k, написанное выше разложение обрывается после члена , то есть k 04 = 0, если k > 4 (или, более общим образом, k n4 = 0, если k > 4).

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

1. Функция  - полином степени, не превосходящей d, тогда и только тогда. когда d+1(n) =0 (или d(n) - постоянная).

2. Если многочлен (n) степени, не превосходящей d, разложен в ряд по базису , 0  k  d, то коэффициенты разложения есть k (0), то есть

.

Еще одна связь комбинаторных объектов с исчислением конечных разностей дается формулой ( 1 .15).