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

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

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

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

Добавлен: 16.10.2024

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

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

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

Теорема 2 (теорема 4.6 в книге Грейбилла). Допустим, что х пред­ ставляет собой вектор случайных независимых переменных, характе­ ризующихся нормальным распределением, нулевым математическим ожиданием и единичной дисперсией; в таком случае квадратичная форма х'Ах характеризуется распределением %2 с г степенями сво­ боды (при условии, что А — идемпотентная симметрическая матрица ранга г).

Доказательство. В разделе г параграфа? главы XII было показано, что квадратичную форму х'Ах {А представляет собой симметрическую

 

 

 

Г

А;г/-, где Агг/?—независи-

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

мые случайные величины,

 

i =

1

распределенные как Агсг2%2, причем все мно­

жители А,г положительны.

В соответствии с условиями теоремы пере­

менные х имеют единичную дисперсию (а2 =

1) и матрица А идемпо-

тентна;

но из этого следует, что

Аг = 1

(см. свойство 5 идемпо-

тентных

матриц). Следовательно,

2Агг/2 — сумма независимых вели­

чин, каждая из которых распределена по закону %?, так что их сумма характеризуется распределением %2.

Мы можем воспользоваться теоремой 2, обратившись к рассмотре­ нию суммы квадратов 55. Допустим, например, что средние значения переменных равны нулю, ковариационная матрица имеет вид о2/, и кроме того, е — вектор нормально распределенных величин; в таком случае в соответствии с теоремой 2 можно прийти к выводу, что вели­ чина 55 распределена как о2% при п — 1 степенях свободы.

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

4. НИЛЬПОТЕНТНЫЕ МАТРИЦЫ

Матрица А иильпотентна в том случае, если А р — 0, а А р~ х Ф 0; р обычно называют индексом нильпотентности.

Пример.

1 2 5 А = 2 4 10 ,

—1 —2 5

о

о

О

А2 = 0

0

0

0

0

0

Матрица А иильпотентна, индекс нильпотентности равен 2. Существуют также матрицы, удовлетворяющие условию А р = /.

357



Например, рассмотрим

 

 

 

 

0

0

6

А

1

о

о

2

 

 

 

о

— о

в таком случае А 3 = I . Возможно,

таким матрицам больше подошло

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

’ Lo — '

Независимо от того, каков вид матрицы В, имеют место следующие соотношения: А 2 = I и А = Л -1.

5. ВЕКТОР ДИФФЕРЕНЦИАЛЬНЫХ ОПЕРАТОРОВ

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

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

Я = 3xi 4х2 “Ь 9х3,

где Я, х1г х2, х 3— скалярные переменные; тогда частные производные функции Я по хх, х 2 и х 3 можно записать в следующем виде:

дХ

д

 

дх1

дхг

3

дХ

д

4

дх2

дх2

9

дХ

д

 

_ дх3_

дх3 _

 

*«uni-potent» ■— сочетание искаженных латинских слов «единица» и

«степень». ■— П р и м , п е р е в .

358


Второй из приведенных векторов показывает, что символ можно

представить вектором дифференциальных операторов. Заметим, что

Xi

X —Зхг + 4 х 2 -f 9х3= [3 4 9] х2 ■а х .

Определив таким образом векторы а и х и подразумевая под ^ век

тор-столбец дифференциальных операторов, можно показать, что ^

равен а, т. е.

■j-(a'x) = -£-(x'a) = a.

(1)

ох дх

Предположим теперь, что такая вычислительная процедура приме­ няется поочередно к каждому элементу вектора у = Ах, имеющего сле­ дующий вид:

1

"а[ х

У2

а'2 х

 

(2)

_Уп_

_а'пх_

где г/; — элементы вектора у и а[,

а2, ..., а'п соответствуют строкам

матрицы А. Теперь вектор операторов ~ будет относиться к соответ­

ствующим элементам вектора у. В результате дифференцирования каж­

дого

элемента получаем

вектор-столбец

частных

производных:

(г =

1, 2, ..., п). Предположим,

что все п значений этого вектора мы

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

матрицу.

Тогда окажется,

что выражение

где у — вектор

вида Ах,

образует матрицу А'.

Действительно,

 

 

 

 

 

 

 

 

 

ду

__ ду г

ду2

дуп

 

 

 

 

дх

дх

дх

 

дх

 

 

 

да[х да2х

да'пх - из соотношения (2)

 

дх

дх

дх

 

 

 

 

 

 

= [а1

а2...ап] из соотношения

(1) =

(3 )

 

= А',

 

 

 

 

 

 

 

так как а[ — это, по определению, строки матрицы А. Таким образом,

(Ах) = Аг.

дх '

359


Рассуждая аналогичным образом, можно прийти к выводу: если а ь а 2,

ап — столбцы матрицы А,

то

х А --\х '

ах х ' а 2... х ' а п\

и

 

(х'А) = А.

Из соотношения (3) следует, что полная запись ~ в случае, если раз­ мер матрицы А, например 3x3, будет выглядеть следующим образом:

 

dyi

дуг

дУз

 

дхх

дхх

дх1

ду

дУг

дуг

дУз

дх

дх2

дх2

(4)

дх2

 

dyi

ду г

ду3

 

дх3

дх3

дх3

Пример. Дан вектор

II =55

!

иначе говоря,

'2

3

6 —Г

1 1

4 7_

--- 1

х2 __Хз_

Ух

2x i A~6x2xs

У2

Зхг2х2-|- 4х3

съ

Злу -f 4x2 4- 7хл_

1

 

 

 

В соответствии с формулой (4)

| |

может быть найден следующим обра­

зом:

 

 

 

ду_

2

3

4 = Л \

6

- 2

дх

 

4

7

—1

Рассмотрим квадратичную форму порядка п:

у х'Ах

V 2 atj xt Xj,

причем суммирование по индексам i и / ведется от 1до п. Дифференци­ руя эту форму по элементам вектора х, получим следующее выраже-

ние дхду :

360