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

Лекции _ Вышка

.pdf
Скачиваний:
21
Добавлен:
05.03.2016
Размер:
905.57 Кб
Скачать

Пусть событие A – деталь стандартная. Условные вероятности этого собы-

тия по условию задачи будут: PB (A) = 0.8,

PB (A) = 0.7 . Тогда по формуле

1

2

полной вероятности имеем: P(A) = 12 ×0,8 + 12 ×0,7 = 12 ×1,5 = 12 × 32 = 34 .

19

P(A)

Лекция 4

Теорема Бейеса. Повторение испытаний. Формула Бернулли. Локальная и интегральная теоремы Муавра - Лапласа.

Английский математик Бейес Томас (1702-1761) поставил и решил одну из основных задач элементарной теории вероятностей, позволяющей

найти значения “апостериорных” вероятностей P (Bi / A) события Bi после наступления события A (“a posteriori – после опыта”) через значения “априорных” вероятностей P (Bi ) (a priori – до опыта). Эта теорема названа теоремой Бейеса и опубликована в 1763 году.

Теорема Бейеса

Т. Пусть событие A может наступить при условии появления одного из несовместных событий B1 , B2 ,..., Bn (гипотез), образующих полную группу

событий. Если событие A уже произошло, то вероятности гипотез могут быть переоценены по формулам Бейеса

PA (Bi ) = P(Bi )× PBi (A) , i = 1, n , где P( A) определяется по формуле

полной вероятности

n

P(A) = åP(Bi )× PBi (A)

i=1

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

Так как событие A может наступить при условии появления одного из несовместных событий B1 , B2 ,..., Bn , образующих полную группу, то вероятность события A определяется по формуле полной вероятности:

P(A) = P(B1 )× PB1 (A) + P(B2 )× PB2 (A) +...+ P(Bn )× PBn (A) .

Предположим, что произведено испытание, в результате которого появилось событие A . Поставим задачу: определить, как при этом изменились вероятности гипотез. То есть найдём вероятности PA (B1 ), PA (B2 ),..., PA (Bn ) .

По формуле умножения зависимых событий:

P(A× B1 ) = P(A)× PA (B1 ) = P(B1 )× PB1 (A)

Следовательно, PA (B1 ) =

P(B1 )× PB1 (A)

=

 

P(B1 )× PB1 (A)

 

P(A)

P(B1 )× PB1

(A) +...+ P(Bn )× PBn

(A)

 

 

Аналогично можно доказать для любого Bi :

19

PA (Bi ) =

 

P(Bi )× PBi (A)

, i = 1, n

 

P(A)

 

 

 

 

– формула Бейеса.

 

n

 

 

 

 

P(A) = åP(Bi )× PBi

(A)

i=1

 

 

Пример: Вероятность попадания в цель для 1-й группы стрелков – 0,8, для второй – 0,9. Один из стрелков выстрелил и попал в цель. Что вероятнее: стрелок был из 1-ой группы или из второй? В 1-ой группе 10 человек, во второй – 8 человек.

Рассмотрим гипотезы

B1 - стрелок был из первой группы.

 

 

 

B2 - стрелок был из первой группы.

 

 

 

P(B ) =

10

, P(B ) =

8

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

18

 

2

 

18

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть событие A - стрелок попал в цель. Условные вероятности этого

события

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

PB

(A) = 0.8, PB2 (A) = 0.9 . По формуле полной вероятности

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P(A) =

10

×0,8 +

8

×0,9 =

1

(8 + 7,2) =

15,2

=

38

.

18

 

 

 

 

 

 

 

 

 

 

18

 

18

 

 

 

18

45

 

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

Бейеса равна:

 

 

 

 

 

 

 

 

 

10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

×0,8×45

 

 

 

 

 

 

 

P (B ) =

P(B1 )× PB1 (A)

=

18

=

10

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

1

 

 

 

 

 

 

P(A)

 

 

 

 

 

 

 

 

 

38

 

19

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вероятность того, что попал стрелок из второй группы:

 

 

8

×0,9×45

 

 

 

7,2

 

 

 

9

 

 

 

 

 

 

 

 

P (B ) =

18

=

 

=

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

2

 

 

 

 

 

 

38

 

 

15,2

 

19

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, вероятнее, что стрелял стрелок из 1-ой группы.

Повторение испытаний. Формула Бернулли.

Если производятся несколько испытаний, причём вероятность A в каждом испытании не зависит от исходов других испытаний, то такие испытания называются независимыми относительно события A .

Предположим, что производиться n независимых испытаний, в каждом из которых событие A может появиться либо не появиться. Будем считать, что вероятность наступления события A в каждом испытании

равна p (0 < p < 1). Вероятность не наступления равна q = 1− p .

20

Последовательность независимых испытаний с двумя исходами впервые была исследована швейцарским математиком Якобом Бернулли (1654 - 1705)

и носит название последовательности испытаний Бернулли.

Поставим задачу: найти вероятность того, что в n испытаниях событие A появится ровно k раз Pn (k) .

Например, если идёт речь о появлении события A в четырёх испытаниях 3 раза, то возможны следующие исходы: AAAA, AAAA, AAAA,

AAAA .

Т. Если вероятность наступления события A в каждом испытании одинакова и равна p(0 < p <1) , то вероятность того, что в n испытаниях

событие A наступит ровно k раз определяется по формуле Бернулли.

 

 

 

 

 

 

 

 

 

 

 

n!

 

 

 

P (k) = Ck × pk ×qnk

 

 

, где Ck =

 

 

 

 

 

 

 

 

 

 

 

 

n

n

 

 

 

 

n

k!(n - k)!

 

 

 

 

 

 

 

 

B

 

A

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

Пусть

событие

– в

n

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

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

события

B

состоит

в появлении одного

из

несовместных событий

B1 , B2 ,..., Bm ,

 

вероятности которых можно найти по

формулам:

 

 

 

 

 

 

pk ×qnk ,

 

 

 

 

P(B ) =

p × p ×...× p ×q ×q ×...×q =

 

 

 

 

1

14243 14243

 

 

 

 

 

 

 

 

 

 

k

 

nk

pk

×qnk , ...

P(B ) = pk ×qnk .

 

P(B ) = q × p × p ×...× p × q ×q =

 

2

14243

{

 

 

 

 

 

m

 

 

 

 

k

nk −1

 

 

 

 

 

 

 

 

Число таких событий m равно числу сочетаний, которые можно составить из n элементов по k элементов. То есть m = Cnk .

P(B1 +...+ Bm ) = P(B1 ) + P(B2 ) +...+ P(Bm ) = Cnk × pk ×qnk

Pn (k) = Cnk × pk ×qnk – формула Бернулли

Следствия:

1)Вероятность того, что в n испытаниях событие наступит не более k раз:

P(1; k) = Pn (1) + Pn (2) +...+ Pn (k) .

2)не менее k раз: Pn (k;n) = Pn (k) + Pn (k +1) + ...+ Pn (n) .

3)более k раз: Pn (k +1; n) = Pn (k +1) +...+ Pn (n) .

4)менее k раз: P(1; k -1) = Pn (1) +...+ Pn (k -1) .

5)от k1 до k2 раз: P(k1 ;k2 ) = Pn (k1 ) + Pn (k1 +1) +...+ Pn (k2 ) .

Пример: Найти вероятность того, что при четырёх бросаниях монеты «герб» выпадет 3 раза.

21

Вероятность выпадения герба в одном испытании p =

1

;

q =

1

. По

2

2

 

 

 

 

 

формуле Бернулли

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P (3) = C3

×

æ

1

ö3 ×

æ

1

ö1 =

4!

×

1

×

1

=

4

=

1

 

 

 

 

 

 

ç

 

ç

 

 

 

 

 

 

 

 

 

 

 

4

4

 

2

÷

2

÷

3!×(4 -3)! 8

2

16 4

 

 

 

 

 

 

 

 

 

è

ø

è

ø

 

 

 

 

 

 

Локальная теорема Муавра-Лапласа.

Пользоваться формулой Бернулли при большом числе испытаний n довольно трудно. Локальная теорема Муавра-Лапласа даёт возможность приближенно вычислить вероятность того, что в n испытаниях событие A появится ровно k раз.

Эту теорему сначала доказал для частного случая при p = 12

английский математик Муавр Абрахам (1667-1754). В 1783 г. Лаплас обобщил ее для произвольного p . Рассмотрим эту теорему без

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

Т. Если вероятность появления события A в каждом испытании постоянна и отлична от нуля и единицы (0 < p < 1) , то вероятность Pn (k) того, что в n

испытаниях

событие

A

появится

 

ровно k

 

раз

 

приближенно можно

вычислить по формуле:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

- np

 

 

 

 

 

 

 

 

x2

 

Pn (k) »

 

1

 

ϕ(x)

 

 

x =

 

 

ϕ(x) =

 

1

 

 

e

 

 

 

,

где

 

,

 

 

 

2

.

 

 

 

 

 

 

 

 

 

 

npq

npq

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Составлена таблица значений ϕ(x) , соответствующая положительным значениям x . При вычислении функция ϕ(x) от отрицательных значений аргумента учитывается то, что функция ϕ(x) чётная, то есть ϕ(-x) = ϕ(x) .

Пример: Найти вероятность того, что событие A наступит 80 раз в 400 испытаниях, если вероятность появления в одном испытании p = 0, 2 .

n = 400, k = 80, q = 0,8

22

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

k - np

 

 

 

 

 

 

 

 

80 - 400×0,2

 

 

 

80 -

400×2

 

 

 

 

P400 (80) =

 

 

ϕ(x),

x =

 

,

 

 

x =

 

 

 

=

 

10

 

= 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

npq

npq

 

 

400×0,2×0,8

400×2×8

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

100

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ(0) =

 

1

 

= 0,3989

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P400

(80) =

 

 

 

 

1

 

 

 

×0,3989 = 0,04986

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

400

×2×8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

100

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если вычислять эту вероятность по формуле Бернулли,

то

получим Pn (k) = 0,0498 , то есть погрешность, как видим, незначительная.

 

 

 

 

 

 

 

 

Интегральная теорема Муавра-Лапласа

 

 

 

 

 

Т. Если вероятность наступления события

A

в каждом испытании

постоянна и отлична от нуля и единицы

 

 

(0 < p < 1) , то вероятность Pn (k1 ,k2 )

того, что событие наступит не менее k1

 

и не более

k2

раз (от

k1 до

k2 )

приближённо определяется по формуле:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x′′

 

 

 

 

 

 

 

 

 

 

1

 

x′′

 

z2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Pn (k1;k2 ) » òϕ

 

 

 

 

 

 

 

 

 

 

 

òe

2

 

¢¢

 

 

 

¢

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(z)dz =

 

 

 

 

 

 

 

dz = F(x ) - F(x )

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

¢

=

k1 - np

 

 

x

¢¢

=

 

 

k2 - np

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

npq ,

 

 

 

 

 

 

npq

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

x

 

 

z2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где F(x) =

 

 

òe

 

dz

 

 

– функция Лапласа.

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При решении задач

пользуются

 

специальными

таблицами

для

 

функции

 

 

F(x) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F(x) = òe

z2

 

Лапласа

 

 

так

как

неопределенный интеграл

 

2

dz

не

выражается через элементарные функции. Доказательство теоремы выходит за рамки данного курса, покажем лишь равенство определенного интеграла разности функции Лапласа в точках x′′ и x′ .

 

 

 

 

 

 

 

 

 

1

 

 

x′′

z2

 

 

 

1

 

0

z2

1

x′′

z2

 

 

 

 

 

 

 

 

 

 

ò

 

 

 

 

 

 

 

ò

 

 

 

 

ò

 

 

 

P (k ;k

2

) =

 

 

e 2 dz =

 

 

e 2 dz +

 

e 2 dz =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

1

 

x

z2

 

 

 

 

1

 

x′′

z2

 

 

 

 

 

 

 

 

 

0

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-

 

 

òe

 

dz

+

 

 

òe

 

 

dz

 

 

¢¢

 

 

¢

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= F(x ) - F(x )

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Функция Лапласа F(x)

нечётная, действительно:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

23

 

 

 

 

 

 

 

 

1

 

x

z2

 

z = -t

 

dz = -dt

 

-1

 

x

z2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F(-x) =

 

 

 

ò e

 

2 dz =

z

0

-x

=

 

 

 

òe

 

2 dt = -F(x) .

 

 

 

 

 

 

 

 

 

0

 

 

 

t

0

x

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример: Вероятность того,

что деталь не прошла проверку ОТК p = 0, 2 .

Найти вероятность того, что среди 400 случайно отобранных деталей окажется непроверенными от 70 до 100 деталей.

 

 

 

 

 

По интегральной теореме Муавра-Лапласа

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

¢¢

 

 

 

¢

 

 

 

 

 

¢¢

 

100 - 400×0,2

 

 

 

 

100 -80

 

 

 

 

 

 

 

P400 (70;100) = F(x ) - F(x ),

 

x

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

400

×0,2

×0,8

 

 

4×100×

 

 

16

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

100

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

20

 

 

10

 

 

 

 

5

 

x

¢

 

 

 

70 - 400×0,2

 

 

 

70 -80

 

-10

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 2×4 = 4 =

 

2 = 2,5,

=

 

400×0,2×0,8 =

2×4

 

 

= 8 = - 4 = 1,25,

 

 

 

 

 

 

 

 

 

P400 (70;100) = F(2,5) - F(-1,25) = F(2,5) + F(1,25) = 0,4938 + 0,3944 = 0,8882

Замечание 1. Пусть m – число появлений события

A в n испытаниях. Если

m изменяется от k1

до

k2

, то дробь

 

m - np

 

изменяется от x¢ =

k1

- np

 

до

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

npq

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

npq

x

¢¢

=

k2 - np

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

. Тогда интегральную теорему Муавра-Лапласа можно записать

 

 

 

 

 

 

 

 

 

 

 

npq

 

в виде:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

 

 

 

m - np

ö

1

 

x′′

e

z2

 

 

 

(

 

 

)

 

 

 

 

(

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

P ç x¢ £

 

£ x¢¢÷ »

 

 

 

2

dz = F

x¢¢

- F

x¢

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ò

 

 

 

 

 

 

 

 

 

 

 

 

ç

 

 

 

 

 

npq

÷

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

 

 

 

 

 

ø

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Замечание 2.

 

 

Таблица

 

F(x)

 

дана для

 

o £ x < 5 . При

 

x ³ 5 полагают

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F(x) = 0,5 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

24

Лекция 5

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

Теорема Пуассона. Найвероятнейшее число появлений события в независимых испытаниях

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

Пусть производится n независимых испытаний, вероятность появления события A в каждом испытании постоянна и равна p . Как

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

Поставим задачу: найти вероятность того, что относительная частота mn появления события A отклонится от p на величину, меньшую ε > 0 , то

æ

 

m

 

ö

 

 

 

 

 

 

 

 

 

 

есть найдем вероятность P ç

 

- p

£ ε ÷

, где m

– число появлений события

n

è

 

 

ø

 

 

 

A . Другими словами, найдем

вероятность

осуществления неравенства

 

m

 

 

 

 

£ ε . Воспользуемся замечанием 1

 

 

 

 

 

 

- p

к интегральной теореме Муавра-

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Лапласа.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

 

m

 

ö

 

æ

 

 

 

m

 

 

 

 

ö

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P ç

 

 

 

 

 

 

- p

÷

£ ε ) = Pç

-ε £

 

 

 

- p £ ε )÷ =

 

 

 

 

 

 

 

 

 

 

n

 

n

 

 

 

 

 

 

 

 

 

è

 

 

 

 

ø

 

è

 

 

 

 

 

 

 

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

-ε £

 

m - np

£ ε

ö

 

æ

 

-εn

 

 

m - np

 

 

εn ö

 

 

 

= P

ç

 

 

 

÷

= P ç

 

 

 

 

£

 

 

 

£

 

 

 

÷

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

ç

 

 

npq

 

 

 

npq

 

 

 

÷

 

 

 

 

 

 

 

è

 

 

 

 

 

 

ø

 

è

 

 

 

 

 

 

 

 

npq ø

 

æ

= P ç -ε

ç

è

n

 

£

m

- np

 

£ ε

 

n

pq

 

 

 

pq

 

npq

 

 

 

 

 

ö

÷ »

÷

ø

æ

ε

ç

è

n

pq

ö

æ

-ε

÷

- Fç

÷

ç

 

ø

è

 

 

 

ö

æ

 

 

 

 

ö

n

 

ε

n

 

÷

= 2

÷

 

 

 

÷

ç

 

 

 

÷

pq ø

è

 

 

pq ø

Таким образом, получаем формулу

æ

 

m

 

 

ö

æ

 

 

n

 

ö

 

 

 

 

 

P ç

 

 

- p

 

£ ε ÷

= 2

ε

 

÷

вероятность отклонения относительно

 

 

 

 

 

 

 

 

ç

 

 

 

÷

 

è

 

n

 

 

ø

è

 

 

pq ø

 

 

 

 

 

 

 

частоты от постоянной вероятности p в n независимых испытаниях.

25

Пример1: Вероятность того, что деталь нестандартна p = 0,1. Найти вероятность того, что среди случайно отобранных 400 деталей относительная

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

p не более чем на

0,03.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение: n = 400; p = 0,1; q = 0,9; ε = 0,03 .

 

 

 

 

 

 

æ

 

m

 

 

 

 

 

ö

æ

 

 

 

ö

æ

 

 

 

 

ö

 

 

 

 

 

 

 

 

400

 

 

 

4

 

 

P ç

 

 

 

 

- 0,1

£ 0,03

÷

» 20,03

 

÷

= 2

0,03×100

 

 

÷

=

 

 

 

 

è

 

400

 

 

 

 

 

ø

ç

 

0,1×0,9

÷

ç

 

 

9

÷

 

 

 

 

 

 

 

è

 

ø

è

 

 

ø

 

 

 

 

 

 

 

 

 

 

æ 0,03

 

 

2

ö

 

 

 

 

 

 

 

 

 

 

 

 

 

= 2Fç

 

 

×100×

 

÷

= 2F(2) » 2×0,4772 = 0,9544

 

 

 

 

 

 

 

10

3

 

 

 

 

 

 

è

 

 

 

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 2: Французский ученый Бюффон (XVII) бросил монету 4040 раз. Причем “герб” появился 2048 раз. Найти вероятность того, что при повторении опыта Бюффона относительная частота появления “герба” отклонится от вероятности появления “герба” по абсолютной величине не более, чем в опыте Бюффона.

Найдем частоту появления “герба” в опыте Бюффона: m = 2048,n = 4040 ,

 

m

=

2048

» 0.506930693 .

Тогда отклонение частоты от вероятности

 

n

4040

 

 

 

составляет ε »

 

 

 

 

выпадения “герба” p = 1/ 2

 

0.5 - 0.506930693

 

» 0.006930693.

 

 

Найдем вероятность того, что при бросании монеты повторно 4040 раз у нас относительная частота отклонится не более, чем на ε .

æ

 

 

m

 

 

 

ö

æ

 

 

 

 

ö

 

 

 

 

 

 

 

 

4040

 

 

P ç

 

 

 

- 0,5

£

0,006930693

÷

» 2

0,006930693

 

 

÷

=

 

 

è

 

4040

 

 

 

ø

ç

 

 

0,5×0,5

÷

 

 

 

 

 

è

 

 

ø

 

 

 

 

 

 

 

 

= 2F(0,881) » 2×0,3106 = 0,6212

Теорема Пуассона.

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

появится ровно k раз в общем случае определяется по формуле Бернулли. Однако, при больших значениях n , а также при малых значения p или q

вычисление по этой формуле становится затруднительным. Французский математик Пуассон Симеон Дени (1781-1840) вывел приближенную асимптотическую формулу расчета этой вероятности, применимую для

больших значений n и малых значений

p

(

p £

)

 

0.1 . Она определяет

вероятность массовых ( n велико) и редких ( p

мало) событий.

26

Т. Если n → ∞; p → 0; np → λ , то Pn (k) = Cnk pk qn-k ® λk

e-λ .

 

 

 

 

 

 

 

 

λn

 

k !

 

 

 

Доказательство: Обозначим λn = npn

® λ; p =

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

P (k) =

n!

pk qn-k =

(n - k +1)(n - k + 2)(n - k + 3)...n

æ

λn ök

æ1-

λn ön-k

=

 

 

 

÷

÷

n

k!(n - k)!

 

k!

 

ç

ç

 

 

 

 

è

n ø

è

n ø

 

 

 

(n - k +1)(n - k + 2)...n λnk

æ

 

λn

ön æ

 

 

 

 

 

λn

 

ö-k

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

ç1

-

 

÷ ç1

-

 

 

 

 

 

÷

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k!

 

 

 

 

 

 

 

 

 

 

 

n

n

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

 

ø è

 

 

 

 

 

ø

 

 

 

 

 

 

 

 

 

 

 

 

æ

 

k

 

 

1

öæ

 

 

 

 

k

 

 

2

ö

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ç1-

 

+

 

 

֍1-

 

 

 

+

 

 

÷...1

 

k

 

 

λn

 

n

 

 

 

 

 

 

 

λn

 

n

 

 

 

 

 

n

 

n

 

 

n

n

æ

 

ö

æ

 

 

 

 

 

 

ö

 

 

 

 

 

è

 

 

 

 

øè

 

 

 

 

 

 

 

ø

 

 

 

λn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ç1-

 

 

÷

 

 

ç1

-

 

 

 

 

 

÷

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k !

 

 

 

 

 

 

 

 

 

 

 

 

1

n

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

 

ø

 

 

è

 

 

 

 

 

 

 

ø

 

 

 

 

 

 

 

 

а) lim

æ

 

-

λn

ök

= 1,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ç1

 

 

n

 

÷

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n®¥

è

 

 

 

 

 

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b) lim

æ

 

 

 

λn

ön

 

 

 

 

 

 

 

æ

 

 

n - λn

 

ön

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

n

 

 

ç1

-

 

 

 

 

÷

 

= lim ç

 

 

 

 

 

 

 

 

÷

= lim(

 

 

 

 

 

 

 

 

 

 

 

)

 

 

=

 

 

 

n

 

n - λ + λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n®¥

è

 

 

 

 

 

ø

 

 

 

n®¥

è

ø

 

n®¥

1+

 

 

λn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n - λn

 

 

 

 

 

 

=

 

 

 

 

1

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

1

 

 

 

 

=

1

æ

 

 

 

 

λn

ö

n

 

 

 

 

 

 

 

 

 

 

 

 

 

n-

λn

 

λn

×n

 

 

 

lim

λn×n

 

 

eλn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n-λn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

 

 

 

 

 

 

ö λn n-λn

 

 

 

 

 

 

 

 

 

 

lim ç1+

 

 

 

÷

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

en®¥

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

n®¥ è

 

 

n - λn ø

 

 

 

 

 

 

 

 

ç

 

 

 

 

 

÷

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lim ç1+

 

 

 

 

÷

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

- λn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n®¥ ç

 

 

 

÷

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ç

 

 

 

 

 

 

÷

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

 

 

 

 

 

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1)

= e-λn

 

 

æ

k

1

öæ

 

k

2

ö

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ç1

- n

+ n

֍1

-

 

+ n

÷

...1

λk æ

 

 

λ

ön æ

 

λ

ök

 

λk -

 

 

 

n

 

 

 

 

 

lim P (k) = lim

è

 

 

øè

 

 

 

ø

 

n

ç

1

-

n

÷

ç

1-

n

÷

=

 

e

λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

k !

 

 

 

1

 

 

n

 

n

 

k!

 

 

n®¥

 

n®¥

 

 

 

 

 

è

 

 

ø

è

 

ø

 

 

 

Таким образом, при приближенных вычислениях вероятности массовых ( n велико) но редких ( p мало) событий можно воспользоваться формулой

Pn (k) » λk e-λ = np – формула Пуассона k!

Пример: Завод отправил на базу 5000 доброкачественных изделий. Вероятность того, что по пути изделие повредится, равна 0,0002. Найти вероятность того, что на базу прибудут 3 негодных изделия.

Решение: n = 5000; p = 0,0002; k = 3;λ = np = 5000×0,0002 = 500010000×2 =1.

27