Материал: Чернова Н.И. - Теория вероятностей - 1999

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

5.6Теорема Пуассона для схемы Бернулли

Предположим, нам нужна вероятность получить не менее десяти успехов в 1000 испытаний схемы Бернулли с вероятностью успеха 0.003. Вероятность этого события равна любому из следующих выражений:

1000

 

 

9

 

X

(0:003)k (0:997)1000 k = 1

 

X

(0:003)k (0:997)1000 k;

Ck

 

Ck

1000

 

1000

 

k=10

 

 

k=0

 

и вычисление даже одного слагаемого в каждом из этих выражений весьма проблематично.

Сформулируем теорему о приближенном вычислении вероятности какого-либо числа успехов в большом числе испытаний схемы Бернулли с маленькой вероятностью успеха. Термин «большое число» должен означать n ! 1. Если при этом p = pn 6!0, то, очевидно, вероятность получить любое конечное число успехов при растущем числе испытаний стремится к нулю. Необходимо чтобы вероятность успеха p = pn ! 0 одновременно с ростом числа испытаний. Но от испытания к испытанию вероятность успеха меняться не может (см. определение схемы Бернулли).

Поэтому рассмотрим «схему серий»: есть

одно испытание

 

с вероятностью успеха p1

два испытания

,

с вероятностью успеха p2

: : :

, : : : ,

: : :

n испытаний

с вероятностью успеха pn

: : :

 

: : :

Вероятность успеха меняется не внутри одной серии испытаний, а от серии к серии, когда меняется общее число испытаний. Обозначим через n число успехов в n-й серии испытаний.

Теорема 16 (Теорема Пуассона).

Пусть n ! 1, pn ! 0 òàê, ÷òî npn ! > 0. Тогда для любого k > 0 вероятность полу- чить k успехов в n испытаниях схемы Бернулли с вероятностью успеха pn стремится

к величине k e : k!

P(

 

= k) = Ck pk

(1

 

p

)n k

!

 

k

e

ïðè n

! 1

; p

 

 

!

0

 

òàê, ÷òî np

n !

> 0:

 

 

k!

 

 

 

n

 

 

 

n n

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

Доказательство.

Положим n = n pn ! > 0. По свойству 4, Cnk

 

nk

при фикси-

k!

рованном k и при n ! 1. Тогда

 

 

n k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

Ck pk

(1

 

p

)n k = Ck

nk

1

 

n

 

 

 

6nk

nk

 

 

1

 

 

 

 

n

 

 

n

1

 

n

 

 

k

e : (8)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

k! nk

 

n

 

 

n

 

 

!

 

k!

n n

 

n

 

 

n nk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

|

 

 

 

 

 

{z

 

 

 

 

}|

 

 

 

{z

 

 

 

}

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

#

 

 

 

 

 

 

 

 

 

#

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 (8) мы использовали свойства nk ! k è 1

n

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

! e . Докажем последнее свой-

n

 

 

 

ñòâî:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

n

 

 

 

 

 

n

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

n2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln 1

 

= n ln 1

 

 

= n

 

 

+ O

 

 

! :

 

 

 

 

 

 

 

 

 

 

n

n

n

 

n2

 

 

 

 

 

Для доказательства теоремы осталось в формуле (8) воспользоваться свойством 5.

30

Определение 23. Пусть > 0 — некоторая постоянная. Набор чисел

k

 

k! e ; k = 0; 1; 2; : : :

называется распределением Пуассона с параметром .

Пользуясь теоремой 16, можно приближенно посчитать вероятность получить не менее десяти успехов в 1000 испытаний схемы Бернулли с вероятностью успеха 0.003, с вычисления которой мы начали. Поскольку n = 1000 «велико», а pn = 0:003 «ìàëî», òî, âçÿâ = npn = 3, можно написать приближенное равенство

 

 

9

 

 

 

9

3k

1

3k

 

 

 

X

 

 

 

X

 

X

 

 

 

 

 

1

 

Ck

(0:003)k (0:997)1000 k

 

1

 

 

e 3 =

 

 

 

e 3

=

 

1000

 

 

k!

k=10

k!

 

 

 

k=0

 

 

 

k=0

 

 

 

 

 

 

 

 

 

 

=

табличное значение

 

3(10) 0; 001: (9)

 

 

 

 

 

Осталось решить, а достаточно ли n = 103 «велико», а pn = 0:003 «мало», чтобы заме-

нить точную вероятность P( n = k) на приближенное значение k e . Для этого нужно уметь оценивать разницу между этими двумя вероятностями. k!

Следующую очень полезную теорему мы докажем в конце курса.

Теорема 17 (Теорема Пуассона с оценкой погрешности).

Пусть A f0; 1; 2; : : : ; ng — произвольное множество целых неотрицательных чисел, n — число успехов в n испытаниях схемы Бернулли с вероятностью успеха p, = n p. Тогда

P( n 2 A)

 

k!

e

=

 

Cnk pk (1 p)n k

 

k!

e

6 np2 = n

:

 

k A

k

 

 

k A

 

k A

k

 

 

2

 

X

 

 

 

 

X

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, теорема 17 предоставляет нам возможность самим решать, достаточно ли n «велико», а p «мало», руководствуясь полученной величиной погрешности. Какова же погрешность в формуле (9)?

 

1

3k

 

 

 

 

 

9

!

1 3k

 

 

 

 

P( 1000 > 10)

 

e 3

=

1

C1000k (0:003)k (0:997)1000 k

 

 

 

e 3

6

k!

k=10

k!

 

k=10

 

 

 

 

 

 

k=0

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

X

 

X

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

np

 

= 0;009

:

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

Погрешность не более 0,009 (при вероятности около 0,001 :-) ). Во всяком случае, можно утверждать, что искомая вероятность никак не больше, чем 0,01=0,001+0,009.

31

Раздел 6. Случайные величины и их распределения

6.1Случайные величины

Мы уже видели, что для очень многих экспериментов нет никаких различий в под- счете вероятностей событий, тогда как элементарные исходы в этих экспериментах очень различаются. Но нас и должны интересовать именно вероятности событий, а не структура пространства элементарных исходов. Поэтому пора во всех таких «похожих» экспериментах вместо самых разных элементарных исходов использовать, например, числа. То есть ввести соответствие (иначе говоря, отображение) между элементарными исходами и вещественными числами (с ними удобно работать).

Пусть имеется случайный эксперимент и задано вероятностное пространство h ; F; Pi.

Определение 24. Функция : ! R называется случайной величиной, если для любого x 2 R множество f < xg = f! : (!) < xg является событием, то есть принадлежит-алгебре событий F.

Замечание 9. Читатель, не желающий забивать себе голову абстракциями, связанными с -алгебрами событий и с измеримостью, может смело считать, что любое множество элементарных исходов есть событие, и, следовательно, случайная величина есть произвольная функция из в R. Никаких неприятностей на практике это обычно не влечет, так что все дальнейшее в этом параграфе можно пропустить. Полезно, тем не менее, помнить: каждая такая «уступка» себе существенно снижает ваши адаптивные способности к жизни.

Определение 25. Будем говорить, что функция : ! R является F-измеримой, если f! : (!) < xg принадлежит F для любого x 2 R.

Итак, случайная величина есть F-измеримая функция, ставящая в соответствие каждому элементарному исходу ! 2 число (!) 2 R.

Пример 23. Подбрасываем 1 раз кубик. Пусть = f1; 2; 3; 4; 5; 6g, и две функции из

в R заданы так: (!) = !, (!) = !2.

Если F есть множество всех подмножеств , то и являются случайными вели- чинами, поскольку любое множество элементарных исходов принадлежит F, в том числе и f! : (!) < xg или f! : (!) < xg. Можно записать соответствие между значениями случайных величин и и вероятностями принимать эти значения в виде «таблицы распределения вероятностей» или, коротко, «таблицы распределения»:

 

1

2

3

4

5

6

 

 

1

4

9

16

25

36

P

1

1

1

1

1

1

 

P

1

1

1

1

1

1

6

6

6

6

6

6

 

6

6

6

6

6

6

 

 

 

Здесь 16 = P( = 1) = : : : = P( = 6) = P( = 1) = : : : = P( = 36).

Пусть -алгебра событий F состоит всего из четырех множеств:

F = ; ?; f1; 3; 5g; f2; 4; 6g ;

то есть событием является, кроме достоверного и невозможного событий, выпадение четного (соответственно, нечетного) числа очков. Убедимся, что при такой «бедной» -алгебре ни , ни не являются случайными величинами, так как эти функции не F-измеримы. Возьмем (например) x = 3;967. Видим, что f! 2 : (!) < 3;967g = f1; 2; 3g 62F и f! 2 : (!) < 3;967g = f1g 62F.

Упражнение. Описать класс всех функций, измеримых относительно -алгебры

F = ; ?; f1; 3; 5g; f2; 4; 6g .

32

Пусть -алгебра событий F есть тривиальная -алгебра : F = f ; ?g.

Доказать, что и не являются случайными величинами, так как эти функции не F-измеримы.

Доказать, что измеримы относительно тривиальной -алгебры только функции вида(!) = c (постоянные).

Теперь попробуем понять, зачем нужна F-измеримость и почему требуется, чтобы f! : (!) < xg являлось событием.

Если задана случайная величина , нам может потребоваться вычислить вероятности типа P( = 5) = Pf! : (!) = 5g, P( 2 [ 3; 7]), P( > 3;2), P( < 0) (и вообще самые разные вероятности попадания в различные множества на прямой). Это возможно только если множества, стоящие под знаком вероятности, являются событиями (напомню, что вероятность есть функция из -алгебры событий в [0,1]).

Но если потребовать, чтобы Ax = f! : (!) < xg было событием при любом x, то мы

из свойств -алгебры сразу получим, что

 

è

 

x = f! : (!) > xg — событие, и f! : x1 6 (!) < x2g = Ax2 nAx1

— событие,

A

 

1

 

 

 

n\

 

è

Bx = f! : (!) 6 xg = Ax+ n1 — событие,

 

 

=1

 

è

f! : (!) = xg = BxnAx — событие,

(10)

и т.д., и т.п. (операции пересечения, объединения, дополнения событий не выводят из класса событий).

Можно потребовать в определении 24 чего-нибудь другого. Например, чтобы событием было попадание в любой интервал: f! : (!) 2 (a; b)g 2 F для любых a < b. Или чтобы f! : (!) > xg было событием для любого x. Любое такое определение эквивалентно исходному.

Замечание 10. Те, кто не поленился прочесть про борелевскую -алгебру в разделе 3.3, могут сформулировать все наши потребности так: мы хотим, чтобы попадание в любое борелевское множество являлось событием. Мы могли это потребовать в определении, но ограничились эквивалентным условием, чтобы попадание в любой открытый интервал ( 1; x) было событием. Эти условия эквивалентны, поскольку борелевская-алгебра порождается интервалами, что мы еще раз показали в формулах (10).

Опишем различные типы распределений случайных величин. Под распределением случайной величины мы будем понимать соответствие

«значение случайной величины $ вероятность принимать это значе-

íèå»,

ëèáî (÷àùå)

«множество на прямой $ вероятность случайной величине попасть в это множество».

6.2Дискретные распределения

Определение 26. Говорят, что случайная величина имеет дискретное распределение, если существует конечный или счетный набор чисел fa1; a2; : : : g такой, что:

 

 

1

à) pi = P( = ai) > 0 äëÿ âñåõ i;

á)

P pi = 1.

 

 

i=1

То есть случайная величина имеет дискретное распределение, если она принимает не более чем счетное число значений.

Определение 27. Если случайная величина имеет дискретное распределение, назовем таблицей распределения соответствие ai $ pi, которое чаще всего рисуют так:

33

a1 a2 a3 : : :

Pp1 p2 p3 : : :

Примеры дискретных распределений

Вырожденное распределение.

Говорят, что случайная величина имеет вырожденное распределение с параметром a, и пишут = Ia, если принимает единственное значение a с вероятностью

1, то есть P( = a) = 1. Таблица распределения имеет вид

 

a

P

1

 

Распределение Бернулли.

Говорят, что случайная величина имеет распределение Бернулли с параметром p, и пишут = Bp, если принимает значения 1 и 0 с вероятностями p и 1 p, соответственно. Случайная величина с таким распределением равна числу

успехов в одном испытании схемы Бернулли с вероятностью успеха p (0 успехов

или 1 успех). Таблица распределения имеет вид

 

0

1

P

1 p

p

 

Биномиальное распределение.

Говорят, что случайная величина имеет биномиальное распределение с параметрами n и p, где 0 6 p 6 1, и пишут = Bn;p, если принимает значения 0; 1; : : : ; n с вероятностями P( = k) = Cnkpk(1 p)n k. Случайная величина с таким распределением имеет смысл числа успехов в n испытаниях схемы Бернулли с вероятностью успеха p.

Таблица распределения имеет вид

 

 

 

 

 

 

 

 

 

 

 

1

 

: : :

 

k

 

: : :

 

n

 

0

 

 

 

 

 

P

(1 p)n

np(1 p)n 1

 

: : :

 

Cnkpk(1 p)n k

 

: : :

 

pn

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

Говорят, что случайная величина имеет геометрическое распределение с пара-

метром p, где 0 6 p

6 1, и пишут = Gp, если принимает значения 1; 2; 3; : : :

с вероятностями P(

= k) = p(1 p)k 1. Случайная величина с таким распре-

делением имеет смысл номера первого успешного испытания в схеме Бернулли с вероятностью успеха p.

Таблица распределения имеет вид

 

1

2

: : :

k

: : :

P

p

p(1 p)

: : :

p(1 p)k 1

: : :

 

Распределение Пуассона.

Говорят, что случайная величина имеет распределение Пуассона с параметром ,

где > 0, и пишут = , если принимает значения 0; 1; 2; : : : с вероятностями

P( = k) = k e . k!

 

 

0

1

: : :

 

k

: : :

Таблица распределения имеет вид

 

 

 

 

k

 

 

P

e

e

: : :

 

 

e

: : :

 

 

k!

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

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

параметрами n, N и K, где K6N, n6N, если принимает целые значения от

Ck Cn k

maxf0; N K ng до minfn; Kg с вероятностями P( = k) = K N K . Случайная

CNn

величина с таким распределением имеет смысл числа белых шаров среди n шаров,

34