Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

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

..pdf
Скачиваний:
10
Добавлен:
05.02.2023
Размер:
3.33 Mб
Скачать

71

X0

x1,0

v0

0

; Fk

x2,0

v1

0

 

 

Кроме того, отметим, что нас интересует лишь решений x1,k vk .

0 0 . fk 1k

первая компонента вектора

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

структуры векторов X 0

и Fk ,

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

системы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

( d k n

d k n )

 

 

d k

k

1

 

d k

 

k 1

v 0

1

2

 

1

 

1

 

n 1

 

2

 

 

n 1

.

 

 

 

 

 

 

 

 

 

 

k

d2

d1

n 1

d2

d1 n 1

d1n

 

d2

d1 n 1 d2n

n 1

 

 

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

 

 

 

 

 

d k

 

 

d k

1

 

 

d k

 

 

 

d k

 

1

 

 

 

 

vk

 

 

 

1

 

 

1

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

d2

d1

 

d k

 

(d 1)

 

d2

d1

 

d k

 

(d

2

1)

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

d k

1

 

 

 

 

 

 

d k

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

2

 

 

 

 

 

;

 

 

 

 

 

 

 

 

(d2

d1) (d1

1)

 

(d2

d1) (d2

1)

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

d k

 

 

 

 

 

 

 

d k

 

hk

vk

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

2

.

 

(1 d1) (1 d2 )

 

(d 2

d1) (1 d1)

 

(d 1 d2) (1 d2)

 

 

 

 

 

 

Полученное выражение совпадает с результатами, операторного метода и метода Лагранжа и описывает переходную характеристику, исследуемой дискретной системы.

Импульсная характеристика дискретной системы. Приступаем к определению импульсной характеристики дискретной системы различными методами.

Как известно, импульсная характеристика представляет собой реакцию дискретной системы, находящейся в исходном состоянии покоя, на входной одиночный единичный - импульс, при k 0 , то есть 10 . Под исходным

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

Отметим, что импульсная характеристика дискретных и цифровых систем определена при k 1.

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

gk hk hk hk 1 hk 1.

Учитывая, что переходная характеристика имеет вид

72

 

 

 

1

 

 

 

 

 

d k

 

 

 

 

 

d k

 

 

 

hk

vk

 

 

 

 

 

 

 

1

 

 

 

 

2

,

 

 

(1 d1) (1 d2 )

 

(d 2 d1) (1 d1)

 

(d 1 d2) (1 d2)

 

 

 

 

 

 

 

 

запишем выражение соответствующей функции, отстающей на один такт

 

 

 

1

 

 

 

 

 

d k

1

 

 

 

d k 1

 

 

 

hk 1

vk 1

 

 

 

 

 

 

 

1

 

 

 

2

 

.

 

 

(1 d1) (1 d2 )

 

 

(d 2 d1) (1 d1)

 

 

(d 1 d2 ) (1 d2 )

 

 

 

 

 

 

 

 

 

 

Применяя уравнение связи, сразу получаем импульсную

характеристику, исследуемой дискретной системы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d k 1

 

d k 1

 

 

 

 

 

 

 

 

 

gk hk

hk

1

1

2

 

,

 

 

 

 

 

 

 

d1

 

d2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при k

1, или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d k

d k

 

 

 

 

 

 

 

 

 

gk

1

 

1

2

 

,

 

 

 

 

 

 

 

 

 

 

d1

d2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при k

0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, по известной переходной характеристике дискретной

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

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

Оригиналу входного воздействия 10 , согласно теории Z -

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

ek 10 Ez 1.

Изображение выходной реакции дискретной системы будет иметь вид

1

Vz (z d1) (z d2 ) .

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

Всоответствии с теоремой о начальном значении функции

 

v0

lim vk

lim Vz ,

 

 

 

k

0

 

z

 

находим

 

 

 

 

 

 

v0

lim Vz

lim

 

1

0 .

 

 

(z d1) (z d2 )

 

z

z

 

 

Далее, используя теорему Z - преобразования об упреждении функции на один такт

v1 z Vz z v0 ,

находим

 

 

 

 

 

 

 

 

 

 

73

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v1

 

 

 

 

 

z

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d1) (z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(z

d2 )

 

 

 

 

 

 

Теперь, используя таблицы обратного

Z -

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

находим

оригинал выходной реакции

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

 

e

 

1 t

 

e

2 t

 

 

 

k

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d1

 

 

 

 

d2

,

 

 

 

 

 

 

k

1

 

 

 

d1

 

d2

 

 

 

 

 

d1

 

 

 

 

d2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при k

0 или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d k 1

 

d k

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

vk

1

 

 

 

 

2

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d1

 

d2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln(d1)

 

 

 

 

 

ln(d2 )

 

 

 

 

 

 

t

 

 

 

 

 

 

1 t

k

2 t

k

при k

1, где

1

T

; 2

 

 

 

T

 

;

k

 

 

 

T

;

e

 

d1 ; e

 

d2 ; T -

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

период входной последовательности.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Полученные выражения описывают импульсную характеристику,

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d k

 

d k

 

 

 

 

 

 

 

 

 

g

k 1

 

v

1

 

 

1

 

 

 

2

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

d1

 

d2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при k

0 или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d k

 

1

 

 

d k

1

 

 

 

 

 

 

 

 

 

 

 

 

gk

vk

1

 

 

 

 

2

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

d1

 

d2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при k

1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отметим, что при k

0 имеем v1

 

g1

0 .

 

 

 

 

 

 

 

 

 

 

 

Построение разностного

 

уравнения

дискретной

системы.

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

системной функции путем

замены

 

 

переменной

z

на k , в левой части

выражения и заменой оператора zn

на оператор сдвига En , в правой части

выражения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Vz

 

 

1

 

 

 

 

 

 

1

 

 

 

 

 

 

Ez (z d1) (z d2)

 

z2

(d d

2

) z d d

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

vk

 

 

 

 

 

1

 

 

 

.

 

 

 

 

10

 

 

E2

(d d

2

) E d d

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

Преобразуя выражение, получаем неоднородное разностное уравнение

второго порядка

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

vk 2

(d1

d2 ) vk 1

d1 d2 vk

10

fk

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

vk 2

 

(d1 d2 ) vk 1

d1 d2 vk

 

10 .

 

 

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

74

Определение начальных условий. Для однозначного определения решения разностного уравнения необходимы дополнительные независимые условия, в качестве которых удобно воспользоваться начальными условиями. Так как исходное разностное уравнение второго порядка и импульсная характеристика определена при k 1, необходимо определить v0 , v1 и v2 .

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

 

v

lim v

lim V

 

lim

 

 

1

 

 

 

 

0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

k 0

k

z

 

z

z

 

(z d ) (z d

2

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

В соответствии

с

теоремой

упреждения,

значение

функции vk 1

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

vk

 

z Vz

 

z

v0

 

 

z

 

 

 

.

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(z

d1) (z

d2 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Применяя повторно теорему о начальном значении функции, получаем

v

lim v

lim (z V

 

z

v )

 

lim

 

z

 

 

 

 

0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

k 0

k 1

z

 

z

 

0

 

z

 

(z d ) (z d

2

)

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

Применяя еще раз теорему упреждения на один такт к последнему результату

vk 2

z Vz

z v1

z2

(z d1) (z d2)

 

 

 

и теорему о начальном значении функции, получим

v

lim v

 

lim (z V

z v )

lim

z2

 

 

 

1.

 

 

 

 

 

2

k 0 k 2

z

z

1

z

(z d ) (z d

2

)

 

 

 

 

 

 

 

 

 

1

 

 

 

С другой стороны, для определения начальных условий можно

воспользоваться

исходным

 

разностным

уравнением,

полагая

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

k

2 , k

1 и k 0, последовательно получаем

 

v0

(d1

d2 ) v 1

d1 d2 v 2

1 2

(d1

d2 ) 0 d1 d2 0 0 0 ;

 

v1

(d1

d2 ) v0

d1 d2 v 1

1 1

(d1

d2 ) 0 d1 d2 0 0 0 ;

 

v2

(d1

d2 ) v1

d1 d2 v0

10

(d1

d2 ) 0 d1 d2 0 1 1 .

Таким образом, получаем, что начальные значения равны v0 0 , v1 0

и v2 1.

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

Метод вариации произвольных постоянных (метод Лагранжа).

Согласно методу Лагранжа, общее решение неоднородного разностного уравнения второго порядка

vk 2 (d1 d2 ) vk 1 d1 d2 vk 10 fk

75

следует искать в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

c

v

c

v

c

d k

c

d k ,

 

 

 

 

 

 

 

k

1,k

1,k

2,k

2,k

1,k

1

2,k

2

 

 

 

где

d ,

d

2

- корни

характеристического

уравнения;

y

d k ,

y

d k -

 

1

 

 

 

 

 

 

 

 

 

1,k

1

2,k

2

фундаментальная система решений соответствующего однородного уравнения; c1,k , c2,k - варьируемые постоянные – неизвестные пока функции.

Варьируемые постоянные находятся из определяющей системы

уравнений Лагранжа

 

 

 

 

 

 

 

 

c

d k 1

c

 

d k

1

0 ;

 

 

1,k 1

2,k

2

 

 

 

 

c

d k 2

c

d k 2

1

f

k

.

1,k

1

2,k

 

2

0

 

 

Напомним, что определяющая система уравнений Лагранжа образуется при подстановке предполагаемого общего решения в исходное разностное уравнение и наложении ограничения на сдвиг функций c1,k и c2,k . Первое

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

Выразим разности

варьируемых

постоянных

 

c1,k

и

 

c2,k из

определяющей системы уравнений, используя правило Крамера

 

 

 

 

 

 

 

d k

1

 

d k

1

 

d k 1

d k 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

1

 

2

 

 

(d

2

d ) ;

 

 

 

 

 

 

 

d k

2

 

d k

2

 

 

1

 

 

2

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d2k

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c1,k

 

fk

d2k 2

 

 

 

 

10 d2k 1

 

 

 

 

 

 

 

 

10

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d k 1

d k 1

 

(d

2

d )

 

d k 1

(d

2

d )

 

 

 

 

 

 

 

1

 

2

 

 

1

 

 

 

 

 

1

 

 

1

 

 

 

 

d k

1

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c2,k

d1k 2

fk

 

 

 

10 d1k 1

 

 

 

 

 

 

 

 

10

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d k 1

d k 1

 

(d

d )

 

 

d k 1

(d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

d )

 

 

 

 

 

 

 

1

 

2

 

 

 

2

1

 

 

 

 

 

2

 

 

 

1

 

 

Для определения

варьируемых постоянных

 

 

c1,k

 

и c2,k

применим

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

1

k

1n 1

 

1

k 1n 1

1

 

 

c1,k

c1,k

 

 

 

 

 

 

 

 

c1

;

 

 

 

d2 d1 n 1 d1n

 

 

 

n 1 d1n (d2

d1)

 

 

d1 (d2 d1)

 

 

 

 

 

 

 

 

 

 

76

 

 

 

 

 

 

 

 

 

 

c2,k

1

c2,k

k

 

1n 1

 

 

 

 

1

 

 

k 1n 1

1

c2 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1 d2n (d2

d1) d2

d1 n 1 d2n

 

d2 (d2 d1)

 

 

 

 

 

где c1, c2 - новые постоянные суммирования.

 

 

 

 

 

 

 

Подставляя найденные значения c1,k

и c2,k в предполагаемое общее

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

 

 

 

 

 

 

d1k

 

 

k

 

 

d2k

 

 

 

 

k

 

 

 

vk

 

 

 

 

c1

d1

 

 

 

 

 

 

 

c2 d2 ;

 

 

 

d1 (d2

d1)

 

d2

(d2

d1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d1k 1

d2k

1

 

 

 

k

 

k

 

 

 

 

 

 

vk

 

 

 

 

 

 

c1 d1

c2 d2 .

 

 

 

 

 

 

d2

d1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для определения постоянных суммирования c1

и c2 воспользуемся,

найденными ранее,

начальными условиями v1

0 и v2

1, так как решение

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

k 1. Так, приравнивая общее

решение,

при k

1 и

k 2 , начальным

условиям, находим

 

 

 

 

 

 

 

 

v1

0 0 c1 d1

c2 d2 ;

 

 

v

1

1

c

d 2

c

d 2

 

 

2

 

 

1

1

2

2

 

 

или

 

 

 

 

 

 

 

 

 

d1

c1

d2

c2

0 ;

 

 

 

 

d 2

c

d 2 c

0 .

 

 

 

 

1

1

2

2

 

 

 

 

Из полученной системы сразу следует, что c1

c2

0 .

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

 

 

 

 

d k 1

d k 1

 

gk

vk

1

 

 

2

,

 

 

d2

d1

 

 

 

 

 

 

 

при k

1 или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d k

d k

 

g

k

v

1

 

1

2

 

,

 

 

 

 

 

 

1 k

 

d2

d1

 

 

 

 

 

 

 

 

 

 

при k

0.

 

 

 

 

 

 

 

 

 

Полученное решение

описывает

импульсную характеристику,

исследуемой дискретной системы, и, как видим, совпадает с выражениями, найденными операторным методом.

Решение в форме Коши (метод Коши). Рассматриваемый нами вариант метода Коши предполагает предварительное преобразование исходного неоднородного разностного уравнения

vk 2 (d1 d2 ) vk 1 d1 d2 vk 10 fk

в эквивалентную систему двух разностных уравнений первого порядка.

77

Так,

вводя новые

переменные

x1,k vk ;

x2,k x1,k 1 vk 1;

x3,k v2,k 1

vk 2 , получаем эквивалентную систему разностных уравнений

первого порядка

 

 

 

 

 

x1,k 1

0

1

x1,k

0

 

x2,k 1

d1 d2

d1 d2

x2,k

10

или

X k 1 A X k Fk .

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

 

 

 

 

 

Ak X

 

k

 

 

 

 

 

X

k

0

Ak n F

1

,

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

n 1

 

 

где X0

x1,0 x2,0

t

v0

v1

t - вектор начальных условий; Ak - степенная

 

функция от матрицы коэффициентов системы.

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

A H H 1,

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

F ( A) H F ( ) H 1,

где - диагональная матрица собственных значений; H - модальная матрица собственных векторов; F ( ) k диагональная матрица указанной функции

от каждого собственного значения.

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

det([ A

])

 

 

1

2

(d1

d2 )

d1 d2 0 .

 

 

 

d1

d2

d1 d2

 

 

 

 

 

 

 

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

1

0

d1

0

.

0

 

0

d2

2

 

Собственные вектора, как столбцы

модальной матрицы H , по

известным собственным значениям матрицы

A, определяются из решения

однородных систем уравнений

 

 

 

 

A

i

hi 0 ,

 

 

где i - диагональная матрица, составленная из i .

78

Доказывается, что модальная матрица H может быть определена алгебраическими дополнениями элементов одной из строк, например первой,

матрицы [ A

 

i ]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H

11(

1)

 

11( 2 )

 

 

d1

d2

d1

d1

d2

d2

 

 

 

d2

 

d1

.

12 ( 1)

 

12 ( 2 )

 

 

d1 d2

 

d1 d2

 

 

d1 d2

 

d1 d2

 

 

 

 

 

 

 

 

 

Определитель модальной матрицы равен

H d1 d2

(d2

d1) .

 

Используя

определитель, выразим матрицу обратную модальной

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

d

d

2

d

 

 

1

 

1

1/ d

 

 

 

 

H 1

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d1 d2 (d2

 

d1)

d1 d2

d2

 

d2

d1

 

 

1 1/ d1

 

 

 

В результате, получаем выражение степенной функции от матрицы

коэффициентов

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

2

d

 

d k

0

 

1

 

 

 

1

1/ d

2

 

Ak H

k H

1

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

d1 d2

d1 d2

0 d2k

d2

 

d1

1 1/ d1

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

d2 d1k

 

d1 d2k

 

 

d1k

d2k

 

.

 

 

 

 

 

 

 

 

 

d2

d1 d1 d2 (d1k

d2k )

d1 d1k

d2 d2k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Заметим, что структура матрицы Ak n аналогична и отличается лишь показателем степени.

Теперь все подготовлено для представления решения в форме Коши. Предварительно отметим, что в нашем случае вектор начальных условий X 0

и вектор правой части эквивалентной системы разностных уравнений Fk имеют вид

X0

x1,0

v0

0

; Fk

0

0

.

x2,0

v1

0

fk

10

 

 

 

Кроме того, отметим, что нас интересует лишь первая компонента вектора решений x1,k vk .

Обратим внимание на тот факт, что, несмотря на определение импульсной характеристики при k 1, в отличие от метода Лагранжа, где в качестве начальных условий используются значения v1 и v2 , в методе Коши

вкачестве начальных условий используются значения v0 и v1.

Всвязи с отмеченными обстоятельствами, используя общее выражение для решения системы разностных уравнений первого порядка, и, учитывая структуры векторов X 0 и Fk , выразим выходное напряжение дискретной

системы

k

d k n

d k n

 

 

d k

k

1

 

d k

 

k 1

v 0

1

2

1

 

1

 

n 1

 

2

 

 

n 1

.

 

 

 

 

 

 

 

 

 

k

d2

d1

n 1

d2

d1 n 1

d1n

 

d2

d1 n 1 d2n

n 1

 

 

Используем для раскрытия сумм формулу геометрической прогрессии. Учитывая тот факт, что входное воздействие в данном случае существует

79

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

 

 

 

 

d k

 

 

d k

 

 

 

d k 1

d k 1

 

 

gk

vk

 

 

1

 

 

2

 

 

 

1

2

,

 

(d2

d1) d1

(d2 d1) d2

d2

d1

 

 

 

 

при k

1 или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d k

d k

 

 

 

 

 

 

 

g

k 1

v

1

1

2

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

d2

d1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при k

0.

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

 

V (z)

 

V

 

z2 b z

b

 

S(z) Sz

 

 

z

 

1

0

,

E(z)

Ez

(z 1) (z d )

 

 

 

 

где Ez - изображение входного воздействия; Vz - изображение выходной

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

Частотная характеристика дискретной системы определяется по системной функции путем замены z e j T

 

V ( )

 

e 2 j T

b e j T

b

 

e 2 j T

b e j T

b

 

 

S( )

 

 

 

1

0

 

 

1

0

 

,

E( )

 

(e j T

1) (e j T

 

 

e 2 j T

(1 d ) e j T

 

 

 

 

d )

 

 

d

где T - период дискретизации по времени.

Амплитудно-частотная характеристика системы соответствует модулю комплексной частотной характеристики

 

S(

 

)

Abs(S(

)) .

 

Фазочастотная характеристика системы соответствует аргументу

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

 

 

 

 

( )

 

Arg(S(

)) 180 / .

Изображение выходной реакции запишется

 

 

 

E

z

 

(z2

b

z

b )

 

Vz

 

 

 

1

 

0

.

 

 

(z 1)

(z

d )

 

 

 

 

 

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

определяет характеристическое уравнение

 

 

 

(z 1) (z d ) z2

(1 d ) z d 0 ,

корни которого, соответственно, равны d1

1;

d2 d .

80

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

Как известно, переходная характеристика представляет собой реакцию дискретной системы, находящейся в исходном состоянии покоя, на входную последовательность 1k (единичных - импульсов при k 0 и периодом T ).

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

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

Оригиналу входного воздействия 1k , согласно теории Z -

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

z ek 1k Ez z 1 .

Изображение выходной реакции дискретной системы будет иметь вид

 

z

(z2 b

 

z

b )

 

V

 

1

 

0

.

 

 

 

 

z

 

(z 1)2

(z

d )

 

 

 

 

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

 

 

 

Vz

 

z (z2

b1 z b0)

 

vk

 

 

A A0 t B e

 

t

A A0

t B d

k

 

 

 

 

 

(z

1)2

(z

d )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 2 d d b b k (1 b b ) (d 2

d b b ) d k

 

 

 

 

 

 

hk

 

vk

 

 

 

 

 

 

1

0

 

 

 

1

0

 

 

 

 

1

0

 

,

 

 

 

 

 

 

 

 

 

(1 d )2

 

 

 

 

 

(1 d )

 

 

 

 

(1 d )2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2 d

d

b

b

 

 

1

b

b

 

 

d 2 d

b

b

 

 

 

 

ln(d )

 

где

 

A

 

 

 

 

 

1

0

;

 

A

 

 

1

0

;

B

 

 

1

0

 

;

 

 

 

;

 

 

 

 

 

 

(1 d )2

 

 

 

0

 

T (1 d )

 

 

 

(1 d )2

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

t

 

; e

 

t

d k ; T - период входной последовательности.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отметим, что при k

0 имеем v0

h0

1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Построение

разностного

уравнения

 

дискретной

системы.

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

выражения и заменой оператора zn на оператор сдвига En , в правой части выражения