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

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

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

140

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

[ГЛ. 3

и

математическое ожидание р2:

 

 

Р2 = $ P2/i (р) Ф = За2.

(3.45)

 

О

 

Последнее равенство позволяет выразить f x (р) через сред­ нее квадратическое модуля скорости:

 

 

Ш

- ' &

 

2р*

(3.46)

 

 

т к

 

 

 

 

 

V я

(р

 

Дисперсия р находится по формуле (2.68),

 

 

ol =

Р2 -

(Р)2 = f3 -

~ ) a2 S 0,454а2,

(3.47)

так что

стандарт — среднее

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

модуля

скорости от

своего

среднего значения — равен

 

 

бр=

 

| / ' з - - | а = 0,674а.

(3.48)

З а д а ч а

63. Найти плотность вероятности

видимой

величины к-й по яркости звезды.

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

находится

в промежутке [тх, тх +

dmx], в

промежутке

[mi, т] звезд нет,

а в промежутке

[т, т Ч~

dm] звезда

есть [см.

задачу

48]:

 

 

e-N(mi) ^'(m i) dmie-^ ^ ) ^ ^ ) ] N' (m) dm =

e-N(m) ftp (m) dm N' (mx) dmx.

Проинтегрировав это выражение по тх от 0 до т, полу­ чим искомую плотность вероятности

/2 (т) — N (т) e_N(m>N' (т).

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

/*(«») в й г гЯ(т)ЛГ'(т )-

§ 37]

МАТЕМАТИЧЕСКОЕ ОЖИДАНИЕ И ДИСПЕРСИЯ

 

141

З а д а ч а 64.

Плотность

вероятности

каждого

из

прямоугольных компонентов X, Y и Z скорости

звезд

есть

нормальная

функция со

средним,

равным

нулю,

и дисперсиями, соответственно равными о?,

 

и о\.

Найти

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

 

62

Р е ш е н и е .

Отличие этой задачи

от

задачи

заключается только в том, что дисперсии компонентов скоростей по трем направлениям различны. Находим

f ( x , у, Z) =

1

(3.49)

310253 (2я)’/г

Распределение (3.49) называется распределением Шварцгиилъда. Оно является эллипсоидальным распределением. Это означает, что в пространстве скоростей плотность

вероятности в точках

поверхности

эллипсоида

 

 

 

 

(3.50)

(при произвольном с)

постоянна.

распределения

таков:

Общий вид эллипсоидального

/ (х, у, z)

= ц

 

(3.51)

З а д а ч а 65. Все

направления

трехмерного

случай­

ного вектора равновероятны. Функция распределения длины проекции X вектора на произвольное направление,

/i (х), известна.

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

/2 (р) модуля р случайного вектора.

Р е ш е н и е .

Если а — угол между случайным век­

тором и заданным направлением, то

X = р cos а,

причем О ^ а ^ я . Рассматривая случайные векторы (р,*Х) и (р, а) с плотностями g и glf можно написать

ё (р, х) ф dx = ё г (р, «) ф da =

= /а’(Р)!Ф • y sin a da = /*(р,) ф . у ^ .

142 СЛУЧАЙНЫЙ ВЕКТОР [ГЛ. 3

Интегрируя по р и учитывая,

что р

х,

находим

о о

 

 

 

А (*) = $ - у А (р)dP-

(3.52)

X

 

 

 

Продифференцировав уравнение (3.52)

по х, окончатель­

но получаем

 

 

 

/а (Р) = -

2рА'(р).

 

(3.53)

З а д а ч а 66. Ротационной скоростью v звезды назы­ вается линейная скорость точек экватора звезды, вызыва­ емая вращением звезды вокруг своей оси. Измерение расширения линий в спектре звезды, обусловленное ее вращением, дает не истинную ротационную скорость, а ее проекцию на луч зрения — видимую ротационную скорость

у = v sin i ,

(3.54)

где г — угол между осью вращения звезды и лучом зре­ ния. Величина v и у = v sin г у различных звезд различ­ ны и при случайном, выборе звезды могут рассматривать­ ся как случайные величины. Из наблюдений можно определить плотность вероятности видимой ротационной скорости f (у). Требуется, считая ее известной, найти плотность вероятности истинной ротационной скорости fy (v) и найти зависимость между математическими ожида­

ниями и дисперсиями у и V. Предполагается,

что все на­

правления осей вращения звезд равновероятны.

Р е ш е н и е .

Как

было определено выше, функция

распределения угла г

есть с sin

i. Поскольку i изменяется

от 0 до я/2, коэффициент с =

1.

 

Рассмотрим взаимно однозначно определяющие друг

друга случайные

векторы (v , у) и (v, i) с плотностями g

и gy. Из физическихсоображений очевидно,

что случай­

ные переменные v и i взаимно

независимы.

Поэтому ра­

венство (3.32) можно

записать

в виде

 

g (v, у) dv dy — gy (v, i) dv di = fy (v) dv sin

i di, (3.55)

где fy (v) — плотность вероятности v. Используя (3.54), перейдем в правой части (3.55) от i к у (при этом v нуж­ но считать фиксированным и все множители брать по

8 37]

МАТЕМАТИЧЕСКОЕ ОЖИДАНИЕ И ДИСПЕРСИЯ

143

абсолютной величине):

 

 

 

8 (v, у) do dy =п U (v) dv

dy.

 

Плотность вероятности / (у) равна интегралу от g (v, у) по всем возможным значениям v. Так как всегда v > у, то

J T - 7 iv-

<3-56)

Уравнение (3.56) определяет зависимость между плотно

стями

вероятностей случайных величин у ж V.

Так как

из наблюдений определяется / (у),

а искомым

является

Л (у),

это уравнение — интегральное. Оно легко приво­

дится к уравнению Абеля и имеет решение

 

 

iL y JL С

/(у)

dy.

(3.57)

 

h И = — я v dv } у

 

Использование (3.57) для вычислений неудобно. На практике основной задачей обычно является нахождение средней истинной ротационной скорости и ее дисперсии для звезд различных классов. Найдем поэтому при помо­ щи (3.56) зависимость между математическими ожидания­ ми v и у, а также их квадратов:

ОО 00 оо

« = $ уШ Н у =

 

=

о

О

V

= 5 пЛ (и) dv \ у- т

у 5 = dy = - - \ v h (V) d v \ = ^ v ,

о

о

о

ОООС

Уг =

- 7 ^ = 7 * “

 

Отсюда следует:

 

 

 

4 -

(3.58)

о» =

— (vf = у у 2 —

(3.59)

144 СЛУЧАЙНЫЙ ВЕКТОР |Г Л . 3

Рассматривая совокупность звезд, у которых измере­ ны видимые ротационные скорости, как статистический

коллектив, и вычисляя в них у

и у 2 по формулам

 

 

П

 

у =

4 -

2

 

 

 

 

г=1

 

 

 

n

 

 

у 5 =

4 "

2

 

»?.

 

 

i=l

 

найдем затем при помощи равенств (3.58) и (3.59) среднее значение и дисперсию истинных ротационных скоростей в этой совокупности звезд.

Аналогично, используя (3.56), можно найти зависи­ мость между моментами любого порядка функций распре­ деления видимой и истин­

 

 

ной ротационной скорости

 

 

звезд.

67.

Найти

 

 

З а д а ч а

 

 

зависимость

между

функ­

 

 

циями распределений квад­

 

 

ратов истинных и видимых

 

 

сферичностей

галактик,

 

 

считая, что галактики яв­

 

 

ляются сжатыми эллипсо­

Рис.

13.

идами вращения и все ори­

 

 

ентации их плоскостей сим­

Р е ш е н и е.

 

метрии равновероятны.

Истинной сферичностью сжатого эллип­

соида вращения называется отношение его малой полуоси к большой. Обозначим квадрат истинной сферичности §. Один и тот же сжатый эллипсоид вращения при наблюде­ нии его по различным направлениям представляется в виде эллипса с различной по величине малой полуосью и постоянной большой полуосью. Назовем видимой сферич­ ностью галактики отношение малой полуоси к большой полуоси ее видимого эллипса. Квадрат видимой сферич­ ности обозначим к). Рис. 13 показывает, как зависит ве­

личина малой полуоси видимого эллипса, равная У ц (большая полуось принята за 1), от угла i между направле­ нием луча зрения и плоскостью симметрии галактики.

i 371

МАТЕМАТИЧЕСКОЕ ОЖИДАВИЕ И ДИСПЕРСИЯ

145

Найдем зависимость между т), | и г. Уравнение касатель­ ной к эллипсу имеет вид

у = х tg г + У 1 + t g 4 .

Находя расстояние этой касательной от начала коорди­ нат, получаем

r 1 У

i + tga <

откуда находим

 

 

sin2i =

j ^ - | .

(3.60)

Рассмотрим два взаимно однозначно определяющих ДРУГ ДРУга случайных вектора (|, ц) и (|, i). Так как | и г взаимно независимы, на основании (3.32) находим

g (£, ri) dl dri = gx (g, i) dl di

= f x (g) d\ cos i di.

(3.61)

Используя (3.60), получаем

 

 

 

g (£. Л) = 2

 

Ml)

(ц- I )

(3.62)

 

 

Интегрируя (3.62) по g (l

rj),

приходим к искомому со­

отношению

 

 

 

 

 

A (О

dl.

(3.63)

 

 

 

Y( 1 - 0 ( л - 0

 

Это интегральное уравнение относительно функции рас­ пределения истинных квадратов сферичностей галактик принадлежит к типу Абеля и разрешается. Однако удобнее использовать решение в моментах. Помножим обе части (3.63) на г]" и проинтегрируем по всем значениям г], т. е. от 0 до 1:

А (I)

dl =

 

V(i —о (л—D

 

 

А(0

: dr\.

 

2 V i - l

 

V Г) —<

146 СЛУЧАЙНЫЙ ВЕКТОР [Г Л . 3

Вычисляя интеграл и используя выражения для моментов

случайных

величин,

получаем

 

 

 

 

 

 

 

п

 

 

 

 

к

 

 

___

 

 

^

=

2

<?«

~

 

2

( - 1 У“

4 г * .

(3.64)

 

 

 

к=о

 

 

г=о

 

 

 

Для п = 1

и

п =

2

равенство

(3.64)

принимает

вид

 

 

-

 

1

, 2

=

 

 

 

 

(3.65)

 

 

4

Т

+

Т

 

 

 

 

 

 

 

 

 

 

 

 

 

ц 2 =

— +

— £ 4- — Р

 

(3.66)

 

 

1

 

15 ^

15 *

'

15 *

 

 

откуда получаем

 

 

3 -

 

 

 

 

 

 

 

 

 

=

 

 

1

 

 

 

(3.67)

 

 

 

^ ~

2

11

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

15 -=■

6 —

1

 

(3.68)

 

 

 

* =

¥ * > ’! - Т т1 - 8 •

 

 

 

 

Находим также

дисперсию

 

 

 

 

 

4

=

i 5- d

) 2

=

¥ v - | - ( n ) 2- 4 -

(3.69)

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

галактик и вычисляя ц и ц2, найдем затем при помощи равенств (3.67) и (3.69) среднюю величину и дисперсию квадратов истинных сферичностей.

3 а д а ч а 68. Определить в скоплении галактик функ­ цию распределения угла между видимым направлением от галактики на центр скопления и видимым направле­ нием большой оси галактики. Рассмотреть два предполо­ жения: 1) все ориентации плоскости симметрии галактик равновероятны; 2) все плоскости симметрии галактик

проходят через

центр скопления.

Р е ш е н и е .

Направление видимой большой оси

галактики есть направление прямой, по которой плоскость симметрии галактики пересекает картинную плоскость (плоскость, перпендикулярную к лучу зрения). Если выполняется предположение 1), то все ориентации види­ мой большой оси на картинной плоскости равновероят­ ны и плотность вероятности острого угла между направле­

1.37] МАТЕМАТИЧЕСКОЕ ОЖИДАНИЕ И ДИСПЕРСИЯ 147

нием на центр скопления и большой осью галактики дает­

ся равенством

 

 

т

= -§.

(3.70)

Пусть теперь выполняется предположение 2). Рассмот­ рим рис. 14, на котором центр сферы Q совпадает с цент­ ром скопления галактик. Сфера проведена через галакти­ ку, находящуюся в точке В. Наблюдатель смотрит в направлении АО, угол между лучом зрения и на­ правлением из центра ско­ пления на галактику ра­ вен г. Картинная плос­ кость проходит через центр скопления перпендикуляр­ но к лучу зрения. Боль­ шие круги, образуемые пе­ ресечением сферы с плос­ костью А ОБ и картинной плоскостью, пересекаются

в точке D. В этой точке картинной плоскости наблюда­ тель видит галактику В. Проходящую через центр скоп­ ления плоскость симметрии галактики пересекает картин­ ную плоскость по прямой ОС. Направление ОС есть направление видимой большой оси галактики и, следова­

тельно, интересующий нас

угол р есть угол COD. Из

сферического треугольника

BCD находим

ч *

<3-7»

Угол ас есть случайная величина с плотностью вероят­ ности

Л И = -I •

(3.72)

Рассмотрим зависимость между функциями распределе­ ния случайных векторов (р, i) и (х, г):

Л

g (Р, i) dp di = g1 (x, i) dxdi = — dx sin i di.

(3.73)

148 СЛУЧАЙНЫЙ ВЕКТО Р- [ГЛ. 3

Заменяя при помощи (3.71) dx в (3.73) и выполняя интег­

рирование по

i, находим

 

 

2

1

71/2

2 In sin p

 

Р cos i sin i di

(3.74)

я

cos2 (E

J cos2 i -f- tg2 P

Я COS2P

 

 

0

 

 

Для того чтобы выяснить, какое

из предположе­

ний 1) и 2) имеет место в скоплениях

галактик,

нужно

сравнивать распределения (3.70) и (3.74) с наблюдае­ мым распределением. Это часто бывает весьма неудобно, в особенности, если число галактик в скоплении не очень велико.

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

Заметим, что в нашей задаче в случае предположения 2) углы р должны ожидаться в среднем меньшими, чем при предположении 1). Поэтому можно рассматривать,

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

cos2 р, тем более,

что оно просто находится.

 

 

В предположении

1)

 

 

cos2 р =

^cos2

=

0,5.

 

о

 

 

В предположении 2)

 

 

[271

 

 

 

= - 5cos2

 

=

In2^0,693.

О

 

 

 

В зависимости от того, какое из чисел, 0,5 и 0,693, ока­ жется ближе к найденному из наблюдений,

 

й 'П

|

c°s2 р =

2

cos2Pi-

можно судить о предпочтительности гипотез 1) или 2).

§ 39]

НЕРАВЕНСТВО ШВАРЦА

149

§ 38. Математическое ожидание функции от случайного вектора

Если f (х1, х2, • • •, хп) есть плотность вероятности случайного вектора (Хг, Х 2, ■■., Х п), а Г] (Хг, Х 2, . . .

. . Хп) — некоторая функция от этого случайного век­ тора, то величина

 

о о

с о

с о

Мг\(Хи Х2, . . Х п) =

[

^

$ Ц{х1, х 2, . . . , х п) х

 

— о о — о о

— оо

X

/ (хх, х2, . .., хп) dxx dx2. . . dxn (3.75)

называется математическим ожиданием т] (Хх, Х 2, ■. ., Хп).

Аналогично тому, как это было сделано для случайной величины в § 20, легко доказать, что математическое ожидание суммы функций случайного вектора равно сумме математических ожиданий этих функций, т. е.

М [тц (Хг, Х 2, . . ., Х п) + 1Ъ (Хг, Х2, . . ., Z J ] =

= Мц1 (Хг, Х 2, . . ., Х п) -f- М у]2 (Xi , Х 2, . . ., Х п). (3.76)

Точно так же очевидно, что

М [сц (Хи Х 2, . . ., Х„)] = сМц (Xlt Х 2, . . ., Х п).

(3.77)

Если X t и Х 2 взаимно независимы, то

оооо

М [тц (Хх) г\2(Х2)] = § \ % 0гх) т]2 (х2) /х (жх) U (х2)dx-г dx2 =

— оо — оо

= Мг\г{Хг)Мх\2(Х2). (3.78)

§ 39. Неравенство Шварца

Докажем неравенство

(ZxZ2)P < М (Х\)М (X?),

(3.79)

называемое неравенством Шварца *).

•) Неравенство Шварца есть не что иное, как известное нера­ венство Коши — Буняковского, если рассматривать как элементы гильбертова пространства со скалярным произведением

( X t , X f) = M X , X t .

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