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

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

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

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

над скалярными

величинами не всегда сохраняются при действиях

с матрицами; эти

причины легче понять после того, как по аналогии

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

6.ЗАКОНЫ АЛГЕБРЫ

Вэтом параграфе содержится формальный анализ следующего вопроса: насколько применимы к операциям сложения и умножения матриц ассоциативный, коммутативный и дистрибутивный законы алгебры?

а) Сложение матриц ассоциативно -в том случае, если матрицы согласованы для сложения. Известно, что если матрицы Л, В и С одного и того же порядка, тогда

+ В) + С = {аи + bi}} + {Cij} =

= {агj + bi} + Cjj} = А + В -f- С.

Точно таким же образом можно показать, что

(ао + Ьи + си } = {аи } + {Ьи + си } =

А + (В + С).

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

+ В) + С = А + + С) = А + В + С.

б) Операция умножения матриц также удовлетворяет требованиям ассоциативного закона, если только матрицы согласованы для умно­ жения. Если Л представляет собой матрицу размером рХд, В —матри­ цу размером q X г, а С — матрицу размером г X s, то

(ЛВ)С - j х ( ьлV )„<:„„} = [ з щД,Ь,„С„ } =

- J i a o ( i 6л ам )}-/Н 8 С ),

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

{АВ) С = Л (ВС) = АВС.

в) Операции с матрицами удовлетворяют также требованиям ди­ стрибутивного закона, согласно которым

Л + С) = АВ + АС

4 0

в тех случаях, когда матрицы В и С согласованы для сложения (т. е. у них обязательно должны быть одинаковые размеры) и матрицы А я В (а следовательно, также А и С) согласованы для умножения. Предпо­ ложим, что А представляет собой матрицу размером р X q, а размер обеих матриц В и С q X г, тогда

А{В-\-С) =

аи (b}k + c}k)! =

(

Я

я

~ |

aij bjk +

^ Яц Cjhj ^- АВ + АС.

U=i

/=1

г) Сложение матриц коммутативно в том случае, когда матрицы согласованы для сложения. Если размер матриц А я В одинаков, тогда

А + В = {аи + Ьи } = {bti + аи } = В + А.

д) Умножение матриц, вообще говоря, не коммутативно, другими словами, АВ не всегда равно ВА. Мы уже видели, что, перемножая матрицы А я В, можно получить два произведения АВ я ВА, я если размеры матрицы А равны г X с, то оба произведения существуют только в том случае, когда размеры матрицы В равны с X г. Тогда про­ изведение АВ образует квадратную матрицу порядка г, а произведение ВА — квадратную матрицу порядка с. Поэтому размеры АВ могут быть равны В А в том случае, когда г — с, другими словами, когда обе мат­ рицы А и В квадратные и имеют один и тот же порядок, равный г. Тогда соответствующие произведения будут составлять

АВ = j 2 aik bhjj и BA = J bth ahjj при t, / = 1, 2.......

r.

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

Из этого, понятно, отнюдь не'следует, что АВ и ВА всегда должны различаться между собой, в отдельных случаях они могут быть равны. Например,

"3

2'

1

2 '

2

3

2

— 1

'1

4 '

 

см

XT'

 

см

-----

 

 

 

1

'3

2 ‘

2 3

Вместе с тем несоблюдение закона коммутативности при умножении матриц можно продемонстрировать с помощью следующего примера:

А 2 “ 0 — Г

2

— 3 ‘

0

— Г

1 2"

— 3

— 4"

 

 

— 7 . Ф

 

 

 

 

3 4 . _1 — 1.

4

J

— 1.

3 4.

— 2

— 2,

41

В двух случаях, имеющих особо важное значение, умножение мат­ риц обладает свойством коммутативности:

1) в случае умножения на нулевую матрицу: если А р представляет собой квадратную матрицу р-го порядка, а Ор — аналогичную матри­ цу, все элементы которой составляют нули, тогда

0 рА р — А р0 р — 0 Р\

2) в случае умножения на диагональную матрицу, все диагональные элементы которой равны единице; такая матрица называется тождест­ венной или в ряде случаев единичной матрицей и обычно обозначается буквой /; иногда это обозначение снабжается индексом, указывающим порядок матрицы, таким образом,

1

0

0'

/ з - 0

1

0

0

0

1

В общем случае, если порядок матрицы А р равен р, то точно так же, как нулевая матрица выполняет роль нуля, единичная матрица того же порядка служит «единицей» в матричной алгебре. Например,

1

0'

' 1

2

—4

1

2

—4'

0 1 _ 9 7

2 _

9

7

2_

Матрицу вида Я/, где Я

 

скалярная величина, называют иногда ска-

лярной матрицей. Так,

матрица 41 =

"4

01

* „

q

4

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

лярную матрицу.

ArXs,

умноженная

слева

или справа на нулевую

Любая матрица

матрицу соответствующего размера, дает нулевую матрицу. Поэтому,

если ОсХг — нулевая

матрица

размера

с X г, то

 

Ocxr Arxs —О,cXs

Аналогичным образом находим

 

 

Arxs Osxp

 

Например,

 

 

 

 

 

1

2

—4"

[0 0 0].

[0

0] 9

7

2

Существует еще один случай, когда при умножении матриц спра­ ведлив коммутативный закон: случай умножения матрицы на скаляр­ ную величину, скажем, на Я:

ЯА = {Яп^} = {<%Я} — ЛЯ.

4 2

7.ВЫВОДЫ

Врезультате рассмотрения законов матричной алгебры становится понятным, насколько важно перед выполнением операции сложения, вычитания и умножения выяснить, согласованы ли матрицы для этих операций.На практике всегда следует так поступать — сначала должен исследоваться вопрос о согласованности матриц, даже в тех случаях, когда это не оговаривается особо. Упоминая произведение АВ, мы в ред­ ких случаях специально оговариваем: «АВ, если А я В согласованы для умножения», однако такая согласованность предполагается, и в каждом случае нужно суметь показать согласованность матриц, умно­ жаемых друг на друга.

Упражнения

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

'17

4

12'

 

 

'2 0

5

10'

 

6

4

13

по магазину 1

,

10

5

15

по магазину 2

11

4

8

20

5

8

 

 

 

- 7

4

6 -

 

 

ДО 5

10-

 

П2 3 4 '

8 3 4 по магазину 3.

10 3 6

. 4 3 7.

а) Покажите, что в каждый сезон 1-й и 3-й магазины, вместе взятые, прода ли больше каждого вида товаров, чем магазин 2 .

б) Покажите, что общие продажи всех трех магазинов составили:

‘49 12 26' 24 12 32

4112 22

,1_2 1 . 12 23.

2.Данные о продукции добывающей промышленности по определенным видам минерального сырья группируются обычно также по странам, являющим­ ся основными поставщиками этого сырья. Пусть по строкам группируются дан­

ные о добыче в странах 1, 2 и 3, а по

столбцам — о минеральном

сырье вида 1,

2 и 3. Тогда добыча минерального

сырья (тыс.

т) в

1960

и

1965 гг. может

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

 

 

 

"

450

780

210'

'

520

910

220'

 

 

А =

1 050

240

90

 

и В =

1 080

580

290

 

 

 

1 500

120

590

 

 

1 460

830

600

 

 

а) Рассчитайте матрицу приростов

добычи

за период с 1960 по 1965 г.

б) Проведя соответствующие арифметические

действия

над матрицами,

43

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

' 485 845 215' 1 065 410 190 . _ 1480 475 595_

3. Покажите, что

г6 31

1

о

СЛ

 

1

 

 

ГЗ

8

 

to

 

-------1

4^

+

 

1

 

1

 

0

'

12

17"

3

6

"

11

2

—13 —4

 

 

3

9

2

12

 

б)

1

+ 5

1 8'

0

—2 5

 

 

4.

В одном из разделов теории игр,

рассматривающем игру одного человека

против природы, используются так называемые платежные матрицы; в заголов­ ках строк такой матрицы указываются возможные действия человека, а в заго­ ловках столбцов —возможные результаты этих действий (степень успеха или не­ успеха), тогда элементы матрицы будут характеризовать размеры соответствую­ щих платежей. При исследовании результатов подобных игр пользуются мини­ максным критерием возможных потерь; такие потери описываются с помощью особой матрицы, исчисляемой на основе платежной. Сначала составляется мат­ рица, в которой каждый элемент /-го столбца равен наибольшему элементу со­ ответствующего столбца платежной матрицы. Затем из этой матрицы вычитается платежная матрица, полученную разность называют матрицей возможных по­

терь (regret matrix). Ее элементы представля'ют

собой максимальные

потери,

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

избрано нелучшее решение.

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

стратегия: мы выбираем действия, соответствующие той строке матрицы

потерь,

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

а) Определите матрицу возможных потерь на основе приведенной платеж­ ной матрицы.

б) Найдите минимаксную стратегию.

~ 0

2

100

0

 

7

8

100

2

 

10

40

1

90

'

_ 1

3

1

1

-

5.В другом разделе теории игр исследуются игры двух лиц, в этом случае

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

вого игрока, а в заголовках

столбцов — действия второго игрока. Например,

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

 

 

10

1

8

4

18

12

Р = 6

7

17 и Q= 8 8 14

22

4

ю

.1

10

13.

а) Как с помощью совместных действий обоих участников игры добиться наибольшего выигрыша, выплачиваемого первому игроку?

б) Как с помощью совместных действий обоих участников добиться наи­ большего выигрыша для второго игрока?

44

в) Допустим, что участники игры вступают в коалицию, как с помощью совместных действий можно максимизировать сумму их общего выигрыша? Как выглядит в этом случае их общая платежная матрица?

г) Какие совместные действия предпочел бы первый участник игры,

если

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

выиг­

рышем второго игрока?

 

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

издержек перевозки (выраженных в долларах)

будет равен 3 x 4 , в заголовках

строк

будут указаны различные заводы,

а

в заголовках столбцов — места на­

значения.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5 =

'10

 

15

9

 

Т

 

 

 

 

 

 

 

 

 

 

 

 

 

14

8

12

 

8 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

14

22

17

 

 

 

 

При

этом

издержки

 

производства

(доллары)

на

первом

заводе равны 40

на втором — 38,

 

на

третьем — 41.

 

 

производства

размером 3 x 4 ,

элементы

а)

Напишите матрицу

Р издержек

которой группируются

по строкам и столбцам так же,

как и в S.

на про­

б)

Определите матрицу

совокупных

издержек,

объединяя

затраты

изводство и транспортирование товаров.

 

 

 

 

 

 

7.

Пусть

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

“3

6 '

 

 

 

1

0

 

3 2'

 

'

Г

Г

 

 

 

А =

,

в =

 

 

 

1

 

 

 

2 1.

0 — 1 — 1 1_

 

 

0

и (/ =

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 .

 

 

 

Покажите, что

—6 3

12

 

 

 

 

 

 

 

 

 

 

 

 

а)

АВ =

3

 

 

 

 

 

 

 

 

 

 

 

 

2 —1 5 5

 

 

' —Г

 

 

 

 

"1

0 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б)

Вх =

 

 

 

В 'В х =

. 2

,

где

В '

-

0 — 1

 

 

 

 

 

 

 

 

 

 

 

 

— 1

 

 

 

 

 

3 — 1

 

 

 

 

 

 

 

 

 

 

 

 

—4-

 

 

 

 

L.2

1J

 

 

 

в)

А у =

 

 

,

А'Ау--

 

—7

,

где

А ’ -

3 2

 

 

 

 

 

 

 

—17

6 1

 

 

 

 

г)

А2 — 4А — 91= 0;

 

 

 

 

 

 

 

 

 

 

 

 

 

д)

{ Л

— 1

 

6

 

1 1

6 '

Л_ 1 0 '

 

 

 

 

 

2 —3 = 9

 

 

Л =

0 1.

 

 

 

 

 

 

 

 

2 -- 3

 

 

 

 

 

 

 

8 .

Убедитесь,

что

 

 

 

 

 

 

 

 

 

 

 

 

 

а)

1 г

2 3

 

2 3'

1 Г

 

 

 

 

 

 

 

 

 

 

1 1 3 2

 

3 2 .1 1.

 

 

 

 

 

 

 

 

 

 

б)

 

 

 

1

1

1

1

,

то А2 =А;

 

 

 

 

 

 

 

если А =

1

1

1

 

 

 

 

 

 

 

 

 

 

 

 

1

1

1

 

 

 

 

 

 

 

 

 

 

 

 

в)

если

С

6

- 4

 

, то С2= 0;

 

 

 

 

 

 

 

 

 

9 — 6

 

 

 

 

 

 

 

 

 

45

 

1 ~4 —5 - - Г

'

1

1

Г

 

 

 

 

 

 

 

 

г) -Г 1

1

2

 

1

0

1 = / 3.

 

 

 

 

 

 

 

9 _4

4 - -1_

 

0

4 - L

 

 

 

 

 

 

 

 

 

 

 

-

 

2

3"

 

 

' 6

 

0

0

 

 

 

9.

Пусть

Х =

 

0 — 1 — 2

и Y =

—3

 

4

0

, найдитеX2, К2, XY, YX

 

 

 

 

 

0

7

 

 

0

—5

2

 

 

 

и покажите,

что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

40

 

5

44"

 

 

(Х + К)2= Х2 + Х К + У Х + К 2=

—28

13 —33 .

 

 

 

 

 

 

 

 

 

 

 

 

1 —62

88

 

 

 

1 0 2 '

 

"1

3

0 '

 

 

'6 5 Т

 

10.

Дано

Л =

0

1

1 ,

в =

0

4

1

,

х

= 2 2 4

 

 

 

 

2

0

2

 

2

3

0

 

 

3 3

6 .

 

покажите, что АХ = ВХ, хотя A=f=B.

11. Пусть А = {aij} при г, / = 1, 2, ... , г,

D =

{d^}

при i, } —

1 ,2 ........ ,

г; d*7- = 0 для

г =£ /';

покажите, что ЛД =

{d;jajj} и

ДЛ = {da

ац] при t, / =

1, 2....... г.

12. Бирман

[1]

исследовал процесс рефинансирования

в условиях, когда

в соответствии с принятыми допущениями вероятные изменения процентных ставок заданы, например, такой матрицей вероятностей перехода:

0

, 6

0 ,2

0 ,2 '

Р = 0,3

0,5

0 ,2

.0

,1

0,3

0 ,6

в заголовках строк и столбцов указываются процентные ставки соответственно

5, 6 , 7%.

а) Текущий вектор вероятностей состояния для каждой из процентных

ставок

в период п х'п определяется в результате умножения справа вектора ве­

роятностей

в период 1 х[ на Р п— Допустим, что начальная процентная ставка

равна 6 %,

тогда х\ будет составлять (0 1 0]. Как будут выглядеть векторы ве­

роятностей для двух последующих периодов?

б)

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

ющих соотношений:

х'а=х'2Р = х(Р2.

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

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

особого вида. Вычислите, чему

равны

х4, х'ъ и Х5, и решите, каким мог бы быть

вид этого вектора.

 

 

 

 

г) Какой вид примет этот вектор,

если порядок матрицы равен 10?

13. Даны следующие две матрицы вероятностей перехода:

Р2 %РЯ

Я2

Р

Я

о

Р = Р2 2р<7

Я2 и Т =

2

2

_р2 2pq

q2_

2

- 0

р

q _

46

П о к а ж и т е ,

что

 

 

а)

Р Т = Р = ТР,

б) Т2 = -—■(Р + Т),

в) Р2= Р,

 

( 1

\ n — l

 

г) Тп = Р +

(Т— Р)

 

 

 

\ 2 /

и если

1

1

1

 

S = — / + — Т + — Р, то

 

 

4

2

4

 

Д) ST = TS = Т2 и е) S2 = — (7+6Т + 9Р). 16

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

Размеры спроса на билеты . . . . .

 

1

2

3

4

Вероятности...........................................

0,3

0,4

0,2

0,1

Пусть матрица условных издержек составляет

 

 

 

 

■ 3

3

3

 

3~

 

 

 

— 2

6

6

 

6

 

 

 

—7

1

9

 

9

 

 

 

_— 12

—4

4

12.

 

 

 

Для того чтобы вычислить ожидаемый уровень издержек при каждом из воз­ можных решений, умножьте эту матрицу на вектор вероятностей. Будут ли раз­ личаться между собой результаты в зависимости от того, как мы будем умножать матрицу на вектор: справа или слева?

15. Рассматривается проект вложения капитала сроком на три года, кото­

рый должен обеспечивать следующую денежную

выручку:

в первый год —

1000 долларов, во второй — 2000 долларов,

в третий

10 000 долларов.

Допустим, что проект будет принят в том случае, если совокупный доход от капи­ таловложений (в пересчете на сегодняшний доход) будет превышать требую­ щиеся затраты, составляющие 5000 долларов, причем дисконтирование ожидае­ мого дохода проводится по годовой ставке 10%. При такой ставке дисконтирова­ ния доход, который будет получен на протяжении первого года, должен быть ум­ ножен на 0,9091, на протяжении второго года — на 0,8264 и на протяжении третьего — на 0,7513. Запишите денежную выручку и дисконтирующие множи­ тели в векторной форме; перемножив их между собой, определите, чему равен дисконтированный совокупный доход от капиталовложений. Будет ли принят рассматриваемый проект?

16. а) Рассмотрим линейную производственную функцию

Qt Ьххц +Ь2х2(,

где хи и X2t обозначают соответствующие затраты за период 7; bt и Ь2 — техно­ логические коэффициенты; qt — выпуск продукции на протяжении периода t. Допустим, что вектор-строка q' = {qt} характеризует выпуск продукции на про­

тяжении соответствующих

лет

t-= 1, ..., п; в

матрице затрат X

—■ {хц}

при i =

1, 2 по строкам

группируются данные

о различных видах

затрат,

а по

столбцам — данные

за

соответствующие

годы; вектор-строка

b’ — Ibt)

характеризует

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

Представьте в

матричной форме зависимость q от X и Ь.

 

47

б ) П р е д п о л о ж и м , ч то п р и п = 3 мы р а с п о л а г а е м с л е д у ю щ и м и д а н н ы м и

320 360 290

и Ь' = [3 5].

700 780 800

Определите вектор выпуска продукции.

 

 

 

 

 

Л И Т Е Р А Т У Р А

 

1. B i c r m a n

Н.

Jr.

(1966). The bond refunding decision as a Markov

process. Management

Science, 12,

545—555.

 

2.

D e m i n g

W.

E.

and

G 1 a s s e r G. J. (1968). A Markovian analy­

sis of the life of newspaper

subscriptions. Management Science, 14, B, 283—293.

 

3.

D e r m a n

C. (1963). Optimal replacement and maintenance under

Markovian deterioration

with probability bounds on failure. Management Science,

9,

478—481.

 

 

 

 

N1

ДРУГИЕ ДЕЙСТВИЯ

ГЛАВА

С МАТРИЦАМИ

1.ЛИНЕЙНЫЕ ПРЕОБРАЗОВАНИЯ

Впредыдущей главе мы уже видели, что произведение матрицы

Аи вектора х само является вектором. Обозначим этот вектор-произве­ дение через у. Тогда у —Ах. Из самого определения процесса умноже­ ния следует, что каждый элемент произведения есть сумма элементов вектора х, каждый из которых умножен на соответствующий элемент Л. Например,

Ух

ац

ап

а13 ]

л2

=у = Ах

а21

 

 

У2 .

а22 ^23 у

дает

 

 

 

л3

 

 

 

 

Ух ~ а11х1 + UX2X2

<213X3

 

и

 

 

 

 

У2 — Я21Х1 &22х2 “Ь <223X3-

Таким образом, если не все члены ап , а12 и а13 равны нулю, то элемент ух равен взвешенной сумме элементов х с весами, равными элементам первой строки матрицы А. Такая взвешенная сумма называется' ли­ нейной комбинацией элементов х\ она линейна, поскольку элементы х не имеют иной степени, кроме первой, и нет членов, содержащих более одного элемента вектора х. Следовательно, в произведении у = Ах элементы у являются линейными комбинациями х-ов; иными словами, вектор х был преобразован в вектор у с помощью умножения. Говорят, что в такой ситуации матрица А представляет собой линейное преоб­ разование х в у. Она и есть то операционное средство, с помощью ко­ торого элементы х преобразуются в элементы у.

Необходимость в линейном преобразовании одного ряда перемен­ ных в другой (или одного вектора в другой) возникает в различных обстоятельствах. В качестве иллюстрации воспользуемся задачей оп­ ределения производственных потребностей.

49

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