Материал: 1395

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

 

 

 

56

 

 

 

 

 

 

 

 

Известно, что система линейных алгебраических уравнений n - го

порядка в векторно-матричной форме записывается в виде

 

 

 

 

 

 

 

A X

Y ,

 

 

 

 

 

 

где

A- матрица

коэффициентов системы; X , Y -

входной

и выходной

векторы. Естественно, что результирующий вектор Y

зависит,

как

от

матрицы коэффициентов A, так и от входного вектора X .

 

 

 

 

 

Возникает вопрос, существуют ли при данной матрице коэффициентов

A

такие векторы xi , преобразование

которых

сводится

лишь

к

их

умножению на

масштабный

коэффициент i ,

то

есть

выполняется

соотношение

 

 

 

 

 

 

 

 

 

 

 

A

xi yi

i

xi .

 

 

 

 

 

Оказывается, что в общем случае невырожденной матрицы имеется ровно n таких векторов, где n - размерность системы векторов. При этом, векторы, удовлетворяющие указанному условию, называются собственными векторами матрицы, а масштабные коэффициенты при векторах называются собственными значениями матрицы. Собственные векторы xi принято

обозначать в литературе как hi .

Записывая данные соотношения, определяющие собственные вектора и значения в векторно-матричной форме, получим

A [hi ] A H H ,

где H - модальная матрица, составленная из собственных векторов; - в случае невырожденной матрицы коэффициентов это диагональная матрица собственных значений.

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

Собственные вектора образуют собственный канонический базис линейно независимых векторов. Каноническое представление матрицы в собственном базисе запишется на основании предыдущего соотношения в виде

H

1 A H .

Наоборот любая невырожденная матрица на основании этого

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

A H

H 1.

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

57

F ( A) H F ( ) H 1,

где F ( ) - понимается как диагональная матрица с элементами f ( i ) .

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

 

A

H 0 ,

где A

- характеристическая матрица.

Далее, учитывая, что в общем

случае модальная матрица H не

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

A 0 .

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

det( A

)

A

0 .

 

 

 

 

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

n

an 1

n 1

a1

a0 (

1) (

2 )

(

n ) 0 .

 

 

Корни

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

полинома

i

соответствуют

собственным

 

 

 

 

 

 

 

 

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

 

 

 

Собственные вектора hi , соответствующие собственным значениям i

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

 

 

 

 

 

 

A

i hi B( i ) hi 0 ,

где

i

[ i ]- диагональная

матрица, содержащая одинаковое i - тое

собственное значение по диагонали.

 

Так как матрица

A

i

вырождена, то есть ее ранг равен, по крайней

 

 

 

 

 

мере,

n

1, получаем, что однородная система, определяющая собственный

вектор hi , должна иметь нетривиальное решение. Аналитические

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

n

 

 

 

 

bij (

i ) hij

0 .

 

 

i 1

 

 

 

 

Далее, раскрывая определитель матрицы

A

i

B( i ) , по какой либо

строке i , и, учитывая, что он равен нулю, получаем

 

n

 

 

 

 

det(B( i ))

bij ( i )

ij (

i )

0 .

i 1

 

 

 

 

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

 

58

 

 

hij k ij (

i ) ,

где j 1, 2,

, n ; k - некоторая константа;

ij ( i ) - алгебраическое дополнение

матрицы A

i .

 

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

C [cij ] Adj(B( i )) [ ij ] .

Из полученного соотношения следует, что собственный вектор hi , с точностью до константы, определяется алгебраическими дополнениями элементов i - той строки матрицы A i , то есть соответствующим

столбцом присоединенной матрицы C .

Разностные уравнения второго порядка. Прежде чем перейти к разностным уравнениям высокого порядка остановимся подробнее на уравнениях второго порядка.

Как уже отмечалось, в изучаемой дисциплине предпочтение отдано универсальным аналитическим методам – операторному, Лагранжа и Коши.

В общем случае нормированное неоднородное разностное уравнение

второго порядка имеет вид

 

 

 

 

 

 

 

 

yk 2

a1 yk 1

a0 yk

 

 

fk

или

 

 

 

 

 

 

 

 

 

(E2

a E

a )

y

f

k

,

 

 

1

 

0

k

 

 

где yk - неизвестная

дискретная

функция;

 

yk n - сдвиг неизвестной

дискретной функции

на n -

тактов;

ai -

коэффициенты уравнения; fk -

произвольная известная дискретная функция, определяемая внешним воздействием.

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

 

zn

 

n 1

m

 

y

Y

y z

,

k n

 

z

m

 

 

 

 

 

m 0

 

 

или частный случай при n

1

 

 

 

 

 

yk 1

z Yz

z y0 .

 

 

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

Yz

Fz

Nz

,

z2 a

z a

 

 

 

1

0

 

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

воздействия; N z - изображение,

обусловленное начальными условиями.

 

 

59

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

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

 

Y

b z2

b

z

b

Sz

z

 

2

1

 

0

,

X z

z2

 

 

 

 

a z

a

 

 

 

 

1

 

0

 

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

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

то, выражая изображение

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

 

 

 

 

 

 

 

 

 

 

 

 

(b z2

b

z

b )

X

z

 

Y

S

z

X

z

2

1

 

0

 

,

 

 

 

 

 

 

z

 

 

z2

a

z

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

0

 

 

 

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

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

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

Рассмотрим аналитическое решение нормированного неоднородного разностного уравнения с постоянными коэффициентами второго порядка

E2 a E a

y

k

y

a y

k 1

a y

k

f

k

1

0

 

 

k 2 1

0

 

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

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

 

 

y

c

 

y

c

y

c

d k

c

d k ,

 

 

k

1,k

1,k

2,k

2,k

1,k

1

2,k

2

где y

d k ,

y

d k

- фундаментальные решения однородного уравнения;

1,k

1

2,k

2

 

 

 

 

 

 

 

 

d1, d2 - корни характеристического уравнения.

С целью определения варьируемых постоянных подставим общее решение в исходное разностное уравнение. Для этого последовательно подействуем на предполагаемое решение оператором сдвига. После первого действия получим

yk 1 c1,k 1 y1,k 1 c2,k 1 y2,k 1.

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

 

60

 

yk 1 c1,k 1 y1,k 1

c1,k y1,k 1 c2,k 1 y2,k 1

c2,k y2,k 1

c1,k y1,k 1 c2,k y2,k 1

 

c1,k y1,k 1

c2,k y2,k 1 c1,k y1,k 1

c2,k y2,k 1 .

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

 

c1,k y1,k 1

c2,k y2,k 1 0 .

 

В результате наложения ограничения имеем

 

 

yk 1

c1,k y1,k 1

c2,k y2,k 1.

 

Теперь повторное действие оператора сдвига дает выражение

yk 2

c1,k 1 y1,k 2

c2,k 1 y2,k 2

 

 

 

 

c1,k y1,k 2

 

c2,k y2,k 2

c1,k y1,k 2

c2,k y2,k 2 .

Подстановка полученных выражений в исходное разностное уравнение

и перегруппировка дают

 

 

 

 

 

 

 

c1,k y1,k 2

 

c2,k y2,k 2

 

 

c1,k [ y1,k 2

a1 y1,k 1

a0 y1,k ]

 

 

c2,k [ y2,k 2

a1 y2,k 1

a0 y2,k ]

fk .

Заметим, что выражения в скобках равны нулю, так как удовлетворяют

однородному

разностному

уравнению

в

силу того,

что y1,k , y2,k - его

фундаментальные решения.

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

c1,k y1,k 2 c2,k y2,k 2 fk .

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

c1,k y1,k 1

c2,k y2,k 1

0 ;

c1,k y1,k 2

c2,k y2,k 2

fk .

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

c1,k

y1,k 1

y2,k 1

0

.

c2,k

y1,k 2

y2,k 2

fk

 

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