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

фарафонов теория вероятностей

.pdf
Скачиваний:
321
Добавлен:
14.03.2016
Размер:
1 Мб
Скачать
☆

2. M[C×ξ]=С× M[ξ];

(7.8)

3. M[ξ+ η]= M[ξ]+ M[η].

(Подробнее о сумме случайных величин см. п. 8). Математическое ожидание случайной величины η, являющейся

функцией случайной величины ξ, т.е. η= ϕ(ξ), определяется формулой

n

 

M[η]= mη = åϕ(xi)× pi.

(7.9)

i=1

 

Определение. Дисперсией Dξ случайной величины ξ называется неслучайная величина, определяемая формулой

D = M[ξ2

 

n

(x -m )2

 

 

]=

å

× p ,

(7.10)

ξ

 

i

ξ

 

i

•

 

i=1

 

 

 

 

 

 

 

 

 

 

 

 

где ξ = ξ-mξ – центрированная случайная величина. Для вычис-

лений дисперсии удобен другой вариант формулы (7.10):

 

•

2

 

 

 

2

 

D = M[(ξ) ]= M[(ξ- M

) ]=

 

ξ

 

 

 

ξ

 

 

 

= M[ξ2]- M2[ξ]= n x2

× p -m2

(7.11)

 

 

 

å i

 

i

ξ

 

i=1

Свойства дисперсии:

1. D[С] = 0, если С – const;

2. D[С×ξ]=С2 × D[ξ]; (7.11а) 3. D[ξ+ η]= D[ξ]+ D[η], если ξ, η – независимые случайные ве-

личины.

Величина σξ = Dξ называется средним квадратическим от-

клонением случайной величины ξ.

Вероятностный смысл Мξ и Dξ заключается в том, что математическое ожидание дает среднее значение случайной величины, а дисперсия (и среднее квадратическое отклонение) характеризует степень ее отклонения от среднего значения.

Начальные моменты αk и центральные моменты βk определя-

ются в соответствии с выражением (7.9) формулами:

α

k

= M[ξk]=

n

xk × p ,

 

(k=1,2,¼);

 

 

 

 

å i

i

 

 

 

 

 

 

 

i=1

 

 

 

 

 

 

 

= M[ξk]=

n

(x -m )k × p ,

 

 

β

k

å

(k=1,2,¼).

(7.12)

 

 

i

ξ

i

 

 

i=1

41

7.3. Непрерывная случайная величина

Определение. Случайная величина ξ называется непрерывной, если ее функция распределения F(x) непрерывна при всех значениях x.

Для любого из возможных значений x непрерывной случайной величины ξ выполняется условие (7.3); поэтому имеет смысл лишь вычисление вероятностей p(ξ B), где B=[a, b) конечный интервал (или объединение таких интервалов).

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

Для вычисления указанной вероятности используется общая формула (7.4).

Определение. Случайная величина ξ называется абсолютно непрерывной (распре­ ­деленной с некоторой плотностью), если существует такая интегрируемая, неотрицательная функция f(x), что при любом x:

x

 

F(x)= ò f(x)dx.

(7.13)

-¥

Функция f(x), называемая плотностью распределения вероятностей, равна

 

 

f(x)= dF(x),

(7.14)

 

 

dx

 

и удовлетворяет условиям, следующим из свойств F(x):

 

1)

f(x)³0, так как F(x) – неубывающая;

 

 

¥

 

 

2)

ò f(x)dx0

=1, так как p(–∞ < ξ < ∞) = 1;

 

 

-¥

 

 

3)

f(±¥)=0.

 

 

Обычно непрерывной называют абсолютно непрерывную случайную величину. В соответствии с формулами (7.4) и (7.13) для вероятности попадания случайной величины в интервал (a£ξ<b) имеет место равенство

b

p(a£ξ<b)= ò f(x)dx. a

42

Определение. Математическим ожиданием непрерывной случайной величины ξ называется неслучайная величина, определяемая равенством

ξ

¥

 

 

ò

 

 

M[ξ]=m =

 

x× f(x)dx.

(7.15)

-¥

Математическое ожидание функции η = ϕ(ξ) находится по аналогичной формуле

 

 

¥

 

 

 

 

 

M[ϕ(ξ)]= ò ϕ(x)× f(x)dx.

(7.16)

 

 

-¥

 

 

 

 

Для дисперсии Dξ имеем

¥

 

 

 

 

 

 

= M[ξ2]=

 

 

 

)2 × f(x)dx

 

D

ò

(x-m

(7.17)

ξ

 

 

 

ξ

 

 

 

 

-¥

 

 

 

 

 

или

 

 

 

 

 

 

 

D = M[ξ2]- M2[ξ]=

¥

x2

× f(x)dx-m2.

 

ò

(7.18)

ξ

 

 

 

 

ξ

-¥

Свойства математического ожидания и дисперсии определяются прежними соотношениями (7.8) и (7.11а).

Начальные αk и центральные моменты βk определяются соответственно выражениями:

¥

 

 

αk = M[ξk]= ò xk × f(x)dx

( k=1,2,¼);

 

-¥

 

 

¥

 

 

βk = M[ξk]= ò (x-mξ)k × f(x)dx

(k=1,2,¼).

(7.19)

-¥

7.4. Функция случайных величин

Величина η = ϕ(ξ), являющаяся функцией случайной величины ξ, также является случайной величиной.

Основные задачи сводятся к определению ее характеристик – законов распределения и числовых характеристик по аналогичным характеристикам случайного аргумента.

Пусть ξ – дискретная случайная величина, заданная законом (рядом) распределения:

43

xi, pi = p(ξ = xi), (i = 1, 2, …, n).

Учитывая, что возможные значения yi случайной величины η

определяютсясоотношениями yi =ϕ(xi) и p(η= yi)= p(ξ = xi)= pi , получаем закон распределения для η:

y =ϕ(x ),

p = p(ξ = x ), (i = 1, 2, ... , n).

(7.20)

i

i

i

i

 

Замечание. При неоднозначности функции ϕ(ξ)нескольким xi может отвечать одно yj (j = 1, 2, …, m, m < n). В этом случае соответствующие вероятности pi следует просуммировать.

Для числовых характеристик справедливы формулы, вытекаю-

щие из (7.9):

 

 

 

 

n

 

 

 

 

 

 

M[η]=m =

å

ϕ(x )× p ;

 

(7.21)

 

 

η

 

 

i

i

 

 

 

 

 

 

i=1

 

 

 

 

 

 

2

n

 

 

 

 

2

 

 

D

å

[ϕ(x )-m

× p .

(7.22)

= M[η ]=

]

η

 

 

 

i

η

 

i

 

i=1

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

Биномиальное распределение (распределение Бернулли)

Распределением Бернулли называется распределение случайной величины x – числа успехов в серии из n независимых одинаковых испытаний. Вероятности отдельных значений x = 0,1,2, ..., n находятся по формуле Бернулли

 

 

p (k)=

kpkqn-k,

 

 

 

 

n

Cn

 

 

 

все вместе они образуют ряд распределения:

 

 

 

 

 

 

 

 

 

 

xk

0

1

 

...

k

....

n

p

qn

npqn–1

 

...

kpkqn-k

...

pn

k

 

 

 

 

Cn

 

 

Условие нормировки для распределения Бернулли следует из формулы (бином Ньютона):

n

n

å pk =åCnkpkqn-k =(p+q)n =1n =1.

k=0

k=0

Числовые характеристики биномиального распределения можно найти разными способами. Можно прямо воспользоваться фор-

мулами (7.7), (7.11).

44

Альтернативный способ вычислений заключается в том, что случайную величину x представляют как сумму независимых случайных величин xk:

n

ξ = åξk. k=1

Пусть каждая случайная величина xk отвечает отдельному испытанию. Очевидно, что xk равна нулю с вероятностью q = 1 – p и равна единице с вероятностью p. Числовые характеристики xk определяются следующим образом:

M[ξk]=0×(1- p)+1× p= p , D[ξk]=02 ×(1- p)+12 × p- p2 = pq .

Математическое ожидание самой случайной величины x равно сумме математических ожиданий xk:

n

n

Mξ = åMξk

= å p= np,

k=1

k=1

а дисперсия (учитывая независимость отдельных xk) равна сумме дисперсий xk:

n

n

Dξ = åDξk

= å pq= npq.

k=1

k=1

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

Распределение Пуассона p(k)= λk e-λ широко применяется в те- k!

ории массового обслуживания. Например, оно определяет вероятность числа k вызовов на АТС в единицу времени при средней интенсивности вызовов l. Другими примерами могут служить поступление вызовов на пункт медицинской помощи, клиентов на предприятие бытового обслуживания, последовательность отказов элементов.

Ряд распределения имеет вид:

xk

0

1

...

 

k

...

 

 

 

 

 

 

 

 

 

pk

e

–l

le

–l

...

λk

-λ

...

 

 

k!

e

 

 

 

 

 

 

 

 

Математическое ожидание и дисперсию пуассоновского распределения можно найти, и они получаются равными l:

45

Mξ = Dξ = λ.

Результат можно проверить прямым счетом по формулам (7.7), (7.11). Например, для математического ожидания имеем

¥

λk -λ

¥

λk-1 -λ

λ -λ

 

Mξ = åk

 

e

= λå

 

 

e

= λe e

= λ.

k!

(k-1)!

k=0

 

k=0

 

 

 

 

 

 

 

 

 

 

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

В теории массового обслуживания время обслуживания одной заявки распределено по показательному закону, плотность вероятности которого имеет вид (μ > 0)

ì

-µx

 

ï

 

x³0,

ïµe ,

fξ(x)=í

,

x<0.

ï0

ï

 

 

î

 

 

График плотности вероятности приведен на рис. 7.1.

fξ(x)

µ–

0x

Рис. 7.1. График плотности вероятности показательного закона

Условие нормировки плотности вероятности легко проверяется:

¥

¥

¥

 

¥ =1,

 

ò

fξ(x)dx= ò

µe-µxdx= ò e-zdz =-e-z

 

-¥

0

0

 

0

 

 

где z = x, dz = dx.

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

x

x

µx

 

µx =1-e-µx,

 

Fξ(x) = ò fξ(y)dy= ò

µe-µydy= ò e-zdz=-e-z

 

-¥

0

0

 

0

 

 

где z = y, dz = dy.

Найдем теперь математическое ожидание и выясним содержательный смысл параметра μ. Математическое ожидание (среднее время обслуживания одной заявки) найдем интегрированием по частям:

46

 

¥

-µx

1¥

-z

1

 

-z

 

¥

¥

-z

1

 

-z

 

¥ 1

 

 

 

 

 

 

 

Mξ =

ò

xµ dx=

µò

ze dz=

µ

(-ze

 

 

0

ò

e dz)=-

µ

e

 

 

0

µ

,

 

 

 

 

 

+

 

 

 

 

=

 

 

0

 

0

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где z = x, dz = dx.

Таким образом, среднее время обслуживания обратно пропорционально μ. В свою очередь μ = 1/Мξ – среднее число заявок, обслуженных в единицу времени, или интенсивность обслуживания.

Дисперсию показательной случайной величины находим, дважды интегрируя по частям:

D = ¥ x2

µe-µxdx-

1

=

1

¥ z2e-zdz-

1

=

2

-

1

=

1

,

 

 

 

 

 

 

ξ

ò

µ2 µ2 ò

µ2 µ2 µ2 µ2

 

 

0

0

 

 

 

 

 

 

 

 

где z = x, dz = dx.

Нормальный закон распределения

Нормальное распределение (распределение Гаусса) занимает в теории вероятностей особо важное место, поскольку согласно центральной предельной теореме достаточно большая сумма сравнительно малых независимых случайных величин имеет распределение, близкое к нормальному. Плотность вероятности нормального закона имеет вид

 

 

 

 

 

 

 

 

 

 

1

 

-

(x-a)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2σ

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (x)=

 

 

 

 

 

 

e

 

 

 

, -¥< x<¥.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ

 

 

σ

 

 

2π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Проверим условие нормировки:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1æx-aö2

 

 

 

 

 

 

 

 

 

 

z

2

 

¥

 

 

 

 

 

¥

 

-

ç

 

 

÷

 

 

 

 

 

 

 

 

¥

 

 

 

 

 

 

 

 

1

 

 

 

ç

 

÷

 

æx-aö

 

 

1

 

-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ç

 

÷

 

 

 

 

 

 

 

 

ò

 

 

 

 

 

ò

 

 

 

2è

σ

ø

 

ç

 

 

÷

 

 

 

ò

 

2

 

 

fξ(x)dx=

 

 

e

 

 

 

 

 

 

 

ç

 

 

÷

=

 

 

e

 

 

dz =1,

 

 

 

 

 

 

 

 

 

 

 

dç

σ

÷

 

 

 

 

 

-¥

 

 

 

2π

-¥

 

 

 

 

 

 

 

 

 

è

 

ø

 

 

2π

-¥

 

 

 

 

где z =

x-a

, dz =

1 dx.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ

 

 

σ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

¥

-

z2

 

 

 

ò

e 2 dz = 2π.

-¥

 

 

 

 

 

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

47

параметр а – точка максимума плотности и центр симметрии; σ определяет крутизну кривой и величину максимума.

fξ(x)

σ1<σ

a<0 a>0

x

Рис. 7.2. Семейство плотностей нормальных кривых

Выясним теперь теоретико-вероятностный смысл параметров. Найдем математическое ожидание случайной величины ξ . С учетом замены переменной z = (x–a)/σ имеем

 

 

 

 

 

 

 

 

 

 

1æx-aö2

 

 

 

 

 

 

 

 

 

¥

 

ç

 

 

 

÷

 

 

 

 

 

 

1

 

 

 

- ç

 

÷

 

 

 

 

 

 

 

 

 

 

ç

 

÷

 

 

M[ξ]=

 

 

 

ò x×e

2è

 

σ ø dx=

 

 

 

 

 

 

σ 2π

 

 

 

 

 

 

 

 

-¥

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ

 

¥

 

 

z2

 

 

 

 

¥

z2

 

 

ò ze-

 

 

 

 

a

 

ò e-

 

 

=

 

2 dz+

 

2 dz =0+a= a.

 

 

 

 

 

 

 

 

 

2π

 

 

 

2π

 

 

-¥

 

 

 

 

 

 

-¥

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Первый интеграл равен нулю, как интеграл от нечетной функции по симметричному относительно нуля интервалу, второй интеграл сводится к интегралу Пуассона.

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

При вычислении дисперсии необходимо проинтегрировать по частям дважды:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1æx-aö2

 

 

 

 

 

 

 

 

 

 

z

2

 

 

 

 

 

¥

 

 

 

 

-

ç

 

÷

 

 

 

 

 

2

¥

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

ç

÷

 

 

 

σ

 

 

 

-

 

 

 

D[ξ]=

 

 

 

 

 

ò (x-a)2

×e 2è

σ ø

dx=

 

 

 

 

ò z2e 2 dz =

 

 

 

 

 

 

 

 

σ 2π

 

 

 

2π

 

-¥

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-¥

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

 

z2

 

 

 

¥

 

z2

 

 

 

ö

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

ç

 

 

 

 

 

¥

 

 

 

 

 

 

÷

 

 

 

 

 

 

 

 

 

 

 

 

σ

 

-

 

 

 

 

 

 

 

-

 

 

 

 

÷

 

 

 

 

 

 

 

 

 

 

 

 

 

ç

 

2

 

 

 

 

 

 

 

 

2

 

 

 

÷

 

 

2

 

 

 

 

 

=

 

 

 

 

 

 

çç-ze

 

 

 

 

-¥

+ ò

e

 

dz

÷÷

= σ

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ç

 

 

 

 

 

 

 

 

 

 

 

 

÷

 

 

 

 

 

 

 

 

 

 

 

 

 

2π ç

 

 

 

 

 

 

 

-¥

 

 

 

 

 

÷

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ø

 

 

 

 

 

 

 

Отсюда σ2 – это дисперсия нормальной случайной величины, а σ – ее среднее квадратическое отклонение.

48

Определение вероятности попадания нормальной случайной величины в интервал требует введения специальной функции – функции Лапласа (см. приложение, табл. 2):

1 x -t2

Φ(x)= 2π ò0 e 2 dt,

поскольку интеграл от нормальной плотности не выражается в известных элементарных функциях.

С использованием этой функции получаем

 

 

 

 

 

β

 

 

β

 

 

1

 

-

(x-a)2

 

 

 

 

 

 

 

 

 

 

 

 

2σ

2

 

p(α <ξ<β)=

ò

f

(x)dx=

ò

 

 

 

 

 

e

 

 

dx=

 

 

 

 

 

 

 

 

 

 

 

 

ξ

 

 

σ

 

2π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

 

 

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

β-a

 

z2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

σ

-

æ

β

 

 

ö

 

æ

 

 

ö

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

ç

-a÷

 

ç

α-a÷

=

 

 

ò

e

 

dz = Φç

 

 

 

÷

-Φç

 

 

÷.

 

 

 

 

 

ç

 

σ

÷

 

ç

σ

÷

 

 

2π α-a

 

 

è

 

ø

 

è

ø

σ

При решении практических задач можно либо заново проделать подобные выкладки, либо пользоваться окончательным результатом:

æ

 

ö

æ

 

ö

ç

β-a÷

ç

α-a÷

p(α <ξ<β)= Φç

 

÷

-Φç

 

÷.

ç

σ

÷

ç

σ

÷

è

ø

è

ø

Если интервал является симметричным относительно математического ожидания, то формула упрощается:

p(а-t×σ<ξ< a+t×σ)= 2Φ(t).

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

стро приближается к своим предельным значениям, например 2Φ(3)»0,9973, то справедливо так называемое правило «трех сигм»: теоретически нормальная плотность вероятности отлична от нуля при любых значениях x, однако практически значения случайной величины ξ сосредоточены на отрезке a±3σ:

p(а-3×σ<ξ< a+3×σ)= 2Φ(3)»0,9973.

Вероятность попадания вне этого отрезка всего 0,0027.

49

Во-вторых, формула позволяет строить интервал с заданной вероятностью попадания – подобную задачу часто приходится решать в математической статистике.

Интервал, симметричный относительно математического ожидания, вероятность попадания в который равна числу p, строится следующим образом.

Вначале, пользуясь таблицей функции Лапласа, находим tp из условия 2Φ(tp)= p. Затем строим интервал, удовлетворяющий заданному условию:

p(а-tp ×σ<ξ< a+tp ×σ)= 2Φ(tp)= p.

7.6. Решение типовых примеров

Пример 7.1. Известна вероятность события А: р(А) = 0,5. Случайная величина ξ – число появлений А в трех опытах. Требуется построить ряд распределения случайной величины ξ в трех опытах; найти ее математическое ожидание и дисперсию.

Решение. Ряд распределения ξ мы построим, если найдем все возможные значения ξи соответствующие вероятности. В трех опытах число успехов может быть 0, 1, 2 или 3. Вероятности определим по формуле Бернулли:

p(ξ =0)=C03(0,5)0(0,5)3-0 =0,125; p(ξ =1)=C13(0,5)1(0,5)3-1 =0,375;

p(ξ = 2)=C32(0,5)2(0,5)3-2 =0,375; p(ξ =3) =C33(0,5)3(0,5)3-3 =0,125.

Получим ряд распределения

Xk

0

1

2

3

pk

0,125

0,375

0,375

0,125

Причем найденные вероятности удовлетворяют условию нормировки:

0,125+ 0,375 + 0,375 +0,125 = 1.

Числовые характеристики величины ξ можно найти непосредственно по формулам (7.7), (7.11):

M[ξ]=0×0,125+1×0,375+2×0,375+3×0,125=1,5;

50