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

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

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

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

Добавлен: 16.10.2024

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

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

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

-

д

~

 

' ^

а ц Х ] +

^ a ilXi

 

дхг

 

 

 

 

 

/

 

1

 

д

 

 

1

a-ijXj - \~ v .a i2Xi

 

дх2

 

 

 

 

 

i

 

i

 

д

i

I

 

 

 

 

дх3

 

 

 

 

 

 

 

i

®njXj “Ь

@in%i

 

 

 

 

 

1

 

 

 

 

_

 

_

 

 

 

~a[ x

a[x

 

 

 

 

a'ix

-!- a'2x

 

 

 

 

 

- •

 

 

 

 

a'n x

Я

 

 

 

 

 

3 1

 

где al и ai представляют

собой соответственно г'-ю строку и г'-й стол­

бец матрицы А.

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

векторы равны Ах и А'х, и в таком

случае

 

 

 

 

 

 

— (х' Ах) = Ах + А' х.

дх

 

Если А — симметрическая матрица, А — А', тогда

 

(х' Ах) = 2Ах.

(5 )

дх

 

Пример 1. Дана квадратичная форма

 

 

 

1

3

5

X i

 

 

у = x' Ax = [Xx

-------l СЛ CO

 

7

(N ____

1-------

= х \ +

 

CO 1

x2 x3]

 

*

 

 

 

 

9

4

 

 

"Е Юх! х3 -(-4x2

~f~ 14х2 Xg -f- 9x1;

i

l

 

 

 

 

 

 

2xx -f 6x2 -f- 10xg

6 х х Х 2 +

ду_

dy

6xx -f 8x2 -f 14Xg

2Ах.

дх

dx%

 

 

 

}

lOxx-p 14x2-f 18Xg

 

 

« *

 

 

 

1

1

 

 

Пример 2. Применение в статистике метода наименьших квадра­ тов предполагает минимизацию суммы квадратов элементов вектора

е, фигурирующего в следующем уравнении:

 

у = ХЬ + е.

(6)

361


Пусть вектор е = {ег} при i = 1, 2, п\ тогда сумма квадратов

S = V е\ ^ е ' е = ( у — ХЬ)'{у — ХЬ) = i = 1

у' у — Ь' X' у — у' Xb -f- b' X' ХЬ .

Поскольку произведение у'ХЬ — скалярная величина, оно не меняется при транспонировании и равно Ь'Ху, следовательно,

S = у'у — 2Ь’Х'у + Ь'Х'ХЪ.

Найдем такие значения элементов вектора Ь, при которых величина S

окажется наименьшей; для этого приравняем к нулю выражение

где под Щ подразумевается описанный вектор дифференциальных опе­

раторов. Тогда с помощью соотношений (1) и (5) можно прийти к сле­ дующему результату:

— = 2 Х 'у + 2Х' ХЬ.

(7)

дЬ

 

Приравнивая это выражение к нулевому вектору, получаем

 

b = { Х ' Х ^ Х ' у

(8)

при условии, что (Х'Х)~1 существует.

 

6. ЯКОБИАНЫ

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

\ ~ \ при t, / —1, 2,..., п I dxj )

может служить матрица, использовавшаяся в уравнении (4). Опреде­ литель такой матрицы называется якобианом; обозначим якобиан бук­ вой J. Описанные в предыдущем параграфе методы дифференцирова­ ния можно применить к любой задаче, в которой переменные у пред­ ставляют собой функции от других переменных х; получаемый в таком случае определитель матрицы первых частных производных называют якобианом функций у%по переменным xj. Допустим, что между пере­ менными у и х существует линейная зависимость, которая может быть выражена формулой у — Тх. В таком случае из уравнения (4) сле­ дует, что J = | Т' |, иначе говоря, якобиан — это определитель тран­ спонированной матрицы Т, или, что то же самое, определитель матри­ цы Т '.

362


Якобианы широко применяются также при интегральном исчисле­ нии в операциях преобразования переменных.

Допустим, что под интегралом стоят дифференциалы dy1> dy2, ....

dyn, тогда переход от переменных уг, у 2, ..., уп к новым переменным xlt х2, ■■■, хп осуществляется следующим образом: вместо у подставля­ ются соответствующие функции от х, а дифференциалы dyt заменяются на | / | dxi, dx2, .... dxn, где J — якобиан функций у по переменным х,

a J — его

абсолютная величина.

 

 

Другой пример матрицы частных производных — это матрица Гесса,

или гессиан. Пусть

у — функция п переменных хх, х2, ..., х„; диф­

ференцируя у по X; (i

= 1, ..., п),

получим матрицу частных производ­

ных второго порядка, которая называется гессианом:

 

н

~ {/ * 1 - ) ПРИ г, /= 1 ,2,..., п.

 

 

{dxi dxj j

 

 

Пример.

Рассматривая задачу

максимизации

потребителем полез­

ности при

наличии

бюджетных

ограничений,

Хендерсон и Квандт

[3] приводят достаточные условия существования максимума функции полезности, причем критерий, указываемый ими, предполагает рас­ смотрение главных миноров окаймленной матрицы Гесса. В книге Самуэльсона [5, глава 6] якобиан применяется в самых различных преобразованиях в моделях потребительского выбора и максимиза­ ции прибыли предпринимателями *.

7.МАТРИЧНЫЕ ФУНКЦИИ

Вряде случаев имеет смысл ввести матричные функции, строящиеся по аналогии со скалярными алгебраическими функциями. Рассмотрим, например, матричную функцию ек :

eK = / + /( + fr_

7+ г2= 1

к г

3!

г\

напоминающую разложение в ряд

скалярной

экспоненциальной

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

представлять собой матрицу того же порядка, что и

К , выраженную

через степенной ряд от К.

зависимостью:

Другая матричная функция может быть задана

tf = loge P ,

(9)

где R и Р обозначают матрицы. Для того чтобы представить R в виде степенного ряда, снова надо обратиться к алгебре скалярных величин:

1о§е х = log [ 1 (1 х)] = -

2

оо

1/(2г — 1).

(1 —х)2г'

 

г

= 1

 

*Описание математического аппарата, применяемого при решении таких

задач, приведено, например, в работе:

Л а н к а с т е р

К- Математическая

экономика. М., «Советское радио», 1972,

гл. 4, — Прим,

перев.

363


Однако это

выражение

справедливо лишь при некоторых значениях

х, поэтому

вместо (9)

предпочтительнее

эквивалентное соотношение

 

 

оо

Rk

 

 

k= i

r -

которое дает возможность представить Р в форме ряда, как ранее ек . Это соотношение позволяет установить определенную зависимость между элементами матриц Р и R:

 

 

Ри ■6„ ■S

ГФ

 

 

и

 

 

k—l

k\

где Р = { р ц } ,

=

1, если i — } и 6;;- =

0 при i Ф /, а гц представляет

собой ij-й элемент

матрицы R k. В некоторых случаях это выражение

можно привести к сравнительно простой форме (см. упражнение 10).

8. ПРЯМАЯ СУММА МАТРИЦ

 

Прямая сумма двух матриц

определяется следующим образом:

Л 0 В =

А 0"

 

0 ВJ

где нули обозначают нулевые матрицы соответствующего порядка.

Такое определение применимо

как

в тех случаях, когда размер мат­

риц А и В одинаков,

так и

тогда,

когда их размеры различны. На­

пример,

 

 

 

 

 

 

 

[1 2

3 ]0

6

7

1

2

30 0

 

0

0

06 7

 

 

8

9

0

0

08 9

 

 

 

 

 

 

Транспонированная матрица прямой суммы совпадает с прямой сум­ мой транспонированных матриц. Ранг прямой суммы равен сумме рангов матриц-слагаемых; это следует из определения ранга матрицы (см. параграф 4 главы VI). Из определения прямой суммы матриц мож­ но заключить, что А 0 (—А) ф 0 (если, конечно, матрица А не пред­ ставляет собой нулевую матрицу). Соотношение

0 В) + 0 D) = (Л + С) 0 + D)

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

(Л 0 В) (C@D) — АС ®BD

364


справедливо лишь в тех случаях, когда выполняется условие согласо­ вания матриц для умножения. Допустим, что матрицы Л и В невырож­ дены; тогда приведенное соотношение позволяет заключить, что

ф В)*1 — А~

| В

Прямая сумма Л ©В

представляет

собой квадратную матрицу,

когда либо обе матрицы А

и В квадратны, либо общая сумма строк двух

матриц равна сумме их столбцов. В том случае, когда А я В — квадрат­ ные матрицы, определитель А ©В равен | А | ] В |, в остальных случаях определитель А ф В равен нулю (можно убедиться в этом, подверг­ нув определитель | А ® В \ разложению Лапласа). Если А я В квад­ ратные матрицы, характеристическими корнями их прямой суммы бу­ дут характеристические корни матриц Л и В, в этом можно убедиться, выписав характеристическое уравнение матрицы Л ©В. Допустим, что и и v соответственно обозначают характеристические векторы матриц А я В, тогда характеристические векторы матрицы Л © В имеют сле­

дующий вид: q и yj. Однако если матрицы Л и В не квадратны, а

матрица Л ф В квадратна, заранее можно утверждать лишь то, что по крайней мере d характеристических корней Л © В будут равны нулю, где d — разность между числом строк и столбцов матрицы Л (или матрицы В). Характеристические векторы прямой суммы матриц рассчитываются тем же способом, что и для любой матрицы.

9. ПРЯМОЕ ПРОИЗВЕДЕНИЕ МАТРИЦ

Прямое произведение матриц ApXq я ВтХп определяется следующим образом:

ап В ... alqB

Лп X Q*BtnХП --

(10)

api В ...

apq В

Из этого определения следует, что прямое произведение можно рас­ членить на подматрицы соответственно числу элементов, содержащих­ ся в матрице Л, причем каждая из этих подматриц будет произведени­ ем матрицы В на соответствующий элемент Л. Элементы прямого произведения состоят из произведений всех возможных пар такого ви­ да: элементы матрицы Л умножаются на матрицу В, а размеры пря­ мого произведения равны pm X qn. Например,

[1 2 3]* '6 г ‘6 7 12 14 18 21" 8 9 8 9 16 18 24 27

Транспонированная матрица прямого произведения совпадает с пря­ мым произведением транспонированных матриц-множителей; в спра­ ведливости этого вывода можно убедиться, транспонировав матрицы, приведенные в определении (10). Заметим, однако, что (Л * В )'= А' *В’

365