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 , |
|
что вторая разность линейной функции равна нулю.
| 05_Холера |
| 1 |
| 10.4. Исследование регистров |
| 1112 |
| 12 |
| 1285 |
| 1560 |
| 1568 |
| 1604248853606027 |
| 17 |