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

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

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

61

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c2,k

1

 

 

c2,k

 

k 1n b 1n 1

 

 

 

 

1

 

k 1n

 

b 1n 1

 

 

 

b

c2 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1 d n (d 1)

 

 

 

d 1 n 1

 

d n

 

 

 

d (d 1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

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

и c2,k

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

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

 

 

 

 

 

 

 

vk

 

 

b 1k

 

c1

1k

b d k

 

 

 

 

c2

d k

 

 

 

 

 

 

 

 

 

 

 

 

d

1

 

d (d

 

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b b d k 1

c

c

 

d k

 

b (d k 1

1)

 

c

 

c

 

d k .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

1

 

d

1

 

1

2

 

 

 

 

 

d

1

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

и

c2

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

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

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

1

и

v2

1

b d , так как

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

 

 

 

 

 

v1

1

0

c1

 

 

c2

d ;

 

 

 

 

 

 

 

 

 

 

 

 

 

v

1

b

d

b

 

 

c

c

d 2

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c1

d

c2

1;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

d 2

c

 

1

d .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

и

c2 ,

используя

правило Крамера для решения линейной системы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

d

 

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

d (d 1);

c1

 

 

1

d

d

2

 

 

1

 

; c2

 

 

1

1

d

 

1

.

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

d 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d 1

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

 

g

 

v

b (d k 1

1)

 

1

 

d k

1 d k

 

 

b (1 d k 1)

,

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

d

1

 

 

d 1

d 1

1

d

1

d

 

 

 

 

 

 

 

 

при k

1 или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

gk

 

vk

 

 

b (1

d k )

1

d k

1

,

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

d

1

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при k

0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

описывает

импульсную

характеристику,

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

62

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

 

vk 2

(1 d ) vk 1

d vk

11

b 10

fk

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

Так,

вводя новые

переменные

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

d 1 d

x2,k

11 b 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

(1 d )

d 0 .

 

 

d

1 d

 

 

 

 

 

 

 

 

 

 

 

 

 

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

1

0

1

0

 

 

d .

0

2

0

 

 

 

 

63

 

 

Собственные вектора, как столбцы модальной

матрицы H , по

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

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

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

 

 

 

A

i hi 0 ,

 

 

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

i .

 

Доказывается, что модальная

матрица

H может

быть определена

алгебраическими дополнениями элементов одной из строк, например первой,

матрицы [ A

i ]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H

11( 1)

11( 2 )

 

1 d 1 1 d d

d 1

12 ( 1)

12 ( 2 )

 

d

 

 

 

 

 

 

d

d d .

Определитель

 

модальной

матрицы равен

 

 

 

H

d (d

1) . Используя

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

H

1

 

1

 

 

 

d

1

 

 

1

 

 

 

1

 

1/ d

.

 

 

 

 

 

 

d

d

 

 

 

 

 

 

1

1

 

 

 

d (d

1)

 

 

 

d

1

 

 

 

 

 

 

 

 

 

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

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ak

H

 

k

H 1

d

 

1

1k

0

 

 

 

1

 

1

1/ d

 

 

d

 

d

0

d

k

 

d

1

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

d 1k 1 d k

1k

 

d k

 

 

 

.

 

 

 

 

 

d

1

k

d

d

k

k

d d

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d 1

 

1

 

 

 

 

 

 

 

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

системы

 

 

 

 

 

 

 

 

 

64

 

 

 

 

 

1

 

 

 

 

 

 

k

 

 

 

 

vk

 

 

( 1k

d k ) 1

( 1k n d k n ) (1n

b 1n 1)

 

 

 

 

d

1

 

 

 

 

n 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

(d k

1

k

(1

b 1

) d k k

1n

b 1n 1

.

 

 

 

 

 

 

 

 

 

 

d

1

 

 

n

n 1

 

 

d n

 

 

 

n 1

 

 

n 1

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

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

 

g

 

v

 

1

 

 

(d k

1) b b d k 1

1 d k

 

b (1 d k 1)

,

 

k

 

 

 

 

 

 

 

 

 

 

 

k

d

1

 

 

 

 

 

 

 

 

1

d

1

d

 

 

 

 

 

 

 

 

 

 

 

 

при k

1 или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g

 

 

v

 

1

d k 1

b (1

d k )

,

 

 

 

 

 

 

 

 

k

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

1

d

1

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при k

0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Полученные

 

выражения

 

совпадают

с результатами,

операторного

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

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

S (z) Sz

V (z)

 

Vz

 

1

,

E(z) Ez

 

(z d1) (z d2 )

 

 

 

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

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

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

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

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

 

T

 

 

 

 

S( )

V ( )

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E( ) (e

j T

d1) (e j T

d2 )

 

 

 

 

 

 

 

 

 

1

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e 2 j T

(d d

2

) e j T

d d

2

 

 

 

 

 

 

 

 

 

1

 

1

 

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

 

 

 

 

 

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

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

 

 

 

 

 

 

 

 

S( )

 

Abs(S( )) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

65

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

(

) Arg(S( )) 180 / .

 

 

 

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

 

 

 

Vz

 

 

Ez

 

 

 

.

 

 

 

 

 

d1) (z

 

 

 

 

 

 

 

 

 

(z

 

d2 )

 

 

 

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

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

 

 

 

 

 

 

 

(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

Vz (z 1) (z d1) (z d2 ) .

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

Vz

 

 

z

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

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

d k

 

 

 

 

d k

 

 

hk

vk

 

 

 

 

1

 

 

 

 

2

 

,

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

d1) (1 d2 ) (d1 d2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

66

 

 

 

 

 

 

 

 

где

 

A

 

 

1

 

;

 

B1

 

 

 

1

 

 

; B2

 

1

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 d1) (d2

d1)

(1 d2 ) (d1 d2 )

 

 

 

(1 d1) (1 d2 )

 

 

 

 

 

 

 

ln(d1)

 

 

ln(d2 )

 

 

t

 

 

1 t

k

 

 

2 t

k

 

 

1

 

T

;

2

T

;

k

T

;

e

 

d1

;

e

 

d2 ;

T - период

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

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

0 имеем v0

 

h0

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

 

 

 

.

 

 

 

 

 

 

 

1k

 

E2

(d d

2

) E d d

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

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

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

vk 2

(d1

d2 ) vk 1

d1 d2 vk

1k

fk

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

vk 2

 

(d1 d2 ) vk 1

d1 d2 vk

 

1k .

 

 

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

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

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

 

v

lim v

lim V

 

lim

 

z

 

 

 

 

 

1

 

 

 

0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

k 0 k

z

 

z

z

 

z 1 (z d ) (z d

2

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

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

значение

функции vk 1

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

vk

 

z Vz

 

z

v0

 

 

 

 

z2

 

 

.

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

v

lim v

 

lim (z V

 

 

z v )

 

lim

 

 

z

z

 

 

 

 

 

0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

k 0 k 1

z

 

z

 

0

 

z

 

 

z 1 (z d ) (z d

2

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

67

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

k

2 и k

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

 

 

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 .

 

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

v0

0 и v1

0 .

 

 

 

 

 

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

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

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

vk 2 (d1 d2 ) vk 1 d1 d2 vk 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

f

k

.

1,k

1

2,k

 

2

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

68

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d2k

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c1,k

 

fk

d2k 2

 

 

 

 

1k d2k 1

 

 

 

 

 

1k

 

 

;

 

 

 

 

 

 

d k 1

d k 1

(d

 

 

d )

d k 1

(d

 

 

 

 

 

 

 

 

 

 

2

2

d )

 

 

 

 

1

2

 

1

 

 

1

 

 

1

 

 

 

 

d k

1

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c2,k

d1k 2

fk

 

 

1k d1k 1

 

 

 

 

 

1k

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d k 1

d k 1

(d

d )

 

d k 1

(d

 

 

 

 

 

 

 

 

 

 

 

 

2

d )

 

 

 

 

1

2

 

 

2

1

 

 

2

 

 

1

 

 

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

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

 

c1,k

и c2,k

применим

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

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

используя

 

 

 

для

 

раскрытия

 

 

 

 

сумм

формулы

арифметической либо геометрической прогрессий

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

1n 1

 

 

 

 

 

 

 

 

 

 

 

 

c1,k

 

 

 

c1,k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1 d1n (d2

d1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

k 1

 

 

 

 

 

 

 

 

(d k

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

c1 ;

 

 

 

 

d2

d1 n 1 d1n

 

 

 

 

(d2

d1) d1k (d1

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

1n 1

 

 

 

 

 

 

 

 

 

 

 

 

 

c2,k

 

 

 

 

c2,k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1 d2n (d2

d1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

k

1

 

 

 

 

 

 

 

 

 

 

 

 

 

d k

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

c2 ,

 

 

 

 

d 2 d1 n 1 d2n

 

 

 

(d 2 d1) d2k (d2

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

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

и c2,k

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

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

 

 

 

 

 

 

 

 

 

d1k

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

d2k

1

 

 

 

 

k

;

vk

 

 

 

 

 

 

 

 

 

 

 

 

c1

 

d1

 

 

 

 

 

 

 

 

 

 

 

 

 

c2 d2

 

(d 2

d1) (d1

1)

 

 

 

(d 2 d1) (d2

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d1k

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d2k

1

 

 

 

 

 

k

k

 

vk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c1

d1

c2 d2 .

 

(d 2

d1) (d1

1)

 

 

(d 2 d1) (d2

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

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

0 и v1

0 . Так, приравнивая

общее решение, при k

0 и k

 

 

1, начальным условиям, находим

 

 

 

 

 

 

 

 

v0

0

 

 

 

 

0

 

0 c1

c2 ;

 

 

 

 

 

 

 

 

 

v1

0

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

1

 

 

c1 d1

c2

d2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d2

 

 

 

 

d1

 

d2

d1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c1

 

 

c2

0 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d1

c1

 

d2

c2

 

0 .

 

 

 

 

 

 

 

 

 

 

 

69

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

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

 

 

 

 

 

d k

1

 

 

d k

1

 

 

 

 

vk

 

 

1

 

 

 

 

 

2

 

 

 

 

 

 

(d 2

d1) (d2

1) (d 2 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)

 

 

 

 

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

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

 

vk 2

(d1

d2 ) vk 1

d1 d2 vk

1k

fk

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

Так,

вводя новые

переменные

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

или

 

 

 

 

 

 

 

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 диагональная матрица указанной функции от каждого собственного значения.

70

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

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 .

Доказывается, что модальная

матрица

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 имеют вид