Материал: 1395

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

86

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

vk

c1,k v1,k c2,k v2,k

c1,k 1k c2,k d k ,

где 1, d - корни

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

уравнения; y1,k 1k , y2,k d k -

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

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

c1,k 1k 1

c2,k d k 1

0;

c1,k 1k 2

c2,k d k 2 10

fk .

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

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

Выразим разности варьируемых постоянных c1,k

 

 

и c2,k из

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

 

 

 

C

 

1k

1

 

 

d k

1

 

1k 1 d k 1

(d 1) ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1k

2

 

 

d k

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d k 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c1,k

 

10

d k 2

 

 

 

 

10 d k 1

 

 

 

10

 

;

 

 

 

 

 

 

 

 

 

1k 1

d k 1 (d 1) 1k 1

 

 

 

 

 

 

 

 

 

 

 

 

(d 1)

 

 

 

1k

1

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c2,k

 

1k 2

10

 

 

 

 

 

10 1k 1

 

 

10

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1k 1

d k 1 (d 1) d k 1

 

 

 

 

 

 

 

 

 

 

 

(d 1)

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

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

 

c1,k

и

c2,k применим

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

1

k

1n 1

 

1 k 1n 1

1

 

 

 

c1,k

c1,k

 

 

 

 

 

 

 

 

c1

;

 

 

d 1 n 1 1n

 

 

 

 

n 11n (d 1)

 

 

d 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

87

 

 

 

 

 

 

 

 

 

 

 

 

 

c2,k

1

 

c2,k

 

k

 

1n 1

 

 

1

 

 

 

k 1n 1

 

 

1

c2 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1 d n (d 1)

 

d 1 n 1 d n

 

d (d 1)

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

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

c1,k

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

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

 

 

 

 

 

vk

 

1k

 

 

c1 1k

 

 

d k

 

c2

d k

 

 

 

 

 

d

1

 

 

d

(d

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d k

1

 

1

c

 

c

d k

1

d k

 

1

 

c

c

d k .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

1

 

 

1

2

 

 

1

 

d

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

и

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

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

 

0 и v2

 

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

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

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

решение при k

1

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

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

 

 

 

 

 

 

 

v1

0

0

c1

c2

d ;

 

 

v

1

1

c

c

d 2

 

 

2

 

 

1

2

 

 

 

или

 

 

 

 

 

 

 

 

c1

d

c2

0 ;

 

 

 

 

c

d 2

c

0.

 

 

 

 

1

 

2

 

 

 

 

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

c2

0 .

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

 

 

1

 

d k

1

,

 

gk

vk

 

 

 

 

 

 

1

d

 

 

 

 

 

 

 

 

при k

1 или

 

 

 

 

 

 

 

 

 

 

 

 

1

d k

 

gk

1 vk 1

 

 

 

,

 

 

1

 

 

 

 

 

 

d

при k

0.

 

 

 

 

 

 

 

 

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

описывает

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

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

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

vk 2 (1 d ) vk 1 d vk 10 fk

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

88

 

Так,

вводя

 

новые

 

переменные

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

 

10

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X k 1

 

A X k

Fk .

 

 

 

 

 

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

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

 

 

 

 

 

 

 

 

Ak X

 

k

n F

 

 

 

 

 

 

 

 

 

X

k

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

(1 d )

d 0 .

 

 

d

1 d

 

 

 

 

 

 

 

 

 

 

 

 

 

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

1

0

1

0

 

 

d .

0

2

0

 

 

 

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

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

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

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

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

 

 

 

A i hi

0,

 

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

89

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

матрицы [ A

i ]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11( 1)

11( 2 )

 

1 d 1 1 d d

d 1

H

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

0

fk

10

 

 

 

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

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

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

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

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

и Fk ,

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

системы

 

 

 

 

 

 

 

 

 

 

 

 

 

v 0

k

1k n

d k n

1

 

1

 

1k k

1n 1

d k

k

1n 1

.

 

 

 

 

 

 

 

 

k

 

d 1

n 1

d 1

n 1 1n

n 1 d n

 

n 1

 

 

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

0, получаем значения сумм равные

90

данном случае существует только при k первым слагаемым.

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

 

 

 

 

 

 

1

 

1 d k

1

1

d k 1

 

g

k

v

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

k

 

d

1

 

 

 

 

1

d

 

 

 

 

 

 

 

 

 

при k

1 или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

gk

 

 

vk 1

1

d k

,

 

 

 

 

 

 

1

 

 

1

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при k

0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Полученные выражения

совпадают

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

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