Теорема 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 |
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» ■— сочетание искаженных латинских слов «единица» и
«степень». ■— П р и м , п е р е в .
Второй из приведенных векторов показывает, что символ можно
представить вектором дифференциальных операторов. Заметим, что
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 ) |
|
= А', |
|
|
|
|
|
|
|
так как а[ — это, по определению, строки матрицы А. Таким образом,
— (Ах) = Аг.
дх '
Рассуждая аналогичным образом, можно прийти к выводу: если а ь а 2,
ап — столбцы матрицы А, |
то |
х А --\х ' |
ах х ' а 2... х ' а п\ |
и |
|
— |
(х'А) = А. |
Из соотношения (3) следует, что полная запись ~ в случае, если раз мер матрицы А, например 3x3, будет выглядеть следующим образом:
|
dyi |
дуг |
дУз |
|
дхх |
дхх |
дх1 |
ду |
дУг |
дуг |
дУз |
дх |
дх2 |
дх2 |
(4) |
дх2 |
|
dyi |
ду г |
ду3 |
|
дх3 |
дх3 |
дх3 |
Пример. Дан вектор
иначе говоря,
|
Ух |
2x i A~6x2—xs |
|
У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до п. Дифференци руя эту форму по элементам вектора х, получим следующее выраже-
ние дхду :