Материал: 1395

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

11

 

Вследствие

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

функции

f (t) , функция

fk,

является

непрерывной по аргументу и удовлетворяет условию

fk 1,1

fk,0 .

 

 

 

Если функция

f (t)

терпит разрыв первого рода в точках tk

 

k

T , то

данное

равенство

не выполняется, так как

lim fk 1,

lim fk,

.

В

этом

 

 

 

 

 

 

 

 

 

 

 

1

 

0

 

 

 

случае,

под

значением

 

функции

fk

принято принимать

предел

справа

fk

lim fk,

,

то

есть

значения

переменной

рассматриваются на

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

полуинтервале

0

 

1.

Функция

fk,

является, вообще говоря,

функцией

двух аргументов k

и

,

однако при

0 второй аргумент, то есть

 

обычно

опускается.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

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

Любая функция, определенная в равноотстоящие отрезки времени, путем выбора масштаба может быть приведена в соответствие обозначению f (k T ) f (k) fk , где k принимает целочисленные значения.

Оператор сдвига. Определим оператор сдвига E (или упреждения), как оператор, сдвигающий отсчет функции на единицу вперед

E

fk

fk 1.

Повторное применение оператора сдвига дает

E2 fk

E (E fk )

fk 2 .

В общем случае для любого целого n имеем

 

En

fk

fk n .

Аналогичным образом может быть определен оператор задержки

E 1, как оператор, сдвигающий отсчет функции на единицу назад

E 1 fk

fk

1.

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

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

12

fk fk 1 fk .

Это так называемый правый разностный оператор или разность вперед, в отличие от используемого иногда левого разностного оператора или разности назад

fk fk fk 1 .

Далее будем использовать в основном правый разностный оператор и будем именовать его просто разностным оператором.

Из определений разностного оператора и оператора сдвига следует их взаимосвязь

 

 

 

fk

(E

1) fk ;

 

 

 

 

E

1.

В общем случае для любого целого k имеем

 

 

 

k

(E 1)k .

 

Разностный оператор или разность второго порядка определяется

выражением

 

 

 

 

2

fk

fk 1

fk ( fk 2

fk 1) ( fk 1 fk ) fk 2 2 fk 1 fk .

 

Разностный оператор или разность третьего порядка, соответственно, запишется

3 fk

2 fk

1

2

fk

 

( fk 3

2 fk 2

fk 1) ( fk 2 2 fk 1 fk )

 

fk 3

3 fk 2

3 fk 1 fk .

В общем случае, разностный оператор или разность k - го порядка определяется выражением

 

 

 

 

 

n fk

n

 

n

 

 

 

 

 

 

 

 

 

( 1)i

i

fk n i ,

 

 

 

 

 

 

 

 

i

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

n

Cni

n!

 

 

n (n

1)

(n

2) (n

i)

- число сочетаний из n по

 

 

i

i! (n i)!

 

 

1 2 3 i

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

Учитывая связь разностного оператора и оператора сдвига, последнее

выражение можно переписать в виде

 

 

 

 

 

 

 

 

n fk

(E 1)n

fk

n

( 1)i

n

En i

fk .

 

 

 

 

i

 

 

 

 

 

 

 

 

i

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Так из выражения первой разности получаем

fk 1 fk

fk .

Расписывая вторую разность в виде

2

fk

fk 2 fk 1

fk fk 2 fk 2

fk ,

 

13

получаем

 

 

 

fk 2

fk

2

 

fk

2 fk .

 

 

Раскрывая выражение разности третьего порядка

 

 

3 fk

2 fk 1

2 fk

fk 3 2 fk 2

fk 1

2 fk

 

 

fk 3

2 fk

4

fk

2

2 fk

fk

fk

2 fk

 

fk 3

fk

3

fk

3

2 fk ,

 

 

 

получаем

 

 

 

 

 

 

 

 

 

 

 

 

fk 3

fk

3

fk

3

2

fk

3 fk .

 

Обобщая данные выражения, получаем общую формулу разложения

дискретной функции в ряд по разностным операторам

 

 

 

 

 

fk n

n

n

i

fk .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

i 0

 

 

 

 

 

В частности при k

0, можем записать

 

 

 

 

 

 

 

 

 

fn

n

n

i

f0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

i

0

 

 

 

 

 

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

Свойства разностных операторов. Отметим, что разностный оператор является линейным оператором. Действительно, рассматривая линейную суперпозицию дискретных функций fi , имеем

k

k

 

k

k

 

ci

fi,n

ci fi,n 1

ci fi,n

ci

fi,n ,

i 1

i 1

i

1

i 1

 

где ci - коэффициенты линейной суперпозиции. Следствием этого выражения являются соотношения

 

 

 

 

c

fk

c

fk ;

 

 

 

 

 

 

 

n [ fk

gk ]

n fk

n gk ;

 

 

 

 

 

m

n

fk

m n

fk .

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, разностный оператор от суммы функций равен сумме

разностных операторов от каждой функции.

 

 

 

Заметим, что эти соотношения справедливы при замене разностного

оператора

на оператор сдвига E . Операторы

и E коммутативны между

собой,

но

не с

функциями,

так,

например

( fk gk ) fk

gk или

2

fk

(

fk ) (

fk ) .

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим разность от произведения двух дискретных функций

 

 

 

14

 

 

( fk gk ) fk 1 gk 1

fk gk

 

 

fk 1

gk 1

fk 1 gk

fk 1 gk fk gk

fk

1

gk

fk gk .

 

Если вместо произведения

fk

1 gk

прибавим и вычтем fk gk 1, то получим

( fk gk )

fk gk

fk

gk .

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

 

2 ( fk gk )

 

fk

2 gk

2

fk 1

 

gk

2 fk gk 2

 

 

 

 

2

fk gk

2

fk

 

gk 1

fk 2

2 gk .

 

В общем случае для разности n - го порядка получим формулу

n

n

n

i

 

n i

 

n

n

n i

 

i gk .

( fk gk )

i

fk

gk i

 

i

fk i

 

i 1

 

 

 

 

i 1

 

 

 

Таким образом, разностные формулы от произведения дискретных функций соответствуют по структуре формулам дифференцирования произведения непрерывных функций.

Рассмотрим разность от отношения двух дискретных функций

 

 

fk

 

fk 1

 

 

fk fk 1 gk

fk gk

 

 

fk gk

fk gk 1

 

 

gk

 

gk 1

 

 

gk

 

 

 

 

 

 

 

gk 1 gk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

fk

gk

 

fk

gk

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

gk 1 gk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Разность второго порядка от отношения функций определяется

аналогично

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 fk

 

 

fk 1

 

 

fk

 

fk 2

 

fk 1

 

 

fk 1

 

fk

 

 

 

 

 

 

gk

 

 

 

gk 1

 

gk

 

 

gk 2

 

 

gk 1

 

 

 

gk 1

 

 

gk

 

 

 

 

 

 

 

 

 

 

 

fk 1 gk 1

fk 1

 

gk 1

 

 

 

fk

gk

 

fk

 

gk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

gk 2 gk 1

 

 

 

 

 

 

 

 

 

gk 1 gk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отметим, что разностный оператор

 

имеет некоторые общие свойства

с оператором

дифференцирования

D

d / dt

непрерывных

функций. Так

производная непрерывной функции определяется в виде

 

df (t)

lim

f (t

T )

f (t)

.

 

dt

 

T

 

 

 

T

0

 

 

 

 

Используя T в качестве периода дискретизации запишем оператор сдвига и

разностный оператор

 

 

 

 

 

 

 

 

 

E f (t)

f (t T ) ,

 

 

f (t)

f (t T )

f (t) .

Из последнего выражения в частности следует, что

 

 

df (t)

lim

f (t)

,

 

 

 

 

 

 

 

 

dt

T

0

T

 

 

 

 

 

15

 

 

или в общем случае

 

 

 

 

d m f (t)

lim

m f (t)

.

 

dtm

T m

 

T 0

 

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

(k)(m) k (k 1) (k 2) (k m 1) ,

где m 1- число сомножителей.

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

 

(k)(m)

m (k)(m 1)

m k (k 1) (k 2)

(k m 2) .

Действительно

 

 

 

 

 

 

(k)(m)

(k 1) k (k 1) (k 2)

(k m 2)

 

k (k 1) (k 2) (k m 1)

 

 

 

 

[k (k 1)

(k 2)

(k

m

2)]

[(k 1) (k

m 1)] m (k)(m 1) .

Разность второго порядка запишется в виде

 

2 (k)(m)

m (k)(m 1)

m (k)(m 1) m (m 1) k(m 2) .

Соответственно разность порядка n

(m

1) запишется в виде

 

n (k)(m)

m (m

1)

[m

(n

1)] k (m n)

m(n) k (m n) .

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

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

Пример 1. Пусть задана постоянная дискретная функция fk

a , где a -

константа, и требуется определить разность первого порядка.

Согласно

определению разностного оператора, получаем

 

 

fk fk 1 fk a a 0 .

Таким образом, первая разность постоянной функции равна нулю.

Пример 2.

Пусть

задана линейная дискретная функция вида

fk a k b , где a

и b -

постоянные коэффициенты. Определим первую

разность этой функции

 

 

fk

[a (k 1) b] [a k b] a .

Вывод – первая разность линейной дискретной функции есть постоянная величина. Вычислив вторую разность, получаем

2

fk

a 0 ,

 

что вторая разность линейной функции равна нулю.