Материал: 1395

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

26

уравнение. Так, если a0 0 , а an

уравнения равен n 1.

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

высшей и низшей степеней оператора сдвига

E .

В связи с этим, при

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

 

 

 

 

 

n

an 1

n 1

 

a1

a0

yk

fk ,

 

 

an

 

 

 

установить

порядок

непосредственно

по

его

структуре

невозможно.

Например,

уравнение

2

3

2

yk

0 эквивалентно

уравнению

 

E2 E yk 0 и, следовательно, это уравнение первого, а не второго

порядка.

Известно, что однородное разностное уравнение n - го порядка

a En

a

1

En 1

a E a y

k

0

,

n

n

 

1

0

 

 

где an 0, a0 0 и ai -

определены для всех целых значений k , имеет n

линейно независимых решений.

Характеристическое уравнение, соответствующее разностному уравнению, соответствует левой части уравнения при замене оператора сдвига переменной p

a pn

a

pn 1

a p a 0 .

n

n 1

 

1

0

В общем случае характеристическое уравнение имеет n - корней pi d i , среди которых могут быть кратные и нулевые.

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

решения можно обозначить y1(k), y2 (k),

, yn (k) или y1,k , y2,k , , yn,k .

Необходимое и достаточное условие линейной независимости решений

имеет вид

 

 

 

 

 

 

y1,k

y2,k

 

yn,k

 

 

 

 

 

 

Ck

E y1,k

E y2,k

E yn,k

0

 

 

 

 

 

 

En 1 y

En 1 y

2,k

En 1 y

n,k

 

 

1,k

 

 

 

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

Общее решение однородного разностного уравнения записывается в

виде

yk c1 y1,k c2 y2,k

cn yn,k ,

27

где ci - произвольные постоянные, не зависящие от k и, определяемые из

независимых дополнительных условий.

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

Элементы Z - преобразования. Важным элементом теории дискретных систем описываемых дискретными или решетчатыми функциями является раздел теории функций комплексного переменного - Z - преобразование дискретных функций в комплексные функции z - плоскости.

Формально,

прямое

Z -

преобразование

дискретной

функции

(последовательности)

{ fk }

( f0, f1, f2,

) , как оригинала, в функцию

комплексной переменной z

определяется суммой ряда по отрицательным

степеням комплексной переменной z

 

 

 

 

 

 

 

 

 

 

F (z) F

f

k

z k .

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

k

1

 

 

 

 

 

 

 

Обратное Z - преобразование функции комплексной переменной z в

дискретную функцию

целочисленной

переменной

k ~ k T ,

где

T

1-

нормированное

значение

периода

дискретной

последовательности

(функции), определяется интегралом

по

замкнутому контуру

C

на

комплексной плоскости z

 

 

 

 

 

 

 

 

 

 

 

 

 

fk

1

 

zk 1 F (z) dz .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

j C

 

 

 

 

 

 

 

 

При вычислении контурного интеграла используется фундаментальное

соотношение Коши из теории функций комплексного переменного

 

 

 

 

 

zk dz

2

 

j,

при

k

1,

 

 

 

 

 

 

0,

 

 

при

k

1.

 

 

 

 

 

 

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для вычисления значения интеграла по контуру C на комплексной плоскости z используются теория вычетов функций комплексных переменных. Вычеты представляют собой значение комплексной функции в особых точках, и значение интеграла определяется суммой вычетов. Оригиналы функций, заданных Z - изображениями, могут быть определены и по таблицам обратного Z - преобразования.

Согласно теории Z - преобразования, оператору сдвига (упреждения) E в области оригиналов соответствует оператор z в комплексной плоскости

изображений, а оператору сдвига (запаздывания) E 1 в области оригиналов

соответствует оператор z 1 в комплексной плоскости изображений.

Дискретные гармонические функции времени. Дискретным гармоническим функциям времени sin( k) и cos( k) , как оригиналам,

соответствуют изображения в комплексной z - плоскости вида

 

 

 

 

 

 

 

28

 

 

 

 

 

 

 

 

 

 

sin(

k)

 

 

 

z

sin(

T )

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

z

cos(

T )

1

 

 

 

 

cos(

k)

 

 

 

z2

z

 

cos(

T )

 

 

.

 

 

 

 

 

z2

 

 

z

 

 

T )

 

 

 

 

 

 

 

 

2

cos(

1

 

Обычно период дискретизации нормируют к единице T 1. Частота

дискретизации

связана

с

периодом

 

дискретизации

соотношением

d

2 /T , то есть fd 1/T .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Последовательность единичных

 

- импульсов. Последовательности

единичных -

импульсов

1k

с периодом следования T

соответствует

изображение в z - плоскости вида

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1k

 

 

z

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

где k

0,1, 2, .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Одиночный единичный

-

импульс.

Одиночному единичному -

импульсу в нулевой момент времени 10

 

соответствует изображение в z -

плоскости вида

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10

 

1,

 

 

 

 

 

 

 

где k

0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Характеристики дискретных систем. По аналогии с непрерывными системами вводят соответствующие характеристики дискретных или цифровых систем.

Системной или (передаточной) характеристикой дискретной или цифровой системы называется отношение изображения реакции к изображению входного воздействия.

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

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

a

n

En y

k

b

m

Em

x .

k

 

k

 

k

n 0

 

 

 

m 0

 

 

 

Заменяя, операторы сдвига Ei

на операторы zi и оригиналы функций

xk , yk на изображения,

в

предположении

нулевых

начальных условий,

приходим к записи

 

 

 

 

 

 

 

29

a

zn Y

b

zm X

z

.

k n

z

k m

 

 

n 0

 

m 0

 

 

 

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

 

 

 

 

b

m

zm

 

 

Yz

 

 

k

 

 

Sz

 

m

0

 

 

,

X z

 

 

a

n

zn

 

 

 

 

 

 

 

 

k

 

 

 

 

 

n

0

 

 

 

где X z , Yz - изображения входного

воздействия и выходной реакции

дискретной системы.

Знаменатель передаточной функции, приравненный нулю, определяет

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

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

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

Начальные условия. Как и в аналоговом случае, начальные условия могут быть определены по системной функции на основании теоремы дискретного операционного исчисления ( Z - преобразования) о начальном значении функции

f0

lim fk

lim

F (z) .

 

k 0

z

 

Можно определить начальные условия и непосредственно по

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

такими, чтобы последовательно

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

начинается при k

0. Так, если исходное разностное уравнение имеет второй

порядок

 

 

 

 

 

 

yk 2

a1 yk 1

a0 yk

fk ,

то последовательно полагая k

2 и k

1

 

 

y0

a1 y 1

a0 y 2

f 2 ,

 

y1

a1 y0

a0 y 1

f 1,

и учитывая, что

y 2 0 и y 1 0 ,

последовательно определяем начальные

значения y0 и y1 .

Частотная характеристика (ЧХ) дискретной или цифровой системы. Частотная характеристика дискретной или цифровой системы,

30

получается, по системной (передаточной) функции при замене оператора z e j T

 

 

 

 

b

m

e j m T

 

Y ( j )

 

 

k

 

 

S ( j )

 

m 0

 

 

.

 

 

 

 

X ( j )

 

 

ak n e j n T

 

 

 

n

0

 

 

 

Амплитудо - частотная (АЧХ) и фазочастотная (ФЧХ) характеристики дискретной или цифровой системы представляют собой,

как и в аналоговом случае, соответственно, модуль и аргумент комплексной ЧХ.

Отметим, что частотные характеристики дискретных или цифровых систем, в отличие от аналоговых систем, имеют периодически

повторяющийся характер, о чем свидетельствуют

множители типа e j k T .

Частотная характеристика дискретной

или цифровой системы

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

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

отсутствие сторонних источников.

 

 

Выполняя

в выражении частотной характеристики замену

вида

j

ln(z) /T ,

перейдем к системной характеристике дискретной

или

цифровой системы.

 

 

Данная связь системной функции и частотной характеристики часто

используются при синтезе дискретных и цифровых фильтров по аналоговому прототипу. При этом, вместо обратной замены p j ln(z) /T ,

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

p

 

j

 

 

2

 

z

1

,

 

 

 

 

T

 

z

1

 

 

 

 

 

 

 

 

 

 

 

что обеспечивает дробно-рациональный вид системной функции.

Тогда, учитывая, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

e

j

d

T

1

 

j

 

 

 

 

 

 

 

 

,

a

 

 

T

 

 

j

d

T

 

 

 

 

 

e

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

получаем зависимость частотных переменных

 

a и d 2 /T аналоговой

и цифровой систем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

2

 

tg

 

 

d

T

 

.

 

 

T

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из приведенного соотношения видно, что связь частот аналогового прототипа и синтезируемого дискретного или цифрового фильтра нелинейная, поэтому при реализации цифрового фильтра с заданной