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

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

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

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

Добавлен: 16.10.2024

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

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

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

дать теми же статистическими свойствами. К тому же величина z = х'Ах, где А — положительно полуопределенная симметрическая матрица r-го ранга, характеризуется распределением о2% с г степенями свобо­ ды. К этим выводам можно прийти в результате следующих рассужде­ ний.

Пусть Е (х) = О, Е (хх') — о2/ и у = U'x, в таком случае матема­ тическое ожидание у имеет следующий вид:

Е(у) = Е (U'x) = U'E (х) = О,

*

аковариационная матрица случайных величин у равна:

Е(уу') = Е (U’xx’U) = U'E (xx')U = a2U'IU = о2/.

Рассмотрим теперь выражение z = х'Ах, где А — положительно полу­ определенная симметрическая матрица ранга г; в таком случае сущест­ вует матрица U, которая удовлетворяет следующим условиям: U'U — = I и U'AU — D, где D — диагональная матрица, содержащая г дей­ ствительных положительных ненулевых диагональных; элементов (обозначим их, к примеру, через %1, ..., Хг). Следовательно, переходя к новым переменным у = U’x, можно записать соотношения:

Г

z х'Ах = y'Dy == У'Ягг/?, i~ l

причем у |— совокупность нормально распределенных случайных вели­ чин с нулевы м математическим ожиданием и ковариационной матрицей, равной о2/. Поэтому, каждая случайная величина у\ характеризует­ ся распределением %2 с одной степенью свободы, т. е. распределением о2хIПоскольку множители Я положительны, z представляет собой сумму г независимых случайных величин Ягу2, каждая из'которых рас­ пределена в соответствии со значениями Ящ2%\. В тех случаях, когда Яга2 = 1 (см. параграф 3 главы XIII), случайная величина х'Ах харак­ теризуется распределением %j.

Упражнения

1. Найдите характеристические корни и характеристические векторы при­ веденных матриц. В каждом случае составьте из характеристических векторов матрицу U и проверьте, справедливо ли равенство U~1AU = D, где D ■— диаго­ нальная матрица, в которой диагональными элементами служат характеристи­ ческие корни исходной матрицы.

1

4

1

—2

3

2

1

0

— 4

 

5

1

3

1

—4

6

 

 

9

15

3

 

 

 

F = 6

10

2

 

 

 

3

5

1

 

. 343


4

— 2

 

о '

 

1

2

1

G= — 2

3 2 ;

н =

1

1 4

0

— 2

 

2

 

2

—4 1

1

 

Г

 

1

4

2'

К = — 2

2

2 ;

L =

4 1 2 у

1

— 2

-

1

 

8

8

1

—9

2

6 "

 

 

“ 2

7

1

М = 2

—9

6

j

 

2

3

8

6

6

7

 

 

9

4

4

2. С помощью значений матриц, приведенных в упражнении 1, покажите, что F = В2 + В. Покажите также, что если X — характеристический корень матрицы В, то X2, + к — характеристический корень матрицы F. Докажите следующее общее утверждение: пусть матрица F представляет собой многочлен от матрицы А, тогда все значения характеристических корней F можно полу­ чить, подставив в этот многочлен соответствующие значения характеристических корней А.

3. Рассмотрим общий вид транспонированной матрицы вероятностей пере­

хода размером 2 X 2 :

 

р>= \ Рп

1~ р22

U — Рп

Р ы . '

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

4. Найдите характеристические корни следующей матрицы:'

 

 

 

1_ _1_

О

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

J_

_1_

j_

 

 

 

 

Т =

2

4

 

 

 

 

 

4

 

 

 

 

 

 

 

_1_

 

 

 

 

 

 

 

2

 

 

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

 

 

 

 

 

 

 

_1_

'

1

\ft+ i

1

_ L

M

\ f e +1

4

,

2

j

2

L4 “

{ 2

)

rpk_

1

 

 

1

 

1

 

4

 

 

2

 

4

 

 

 

 

 

 

 

'

1

’\ft-H

1

 

 

 

4 ■'

,

2

j

2

 

 

 

5. Льюкэс [6 ] предложил теоретическую модель спроса на капитальны блага, в которой предусмотрен механизм приспособления к новым условиям, известный под названием «гибкого акселератора». Рассмотрим упрощенный

244


вариант этой модели. Обозначим через xf желательные размеры капитала, а че­ рез х'% — желательные масштабы применения труда. В таком случае вектор х* характеризует желательные масштабы применения капитала и труда, тогда как вектор хп — размеры действительно используемого в период п капитала и тру­ да. Введем следующее уравнение, характеризующее механизм приспособления производственных ресурсов к новым условиям:

*п+1 — *п = В (х* — хп) или

=

(/ — В) хп +

Вх*,

где В — матрица размером 2 X 2 ,

называемая матрицей

коэффициентов

приспособления. Обратимся теперь к числовому

примеру:

 

0,1

—0,5"

 

 

0,3

0,9J *

 

 

Вычислите характеристические корни матрицы В и проверьте, действительно ли эти корни по своей абсолютной величине меньше 1 (условие, обеспечивающее

сходимость

рассмотренного выше

процесса приспособления к х* при п

оо).

6 .

Розничный

торговец,

продающий два товара, следующим образом регу­

лирует размер своих запасов: он вновь заказывает столько же товара, сколько

в прошлом месяце, плюс некоторое дополнительное количество, которое, по пред­

положению,

пропорционально разности (д.Т — хп< ,),

где xt

— желательный,

с точки зрения торговца, уровень запаса г'-го товара, а хп ,

i — размеры

сде­

ланных в предыдущем месяце заказов на i-й товар. Предположим,что нам заданы

следующие величины:

 

 

 

 

 

 

 

*

Разм еры п оследнего

К оэффициент п ропор-

 

Т о в а р

xi

м есячного за к а за

диональности (коэф ф и ­

 

 

 

 

 

циент п риспособления)

 

 

1

20

15

0,5

 

 

2

50

40

0,1

 

а) Запишите в векторной форме разностное уравнение, которое характе­ ризовало бы размеры заказа (на оба товара) в п-й месяц как функцию заказов, сделанных в предыдущем месяце.

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

в) Предположим, что коэффициент приспособления для первого товара ра­ вен не 0,5, а 1,5. Каковы теперь значения характеристических корней? Будут ли оба корня меньше 1 по своей абсолютной величине? Опишите, что будет проис­ ходить с первым элементом вектора хп по мере того, как п будет неограниченно увеличиваться.

7. Езекиэл [2] при построении паутинообразной модели " использовал простейший принцип рыночного равновесия, который можно представить с по­ мощью системы конечно-разностных уравнений. Рассмотрим, например, сле­ дующую модель (все символы в ней представляют собой скалярные величины):

Dt = — md pt + bd

при

md > 0

и

bd > 0;

S t+1= m s pt + bs

при

ms > 0

и

bs > 0 и Dt+i = St+i,

где Dt — размеры спроса на товар; S; — размеры предложения; pt — уровень рыночной цены. (Берутся значения всех этих переменных за период t, a md, ms, bd и bs представляют собой параметры.) Поскольку в нашей модели

13 Зак. 425

345


S;+i — функция от pt, тем самым задается динамический рыночный меха­ низм, функционирование которого приводит к изменению цен во времени.

а) Составьте конечно-разностное уравнение первого порядка, характери­

зующее зависимость уровня цены от времени.

б)

Представьте pt в виде функции от р0 — начальной рыночной цены.

в)

При каком условии pt будет стремиться к конечному пределу при t ос?

Допустим,

что это условие выполнено. Чему в этом случае равен предел после­

довательности pt?

 

 

 

 

г) Будут ли наблюдаться циклические колебания цен?

д) Пусть даны

= 2,

т&— 1, bd = 4, bs

= 1

и р0 = 2. Представьте

в общем виде pt как функцию времени.

 

 

8 .

В модели

взаимодействия мультипликатора—акселератора Самуэльс

на (см. упражнение

13 главы VI) выводится следующее конечно-разностное

уравнение

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

описывающее изменения

национального дохода

во времени:

 

 

 

 

 

 

yt = (a + аР) y t- i — a$yt- 2

+ gt,

 

где а — предельная склонность к потреблению; Р — коэффициент акселерации;

gt

— государственные расходы

на

протяжении

периода

t.

Предположим, что

а

= 0,92 и Р =

12

g0 =

160

и г/ _ 2 = у~х = г/0 =

2000, то достигается

2 3 ' Если

равновесный уровень национального дохода.

 

 

 

 

 

Теперь предположим,

что

=

176.

 

 

 

 

 

а) Представьте конечно-разностное уравнение второго порядка в виде сис­

темы конечно-разностных

уравнений

первого

порядка,

подставив вместо а ,

Р и g0

приведенные числа; t принимает значения 0,1. (Указание. оф =

0,48.)

 

б)

Достигнет ли yt некоторого предела при t

со ? Если да, то чему равен

этот предел?

 

 

 

 

 

 

как функции

 

 

в)

Решите эту систему уравнений для определения y t

t.

 

9.

Покажите,

что если а =

d,

то

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

векторы матрицы

[ с d ] имеют следующий вид:

~

^ J

h

где

и Х2—характеристи­

ческие

корни. Проверьте, справедливо

ли равенство 2а

= ^-2-

 

 

10.

С помощью итерационной вычислительной процедуры определите (с точ­

ностью

до второго десятичного

знака)

главный характеристический

корень

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

11. Разложите на множители характеристическое уравнение матрицы

a b b

Л = b а b

b b а

взяв значение характеристического корня X = а + 26. Затем проделайте ту же процедуру при значении корня X = а Ь.

" 1 0 0

12. Выразите А = и V 0

в форме произведения А = UDU*1

X У г

1 0 0

0 V 0

0 0 г

346


У к а з а н и е . Н а п о м н и м ,

что

 

 

 

 

1

0

O' 1

1

0

o'

 

a

1

0

— a

1

0

 

b

c

1

ac— b

—c

1

13.

С помощью

решения

упражнения

12

покажите, что матрица вероят

ностей

перехода

 

 

 

 

 

 

 

 

1

О

 

О

 

 

 

1

1

 

О

 

Р =

 

 

 

2

2

 

 

 

1

 

 

 

i

 

 

— —с (1 —с) 2с (1 —с)с (1 с)

возведенная в k-ю степень, равна:

 

 

 

О

О

 

Р* =

1 -

- F

°

 

2

 

1 +

1 N*1 / 1 \А 1

 

 

г6 -

2zk zk

где z =

— + с(с — 1).

 

 

 

14.

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

равенство

A kU i= kki ui будет

выполняться при любом

целом отрицательном k, если существует А -1 .

15. а) В упражнении 3 была приведена транспонированная матрица ве­ роятностей перехода размером 2 X 2 .

1 — P2Z

Ра .

С помощью решения упражнения 3 найдите диагональную матрицу D, в которой диагональными элементами являются характеристические корни Р ', а также соответствующую матрицу U, удовлетворяющую условию Р' = UDU~х.

б) Полагая, что все элементы матрицы Р положительны и используя вы­

полненное требование а, покажите, что при п

оо Рп стремится к конечному

пределу. Какими

свойствами обладает предельное значение матрицы?

 

 

Л И Т Е Р А Т У Р А

 

1. C h o w G.

С.

(1957). Demand for Automobiles in the United States: A Stu­

dy in Consumer Durables. North-Holland Publishing Company, Amsterdam.

2. E z e k i e l

M.

(1938). The cobweb theorem. Quarterly Journal of. Econo­

mics, 52, 255—280.

 

 

13*

3 4 7