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

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

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

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

Добавлен: 16.10.2024

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

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

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

5)

умножим первый

столбец на а второй — на В результате,

как и следовало ожидать, мы получим

 

 

 

 

А ^ С =

А

0

0

0

12

0

'

 

0

1

0

0

 

0

0

 

 

 

-0

0

0

0 _

 

 

 

 

 

 

поскольку ранг матрицы А равен 2.

Аналогичные операции с матрицей размером т X я и ранга г при­ водят к следующим результатам:

1г Огхп—г -

АтХп

От— гХ г От—г Х п —г

где нули обозначают нулевые матрицы, размеры которых указаны со­

ответствующими индексами. В случае, если г = т <

п, матрица А

приводится к виду 10], а если г — п<.т, то — к виду

; если же

г= т = п, то матрица приводится к виду 1п.

Внашем примере предполагалось, что существуют такие матрицы

Ри Q, при которых PAQ = С, где Р и Q представляют собой произве­

дения

элементарных операторов. Если

А — это матрица размером

т х п ,

Р —квадратная матрица порядка т,

a Q—- квадратная матри­

ца порядка п, тогда

 

 

 

Рт X п Qn С

 

0'

 

о

о

Матрицу С обычно называют эквивалентной канонической формой мат­ рицы А, или канонической формой, полученной с помощью эквивалент­ ного преобразования*. Процедуру перехода от матрицы А к матрице С часто называют приведением к эквивалентной канонической форме.

Если требуется найти только общий вид С, можно поступить так, как было сделано выше, т. е. выполнить необходимые элементарные

операции с матрицей А,

не выписывая при этом в явном виде матриц

Р и Q. А в тех случаях,

когда требуется полностью выписать эти опе­

раторы, можно воспользоваться следующим соотношением: PAQ = = (PI)A (1Q). Поэтому, преобразовав единичную матрицу с помощью тех же действий, которые проделывались над матрицей А для того, чтобы привести ее к канонической форме, мы получим матрицы Р и Q, при этом Р объединяет все операции со строками, a Q — все операции

*В нашей литературе в таких случаях обычно употребляются термины «ка­ ноническая матрица» или «каноническая диагональная матрица» (см. например, Г а н т м а х е р Ф. Р. Теория матриц. М., «Наука», 1967). — Прим, перев.

153


со столбцами. Так, в нашем примере можно получить Р, проделав с / 3 операции 1 и 2. В этом случае:

"1 0

0 '

-

1 0

0 “

'

1

0

0 “

0 1

0

 

— 2 1 0

 

— 2

1

0

о 0

1

 

— 1 0

1_

 

1 — 1

1_

Кроме того, проделав с /4 операции 3, 4 и 5, мы получим Q:

 

-1

0

0

о-

1 —3 —2 1-

h

0

1

0

0

0

1

0

0

0

0

1

0

0

0

1

0

 

 

^0

0

0

1-

0

0

0

п

1 —3

4

 

2-

1

1

4

2

 

2

0

1

2

 

- 1

0

1

2

1

0

0

1

 

0

3

.0

0

0

 

1_

0

0

1

' 0

 

-0

0

0

1 .

 

 

 

 

 

Проверим этот результат, проведя следующие вычисления:

'

1

0

0

“2

6

4

2“

=

РА =

2

1

0

4

15

14

7

 

1

— 1

1

2

9

10

5

 

и

“2

6

4

2“

0

3

6

2

О

О

О О

“2

6 4 2“

- 1

4

2

“ 1 0 0 0

PAQ = 0

3

6 3

0

4 -

- 2 - 1

 

0

1 0

0

О

О

О О

о

о

«—

О

о

о

О о

 

 

 

о

о

О

1 гН

 

 

 

 

 

 

1

 

 

 

что и требовалось доказать.

В вычислениях такого рода матрицы Р и Q определяются неодно­

значно. К примеру,

мы

могли в

наших расчетах операцию, эквива­

лентную операции 5, применить

к строкам, а

не к столбцам.

Тогда

по-прежнему

сохранилось бы равенство С = PAQ, но Р и Q

теперь

имели бы иной вид:

 

 

 

 

 

 

 

 

 

1

0

0

 

“ 1

—3

4

2

 

 

2

 

0

1

2

1

 

 

 

 

 

 

р =

2

1

 

И

 

0

Q =

0

1

0

 

 

Т

3

 

0

 

 

1

1

ь

 

_0

0

0

1 _

 

154


Конечный вид операторов Р и Q зависит от того, в каком порядке проводятся отдельные элементарные операции, однако матрица С всегда будет иметь один и тот же вид независимо от того, как выглядит исходная матрица. У нее те же размеры, что и у матрицы А, и все мат­

рицы ранга г могут быть приведены к такому виду, когда С =

®

(нули здесь обозначают нулевые матрицы соответствующих размеров).

10. КОНГРУЭНТНОЕ ПРИВЕДЕНИЕ СИММЕТРИЧЕСКИХ МАТРИЦ

Процесс приведения симметрических матриц к каноническому виду более прост; этим обстоятельством часто пользуются в тех слу­ чаях, когда особенно удобна каноническая форма матриц, и прежде всего при упрощении вида квадратичных форм.

В уравнении (10) матрица приведена к так называемой диагональ­ ной форме:

где Dr — диагональная матрица, содержащая г отличных от нуля элементов. Рассмотрим процесс приведения симметрической матрицы к форме А. В симметрической матрице элементы каждого столбца сов­ падают с элементами соответствующей строки; поэтому, если осущест­ вить с соответствующими столбцами те операции, которые необходимо проделать со строками матрицы Л для того, чтобы обратить в нуль эле­ менты, расположенные ниже диагонали, то при этом обратятся в нуль также и те элементы, которые расположены выше диагонали. Следо­ вательно, если матрица А симметрична, проводя со столбцами и стро­ ками одни и те же операции, можно привести А к диагональной форме Д. Но это означает, что в произведении PAQ матрица Q равна транс­ понированной матрице Р, т. е. Q Р '. Следовательно, для симметри­ ческой матрицы РАР' = А.

Предположим, что все диагональные элементы положительны; тогда можно построить диагональную матрицу RT, элементы которой будут представлять собой квадратные корни из единицы, деленной на соответствующий диагональный элемент матрицы D, т. е. Rf =■■D p \ Допустим, что порядок матрицы А равен п и что

 

F =

R r

0 '

 

 

 

0

1 п - Г _

 

 

 

 

 

 

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

(FP)A (FP)r =

FPAP'F = С

U r

01

о

о , потому что

F' = F. Это означает, что в том случае, когда все элементы матрицы D,

положительны,

симметрическая

матрица может быть приведена к ка­

нонической форме с помощью

эквивалентного

преобразования: для

155


этого нужно умножить А слева на матрицу FP и справа на ее транспо­ нированную матрицу P'F. Такие вычисления называются конгруэнт­ ным приведением симметрической матрицы, а С называется канони­ ческой формой, полученной с помощью конгруэнтного приведения.

Теперь предположим, что не все диагональные элементы матрицы Dr положительны, пусть q из них отрицательны. Тогда произведе­ ние РАР' умножим слева и справа на такой элементарный ^-оператор, чтобы первые r—q диагональных элементов матрицы Dr оказались по­ ложительными, а последние q элементов — отрицательными. Новое произведение будем обозначать Р*АР*'. Под F меч будем понимать ту же матрицу, что и раньше, только R T теперь образуют элементы, равные квадратному корню из единицы, деленной на соответствующие элементы матрицы Dr, независимо от того, какой у них знак. В резуль­ тате умножения получаем

-

F-q

 

0

0

(FP*) A (FP*)' = FP* АР*’ F =

0

-

v

°

 

 

.

0

 

0

0_

Мы рассмотрели общую методику конгруэнтного приведения симмет­ рической матрицы А к канонической форме. Если q = 0, тогда канони­ ческая форма сводится к С. Разность между порядками Ir_q и I q назы­ вается сигнатурой Л; другими словами, сигнатура равна г—2q. На­ личие на диагонали символов со знаком минус означает, что для при­ ведения матрицы к этой форме необходимы только действительные чис­ ла. Если читатель настолько подготовлен, что может оперировать и с

мнимыми числами, включающими i = 1 , тогда в процессе приведе­ ния матрицы к форме С можно воспользоваться сомножителем F, со­ держащим мнимые числа.

Пример. С помощью следующих операций можно привести симмет­

рическую матрицу

 

 

 

 

' 4

12

 

 

 

12

27

 

 

 

к диагональной форме.

 

Вид матрицы А

Операции

 

 

 

после преобразований

1) строка 2—3 X (элементы строки 1)

 

у

4

12

 

0 —9

 

 

 

2) столбец 2—3 X (элементы столбца

1)

/

4

0

0 —9

Р можно получить, проделав операцию

с матрицей

 

 

РАР’

1 с о

' 1

0

'

0

1

 

О

■ 4

12'

 

12

27

10'

3 1

 

СО

О

-------

 

1

'4

0'

0

—9

156


В этом случае

" 1

 

- 1

 

 

 

 

 

0 "

 

 

R

п

0

2

1

 

 

 

 

 

0

 

 

 

L

У9 J

 

 

 

 

3_

 

 

и поскольку здесь F = R, то

 

 

 

 

 

 

2

0

 

 

" —

0

 

FP>= RP =

'

1

0

2

 

 

 

 

3

1

 

 

 

 

0

- 1

- L

 

 

 

 

 

 

 

3

 

 

 

3

 

Тогда

 

 

 

 

 

 

 

 

 

4

12

Г -

- О

1

0~

FPA(FP)'

 

2

 

1

12

27

0

0

1 _

 

 

 

 

 

 

3

 

 

что и требовалось доказать.

Пример. Когда приходится сталкиваться с квадратичными форма­ ми (см. параграф 3 главы III), во многих случаях полезно выяснить, можно ли представить х'Ах в виде суммы квадратов линейных комби­ наций из элементов вектора х; есди это возможно, тогда форма х'Ах положительно полуопределена. Пусть даны, например,

'хГ

"

4

12

20~

х2

и А =

12

45

78

,х3_

20

78

136_

Тогда (см. равенство (10) из параграфа 3 главы III)

х'Ах — 4х\ + 45x1 + 136x1 + 24x^2 + 40x^3 +

+1 5 6 х 2х 3.

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

Ч

0

4)

 

' 1

0

0

Р==

— 1

4 -

0

, то РАР'

0

1

0

 

 

3

 

 

-0

0

0-

_

1

— 2

1 _

 

 

 

 

 

Допустим, что мы теперь можем прибегнуть к линейному преобразо­ ванию у = Р'-Н или (что равносильно приведенному выражению)

157