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

Лекции_тервер

.pdf
Скачиваний:
28
Добавлен:
11.03.2016
Размер:
2.17 Mб
Скачать

Определение 24.1. Пусть случайные величины ξ1, . . . , ξn независимы и одинаково распределены по нормальному закону распределения

1

 

2

2

fξi(xi) =

 

σ

e−xi

/2σ ,

2π

то есть ai = 0 и σi = σ для всех i = 1, . . . , n.

Тогда распределение случай-

ного вектора ξ = (ξ1, . . . , ξn) называется нормальным сферическим распределением или проще сферическим распределением.

Лемма 24.2. Плотность вероятности сферического распределения есть

fξ(x1, . . . , xn) =

1

 

e

1

(x12+ . . . +xn2 ).

 

2σ2

(2π)n/2

σn

 

 

 

 

Доказательство. Это утверждение очевидно следует из 24.1 и 13.6.2).

Замечание 24.3. Так как поверхности равной вероятности, то есть fξ(x1, . . . , xn) = const, являются концентрическими (n 1)-мерными сферами x21 + . . . + x2n = R2 с центром в начале координат, то сферическое распределение инвариантно относительно любого ортогонального преобразования

A : Rn Rn. То есть, если η = A · ξ, то fη(x1, . . . , xn) = fξ(x1, . . . , xn), где через A обозначена ортогональная матрица.

Положим в сферическом распределении σ = 1, получим плотность вероятности

fξ(x1, . . . , xn) =

1

e21 (x12+ . . . +xn2 )

n/2

 

(2π)

и назовём её плотностью вероятности стандартного сферического распределения.

χ- и χ2-распределения

Определение 24.4. Пусть случайные величины ξ1, . . . , ξn независимы и одинаково распределены по стандартному нормальному закону распределения. Случайные величины

χ = ξ12 + . . . + ξn2

и

χ2 = ξ12 + . . . + ξn2

называются соответственно χ- и χ2-распределениями с n степенями свободы. Наша ближайшая задача найти их плотности вероятности fχ(x) и fχ2 (x).

91

Рис. 30: К вычислению плотности χ-распределения.

Теорема 24.5. Плотность вероятности χ-распределения задаётся по

формуле

xn−1

 

 

 

1 2

 

fχ(x) =

 

e

2 x ,

x ≥ 0,

2n2 1Γ (n2 )

где Γ (x)— гамма-функция.

Доказательство. Вероятность события {x < χ < x + dx} можно вычислить, интегрируя плотность вероятности fξ(x1, . . . , xn) стандартного сферического распределения по n-мерному сферическому слою радиуса x и толщины dx (см. рис. 30).

P(x < χ < x + dx) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

. . .

 

 

fξ(x1, . . . , xn) dx1 . . . dxn.

 

 

 

 

 

 

 

 

 

 

x <

 

 

 

< x+dx

 

 

 

 

 

 

 

 

 

 

 

 

x12+ . . . +xn2

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

< x

+

dx

)

 

 

 

. . .

 

 

fχ(x) =

P(x < χ

 

 

=

 

fξ(x1, . . . , xn) dx1 . . . dxn =

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x12+ . . . +xn2

 

 

 

1

 

 

 

 

. . .

1

2

 

2

 

 

 

=

 

 

 

 

 

 

e2

(x1+ . . . +xn) dx1 . . . dxn =

 

 

 

(2π)n/2

 

 

 

 

 

 

 

 

 

x12+ . . . +xn2

= x

 

 

 

 

 

 

1

 

 

1

 

2

 

. . .

 

 

1

2

=

 

 

e

2 x

 

 

dx1 . . . dxn = Cxn−1e2 x ,

 

(2π)n/2

 

 

 

 

 

 

 

 

 

 

 

x12+ . . . +xn2

= x

 

 

 

 

92

где последнее равенство справедливо потому, что площадь (n − 1)-мерной

сферы радиуса x пропорциональна xn−1. Чтобы найти константу C, восполь-

зуемся свойством плотности fχ(x) dx = 1, откуда получаем

0

Cxn−1e12 x2 dx = 2n2 1Γ (n2 ) C = 1,

0

откуда следует требуемый результат.

Теорема 24.6. Плотность fχ2 (x) случайной величины χ2 есть

f 2 (x) =

 

xn2

1

ex2 ,

x

 

0.

2

n

(n2 )

χ

2 Γ

 

 

 

 

 

 

 

 

 

 

Доказательство. Ясно, что случайные величины χ и χ2 связаны функцией g(x) = x2. Поэтому по теор. 12.2, находим, что

f 2 (x) =

 

g1

(x)

 

f (g1

(x)) =

1 f (x) = nxn2 1

 

ex2 , x 0.

 

(

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

)

 

χ

 

 

 

 

· χ

 

2x χ

 

 

2 2 Γ

n2

 

Замечание 24.7. Заметим, что χ2-распределение с двумя степенями свободы (n = 2) является показательным распределением 10.7 с параметром λ = 1/2 и с плотностью вероятности 12 ex2 , x ≥ 0. А χ-распределение с тремя степенями свободы (n = 3) называется распределением Максвелла и описывает в кинетической теории газов распределение модуля скорости частиц в состоянии термодинамического равновесия.

О распределении Стьюдента

Замечание 24.8. Этот пункт имеет информативный характер. Распределение Стьюдента26 часто возникает в прикладных задачах матстатистики.

26Стьюдент — псевдоним английского химика и статистика Вильяма Госсета (William Sealy Gosset, 1876

— 1937). По окончании университета был (с 1899) младшим управляющим пивоваренной компании Arthur Guinness & Son в Дублине. В то время в совете директоров компании уже понимали важность методов теории вероятностей и статистики для массового производства и направили Госсета на стажировку в Лондонский университет для изучения статистики. Вскоре Госсет получил теоретические результаты в статистике, которые оказались важными для пивоварения и которые он хотел опубликовать. В компании не возражали против его публикаций, однако совет директоров хотел, чтобы конкуренты как можно дольше не знали о применении статистики на пивоварне Arthur Guinness & Son. Это объясняет появления псевдонима Стьюдент.

93

Пусть случайные величины ξ0, ξ1, . . . , ξn независимы и одинаково распределены по стандартному нормальному закону распределения

fξi (xi) = 1 e−x2i /2,

2π

для всех i = 0, 1, . . . , n. Тогда случайная величина

τn =

 

ξ0

 

 

n1 (ξ12 + . . . + ξn2)

называется отношением Стьюдента, её распределение называется распределением Стьюдента с n степенями свободы. Формула плотности вероятности

τn имеет вид

(

Γ

 

)2

 

(

 

)

n+12

Γ

n+1

 

 

 

 

x2

 

fτn (x) =

 

2

(

 

)

 

 

 

, −∞ < x < ∞.

 

 

 

n

 

1 + n

 

Замечание 24.9. Легко проверить, что последовательность отношений Стьюдента n}n=1 слабо сходится к стандартному нормальному распределению, то есть

 

1

e

x2

nlim fτn

(x) =

 

2

.

2π

→∞

 

 

 

 

 

§25. Цепи Маркова

Непосредственным обобщением схемы независимых испытаний (см. опред. 9.3) является схема цепей Маркова.

Определение 25.1. Пусть эксперимент состоит в том, что

1)проводят последовательность испытаний с номерами s = 1, 2, . . . ;

2)в каждом испытании появляется одно из k несовместных событий A1,

A2, . . . , Ak, при этом то, что в s-ом испытании произошло событие Ai, будем обозначать Asi ;

3)вероятность появления события Asj+1 (то есть j-го события в (s + 1)-ом испытании) зависит только от того, какое событие произошло в предыдущем s-ом испытании и не зависит от того, какие события произошли в испытаниях

сномерами s − 1, s − 2, . . . .

Тогда такая последовательность событий называется (простой) цепью Маркова.

94

Пример. Согласно предложенной Н. Бором27 модели атома водорода электрон может находится только на одной из допустимых орбит. Обозначим через Ai событие, состоящее в том, что электрон находится на i-ой орбите. Эти события Бор назвал стационарными состояниями атома. Вероятность перехода электрона с i-ой на j-ю орбиту зависит только от i и j, потому что разность j − i зависит от количества испускаемой или поглощаемой атомом энергии и не зависит от того, на каких орбитах находился электрон в прошлом. Этот пример доставляет цепь Маркова теоретически с бесконечным числом состояний атома.

Теория цепей Маркова довольно обширна, поэтому мы ограничимся изу-

чением так называемых однородных цепей Маркова.

( )

Рассмотрим условную вероятность P Asj+1 Asi , т.е. вероятность появления события Asj+1 при условии, что событие Asi произошло.

 

Определение 25.2. Цепь Маркова называется однородной, если услов-

ная вероятность P Ajs+1

 

Ais

не зависит от номера s, при этом вероятность

P As+1 As

называется вероятностью перехода от события As к событию

 

 

j

 

i

(

 

)

i

As+1

и

обозначается

 

 

 

j(

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

( )

pij = P Asj+1 Asi .

Замечание 25.3. Ясно, что

1)0 ≤ pij 1 и

2)pi1 + pi2 + . . . + pik = 1 для всех i = 1, 2, . . . , k.

Определение 25.4. Полная информация о вероятностях перехода в цепи Маркова содержится в таблице

π1

=

p11

p12 . . .

p1k

 

,

p21

p22 . . .

p2k

 

 

 

 

 

 

 

 

 

 

.p.k.2 .. .. ..

 

 

 

 

 

.p.k.1

p. kk. .

 

которая называется матрицей перехода.

Замечание 25.5. Главной задачей теории цепей Маркова является нахождение вероятности перехода от события Asi к событию Asj+n, произошедшему через n испытаний. Обозначим эту вероятность

( )

pij(n) = P Asj+n Asi ,

27Нильс Хенрик Давид Бор (Niels Henrik David Bohr, 1865 — 1962), датский физик, один из создателей квантовой физики. Создал первую квантовую модель атома, участвовал в разработке основ квантовой механики, теории атомного ядра, ядерных реакций и взаимодействия элементарных частиц со средой.

95

а искомую матрицу перехода через n испытаний через

 

p11(n)

 

p21(n)

 

 

πn = pk.1.(.n)

p12(n) . . . p1k(n)

p22(n) . . .

p2k(n)

. . . . . . . . .

pk2(n) . . .

pkk(n)

.

Теорема 25.6. Если 0 < m < n, то πn = πm · πn−m.

Доказательство. Рассмотрим какое-нибудь промежуточное испытание с номером s + m, то есть s < s + m < s + n. Пусть в этом испытании появится

s+m

 

 

 

 

 

какое-то событие Ar

s, где 1 ≤ r ≤ s+.mВ наших обозначениях вероятность

 

 

 

k

 

 

 

перехода от события Ai к событию Ar

равна

 

 

 

s+m

 

s

 

 

 

 

P (Ar

Asi+)m= pir(m),

s+n

равна

а вероятность перехода от события Ar

к событию Aj

 

P (Ajs+n Ars+m) = prj(n − m).

 

 

По формуле полной вероятности

7.2 имеем

 

 

 

k

 

 

 

 

 

 

r

pir(m) · prj(n − m).

 

 

 

pij(n) =

 

 

 

=1

 

 

 

 

 

Последняя формула в точности совпадает с формулой произведения матриц, поэтому получаем πn = πm · πn−m.

Теорема 25.7. Для любого n ≥ 1, имеет место формула πn = π1n.

Доказательство. Доказательство проведём методом математической индукции.

1)Проверяем тождество при n = 1: π1 = π11 = π1.

2)Пусть при n = q выполнено тождество πq = π1q.

3)По теор. 25.6 имеем πq+1 = π1 · πq = π1 · π1q = π1q+1.

Что и требовалось доказать.

§26. Последовательности случайных величин

Замечание 26.1. Напомним, что по опред. 9.1 случайная величина есть функция ξ : Ω R. Поэтому последовательность случайных величин

короче

ξ1, ξ2, . . . , ξn, · · · = n}n=1 = n}

96

есть на самом деле последовательность функций {ξn : Ω R}n=1, определённая на одном и том же пространстве элементарных событий Ω.

В матанализе мы изучили следующие две сходимости последовательности функций.

Определение 26.2. Последовательность функций ξ1, ξ2, . . . , ξn, . . . сходится к функции ξ поточечно, если для любой точки ω Ω последовательность чисел (значений) ξ1(ω), ξ2(ω), . . . , ξn(ω), . . . сходится к значению ξ(ω)

функции ξ; или короче28 lim ξn(ω) = ξ(ω) для любой точки ω Ω.

n→∞

Определение 26.3. Последовательность функций ξ1, ξ2, . . . , ξn, . . . сходится к функции ξ поточечно почти всюду в Ω, если подмножество A Ω, в которых поточечная сходимость не выполняется, имеет меру 0, т.е. µ(A) = 0.

В теории вероятностей объекты ω, A, Ω являются событиями, а мерой их наступления является вероятность, поэтому в теории вероятностей сходимости «почти всюду» соответствует так называемая сходимость «почти наверное».

Определение 26.4. Говорят, что последовательность случайных величин n} сходится почти наверное к случайной величине ξ, если имеют место

эквивалентные друг другу равенства

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)} = 1

 

P {

 

 

 

 

 

ω

 

 

lim ξ

n(

ω

) =

ξ

ω

и

ω

 

lim ξ

(ω) = ξ ω = 0.

P {

 

 

n→∞

 

(

 

 

 

n→∞

n

̸ ( )}

Замечание 26.5. Чтобы пользоваться на практике сходимостью почти наверное необходимо знать, как устроены отображения ω 7→ξn(ω). Как пра-

вило в задачах теории вероятностей известны не сами случайные величины

x

ξn, а их функции распределения, скажем P(ξn ≤ x) = Fξn (x) = fξn (t)dt.

−∞

Можно ли в таком случае, обладая информацией только о функциях распределения, каким-нибудь образом исследовать сходимость последовательности случайных величин n}? Ответ: да, можно, если потребовать, чтобы вероятность тех элементарных исходов ω, для которых значение ξn(ω) не попадает в ε-окрестность числа ξ(ω), сходилась к нулю при n → ∞. В теории вероятностей эту идею реализуют с помощью так называемой сходимости по вероятности.

Определение 26.6. Говорят, что последовательность случайных величин n} сходится по вероятности к случайной величине ξ, если для любого ε > 0 имеют место эквивалентные равенства

nlim P ( n − ξ| ≥ ε) = 0

и

nlim P ( n − ξ| < ε) = 1. ()

→∞

 

→∞

28Используя определение сходимости числовой последовательности.

97

Если n} сходится по вероятности к ξ, то вместо пределов () коротко

P P

пишут: ξn → ξ при n → ∞, а т.к. n → ∞ всегда, то ещё короче: ξn → ξ. Хотя на практике для проверки такой сходимости проверяют выполнение одного из равенств ().

Пример. Рассмотрим последовательность случайных величин n}n=1, в которой для каждого n случайная величина задана следующим рядом распределения

ξn

0

n3

.

P

1 1/n

1/n

Докажем, что эта последовательность сходится по вероятности к вырожденной случайной (детерминированной) величине ξ, имеющей ряд распреде-

ления

ξ

0

, а проще говоря, сходится по вероятности к нулю.

P

1

 

 

 

 

 

 

Действительно, зафиксируем произвольное ε > 0. Для всех n, начиная с некоторого n0 такого, что n30 > ε, выполнено равенство ниже

( )

P (

ξn = n3)

1

 

P ( n 0| ≥ ε) = P ( ξn ≥ ε) =

=

 

.

n

Применяя теперь опред. 26.6 получаем

 

 

 

 

 

 

nlim P ( n 0| ≥ ε) = nlim

1

= 0,

 

 

 

 

 

 

 

n

 

 

 

→∞

→∞

 

 

 

 

 

т.е. последовательность случайных величин ξ1, ξ2, . . . , ξn, . . . сходится по ве-

P

роятности к ξ, т.е. ξn → ξ.

В нашем примере случайные величины ξn с ростом n могут принимать всё б´oльшие и б´oльшие значения, но с всё меньшей и меньшей вероятностью.

Замечание 26.7. Сходимость по вероятности может не сопровождаться сходимостью матожиданий (или других начальных и центральных момен-

P

тов): из ξn → ξ не следует, что Mξn Mξ. Например, в предыдущем примере

P

= n2 и Mξ = 0, и ясно, что последовательность 1, 4, 9, ...

ξn → ξ, однако Mξn

не сходится к нулю.

 

Сходимость по вероятности обладает теми же свойствами, как и поточечная сходимость.

Теорема 26.8. Если

P

1) ξn + ηn → ξ + η,

P

2) ξn · ηn → ξ · η.

P P

ξn → ξ и ηn → η, то

98

P

3) Если g — непрерывная функция, то g(ξn) → g(ξ).

P P

4) Если ξn → C и функция g непрерывна в точке C, то g(ξn) → g(C).

(Без доказательства.)

Другой тип сходимости случайных величин, который нам понадобится в дальнейшем, определяется через поточечную сходимость функций распределений. Это так называемая слабая сходимость. Пусть задана последовательность случайных величин n}n=1 с функциями распределения {Fξn (x)}n=1, и задана случайная величина ξ с функцией распределения Fξ(x).

Определение 26.9. Говорят, что последовательность случайных величин n}n=1 слабо сходится к случайной величине ξ, если последовательность их функций распределения {Fξn (x)}n=1 сходится поточечно к функции распре-

деления Fξ(x), т.е. lim Fξn (x) = Fξ(x).

n→∞

Если n} слабо сходится к ξ, то коротко пишут ξn = ξ. На практике для

проверки

слабой

сходимости проверяют

выполнение равенства

nlim Fξn (x) = Fξ(x) для любого x R.

 

→∞

 

P

ξ,

Теорема 26.10. 1) Если ξn → ξ, то ξn =

2) Если

 

P

 

ξn = C = const, то ξn → C.

 

3) Если

P

и ηn = η, то ξn · ηn = Cη.

ξn → C

4) Если

P

и ηn = η, то ξn + ηn = C + η.

ξn → C

(Без доказательства.)

99

100