Используя связь между плотностью распределения и функцией распределения, получим
|
|
|
x |
|
|
x |
|
|
|
F(x) = f (x)dx = f (x)dx. |
|||||||||
Тогда |
|
|
−∞ |
|
|
α |
|
|
|
|
|
|
|
|
|
|
|
||
x |
|
|
|
x − α |
|
|
|||
|
|
f |
(x)dx = |
, |
α < x ≤ β; |
||||
β − α |
|||||||||
α |
|
|
|
|
(6.2) |
||||
F(x) = |
β |
|
|
|
|
|
|
||
|
|
1 |
|
|
|
|
|
||
|
|
dx = 1, |
x > β. |
||||||
|
|
|
|||||||
β − α |
|||||||||
α |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|||
Числовые характеристики равномерно распределенной случайной величины
Математическое ожидание
|
|
|
|
|
|
+∞ |
|
β |
|
β − α |
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
mx = |
|
xf (x)dx = |
|
x |
1 |
dx = β + α . |
(6.3) |
||||
|
|
|
−∞ |
α |
|
|||||||||
Дисперсия |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
+∞ |
|
|
)2 |
|
β |
|
|
β + α )2 |
1 |
dx = (β − α)2 . |
|
|
D |
= |
|
(x − m |
x |
f (x)dx = (x − |
(6.4) |
||||||||
|
|
|||||||||||||
|
x |
|
|
|
|
|
|
|
2 |
|
β − α |
12 |
|
|
|
|
−∞ |
|
|
|
|
α |
|
|
|
|
|
|
|
В силу симметрии медиана совпадает с математическим ожиданием, а так как функция плотности распределения не имеет максимума, мода отсутствует.
Для получения вероятности попадания случайной величины в интервал (a, b) можно воспользоваться явным видом функции рас-
пределения (6.2):
P(a < x < b) = F(b) − F(a) = |
b − a |
. |
(6.5) |
|
|||
|
β − α |
|
|
На рис. 6.1. показано, что вероятность попадания случайной величины в заданный интервал есть площадь прямоугольника, огра-
ниченного прямыми x = a и x =b . Из выражения (6.5) видно, что вероятность попадания в интервал (a, b) равна отношению длины
интервала ко всей длине участка (α, β) и не зависит от его месторасположения на этом участке.
51
Закон Пуассона
Закон Пуассона, как мы знаем, выражается формулой
P |
(m) = |
λ me−λ |
. |
(6.6) |
n |
|
m! |
|
|
|
|
|
|
При этом параметр λ мы называли средним значением числа появившихся событий. Математическое ожидание случайной величины, распределенной по закон Пуассона, в соответствии с общим правилом (5.3) нужно определять так:
N
M [X ] = mX = xi P[X = xi ].
i=1
Но случайная величина X в данном случае может принимать только значения x = m = 0,1, 2... , тогда математическое ожидание
есть
∞ |
∞ |
λ |
m |
e |
−λ |
|
|
mx = mPn (m) = m |
|
|
. |
(6.7) |
|||
m=0 |
m=0 |
|
m! |
|
|
||
Можно показать, что mx = λ. |
Таким образом, параметр λ в за- |
||||||
коне Пуассона есть математическое ожидание числа появлений редких (т.е. вероятность каждого отдельного события мала), но массовых (т.е. опытов может быть много) случайных событий. Оказывается, что дисперсия случайной величины тоже равна λ :
∞ |
∞ |
|
m |
−λ |
|
|
D = |
(m − λ)2 Pn (m) = (m − λ )2 |
λ |
e |
|
= λ. |
(6.8) |
m=0 |
m=0 |
|
m! |
|
|
|
Рассмотрим более подробно смысл и условия применимости закона Пуассона на таком примере.
Пример. Пусть в качестве события рассматривается дорож- но-транспортное происшествие в городе Москве. Чтобы можно было применить закон Пуассона для определения вероятности заданного количества ДТП за определенный период времени, должны соблюдаться следующие условия:
вероятность количества ДТП за время τ зависит только от величины этого промежутка времени τ и не зависит от того, в каком месте временнóй оси этот промежуток времени находится. То есть вероятность того или иного числа ДТП в течение часа не
52
зависит от того, будет это в интервале, например, 8.00−9.00 или
20.00−21.00;
число ДТП за интервал τ не зависит от того, сколько ДТП произошло за предыдущий интервал τ;
вероятность двух и более ДТП за практически малый промежуток времени пренебрежимо мала по сравнению с вероятностью одного ДТП в этот промежуток. (Опытные водители сразу скажут, что это не совсем так и приведут массу примеров из своего опыта, но если под ДТП мы понимаем не количество столкнувшихся машин, сам факт происшествия, то наше приближение справедливо.)
Пусть в среднем в Москве за сутки происходит 240 ДТП. Понятно, что за один час это число будет равно 10. Вероятность этого события легко подсчитать с помощью функции Пуассона в Exel. Она равна P(10) = 0,125. Определим вероятность того, что
за час произойдет ровно в два раза большее число |
ДТП, т.е. |
20 ДТП. Значит, в условиях этого примера λ = 10, а |
m = 20. Ис- |
пользуя функцию Пуассона в Excel, получим P(20) = 0,0018. Отме-
тим, что по сводкам ГИБДД за 29.01.07 отмечается, что число ДТП было равно 641, т.е. более чем в два раза превосходило среднее число за сутки.
Вероятность события существенно увеличивается, если мы хотим знать вероятность появления числа ДТП из некоторого интервала. Например, вероятность того, что за час произойдет от
9 до 11 ДТП, есть P(9 ≤ m ≤ 11) = P(9) + P(10) + P(11) = 0,363.
Показательный закон распределения
Показательным называется закон распределения случайной величины, в данном случае обозначим ее t, имеющий плотность вероятности вида:
|
−λ t |
, t ≥ 0; |
|
λe |
|
(6.9) |
|
f (t) = |
t < 0. |
||
0, |
|
||
|
|
|
|
Показательным распределением описывается довольно большой круг задач: от задач надежности до задач ядерной физики.
Применительно к бизнесу можно привести такой пример. Понятно, что ни одна из компаний не застрахована от банкротства (не
53
имеется в виду страхование страховой кампанией). Пусть вероятность банкротства компании в единицу времени не зависит от вре-
мени существования компании и равна λ. Тогда вероятность банкротства за время dt есть λ dt. Пусть P(t) − вероятность того, что
компания не обанкротилась за время t своего существования. Тогда событие, состоящее в том, что компания обанкротится именно в интервал времени от t до t + dt есть произведение событий C = A B. Событие A заключается в том, что компания не обанкротится к моменту времени t, а событие B заключается в том, что за время dt компания обанкротится. По теореме умножения вероятностей для независимых событий получим, что вероятность того, что компания обанкротится в интервале времени dt вслед за моментом времени t есть P(t)λdt. Понятно, что изменение вероятности не
обанкротится за время от t до t + dt есть P(t + dt) − P(t) = dP. И эта
величина, в свою очередь, равна вероятности банкротства фирмы в интервал времени (t, t + dt). То есть
dP = −λ P(t)dt. |
(6.10) |
Знак «минус» в выражении (6.10) стоит потому, что dP − отри-
цательная величина. |
|
Если обе части выражения (6.10) разделить на dt ≠ 0, |
то получим |
дифференциальное уравнение первого порядка |
|
dP = −λ P(t) |
(6.11) |
dt |
|
с начальным условием P(0) = 1, поскольку обанкротится в самый
момент начала работы невозможно.
Легко проверить, что решением данного уравнения с учетом начального условия является функция P(t) = e−λ t . Таким образом, вероятность избежать банкротства за время t со дня основания есть
P(t) = e−λ t . |
(6.12) |
Вероятность обанкротится, соответственно, будет |
F(t) = |
= 1− P(t) = 1− e−λ t . Функция F(t) в данном случае есть функция распределения вероятности банкротства. Тогда плотность вероятно-
54
сти банкротства есть f (t) = dFdt = λ e−λ t . Таким образом, примени-
тельно к данному примеру выражение (6.9) описывает плотность вероятности банкротства.
Математическое ожидание случайной величины T, |
распреде- |
||||||||
ленной по показательному закону есть |
|
|
|
||||||
|
+∞ |
+∞ |
1 |
|
|
||||
M [T ]= mt = t f (t)dt = t λ e−λ t dt = |
. |
(6.13) |
|||||||
|
|||||||||
|
−∞ |
0 |
|
|
λ |
|
|||
|
|
|
|
|
|
||||
По своему смыслу m = |
1 |
есть среднее время жизни кампании. |
|||||||
|
|||||||||
t |
|
λ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Легко показать, что дисперсия случайной величины T есть |
|
||||||||
D[T ] |
∞ |
|
1 |
|
|
|
|
||
= (t − mt )2 f (t)dt = |
. |
|
|
(6.14) |
|||||
2 |
|
|
|||||||
|
0 |
|
λ |
|
|
|
|||
|
|
|
|
|
|
|
|||
Геометрическая интерпретация показательного закона распределения приведена на рис. 6.3.
Рис. 6.3. Плотность распределения показательно распределения
Площадь заштрихованной области есть вероятность компании распасться за время от t1 до t2 , то есть F(t1) − F(t2 ) = e−λ t1
Лекция 7 . НОРМАЛЬНЫЙ ЗАКОН РАСПРЕДЕЛЕНИЯ
Нормальный (гауссов) закон распределения случайных величин играет особую роль. Во-первых, это наиболее часто встречающийся на практике закон распределения. Данное обстоятельство объясня-
55