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

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

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

111

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 d2k

 

 

 

 

 

 

 

 

 

 

 

d1 d2

d1 d2

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

1

fk

11 b 10

 

 

 

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

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

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

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

системы

 

 

 

 

1

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

v

 

 

 

( d k

d k ) 1

 

( d k n

d k n ) (1 b 1

)

;

 

 

 

 

 

 

 

 

 

 

 

 

k

 

d2

 

d1

1

 

2

 

 

1

 

2

n

n 1

 

 

 

 

 

 

 

 

 

 

n 1

 

 

 

 

 

 

 

 

 

v

 

1

 

 

( d k

d k ) d k

k

1n

b 1n 1

d k

k

1n

b 1n 1

.

 

 

 

 

 

 

 

 

 

k

 

d2

d1

1

2

1

 

 

d1n

 

2

 

 

d2n

 

 

 

 

 

 

 

 

n 1

 

 

n 1

 

 

 

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

112

 

 

 

 

 

d k

d k

 

 

 

b d k

 

 

 

 

 

b d k

 

 

 

gk

vk

1

2

 

 

 

1

 

 

 

 

 

2

 

 

 

 

d2

d1

 

 

d1 (d2

d1)

 

d2 (d2

d1)

 

 

 

 

 

 

 

 

 

 

 

 

 

d k

d k

 

b (d k 1

d k 1)

 

 

 

 

 

 

 

 

 

1

2

 

 

 

 

 

1

2

 

 

,

 

 

 

 

 

 

 

 

d1

d2

 

 

 

d1 d2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при k

1 или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d k 1

d k 1

 

b (d k

d k )

 

 

 

 

g

k 1

 

v

1

 

 

1

 

2

 

 

 

 

1

2

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

d1

d2

 

 

d1

d2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при k

0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Полученные

выражения

совпадают

с

результатами, операторного

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

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

 

V (z)

 

V

 

z2

b

z

b

 

S(z) Sz

 

 

z

 

 

1

0

,

E(z) Ez

 

 

 

 

 

 

(z d1)

(z d2)

 

где Ez - изображение входного

воздействия;

Vz -

изображение выходной

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

Частотная характеристика дискретной системы определяется по

системной функции путем замены z e j

 

 

T

 

 

 

 

 

 

 

V ( )

 

 

e 2 j T

 

 

b e j T

b

 

 

 

S ( )

 

 

 

 

 

 

 

 

 

 

 

1

 

 

0

 

 

 

E( )

 

 

(e

j T

 

d1) (e j T

d2 )

 

 

 

 

 

 

 

e 2 j

 

T

b

 

 

e

j T b

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

0

 

 

,

 

 

e 2 j T

 

(d d

2

) e j T

d d

2

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

1

 

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

 

 

 

 

 

 

 

 

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

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

 

 

 

 

 

 

 

 

 

S( )

 

Abs(S(

)) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

 

 

 

 

 

 

 

 

 

 

 

(

 

) Arg(S(

 

)) 180 / .

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

E

z

 

(z2

b z

b )

 

 

 

 

 

 

 

V

 

 

 

 

 

 

1

0

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

(z

 

d1)

 

(z

d2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

113

(z d ) (z d

2

) z2

(d d

2

) z d d

2

0 ,

1

 

1

 

 

1

 

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

d1; d2

d2 .

 

 

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

дискретной

системы. Приступаем к

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

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

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

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

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

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

z ek 1k Ez z 1 .

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

 

z (z2 b z

b )

 

Vz

1

0

.

(z 1) (z d1) (z d2 )

 

 

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

 

 

Vz

 

 

 

z (z2

 

b1 z b0 )

 

 

 

vk

A B1 e

1 t

B2 e

 

2 t

 

 

 

 

 

 

 

(z 1) (z d1) (z d2 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

B

 

d k B

 

d k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

 

2

 

2

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h v

 

 

 

1 b b

 

(d

2

 

 

b d b ) d k

 

 

 

(d 2

 

b d

2

b ) d k

 

 

 

1

0

 

 

1

 

 

1 1

0

1

 

 

 

 

2

 

1

 

0

2

,

 

k

 

k

 

 

 

(1 d1) (1 d2 )

 

(1 d1) (d2

d1)

 

 

 

 

(1 d2 ) (d1

d2 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 b b

 

 

 

 

d 2

b d b

 

 

 

 

 

 

 

d 2

 

 

b d

 

b

где

 

A

 

 

 

1

0

;

 

B

 

 

 

1

1

1

0

 

;

 

B

 

 

2

 

 

1

2

0

;

 

 

 

(1 d1) (1 d2 )

1

 

(1 d1) (d2

d1)

 

2

 

 

(1 d2 ) (d1

d2 )

 

 

 

 

 

 

 

 

 

 

 

 

ln(d1)

;

 

 

ln(d2 )

 

; k

 

t

;

 

1 t

k

;

 

 

2 t

 

k

;

 

T -

период

1

 

T

 

2

 

T

 

T

e

 

d1

 

e

 

 

 

d2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

0 имеем v0

h0

1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Построение

 

разностного

уравнения

 

 

 

дискретной

 

системы.

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

114

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

zn

 

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

выражения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

z

 

 

z2

b z b

 

 

 

 

 

z2

b z b

 

 

 

 

 

 

 

 

1

0

 

 

 

 

 

 

1

0

 

 

 

 

Ez

 

(z d1) (z d2 )

 

z2

(d d

2

) z d d

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

 

 

 

v

 

 

E2

 

b

 

 

E

b

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

1

 

 

0

 

 

.

 

 

 

 

 

 

 

1k

 

E2

(d d

2

) E d d

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

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

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

vk 2 (d1 d2 ) vk 1

d1 d2 vk

1k 2

 

b1 1k 1

b0 1k fk

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

vk 2

 

(d1

d2 ) vk 1

d1 d2 vk 1k 2

b1 1k 1

 

b0 1k .

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

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

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

z2

 

b

z b

 

 

 

 

 

 

 

 

v

 

lim v

 

lim V

lim

 

 

 

 

 

1

 

0

 

1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

k 0

k

z

 

z

z

 

z 1 (z d ) (z d

2

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

В соответствии с теоремой упреждения, значение

функции vk 1

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z2 (z2

b z

b )

 

 

 

 

 

 

 

 

 

 

 

vk 1

 

z Vz

 

z v0

 

 

 

 

 

1

0

 

z 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(z 1) (z d1) (z d2)

 

 

 

 

 

 

 

 

 

 

 

 

(1 b d d

2

) z3

(b d d

2

d d

2

) z2

 

d d

2

z

 

 

 

 

 

 

1

 

1

 

 

 

0

1

 

 

1

 

 

 

1

 

.

 

 

 

 

 

 

 

 

 

 

 

(z 1) (z d1) (z d2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

v1

lim vk 1

lim (z Vz z

v0 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

0

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 b d d

2

) z3

(b d d

2

d d

2

) z2

 

d d

2

z

 

 

lim

 

 

1

 

1

 

 

 

0

1

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

(z 1) (z d1) (z d2 )

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

b1

 

d1

d2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

 

исходным

разностным

 

уравнением,

 

 

полагая

115

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

k

2 и k

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

 

 

 

 

 

 

v0 (d1

d2 ) v 1 d1 d2 v 2

10 b1 1 1 b0 1 2

 

 

 

(d1

d2 ) 0 d1 d2 0 1 b1 0 b0 0 1;

 

 

 

v1

(d1

d2 ) v0 d1 d2 v 1

11

b1 10

b0 1 1

 

 

 

 

(d1

d2 ) 1 d1 d2 0 1 b1 1 b0 0 1 b1

d1

d2 .

 

Таким

образом, получаем, что

начальные

значения

равны v0 1 и

v1 1 b1 d1 d2 .

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

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

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

vk 2 (d1 d2 ) vk 1 d1 d2 vk 1k 2 b1 1k 1 b0 1k fk

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

 

 

 

 

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

b 1

 

b 1

f

k

.

1,k 1

2,k

2

k 2

1 k 1

0 k

 

 

Напомним, что определяющая система уравнений Лагранжа образуется при подстановке предполагаемого общего решения в исходное разностное уравнение и наложении ограничения на сдвиг функций 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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

116

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d2k

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c1,k

 

 

fk

d2k 2

 

 

 

 

(1k 2

b1 1k 1

b0 1k ) d2k 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d k

1

d k

 

1

(d

2

d )

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1k 2

b1 1k 1

b0 1k )

;

 

 

 

 

 

 

 

 

 

 

 

 

 

d k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

(d

2

 

d )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

d k

1

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c2,k

 

d1k 2

fk

 

 

 

(1k 2

b1 1k 1

b0 1k ) d1k 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d k

1

d k

1

(d

2

d )

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

 

 

1

 

 

 

 

 

 

 

1k 2

b1 1k 1 b0 1k

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d k

1 (d

2

d )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

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

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

 

c1,k

1

c1,k

 

 

 

k 1n 1

b1 1n

b0 1n 1

 

 

 

 

 

 

 

 

 

 

 

 

n 1

 

 

 

d1n (d2

d1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 b1

b0 )

 

 

k

 

1

 

 

 

(1 b1

 

 

 

k

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b0 ) (d1

 

 

c1 ;

 

 

 

 

 

d2

d1

 

n 1 d1n

(d2

d1) d1k (d1

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c2,k

1

c2,k

 

 

 

 

k 1n 1

b1 1n

b0 1n 1

 

 

 

 

 

 

 

 

 

 

 

 

n 1

 

 

 

 

d2n (d2

d1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 b1

b0

 

k

 

 

 

1

 

 

 

(1 b1

 

 

k

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b0 ) (d2

 

 

c2

,

 

 

 

 

 

d 2 d1 n 1 d2n

 

 

(d 2 d1) d2k (d2

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

и c2,k

в предполагаемое общее

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

 

 

 

 

 

 

 

 

 

 

 

 

(1 b1

b0) (d1k

 

1)

 

 

 

 

 

 

 

 

k

 

(1 b1

b0) (d2k

1)

 

 

k

vk

 

 

 

 

 

 

 

 

 

 

c1 d1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c2

d2 ;

 

(d 2

d1) (d1

1)

 

 

 

 

(d 2 d1) (d2

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 b1

b0) (d1k

 

1) (1 b1

b0) (d2k

1)

 

 

 

 

 

k

 

k

vk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c1

 

d1

c2

d2 .

 

(d 2

d1) (d1

1)

 

 

 

 

 

 

 

(d 2

d1) (d2

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

и

 

c2

 

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

найденными ранее, начальными условиями v0

1 и v1

1

 

b1

d1

d2 . Так,

приравнивая общее решение, при k

0

и k

1,

 

начальным

условиям,

находим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v0

 

1

 

 

 

 

 

0 0 c1

c2 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

117

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

b1

 

b0

1

 

b1

b0

 

 

 

 

 

 

v1

1 b1

d1

d2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c1

d1

c2 d2

 

 

d2

d1

 

 

 

d2

d1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c1

c2

1 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d1 c1 d2 c2

1 b1

d1

d2 Q .

 

 

 

Постоянные суммирования c1

и

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

используя правило Крамера

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

 

 

d2

 

 

d1;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d1

d2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

Q d2

 

(1 b1

d1) ; c

 

 

d1

Q

1 b1

d2 .

1

 

 

 

 

d2

d1

 

 

 

2

 

 

 

 

 

 

 

d2

d1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

(1 b b ) (d k

 

1) (1 b b ) (d k

 

 

1)

 

 

 

 

 

 

 

 

 

 

vk

 

 

 

1

 

0

 

1

 

 

 

 

 

 

 

 

1

 

0

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(d 2 d1) (d1

1)

 

 

 

 

(d 2 d1) (d2

 

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 b d ) d k

 

(1 b d

2

) d k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

1

 

 

 

 

 

1

 

 

2

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d 2

d1

 

 

 

 

 

 

d 2

d1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 b b ) (d k

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

 

 

 

 

 

 

 

 

vk

 

1

0

 

 

1

 

 

 

 

 

1

 

 

 

 

1

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(d 2

d1) (d1

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 b b ) (d k

1) (d

2

1) (1 b d

2

) d k

 

 

 

 

 

 

 

 

 

 

 

1

0

 

2

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

;

 

 

 

 

 

 

 

 

 

 

 

 

 

(d 2 d1)

 

 

(d2

 

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 b b ) (d 2

b d b ) d k

 

 

(1 b b ) (d 2

b d

2

b ) d k

vk

1

0

 

1

1

 

1

 

0

 

1

 

 

 

 

 

 

1

0

 

 

 

2

 

1

 

0

2

 

 

(d 2

 

d1) (d1

1)

 

 

 

 

 

 

 

 

 

 

 

 

(d 2 d1) (d2

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h v

 

1 b b

 

 

 

(d 2

b d b ) d k

 

 

 

(d 2

 

b d

 

b ) d k

 

 

 

1 0

 

 

 

1

 

1 1

 

 

0

 

 

1

 

 

2

 

 

 

1

 

2

 

 

0

2 .

k

k

 

(1 d1) (1 d2 )

 

 

 

(1 d1) (d2 d1)

 

 

 

(1 d2 ) (d1

d2 )

 

 

 

 

 

 

 

 

 

 

 

 

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

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

vk 2 (d1 d2 ) vk 1 d1 d2 vk 1k 2 b1 1k 1 b0 1k fk

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

118

 

Так,

вводя

 

новые

 

переменные

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

1k 2

 

b1 1k 1

b0 1k

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X k

1

 

A X k

Fk .

 

 

 

 

 

Согласно методу Коши, частное решение неоднородной системы

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

 

 

 

 

 

 

 

 

Ak

 

 

k

n F

 

 

 

 

 

 

 

 

 

X

k

X

0

Ak

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 .

119

Доказывается, что модальная матрица 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

1

; Fk

0

0

.

x2,0

v1

1 b1 d1 d2

fk

1k 2 b1 1k 1 b0 1k

 

 

 

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

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

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

и

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

системы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(d

2

d k

d d k ) 1 ( d k

 

d k ) (1 b d d

2

)

 

 

1

 

 

1

1

2

 

1

 

2

 

1

1

 

 

 

vk

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

( d k n

d k n ) (1

 

 

 

 

 

 

 

d2

d1

 

b 1

b 1

 

)

 

 

 

 

 

 

 

 

 

1

2

 

n 1

1 n

0 n 1

 

 

 

 

 

 

 

 

 

n 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 b ) ( d k

d k ) d d k

d

2

d k

 

 

 

 

 

 

 

 

1

 

 

 

1

1

 

2

1

 

1

 

2

 

 

 

 

 

vk

 

 

 

 

 

 

 

 

k

1

 

 

 

 

 

 

k

1 .

 

 

 

 

 

 

 

b0 ) d1k

 

 

 

 

b0 ) d2k

 

 

d2

d1

(1 b1

(1 b1

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1 d1n

n 1 d2n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

120

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

 

 

 

 

 

 

 

 

(1 b d ) d k

(1 b d

2

) d k

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

1

1

1

 

 

 

 

1

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

vk

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

(1 b1

 

 

 

1)

 

 

(1 b1

 

 

 

 

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d2

d1

 

b0 ) (d1

 

 

b0 ) (d2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d1

1

 

 

 

 

 

 

 

 

 

 

d2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 b b ) (d k

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

 

 

 

 

 

 

 

 

 

vk

 

1

 

 

0

 

1

 

 

1

 

 

 

 

 

1

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(d 2

d1)

(d1

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 b b ) (d k

1) (d

2

 

1) (1 b d

2

) d k

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

0

 

2

 

 

 

 

 

 

 

 

1

 

2

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(d 2 d1)

 

(d2

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 b b ) (d 2

b d b ) d k

 

(1 b b ) (d 2

b d

2

b ) d k

vk

1

0

 

1

1

1

 

0

1

 

 

 

 

1

 

0

 

 

2

 

1

 

0

2

 

 

(d 2

 

d1) (d1

1)

 

 

 

 

 

 

 

 

 

(d 2 d1) (d2

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h v

 

1 b b

 

 

 

 

(d 2

b d b ) d k

 

 

(d 2

b d

 

b ) d k

 

 

 

1 0

 

 

 

 

1

1 1

 

0

1

 

 

 

 

2

 

1

2

 

 

 

0

2 .

k

k

(1 d1) (1 d2 )

 

 

(1 d1) (d2 d1)

 

 

 

 

(1 d2 ) (d1

d2 )

 

 

 

 

 

 

 

 

 

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

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

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

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

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

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

 

 

gk

hk

 

hk

hk 1

hk 1.

 

 

 

 

 

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

 

 

 

 

h

1 b b

(d 2

b d b ) d k

 

(d 2

b d

2

b ) d k

1

0

 

1

1

1

0

1

 

2

1

0

2

,

k

(1 d1) (1 d2 )

(1 d1) (d2

d1)

 

(1 d2 ) (d1

d2)

 

 

 

 

 

 

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