Файл: Сирл, С. Матричная алгебра в экономике.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 16.10.2024

Просмотров: 139

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Подставляя в них ^ = 1 и произвольно приписывая аг значение 1,

Г

. Подставляя в уравне­

находим характеристический вектор и г = 1

1

 

ние Я,2 = — и полагая а 2=^0, вычисляем значение другого характери­

О

стического вектора: и2 3 ; наконец, при Х3= и а3 = 0 находим

—4

третий характеристический вектор:

и3=

 

 

 

. Следовательно,

 

'

1

 

0

0 '

я

 

 

 

 

"

24

0

 

 

 

 

 

 

 

 

 

c f

 

 

U --

1

 

3

3

И и~1= —

 

—1

4

—3

 

 

 

1

—4

4

 

 

 

24

 

—7

4

3

 

 

 

 

 

 

 

 

 

 

Таким образом,

k-ю степень матрицы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

0

 

 

0

 

 

 

 

Pk = UDk U~l

 

 

1

О

О

0

/

1

'

 

 

 

24

О

О

 

 

1

3

3

(

т .

)

 

0

 

—1

4

—3

24

 

 

 

 

 

 

1

—4

4

 

 

 

 

/

5

 

—7

4

3

 

 

 

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т

 

 

 

 

можно представить

в следующем

виде:

 

 

 

 

 

 

 

 

Pk =

 

 

 

1

 

 

0

 

 

 

0

 

 

 

 

1 rk/2 3sh/8

rJ2

 

 

3sh/8

 

 

 

 

 

 

_ l ~ r h/ 2 - 2 s h/3

2sk/3

 

rh/2

 

 

 

где гЛ=1 - ‘ H i

 

\k

 

j 5 \k

H i f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С помощью этих выражений мы можем легко рассчитать значения элементов Рк при любом k. Пусть, например, k = 2,

/ 5 \ 2

+

1

у

__ 26

и

S2с

 

 

 

24

_

сле-

тогда г2 == ( — J

 

 

 

( 1

у

/ i y _

 

 

довательно,

 

6

J

~ 3 6

 

 

\ 6

j

e 1J

36

 

 

 

 

 

0

0

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

26

 

1

26

1

 

"1

0

0

 

 

Р2 --=

 

72

 

4

72

4

=

0,39

0,36

0,25

 

 

 

26

 

4

4

26

 

0,20

0,44

0,36

 

 

 

72

 

9

~9

72

 

 

 

 

 

 

320


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

k значения г* и sk стремятся к нулю, а следовательно, и матрица Рк

~1 О (Г

. Такое предельное

в пределе должна иметь следующий вид: 1 0 0

1 0 0

от начального со­

значение Рк свидетельствует о том, что независимо

стояния машина в конечном счете придет в

состояние 1 (непоправи­

мо изношена). Такое состояние, которое в

конце концов распро­

страняется на всю систему, называется состоянием ловушки, или поглощающим состоянием (см. параграф 3 главы VIII).

в) СЛУЧАЙ, КОГДА ХАРАКТЕРИСТИЧЕСКИЕ КОРНИ КРАТНЫЕ

До сих пор мы предполагали, что все характеристические корни различны, в таком случае U представляет собой матрицу полного ранга; существует обратная матрица {/-1 и А = UDU-1. Если не все характеристические корни различаются между собой, тогда может не существовать обратной матрицы.

Если обратной матрицы не существует, матрицу А невозможно преобразовать в UDU_1. В приводимой теореме рассматриваются условия, позволяющие представить матрицу А в виде произведения

UDU~\

Обозначим через гп* показатель кратности корня Хк в характерис­ тическом уравнении | А — XI | = 0. В таких случаях принято говорить, что Хк — это корень кратности лг* (при k = 1, ..., s). Подобное вы­ ражение означает, что при разложении уравнения | А XI | на мно­

жители, можно выделить в качестве сомножителя (X Хк ) ть. А так

как характеристическое уравнение

представляет собой

полином

л-й

 

S

тк =

 

 

 

 

 

степени,

то ^

л.

Приведенные утверждения включают и

тот

случай,

k =

1

1

для всех k.

Подобная ситуация

анализирова­

когда тк =

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

Теорема 1. Дана квадратная матрица л-го порядка, имеющая ха­ рактеристические корни Яь Я2, ..., Xs кратности ти тг, ..., та, причем

S

для того чтобы А была представима в форме произ-

2,тк ^=п, тогда,

k = 1

где D — диагональная матрица, содержащая все п

ведения UDU-1,

характеристических корней, необходимо и достаточно, чтобы ранг матрицы Хк1) был равен л — тк при каждом значении k —

= 1, 2 , . . . , s.

Выводы из этой теоремы имеют большое значение: при любом зна­ чении характеристического корня X кратности т матрицу можно представить в форме произведения UDU~l тогда и только тогда, когда

г (А XI) = л — т.

Доказательство этой теоремы см. в работе.Сирла [9, с. 175].

11 Зак. 4 25

321


Это утверждение включает также случай, когда т = 1 и кратность любого корня равна единице.

Пример. Пусть дана матрица

— 1 — 2 — 2

1 2 1

— 1 — 1 0

Ее характеристическое уравнение можно привести к следующему виду: (А, — I)2 (А, + 1) = 0; оно имеет корни 1,1, —1. Следовательно, А = 1 — это корень кратности 2. Теперь надо исследовать ранг мат­ рицы А — А/. Для значения А = 1

 

 

А — II

—2

—2

—2~

 

 

 

 

 

 

1

1

1

 

 

 

 

 

 

 

—1 —1 —1

 

 

 

 

 

 

 

3 — 2 — п т

Следовательно, мат­

рицу А можно представить в виде произведения

U D U - 1.

 

 

 

Для того чтобы составить матрицу U ,

нужно отыскать такие ли­

нейно-независимые векторы и, которые при А = 1

удовлетворяют урав-

нению — А/) и =

0ЧТакими

 

 

 

Г

и и„ =

векторами будут иг

 

При А3 =

—1 решением

уравнения (Л — А3/) и 3 =

0

слу-

2"

 

 

 

Г

1

1

21

жит вектор и 3=

—1 , а следовательно, U =- [wi и2 и 3]— —1

0—1

.

L

 

 

 

.

0—1

1_

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

что произведение и~гА и

равно соответствующей диагональной матри­

це D и, следовательно, А =

UDU-1.

 

Пример. Дана матрица

 

 

 

 

 

2

—1

1

Л =

3

3

—2

 

 

4

1

0

Ее характеристические

корни

имеют

следующие значения: А = 1;

1 и 3. Таким образом, А =

1 — это кратный корень, т — 2. Подставив

его в выражение А — А/, получаем Матрицу

1 —1 1 ( Л - / ) 3 2 —2 4 1 —1

3 2 2


Ранг этой матрицы равен 2 (элементы строки 3 равны сумме соответст­

вующих элементов строк 1 и 2). Таким образом, при я = 3

и т — 2

г (А — /) = 2 Ф 3 — 2, так что матрица А неприводима

к виду

D = U-'AU.

 

3. НЕКОТОРЫЕ СВОЙСТВА ХАРАКТЕРИСТИЧЕСКИХ КОРНЕЙ

Рассмотрим некоторые свойства характеристических корней и век­ торов, поскольку ими часто приходится пользоваться при решении различных задач.

а) СТЕПЕНЬ ХАРАКТЕРИСТИЧЕСКИХ КОРНЕЙ

Пусть Яг представляет собой решение характеристического уравне­ ния матрицы А и Ui —- соответствующий характеристический вектор. Тогда, следуя (4), можно записать Лм;= Ящг. Умножив это выражение слева на А, получим

A 2Ui = AliUi = Я;Ляг,

(10)

Яг — это скалярная величина, поэтому вновь

воспользуемся соотно­

шением (4); тогда выражение (10) приобретает следующий вид: Л2«г = = Я-цг. Этот результат легко обобщить, показав, что

A kUi = %ktUi,

(11)

где k — любое положительное целое число. Приведенное соотношение выполняется и для отрицательных целых чисел: так, А~к = ( Л -1) 6, если существует обратная матрица Л -1 (см. упражнение 14). Но Я и и удовлетворяют условию Аи = Яя лишь в том случае, если Я пред­ ставляет собой характеристический корень матрицы Л, а и — соот­ ветствующий характеристический вектор. Сравнивая последнее вы­ ражение с (11), легко прийти к следующему выводу: Я* в выражении А ки = Яки является характеристическим корнем А к. Таким образом, отсюда следует, что если Я представляет собой характеристический корень Л, а Л — вырожденная матрица, Я* окажется характеристи­ ческим корнем матрицы А к при всех положительных целых k. Если же

Л — невырожденная

матрица, Я* — характеристический корень мат­

рицы А к при всех целых значениях k (положительных

или отрица­

тельных).

Допустим,

например, что

матрица Л — невырожденная

и k = —1,

а Я представляет собой

характеристический

корень Л,

в таком случае -j- — это характеристический корень матрицы Л -1. Пример. Дана матрица

 

 

 

Г 4 31

 

 

В таком случае

 

 

 

 

 

"28

36"

и

"

8

—3

48

76

Л -1 = —

—4

4

 

20

и *

3 2 3


Характеристические уравнения этих матриц имеют следующий вид:

для А: V — 12^ + 20 =

0, его корни равны 2 и 10;

для А 2:

— 104Х + 400 =

0, его корни равны 4 и 100;

для А -1:

20Х2— 12А, +

1 =

0, его корни равны -j и ^ .

Характеристические корни матрицы А 2 равны возведенным в квадрат характеристическим корням матрицы А, а характеристические корни матрицы Л -1 представляют собой обратные величины по отношению к характеристическим корням А.

6) СЛОЖЕНИЕ И УМНОЖЕНИЕ ХАРАКТЕРИСТИЧЕСКИХ КОРНЕЙ

! .

Рассмотрим следующее уравнение:

 

 

 

\А — М\ =

й \ \

^

а 12

а 13

= 0.

 

^ 2 1

^ 2 2 ^

^ 2 3

 

Разлагая определитель

a3l

®32

^33

'

записать

по

элементам

диагонали, можно

это выражение:

 

 

 

 

 

 

-А.® + (_я)* tri (А) +

(—Л) tr2 (Л) +

| А | = 0,

 

где символом tr* (А) обозначены суммы главных миноров i-го порядка матрицы А (см. параграф 6 главы IV). Если порядок матрицы А равен п, тогда это уравнение выглядит следующим образом:

(—^)п + (—X)"-1 ^гх(Л) + (—К)п~ 2 tr2(Л) + ...+

+ (—A) tr„_ 1(Л) 4 -1А | = 0.

/

[ Используя теорию уравнений, можно показать, что

 

2 ^ = *Г!(Л) и

П Я, = |Л |.

г=1

i=i

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

Пример. Характеристическими корнями матрицы Л

1

4

Г

2

1

0

будут следующие числа: 0, —1 и 4. Их сумма равна 3, и след мат­ рицы Л равен 1 + 1 + 1 = 3, произведение корней равно нулю, нулю равен и определитель этой матрицы. (Элементы строки 1 рав­ ны сумме соответствующих элементов строк 2 и 3.)

в) НЕНУЛЕВЫЕ ХАРАКТЕРИСТИЧЕСКИЕ КОРНИ

Мы рассмотрели условия, при которых матрицу Л, имеющую крат­ ные характеристические корни, можно представить в виде произведе­ ния и~хА и = 0. Если такое преобразование возможно, U представляет собой невырожденную матрицу и, следовательно,

ранг (D) = ранг (Л) = г.

324