книги из ГПНТБ / Агекян Т.А. Теория вероятностей для астрономов и физиков учеб. пособие
.pdf140 |
СЛУЧАЙНЫЙ ВЕКТОР |
[ГЛ. 3 |
и |
математическое ожидание р2: |
|
|
Р2 = $ P2/i (р) Ф = За2. |
(3.45) |
|
О |
|
Последнее равенство позволяет выразить f x (р) через сред нее квадратическое модуля скорости:
|
|
Ш |
- ' & |
|
2р* |
(3.46) |
|
|
|
т к |
|||||
|
|
|
|
|
V я |
(р2Ф |
|
Дисперсия р находится по формуле (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 |
“ |
Т |
+ |
Т |
5» |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
ц 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. Неравенство Шварца
Докажем неравенство
1М (ZxZ2)P < М (Х\)М (X?), |
(3.79) |
называемое неравенством Шварца *).
•) Неравенство Шварца есть не что иное, как известное нера венство Коши — Буняковского, если рассматривать как элементы гильбертова пространства со скалярным произведением
( X t , X f) = M X , X t .
