Анализ временных характеристик дискретных и цифровых устройств
..pdf
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 , в правой части выражения
