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

книги из ГПНТБ / Агекян Т.А. Теория вероятностей для астрономов и физиков учеб. пособие

.pdf
Скачиваний:
63
Добавлен:
23.10.2023
Размер:
8.81 Mб
Скачать

170

СЛУЧАЙНЫЙ

ВЕКТОР

ЕГЛ. 3

Из

(3.142) следует, что

V2, . .

Vh-i являются

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

На основании результатов § 46 заключаем, что случай­

ная величина

 

 

 

 

 

Н-1

 

= v = 2 т? = 2

пр{

= 2 Vi

(3.143)

г=1

 

1=1

 

при

п

со,

гг'!* (mi npi)3(npi)-’> 0,

i = l , 2 , . . .

..., к, имеет асимптотическое раснределение

 

 

/(»)•=

2

 

 

 

 

(3.144)

 

 

§ 48. Моменты случайного вектора.

 

 

 

Коэффициент корреляции

 

Математическое ожидание

функции

 

 

 

 

 

 

П

 

 

(3.145)

 

 

 

 

 

1=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П

 

 

называется моментом

порядка 2

относительно начал

{аи а2, . . ., ап)

случайного

г—1

Х 2, . ■., Х п).

вектора (-Х4,

 

Если все аг =

0, то моменты называются начальными,

если все а; =

MXi, то моменты называются центральными.

Математическое ожидание

X t

равно

 

 

■^i ~

\

\

\

ХЛ (ж1>х2, ■■- , хп) dxx dx2. . . dxn

а дисперсия

X t—

 

 

 

 

 

оо

ос

со

 

 

 

 

 

Oi

^

^

^

(*^г

-^г) / (^l) *^2) • **?%п) dX\ dx<2, • . . dxn.

— оо — оо

— оо

§ 48] МОМЕНТЫ СЛУЧАЙНОГО ВЕКТОРА 171

Эти определения совпадают, очевидно, с определениями, данными ранее, так как

СО ОО ОО

^

^ ...

|, Х%, • • • , Хп) (lxj(Ix^...dj'i -jdXi_j. . . dxn= f (*T;)•

— 00 — 00

— 00

ми

Моменты являются информативными характеристика­

случайного вектора.

Для двумерного случайного вектора (X , Y) важной характеристикой является смешанный центральный мо­ мент второго порядка:

оо

оо

 

[H,i = I

$ (х — Z ) ( y — Y)f{x,y)dxdy.

(3.146)

00—00

Рассмотрим смысл этой характеристики. Подынтеграль­ ное выражение в (3.146) положительно, когда отклонения х и у от математических ожиданий имеют одинаковый знак, и отрицательно, когда они имеют разный знак. Если отклонениям одного знака соответствуют, в общем, большие значения функции / (х , у), чем отклонениям разного знака, то интеграл (3.146) оказывается положи­ тельной величиной. В противоположном случае он отри­ цателен.

Если X viY взаимно независимы, то

ОО оо

P i,i= $

%) h (х) dx § — У )/2 {у) dy — 0, (3.147)

— о о

— оо

так как центральные моменты первого порядка всегда равны нулю. Обратное утверждение неверно. Из равен­ ства нулю р.1,1 не следует, что X и Y взаимно независимы. Оно указывает лишь, что в среднем положительные от­ клонения одной из случайных величин X, Y компенси­ руются отклонениями определенного знака у другой из них. Можно сказать, что равенство нулю p.i,i означает

отсутствие линейной статистической

зависимости между

X и Y.

р.1,1

0, то между X и Y

существует положи­

Если

тельная,

а при

р.!,! < 0 — отрицательная линейная ста­

тистическая зависимость.

pi,i имеет размерность произведения размерностей случайных величин X и Y. Удобно ввести в рассмотрение

172

СЛУЧАЙНЫЙ ВЕКТОР

[ГЛ. 3

безразмерную

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

 

 

pi, 1

 

 

бхС5у

 

где Ox, cry соответственно стандарты X и Y. Характеристику г называют коэффициентом линей­

ной корреляции или просто коэффициентом корреляции

случайных величин X и Y.

Шварца, всегда

Как это следует из неравенства

- 1 < г < 1 .

(3.148)

Если X и Y взаимно независимы,

то согласно (3.147)

их коэффициент корреляции равен нулю.

Определим коэффициент корреляции в том случае,

когда Y является

линейной функцией X:

 

 

 

 

Y =

аХ +

Ь.

(3.149)

В этом

случае Y

= аХ -j- b,

у Y = а (X — X) и

/ (х , y)dx dy= б ах

b)f1 (x)dx

 

Поэтому

находим

 

 

 

= 6 (у — ах — b)f, (y)dy.

 

 

 

 

 

Рт, 1 =

о о

СО

 

 

 

 

 

\

$ (x — X)(y — Y)f(x,y)dxdy =

 

 

 

 

 

 

оо

 

 

 

 

 

 

=

§

а (х — X)2(х) dx = ав\,

Оу =

$

$ Y)2f(x, y)dxdy =

 

 

 

 

 

 

ОС

 

 

 

 

 

=

jj а2(х X)2fi(x)dx

а2Ох

 

 

 

 

 

— СО

 

 

 

 

ау =

| л |

 

(3.150)

так как стандарт всегда положителен.

 

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

когда Y есть линейная функция X, то

§ 48]

МОМЕНТЫ СЛУЧАЙНОГО ВЕКТОРА

173

Коэффициент корреляции равен 1, если а

0, и ра­

вен

—1, если

а < 0 .

 

Отсутствие

линейной статистической зависимости и

функциональная линейная зависимость между X и Y — это крайние частные случаи. В общем случае между X и Y существует линейная статистическая зависимость и их коэффициент корреляции отличен от 0 и от + 1 : он за­ ключен в промежутке (—1, 1). Чем больше | г |, тем силь­ нее линейная статистическая зависимость между X и Y. Она положительная, если г 0, и отрицательная, если

г< 0 .

За д а ч а 70. Найти параметры двумерного нормаль­ ного распределения:

 

- Q ( i , V).

(3.152)

/ (я, У) =

 

 

где

2 я а 1 3 2 1 ^ 1 — р 2

 

 

 

 

 

( » - 1 )(у -У ) I

Q (*. у) = -рг-;

г , (у-У )2

— 2р

3 l3 2

 

Р е ш е н и е .

Находим

 

 

(3.153)

 

(х~ХР

 

 

 

 

 

 

 

 

„ 2

 

 

 

 

2я1

(3.154)

 

 

 

 

М У )=

= — L = r '.

.

(3.155)

 

-- 50

 

.

*

Таким образом, каждая из случайных величин X и Y

распределена нормально,

 

 

 

 

(Ti — Ох, й 2 =

0 у .

 

 

Найдем смешанный центральный момент второго порядка:

ооС»

Hi, i = 5 5 — (г/ — у ) / (а>г / ) =p<3i<32- (3.156)

—оо —оо

Равенство (3.156) показывает, что р есть коэффициент

корреляции X и Y.

Глава 4

ОЦЕНИВАНИЕ ПАРАМЕТРОВ РАСПРЕДЕЛЕНИЙ

ИСТАТИСТИЧЕСКИЕ ГИПОТЕЗЫ

§49. Статистические коллективы

Пусть рассматриваются объекты, которые могут отли­ чаться друг от друга значением некоторой определенной характеристики, Всякое каким-то образом выделенное множество таких объектов называется статистическим коллективом (или генеральной совокупностью), а объекты,

в него входящие,— членами этого статистического кол­ лектива. Число членов статистического коллектива т называется объемом статистического коллектива. Харак­ теристика X, которая может принимать различные зна­ чения у различных членов коллектива, называется аргу­ ментом статистического коллектива.

Можно, например, рассматривать как статистический коллектив звездное скопление, считая аргументом массу (или температуру, или показатель цвета и т. д.) его чле­ нов-звезд. Можно также мысленно выделить как стати­ стический коллектив множество всех звезд спектрального класса В, входящих в состав Галактики, считая их све­ тимость (количество энергии, излучаемой в единицу времени) аргументом.

Примерами статистических коллективов будут мно­ жество частиц плазмы в некотором объеме с аргумен­

том,— величиной заряда

у частицы,— или множество

всех молекул кислорода,

находящихся в воздухе внутри

данного помещения, с аргументом,— модулем скорости молекулы.

Пусть в статистическом коллективе значение аргу­ мента хг имеют mi членов, хг — т 2 членов и т. д., x h mk

членов. Числа

 

ти тг, . . ., mh

(4.1)

называются частотами соответственно

значений

хи хг, . . ., x h

. (4.2)

аргумента X.

 

§ 49]

СТАТИСТИЧЕСКИЕ КОЛЛЕКТИВЫ

175

Если т — объем статистического коллектива, то, оче­

видно,

что

н

 

 

 

 

 

 

 

 

2 т{ =*т.

 

(4.3)

 

 

i—1

 

 

Величины

 

 

 

 

mi

m2

т /с

(4.4)

 

т

’ • ’

т

 

 

называются относительными частотами значении аргу­ мента (4.2).

Очевидно, что если случайным образом извлечь из статистического коллектива один член, затем, после возвращения его обратно, снова случайно извлечь какойто член и т. д., то значения аргумента извлекаемых членов можно рассматривать как значения случайной величины. Если объем статистического коллектива ограничен, то эта случайная величина может быть только дискретной. Согласно классическому определению вероятность p t того, что случайная величина примет значение х ;, равна

относительной частоте аргумента т г .

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

Статистический коллектив определяется интеграль­ ным законом распределения аргумента, рассматриваемого

как случайная величина:

 

 

F(*) = P ( X O ) = 2

2 гп,.

(4.5)

■Ti<X

X.J<X

 

Произведение

mF (х)

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

S TYL

\

^ - i

Xi)

(4-6)

б X i) =

2

г

 

i

 

 

176 ОЦЕНИВАНИЕ ПАРАМЕТРОВ РАСПРЕДЕЛЕНИЙ [ГЛ. 4

есть плотность вероятности. Функцию

 

X (х) = mf (X) = 2

— ^i)

(4.7)

i

 

 

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

 

ь

 

 

^х (х) dx

(4.8)

 

а

 

дает число членов

коллектива, аргумент

которых за-

ключей между а и

Ь.

 

Статистический коллектив характеризуется моментами распределения. В частности,

yr\ m.

1 VI

(4-9)

=

— 2l mixi

i

i

 

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

°о = 2

(xi — х0)2 =

2 mi (*i — хо)2

(4-10)

i

 

i

 

есть дисперсия аргумента. Часто их называют также сред­ ним и дисперсией статистического коллектива. Равенства (4.9) и (4.10) показывают, что если каждому члену стати­ стического коллектива (а не каждому значению аргумен­ та) присваивать свой номер, так что все mi = 1, то эти равенства должны быть записаны в виде

т

 

=

(4-И)

i= l

 

т

 

ol = - ^ ^ ( X i - x 0)\

(4.12)

i= l

 

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

Статистический коллектив может иметь и бесконечно большой объем. Такой статистический коллектив счита­

№9]

СТАТИСТИЧЕСКИЕ КОЛЛЕКТИВЫ

177

ется заданным, если заданы функции распределения аргумента — плотность вероятности или интегральный закон распределения. Рассмотрим, например, статисти­ ческий коллектив — множество всех прямых, различным образом ориентированных в пространстве, в котором аргу­ ментом является угол между прямой и некоторой фикси­ рованной плоскостью. Объем этого статистического кол­ лектива бесконечен и даже несчетен. Аргумент в нем

может принимать любые значения в промежутке

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

/ (a)da = cos a da .

(4.13)

(4.13) имеет смысл доли членов коллектива, аргумент которых заключен между а и а + da.

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

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

Если Ах мало, то / (х)Ах приближенно (тем точнее, чем меньше Ах) равно вероятности попадания случайной

178

ОЦЕНИВАНИЕ ПАРАМЕТРОВ

РАСПРЕДЕЛЕНИЙ

1ГЛ. 4

величины в промежуток

 

 

 

X---- |-Д х, x +

- |- A x j.

(4Д4)

Когда Аж конечно, приближение в общем случае будет лучшим, когда ж совпадает с серединой промежутка. В таком случае плотность вероятности определяется из приближенного равенства

/ (ж) Аж =

,

(4.15)

где mi — частота значении агрумента в данном интервале. Положив в равенстве (4.15) Ах = 1, мы видим, что плотность вероятности / (ж) приближенно равна вероят­ ности того, что случайная величина примет значение

в интервале

(4.16)

Для примера рассмотрим статистический коллектив звезд в окрестности Солнца, если аргументом служит абсолютная величина звезды М. Так как абсолютная

величина звезды однозначно

определяет

ее светимость,

то функцию распределения

абсолютных

величин звезд

кратко называют функцией светимости.

Из многих определений функции светимости приведем те, которые были получены П. П. Паренаго (табл. 4) и

В. Ж. Лейтеном (табл.

5).

 

 

 

Т а б л и ц а 4

 

Функция светимости Паренаго

 

 

м

НМ)

м

 

НМ)

м

НМ)

- 6

0,00000013

+

4

0,017

+ 13

0,107

—5

0,00000126

+

5

0,024

+14

0,118

—4

0,0000048

+

6

0,033

+15

0,018

- 3

0,000019

+

7

0,036

+16

0,102

—2

0,000060

4- 8

0,031

+17

0,079

—1

0,000126

+

9

0,033

+18

0,049

+ 0

0,00051

+10

0,047

+ 19

0,020

+ 1

0,0023

+11

0,069

+ 20

0,0041

+ 2

0,0063

+12

0,100

+21

0,0008

+ 3

0,0078

 

 

 

 

 

S 49] СТАТИСТИЧЕСКИЕ КОЛЛЕКТИВЫ 179

Т а б л и ц а 5

Функция светимости Лейтена

м

Д М ) - 100

м

Д М ) -100

м

Д М ) .100

- 0 ,5

0,009

8.5

3.6

17.5

6,7

+ 0 ,5

0,094

9.5

4.4

18.5

4,1

1.5

0,23

10.5

5.4

19.5

2.3

2.5

0,58

11.5

6,3

20.5

1.4

3.5

0,86

12.5

7.6

21.5

0,7

4.5

1,26

13.5

9.7

22.5

0,4

5.5

1,7

14.5

12,1

23.5

0,05

6.5

2,3

15.5

14,4

24.5

0,005

7.5

2,9

16.5

10,8

 

 

Как отмечалось выше, в таблицах4и5 плотность веро­ ятности / (М) приближенно равна доле звезд с абсолютной

величиной, заключенной в промежутке м - i м

Я -

Допустим, что табл. 4 и 5 дают истинное распределение по абсолютным величинам в двух каких-то статистических коллективах звезд. (На самом деле они являются двумя приближенными распределениями абсолютных величин звезд в окрестности Солнца.) Тогда имеется возможность при помощи моментов выполнить сравнения этих двух статистических коллективов. Нужно вычислить средние (Af0), стандарты (о0), асимметрии (As) и эксцессы (Ех).

Моменты относительно начала а (а удобно принять рав­ ными 14 или 15) вычисляются по формуле

>ч,а = 2 ( М - я )* /( М ),

(4.17)

м

 

которая отвечает и равенству (2.48) и (2.49). Затем по фор­ мулам (2.63) и (2.65) — (2.67) определяются М 0, ц2, р,3

и р,4, и при помощи (2.89) и (2.92) асимметрия и эксцесс. Результаты вычислений для функций светимости Па-

ренаго и Лейтена имеют вид

 

Функция

Функция

 

Паренаго

Лейтена

Мо

+ 12,73

+13,55

бо

3,850

3,873

As

— 0,708

— 0,610

Ех

0,0059

+ 0,230

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