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

книги из ГПНТБ / Сирл, С. Матричная алгебра в экономике

.pdf
Скачиваний:
20
Добавлен:
22.10.2023
Размер:
13.54 Mб
Скачать

(Далее мы покажем, как подсчитать Л -1.) После умножения это урав­ нение примет вид:

1

0

*1

4 200

0

1

_Х2_

39 000

или

и' 4 200

 

х2

39 000

 

 

 

Это означает, что решение имеет вид: хг

4 200;

х2 == 39 000. Таким

образом, мы нашли, что уравнение

 

 

 

 

1

0,2

"*i'

12 000

 

 

 

1 0 — 1

х2

3 000

 

 

имеет решение

 

 

 

 

 

 

1

' — 1

0,2

" 12 000

'

4 200

Х2 * —3

10

1

3 000

 

39 000

Взятая для получения данного результата матрица Л ^1 была так скон­ струирована, что ее произведение о исходной матрицей равно единичной матрице, т. е.

1

1 0,2

1

0 ,2

1 0 ~

3

— 10

1

10

1

0

1

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

1

0,2

12 000

10

1

и b

3 000

Решение уравнения

Ах — b

получаем, умножая его слева на А~х, что дает

А - 1Ах = А Ч .

Поскольку

л -м = /,

получим

х = А-Ч.

(3 )

(4)

(5)

(6)

В матричной алгебре матрица / служит единицей, поэтому матрица Л -1 называется обратной к Л. Определим ее теперь (см., например, урав­ нение (5)) как матрицу, которая в результате умножения справа на Л равна единичной матрице /.

Заметим, что это решение было получено из уравнения (3) без деле­ ния обеих частей выражения (3) на Л; иначе говоря, в процессе решения

100

отсутствовал этап деления на Л, поскольку в матричной алгебре деле­ ние не определяемо и не имеет смысла. Вместо этого решение было по­ лучено умножением слева обеих частей (3) на матрицу Л ~\ что дало (4). Затем, воспользовавшись определением А _1, содержащимся в (5), можно упростить уравнение (4) и получить решение (6).

Оставшиеся разделы этой главы посвящены прежде всего уточнению определения матрицы Л -1, заданной уравнением А*1А = I, и детализа­ ции процесса построения Л -1 на основе Л. Когда матрица Л -1 известна, то любая совместная система линейных уравнений, имеющая единст­ венное решение, может быть представлена в виде Ах = b и решена как х = А ^ Ь вне зависимости от того, сколь велико число уравнений. Иногда могут возникнуть затруднения вычислительного порядка, однако форма решения и процедура его получения остаются теми же самыми.

2. ПРОИЗВЕДЕНИЯ, РАВНЫЕ /

Понятие обратной матрицы было введено нами при изучении реше­ ния системы линейных уравнений. Однако до сих пор мы не ответили на основной вопрос: существует ли такая матрица R, что RA = /? На этот вопрос можно получить три ответа, зависящие от характеристи­ ки Л:

1) в некоторых случаях R

существует и

является

единственной

для данной матрицы Л;

 

 

 

 

2) иногда для конкретной матрицы Л можно получить различные

матрицы R; иными словами,

R существует, но она не единственна;

3) в ряде случаев R просто не существует.

 

 

Приведем теперь примеры этих трех ситуаций.

 

1. Если Лх

2

8

—5-

то RxAx =

/ и, как бу­

3

—3

2

дет показано, Rx определяется единственным образом для заданной матрицы Л.

1

. Тогда имеется беско­

2. Пусть нам дана матрица Л 2-= —1

3

 

нечное множество матриц R, для которых R A 2 /; например,

II я

3. Когда

 

СО

 

и Rz =

4

15

4

2

5

1

7

25

6

 

 

 

л 8

 

"0

3

7"

 

 

 

= 0

2

5 >

 

 

то не существует матрицы R, такой, что R A 3 = /, поскольку первый элемент RA 3 всегда будет равен нулю и, таким образом, произведение матриц не может быть равным /.

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

101

A XRX=

/,

также

как

R XA X= /;

 

 

 

 

A%R2 ф /,

хотя Т?2Л 2 =

/;

 

 

 

 

если 5

=

4

 

8

]

то A 3S

/, хотя не

существует

матри­

5

—7

,

 

 

-2

 

3

 

 

 

 

 

 

цы R , такой,

что R A . a -

 

I.

 

 

 

 

 

В рассмотренных примерах решения системы линейных уравнений

матрица А

определялась из следующего условия:

А " 1А

I.

Однако

мы показали,

что соотношение R A =

/ не обязательно предполагает

единственную матрицу R

при заданной матрице Л,

и произведение A R ,

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

не всег­

да будет равно I даже,

если R A равно I . Только для А х матрица R

обладает этими двумя свойствами, а именно: она определяется единст­ венным образом для данной матрицы Л и умножение ее на Л как справа, так и слева дает единичную матрицу. Будучи аналогами свойств знако­ мых нам обратных величин в скалярной алгебре, эти свойства пригод­ ны и для характеристики обратных матриц. В соответствии с этим выве­ дем для Л матрицу, которая обладает этими свойствами. Полученная матрица будет обратной к матрице Л, и вместо того, чтобы удовлетво­

рить только условию Л _1Л =

/, она будет обладать двумя свойствами:

а) Л _1Л = ЛЛ -1 =/ и б) Л -1

является единственной для данной матри­

цы Л.

 

ЛЛ-1, предполагает, что

Свойство (а), согласно которому А - 1 А =

существуют оба произведения Л _1Л и ЛЛ-1.

Как было показано в па­

раграфе 5 главы II, это может быть только когда и Л и Л л — квадрат­ ные матрицы одного и того же порядка. Это означает, что матрица Л -1, удовлетворяющая условию (а), существует, если Л — квадратная мат­ рица. Речь здесь идет о следующем: из сказанного вытекает, что обрат­ ная матрица может существовать в том случае, если исходная матрица квадратная (обратная матрица также будет квадратной и того же по­ рядка). В противоположность этому прямоугольные матрицы не имеют обратных матриц, хотя для некоторых из них можно получить более узкий класс обратных матриц; такие матрицы рассматриваются нами при изучении случая (б) (см. приложение к данной главе).

Ранее было установлено, что обратная матрица Л *1 может су­ ществовать только тогда, когда Л — квадратная матрица; определи­ тель Л также существует только в том случае, если матрица Л квадрат­ ная. Следовательно, обратная матрица может существовать только тог­ да, когда Л имеет определитель. Теперь введем такое определение об­ ратной матрицы, которое использует этот определитель. Такой способ определения матрицы избран в связи с тем, что так легче представить процесс определения элементов обратной матрицы непосредственно из элементов самой матрицы (хотя подобный способ вычисления обратной матрицы и является трудоемким). В практических ситуациях обратную матрицу редко определяйэт таким путем. Тем не менее этот метод поз­ воляет ясно представить общую взаимосвязь между элементами обрат­ ной матрицы. Кроме того, с его помощью легко проверить только что рассмотренные свойства (а) и (б).

102

3. В Ы В О Д В Ы Р А Ж Е Н И Я Д Л Я О Б Р А Т Н О Й М А Т Р И Ц Ы

Уравнение (1) главы IV показывает, как можно разложить опреде­

литель матрицы А =

{atj}, i, j = 1,2,

..., п в виде следующей суммы:

П

 

 

| А | 2

ац ( —1)^+' | M tj |

для любого i,

!=1

 

 

где | Л4jj | представляет собой минор элемента а^,'т. е. определитель, полученный из | Л | путем вычеркивания i-й строки и /-го столбца. Про­ изведение (— называют алгебраическим дополнением аи .

Оно является минором элемента atj, взятым с соответствующим зна­ ком, поэтому алгебраическое дополнение часто называют знаковым ми­

нором. Обозначим через

 

алгебраическое дополнение аи и получим

 

 

 

 

Р» = (— \)i+i\Mu\-

 

 

 

Пример. Рассмотрим снова матрицу, которая фигурировала в при­

мере рыночного равновесия в начале этой главы:

 

 

 

 

 

 

 

1

0,2

 

 

 

 

 

 

 

 

 

 

А =

ТО

-1

 

 

 

 

(7)

 

 

 

 

 

 

 

 

 

 

 

Ее определитель равен: | =

— I

- 2

=

 

-3, а

алгебраические до-

полнения элементов имеют вид:

 

 

 

 

 

 

Рп =

(- 1 )1+1(-1 ) =

— и

 

Р» =

( - 1 )1+2(10) =

-1 0 ;

Р21 =

(—1)2+1(0,2) =

-0 ,2 ;

 

р22 =

(-1 )2+2(1) = 1.

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

 

 

 

 

Р п

P l 2

 

— 1

 

- 1 0 '

 

(8)

 

 

 

Р 21

И-22.

 

— 0 , 2

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и результат

ее транспонирования

 

 

 

 

 

 

 

 

 

Р и

 

Р 2 1

 

1

0,2

 

 

 

 

 

Р 12

Р-22 .

 

10

1

 

 

 

Произведение транспонированной матрицы на матрицу А дает

Пи

Р21

«11

«12 "

 

- 1

0 , 2 '

'

1

0,2

' — 3

0

P'12

р22

«21

«22 .

 

— 10

1

.10 --1

0

- 3

и поскольку

| А |

= —3,

то

 

 

 

 

 

 

 

 

Р и « 1 1 ~Т Р 2 1 « 2 1

Р и « 1 2 “Ь Р 2 1 « 2 2

 

1

О

 

 

 

 

 

 

 

Р 12 fljx ~Т Р 22 « 2 1

P l2

« 1 2

Л" Р 2 2

« 2 2 .

0

1 Л 1

 

 

 

 

 

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

103

ражения (9) на скаляр 1/| А | (обратную величину определителя А), мы получим единичную матрицу:

(М И'гг Й11 а 12

1 'M i

0 -

Р-12 Р22. Й21 ^22

~~мТ О

1---

Таким образом, на основе предварительного определения Л -1, а имен­ но А - 1А = /, получим

1

Вп

Р 21

1

rB u

P l 2

А - 1*

 

 

 

 

 

M l

9 j2

^22 j

M l

И-21

Р 22

т. е. матрица

 

 

 

 

 

 

 

- 1

0,2

( 10)

 

 

10

1

 

 

 

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

Соотношение (9) содержит общий результат, согласно которому оп­ ределитель | А | представляет собой сумму произведений элементов строки (столбца) на их алгебраические дополнения:

П

П

 

 

1А | - 2 аи \хи = 2 аи \х,и .

(11)

/ = 1

i =

1

 

Этот результат объясняет происхождение диагональных элементов | А |

в уравнении (9). Подобно этому нули в (9)

являются следствием другого

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

любого определителя сумма

произведений элементов одной строки (или столбца) с алгебраически­

ми дополнениями элементов другой строки (или столбца) равна нулю1:

П

аи цк}^-0 при i^=h

/= 1

и

П

 

 

2

auVih = 0 при i=h k.

 

(12)

г=1

 

 

1 Выражение (11) справедливо, поскольку, по определению

оно идентич­

но разложению | А [ на миноры, приведенному в параграфе 2 главы IV. Для

того

 

 

П

при

чтобы убедиться в справедливости (12), заметим, что выражение 2

 

 

/=1

 

i Ф h представляет собой аналог разложения определителя, содержащего стро­

ку элементов а^, а,2, •••,

Щп> 3 в остальных п — 1 строках те элементы

| А |,

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

p/j2, ....

р

104

Общую технику определения Л -1 можно более четко показать на матрицах размером 3x3. Возьмем матрицу

1

2

3"

 

А = 4

5

6

(13)

7

8

10

 

Для нахождения обратной к ней матрицы необходимы алгебраические дополнения всех элементов этой матрицы. Начнем с первого столбца. Алгебраические дополнения его элементов составят соответственно

5

6

 

2

3

( - D 1+1 8

10 2

( - 1)2+1 8

10 = 4

и

2

3

 

 

 

= —3.

 

Г_

1)3+1

6

(И )

 

5

 

 

Аналогично для второго столбца алгебраические дополнения элементов равны:

4

 

6

2,

 

1

3

-

— 11 и (— 1)

1

 

(15)

( - 1) 7

 

10

( +

1) 7

10

4

 

а для третьего столбца:

 

 

 

 

 

 

 

 

 

4

5

=

—3,

 

1

2

 

6 и ( + 1)

1

2

— 3.

(16)

( 1- 1) 7

8

( -

1

 

 

4

5

 

 

) 7 8

 

 

 

 

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

~2 ■2_ —3

4

— 11 6 ,

—3

6 —3

полученную путем замены в матрице А каждого элемента на его алгеб­ раическое дополнение, как это было сделано в (8); таким образом, столбцы новой матрицы состоят из величин, полученных в результате операций (14), (15) и (16). Транспонируем новую матрицу и умножим

Однако,

поскольку эти алгебраические дополнения относятся

к элементам а/ц,

йдг, •••.

при этом эти строки охватят все

строки А

за

исключением (а/ц,

вд2........

aim). Поэтому среди

них окажется и строка (ощ, аг-2........

ain), другими

словами,

п

представленное

в форме

определителя, будет

выражение 2

1=1

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

105

ее на скаляр, равный обратной величине определителя. Как и в (10), результирующая матрица представляет собой матрицу, обратную к А:

 

 

Г

2

4

- 3

 

(17)

 

Л- 1- —

 

2— 11

б

 

 

 

" 3 [ - 3

6

—3

 

 

Умножая слева (13) на (17), покажем,

что А - 1А =

/:

 

2 4 - 3

1 2

3

 

- 3 0 о

2 — 11

6 4

5 6

 

0—3 о

 

 

7

8

10

 

о

о

3

 

 

1

0

0

 

 

 

(18)

 

 

0

1

0

 

 

 

 

 

0

0

1

 

 

 

 

Отсюда (17) есть обратная

к А матрица.

 

 

2x2, диа­

По тем же самым причинам, что и в примере с матрицей

гональная матрица выражения (18) содержит в качестве ненулевых элементов определитель | А | (в данном примере | А | равен —3). Зна­ чения 0 и —3 были получены из (11) и (12). Например, суммируя про­

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

первого столбца

с этими элементами, получим определитель в соответствии с (11):

2 (1) + 4 (4) — 3 (7) = —3 = | А | .

(19)

Умножая алгебраические дополнения (14) на элементы второго и треть­ его столбцов Л и суммируя результаты, получим в соответствии с (12)

нули:

(2) + 4 (5) — 3 (8) = О

2

и

(20)

2(3) + 4(6) — 3(10) = 0.

Полученные суммы произведений (19) и (20) являются элементами пер­ вой строки матрицы-произведения (18). Взяв алгебраические дополне­ ния элементов других столбцов Л, как показано в (15) и (16), получим сходные результаты для второй и третьей строк матрицы-произведения

(18).

Как показано в (18), матрица Л -1 (см. (17)) такова, что А - 1А = /. До сих пор мы исследовали лишь одно свойство матрицы Л-1. Не было показано, ни что ЛЛ-1 также равно /, ни что Л -1 — единственная матрица, для которой А~1А = АА~1= /. Мы уже упоминали эти свой­ ства (свойства (а) и (б)). Далее мы немного остановимся на них, однако сейчас пересмотрим метод получения Л -1.

Начиная изложение, мы рассматривали матрицу Л:

аи

а12

а13

Л = аи

а21

0-23

а31

а32

азз

106

затем сформировали новую матрицу, заменив каждый элемент А его алгебраическим дополнением:

Р а

Pl2

P'13

Р 21

P'22

Раз

Р 31

P'32

P'33

После транспонирования она равна

P'11 P

Pl2 P

P'13 P

2 1

P 3 1

(21)

22

P.32

2 3

Рзз

 

Умножив на нее скаляр ^ , получили матрицу, обратную к А:

P ll

P 2 1

P 31

Pl2

P 2 2

Р з2

Pl3

P 2 3 '

Рзз

. Умножение А слева на эту матрицу дает единичную матрицу, посколь­ ку в матричном произведении каждый диагональный элемент равен сумме произведений элементов столбца А на их алгебраические допол­ нения и, следовательно, равен | А | (см. соотношение (11)). Умножение

на скалярную величину т-^-. делает результат равным единице. Кроме

того, недиагональные элементы матричного произведения представляют собой суммы произведений элементов столбцов А на алгебраические дополнения элементов других столбцов, приводимые в определителе | Л |, и поэтому равны нулю (см. соотношение (12)). Отсюда Л _1Л == / (произведение А А ~1 будет рассмотрено далее).

Матрица (21), а именно матрица Л, в которой элементы замещены их алгебраическими дополнениями, а затем произведено транспонирова­ ние, называется присоединенной или иногда взаимной с Л матрицей. Таким образом, обратная матрица Л^1 может быть представлена как

присоединенная к Л матрица,

умноженная на скаляр Д -,.

Пример. Определитель матрицы

 

IА I

 

 

2

5

 

 

2

5

Л = 3

9

равен

3

9 - 18— 15 = 3;

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

имеет вид

 

 

Таким образом, обратная к Л матрица равна:

А~* =

1

[

9

—5'

 

3

3

2 '

107

4.УСЛОВИЯ СУЩЕСТВОВАНИЯ ОБРАТНОЙ МАТРИЦЫ

Впредыдущем параграфе матрица Л -1, удовлетворяющая условию

А"1А = /, была получена путем умножения присоединенной к А мат­

рицы на скалярную величину

Для того чтобы существовала об­

ратная матрица Л -1, матрица

Л должна удовлетворять следующим

двум условиям:

1) Л-1 может существовать только тогда, когда Л — квадратная матрица (см. параграф 2 главы V);

2) Л -1 существует только в том случае, если определитель | Л | не равен нулю. (Если | Л | равен нулю, то скалярный множитель | Л | в вы­ ражении для Л -1 не определен и Л -1 не существует. Следовательно, для того чтобы существовала обратная матрица Л -1, определитель | Л | не должен быть равен нулю.)

Квадратную матрицу называют вырожденной, когда ее определи­ тель равен нулю, и невырожденной, когда ее определитель отличен от нуля. Вырожденность, таким образом, может быть присуща только квадратной матрице, но не может относиться к прямоугольной; кроме того, только невырожденные матрицы имеют обратные к ним матрицы. Так же как для существования произведения матриц необходимо их соответствие, так и для существования обратной матрицы необходимое условие —- ее невырожденность. В обоих случаях необходимое условие не всегда записывается, однако оно всегда должно быть удовлетворено.

5. СВОЙСТВА ОБРАТНОЙ МАТРИЦЫ

Если Л есть квадратная невырожденная матрица, то обратная к ней матрица Л -1 обладает следующими свойствами.

Свойства

1. Обратная матрица перестановочна с Л; оба произведения дают единичную матрицу Л _1Л = А А - 1 = I.

2. Обратная к Л матрица является единственной. 5Л = AS = I тогда и только тогда, когда S = Л -1.

3.Определитель обратной к Л матрицы равен обратной величине определителя матрицы Л: ( Л_1| = j-^-j.

4.Обратная матрица является невырожденной.

5. Обратной матрицей к Л 1 будет матрица Л: (Л 'Щ 1 = Л.

6. Обратная к транспонированной матрица равна транспонирован­ ной обратной матрице:

')-1 = (Л-1)'.

7. Если матрица Л симметрическая, то такой же будет обратная матрица: Л' = Л, (Л-1)' = Л^1.

8. Матрица, обратная к произведению матриц, равна произведению обратных матриц, взятых в обратном порядке при условии, что обрат­ ные матрицы существуют. Если существуют Л ”1 и В ' 1, то

(.АВ)~1 = В-1А-1.

108

о

9.Если А такова, что обратная к ней матрица равна транспони­

рованной матрице А, то

говорят, что А — ортогональная матрица

и А А' = /.

раньше (см. стр. 102) мы назвали свойства­

Первые два свойства

ми (а) и (б). Доказательство этих свойств приведено далее, а подбор

примеров предоставляется самому читателю.

Доказательства

 

1. Как было показано, для А 1 справедливо отношение А - 1А = I.

Это свойство объясняется тем, что элементы матричного произведения А - 1А представляют собой суммы произведений элементов столбцов А на алгебраические дополнения элементов тех же самых и других столб­

цов А.

Точно так же произведение АА~г содержит элементы, которые

равны суммам произведений элементов строк А, умноженных на ал­

гебраические дополнения элементов этих и других строк А. Следова­

тельно,

А А - 1 = I.

2.

Предположим, что А 1 не является единственной обратной к А

матрицей и что существует другая, отличная от Л -1 обратная матрица,

такая,

что ВЛ = I. Тогда, умножая справа обе стороны выражения

5Л =

/

на Л -1, получим

 

 

5ЛЛ' 1 = /Л -1 = Л -1,

а поскольку ЛЛ 1 = /, то SI = А 1. Следовательно, S = Л-1. Таким

образом, Л -1 — единственная обратная к Л матрица.

3

 

и 4. В параграфе 5 главы IV было показано, что если две квадрат­

ные матрицы А и В имеют один и тот же порядок, то | Л | | В \ = | АВ |,

поэтому

 

 

 

 

I л||Л -1|=-"|ЛЛ-1|н/|

---=1

 

и,

таким образом, | Л- 1 1

^ .

 

 

Из

этого

непосредственно

следует невырожденность матрицы

Л-1.

 

5,

6 и 7. Рассмотрим тождество I =

Л -1Л. Умножение его слева

на | Л _1| -1дает результат (Л-1)-1 = Л. Транспонируя затем это

тож­

дество и умножая обе его части слева на (Л')-1, мы получим, что (Л')-1 =

=

(Л-1)'

и если Л'

Л, то А~х = (Л-1)'.

 

 

 

8.

Предположим, что Л и В — квадратные невырожденные матри­

цы, имеющие один и тот же порядок. Тогда можно написать

 

 

В-М -М В - В-1 (Л-М)В = в - Ч в = В 'В = / ,

 

и,

следовательно, умножив это выражение справа на (АВ)-1, получим

(АВ)-1 - В 1А - 1,

поэтому правило инверсии применяется не только при транспониро­ вании произведений, но и при их обращении.

109

I

Соседние файлы в папке книги из ГПНТБ