Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Литература / книга.pdf
Скачиваний:
147
Добавлен:
17.04.2013
Размер:
1.36 Mб
Скачать

Глава 2 Случайные величины

2.1. Основные понятия, связанные со случайной величиной

Пусть опыт полностью формализован, т.е. Э {, F, P}. Рассмотрим какую-либо функцию X(ω), ω Ω, отображающую в R (R). При некоторых дополнительных условиях X(ω) называется случайной величиной.

Пример 1. Монета подбрасывается три раза. Определим случайную величину Х - число выпавших гербов. Построим явно эту функцию.

Записывая исходы эксперимента в виде слов из нулей и единиц, получим:

= {111, 011, 101, 110, 100, 010, 001, 000}.

Всего исходов 23 = 8. Закон соответствия ω →X(ω) описывается следующей таблицей:

ω

111

011

101

110

100

010

001

000

X(ω)

3

2

2

2

1

1

1

0

Обозначим ЕX = {0,1,2,3} множество возможных значений случай-

ной величины X(ω). Множество EX часто называют спектром случайной величины Х.

Пример 2. Наудачу выбирается точка в единичном квадрате. Определим случайную величину Z(ω) - расстояние от точки до начала координат. Выразить явно Z через исходы.

= {(x,y)}| 0 x 1, 0 y 1}; Z(ω) = x2 + y2 , x [0,1], y [0,1];

Ez = [ 0, 2 ] - спектр Z(ω); Z - случайная величина непрерывного типа (в примере 1 X - случайная величина дискретного типа).

Определение. Случайной величиной называется числовая функция X(ω), определенная на множестве , такая, что для любого х из множества действительных чисел R (x R) множество тех ω, для которых выполняется условие X(ω) < x, принадлежит алгебре событий F для

28

данного опыта. (Краткая запись: {ω | X(ω) < x, x R} F; в дальнейшем будем писать сокращенно: {X < x}.)

Замечание 1. Дополнительное условие, фигурирующее в определениислучайнойвеличины, обеспечиваетизмеримость события{X(ω) < x дляx R}. При этом измеримость других интервальных событий будет иметь место автоматически, поскольку F - алгебра.

Действительно, по определению алгебры {X x} = {X < x} F.

Событие {x1 X < x2} = {X x1} {X < x2} F, поскольку представляет произведение и поэтому также является элементом алгебры F.

Замечание 2. Рассматриваемая в примере 1 функция X(ω) является случайной величиной, так как - конечно, а множество всех подмножеств конечного множества образует алгебру. Рассматриваемая в примере 2 функция Z(ω) - случайная величина, если в качестве алгебры событий F рассматривается совокупность всех квадрируемых (по Лебегу) подмножеств из .

Действительно, множество {Z(ω) <z, z R+} представляет собой круг радиусом z > 0, поэтому является квадрируемым множеством.

2.2.Функция распределения случайной величины

иее свойства

Определение 1. Функцией распределения случайной величины X

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

 

FX(x) = P{X < x}.

(2.1)

Свойства FX(x):

 

1.

FX(x) определена для любого х R, причем 0 FX(x) 1;

 

2.

FX(– ) = 0, FX(+ ) = 1, так как события {X < – }

= ,

{X < +

} = ;

 

3.

FX (x) - неубывающая функция переменной х;

 

(Действительно, пусть х1 < x2 {X < x1} {X < x2} {по следст-

вию из аксиом} FX(x1) FX(x2).)

 

4.

FX(x) непрерывна слева для х R;

 

5.

вероятность попадания на полуинтервал

 

 

P{x1 X < x2} = FX(x2) – FX(x1).

(2.2)

Ясно, что {X < x2} = {X < x1} + {x1 X x2} {по аксиоме 3 в

силу несовместимости слагаемых} P{x1 X x2} = FX(x2) – FX(x1).

29

Эти свойства универсальны, т.е. не зависят от того, дискретна или непрерывна функция X(ω).

Определение 2. Случайная величина X называется случайной величиной дискретного типа (СВДТ), если множество ее значений EX конечно или счетно.

Определение 3. Случайная величина Х называется случайной величиной непрерывного типа (СВНТ), если множество ее значений EX образует интервалы на действительной оси.

2.3. Закон распределения случайной величины дискретного типа

Определение 1. Законом распределения СВДТ называется таблица,

состоящая из двух строк: в первой строке перечисляются все возможные значения случайной величины, а во второй строке указываются соответствующие вероятности их реализаций. При этом должно выпол-

няться обязательное условие нормировки:

P{X = xk } =1 , хk ЕX. k

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

 

Х

0

1

2

 

 

 

 

3

 

Pk =1

 

P

1/8

3/8

3/8

 

 

 

1/8

 

 

 

 

 

 

k

 

 

 

 

1

2

1

 

 

 

 

 

Например, P{X = 0} = P{(000)} =

 

=

 

;

 

 

 

 

 

8

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

P{X = 1} = P{(100)} + P{(010)} + P{(001)} =

 

3

и т.д.

8

 

 

 

 

 

 

 

 

 

 

 

 

Пример 4. Пусть Э - последовательность испытаний по схеме Бернулли. Определить случайную величину Х - число успехов в n опытах посхемеБернулли. Описатьзаконраспределения случайнойвеличиныХ.

Очевидно, что ЕX = {0,1,2,…,n}. Ясно, что событие {X = k} = Bn,k

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

 

P{X = k} = Cnk pk qnk ,

(2.3)

30

 

где p - вероятность успеха в одном опыте; q = (1 – p); k = 0,1,…,n. Определение 2. Распределение, описываемое формулой (2.3),

называется биномиальным распределением с параметрами n и p (кратко:

Х B(n,p)).

Пример 5. Убедиться, что распределение, полученное в примере 3,

является биномиальным с параметрами n = 3, p = 12 .

Например, как следует из таблицы распределения, P{X = 2} = 83 . С

 

1

), тоP{X = 2} = C2

 

1

2

 

1

1

3

 

другойстороны, еслиХ B(3,

 

 

 

 

 

 

 

=

 

ит.д.

 

2

3

 

2

 

 

2

 

 

8

 

Пример 6. Описан закон распределения СВДТ, т.е. задана таблица. Построить функцию распределения FX(x).

Согласно определению (2.1)

FХ(x) = P{X < x} = P{X = xk } .

(2.4)

xkk<x

 

Формула (2.4) иллюстрируется рис.2.1: складываются вероятности тех дискретов из спектра, которые расположены левее пробной точки x.

Формула (2.4) свидетельствует о том, что FХ(x) является функцией накопленных вероятностей.

График FХ(x) для Х из примера 3 приведен на рис.2.2.

2.4. Основные дискретные распределения и их

x1 x2 x3

xn

Пробная точка

Рис.2.1.

FX(x)

 

 

 

 

1

 

 

 

 

5

 

 

 

 

8

 

 

 

 

1

 

 

 

 

2

 

 

 

 

1

 

 

 

 

8

 

 

 

 

0

1

2

3

x

 

 

Рис.2.2.

 

 

31

числовые характеристики

Ключевым понятием является "момент распределения".

Момент распределения

Начальный момент s-го порядка

Центральный момент s-го порядка

αs = xks pk ;

µs = (xk −α1)s pk ;

k

k

pk = P{X = xk}, xk Ex;

pk = p{X = xk}, xk Ex;

s = 0,1,…;

s = 0,1,…;

α0 = 1 (для распределения).

µ0 = 1 (для распределения).

Определение 1. mХ = α1 = x1k pk называется математическим

k

ожиданием. Это среднее значение случайной величины Х по распределению (аналогия с формулой центра тяжести в механике). Часто mХ

называют "центром распределения".

Определение 2.

DХ = µ2 = (xk mX )2 pk

(2.5)

k

 

называется дисперсией. Дисперсия - показатель степени разброса случайной величины.

Определение 3. Величина aХ = µ3 называется асимметрией, где

σ3x

σX = DX - среднеквадратическое отклонение ( σX - эталон разброса).

Если ax = 0, то распределение случайной величины X симметрично относительно математического ожидания.

Определение 4. eХ = µ4 – 3 называется эксцессом. Эксцесс -

σ4X

характеристика островершинности распределения.

32

Преобразуем формулу (2.5):

DX = (xk mX )2 pk = (xk2 2mX xk + m2X ) pk = xk2 pk

k

 

k

 

 

 

k

 

 

 

 

xk2 pk

 

 

 

 

 

 

 

 

= α2

 

 

 

 

 

 

k

 

 

 

 

 

 

2

 

 

 

xk pk = mX

 

= α2

2

2mX xk pk + mX pk

=

 

 

 

mX ;

k

k

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pk =1(нормировка)

 

 

 

 

 

 

k

 

 

 

 

 

 

 

DX = α2 m2X .

 

 

(2.6)

1.Равномерное распределение (дискретное).

ХR(n), n - число дискретов. Таблица распределения в этом случае имеет вид:

Х

х1

х2

хn

 

 

 

 

 

 

 

 

P

 

1

 

 

1

 

 

1

 

 

n

 

n

 

n

 

 

 

 

 

mX = 1 n xk - среднее арифметическое. n k =1

 

 

1

n

1

 

n

 

2

 

1

 

n

 

1

 

n

 

2

 

 

 

 

 

 

 

 

 

DX

=

xk 2

 

xk

 

=

xk 2

 

xk

 

.

 

2

 

n

 

 

 

 

n k =1

n

k =1

 

 

 

k =1

 

n k =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Пусть Х B(n,p). Это значит, что спектр возможных значений

EX = {0,1,2,…,n} и

P{X = k} = Pn,k(p) = Cnk pk qnk , k EX.

33

Найдем mX и DX:

n

mX = k Cnk pk qnk .

k =0

Эту сумму можно свернуть двумя способами: 1) с использованием бинома Ньютона.

Формула бинома Ньютона:

n

(q + p)n = Cnk pk qnk , q иp - произвольныедействительныечисла. (2.7)

k =0

Если q + p = 1, то бином Ньютона дает нормировку и возможность использовать различные операции.

Чтобы вычислить mX, надо (2.7) продифференцировать по p и

умножить на p:

 

p[(q + p)n]' = mX = p(n(q + p)n–1) (q + p) = 1 = np

 

mX = np.

(2.8)

Для вычисления дисперсии по формуле (2.6) находим α2:

n

α2 = k 2 Cnk pk qnk . k =0

Продифференцируем (2.7) дважды по р:

n

[(q + p)n]' = [n(q + p)n–1] = k Cnk pk 1 qnk ;

k =0

p[np(q + p)n–1]'p = α2 = pn((q + p)n–1 + (n – 1)p(q + p)n–2) =

= (таккакq + p = 1) = np(1 + np p) = np(q + np) = {q = 1 – p} = n2p2 + npq;

DX = α2 m2X = n2p2 + npq n2p2 = npq; DX = npq;

2) c помощью производящей функции.

Производящая функция для биномиального распределения:

n

Ψ(x) = (q + px)n = Pn,k ( p)xk , где Pn,k - биномиальные вероятности.

k =0

mX = Ψ '(x)|x=1 = n(q + px)n–1 p = np(q + 1 – q)n–1 = np;

34

α2 = [xΨ '(x)]'x=1 = np + x np(n – 1)p(q + px)n–2 = np(1 + np p) = np(np + q); DX = (np)2 + npq – (np)2 = npq;

DX = npq.

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

Определение 5. Говорят, что Х подчиняется распределению Пуас-

сона с параметром λ > 0 (кратко пишут: X Pu(λ)), если множество

0

возможных значений EX = N = {0,1,2,…}, а соответствующие вероятности вычисляются по формуле:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pk = P{X

= k} =

λk

e

−λ

.

 

 

 

 

 

(2.9)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычислим mX и DX .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

e− λ

 

 

 

 

k

 

= e−λ eλ

 

 

 

 

 

 

Нормировка:

λ

= e− λ

λ

 

= 1;

 

 

 

 

 

 

 

 

 

 

 

 

 

k =0

k!

 

 

 

 

 

 

k =0

k!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λk

 

 

 

 

 

λk

 

 

 

 

λk

 

 

 

 

 

 

 

 

mX = k

e

− λ

= k

e− λ

=

 

 

e− λ = {m = k 1} =

 

 

 

 

 

(k 1)!

 

k =0

 

 

 

k!

 

 

 

k =1

 

k!

 

 

 

 

k =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

− λ = λ eλ e−λ = λ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= λ

λ

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m=0

m!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

2

 

 

λk

 

− λ

 

k

 

 

 

k

 

λ

 

 

 

 

(k 1)

+1

k

−λ

 

α2

 

k

 

 

 

 

 

e

 

 

=

 

 

 

 

 

λ e

 

=

 

 

 

 

 

 

 

λ e

 

=

 

 

 

k!

 

 

(k 1)!

 

(k 1)!

 

 

 

k =1

 

 

 

 

 

 

 

 

 

 

k =1

 

 

 

 

 

 

 

 

k =1

 

 

 

 

 

 

= (k 1) λk k =1 (k 1)!

λ k 1

+λ (k 1)!

k =1

 

1

 

 

(k 1)

 

 

 

e−λ

+

 

λk e−λ

= λ

λk 1 e−λ

+

(k 1)!

(k 1)!

 

k =1

 

k =1

 

 

 

 

 

 

 

 

 

 

e−λ = λ2 + λ ;

DX = m2X + α2 = λ2 + λλ2 = λ. Такимобразом, дляраспределенияПуассона

DХ = mХ = λ.

35

Распределение Пуассона описывает:

-число атомов, распавшихся в единицу времени (радиоактивный распад);

-число независимых вызовов на АТС в единицу времени и т.д.

 

 

Теорема 2.1. Пусть X ~

B(n,p). Положим λ =

np, т.е. p =

λ

, и

устремим

n

→ ∞ в формуле

 

 

 

 

 

 

 

 

 

 

λ =

 

 

 

n

 

(2.3). Если при этом

const,

P

 

λ

λk

e

−λ

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n,k n

 

n→∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ

 

 

k

 

λ

k

 

 

λ

nk

n!

 

 

 

k

(1

λ)n

 

 

 

 

 

 

Pn,k(

 

 

 

 

 

 

 

 

 

 

 

λ

 

 

 

 

n

 

 

 

 

 

 

 

 

 

) =

Cn

 

 

 

 

 

1

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

n

 

 

k!(n k)!

nk

 

 

 

λ

k

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1

n )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ)n

 

n(n 1)...(n k +1)

 

1

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

(1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n(n 1)...(n

k

+1)

 

 

 

 

k

 

 

 

 

 

 

 

λ

 

 

 

n

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

nk

 

 

 

 

 

(1

 

 

 

 

 

 

 

 

λ

 

k!

 

 

 

 

 

 

 

 

 

λ)k

 

(1

λ) 1, так как

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

n

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λk e−λ при n → ∞, т.е. получаем формулу (2.9). k!

Из теоремы следует, что при выполнении условий n велико, p мало, λ - фиксированная величина, биномиальные вероятности можно аппроксимировать пуассоновскими. Качество аппроксимации гарантируется оценкой [3]:

P

 

λ

λk

e

−λ

λ2

 

 

k!

 

0 , для любых k = 0,1,…,n.

n,k

n

 

 

 

 

n n→∞

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

Пример 7. Система коммутаций на АТС содержит n = 50000 независимо работающих элементов. Вероятность отказа каждого за время Т равна p = 2 10–5.

1. Найти среднее (по распределению) число отказавших за время Т элементов и наиболее вероятное число отказавших элементов.

36

2. Какова вероятность, что за время Т откажет не менее 2-х элементов?

Обозначим Х число отказавших за время Т элементов. Согласно условию Х ~ В(50000, p = 2 10–5). При этом np = λ = 5 104 2 10–5 = 1.

1.Согласно (2.8) mХ = np = 1 - среднее число отказавших элементов

математическое ожидание. Наиболее вероятное число отказавших элементов (мода dX) удовлетворяет неравенству:

np q dX np + p,

отсюда следует, что dX = 1.

2. Используем аппроксимацию Пуассона:

P{X 2} = 1 – P{X < 2} = 1 – p0 p1 = 1 – e–1 e–1 = 1 – 2 e–1 = 1 – 0,7358 = 0,2642.

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

Определение 6. Говорят, что Х распределена по геометрическому за-

конуспараметромр> 0 (запись: Х~ Гео(р)), еслиEX = N, P{X = xk} = qk–1 p.

Это распределение встречается в опытах до 1-го успеха по схеме Бернулли, где X - число проведенных опытов до 1-го успеха (включительно).

Нормировка:

 

 

1

 

qk 1 p = pqk 1 =

p

 

= 1.

1

q

k =1

k =1

 

 

 

 

 

 

Здесьиспользованаформуласуммычленовгеометрическойпрогрессии:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

qk =

 

 

.

 

(2.10)

 

 

 

 

 

 

 

 

 

 

 

 

1

q

 

 

 

 

 

 

 

 

 

 

 

 

k =0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычислим математическое ожидание:

 

 

 

 

 

 

 

k 1

 

d

k

 

 

 

 

 

 

 

 

'

 

mX = pkq

= p

q

=

{дифференцируем(2.10)} =

 

1

 

 

 

 

 

 

p

 

 

=

 

 

 

dq

 

 

 

 

k =1

 

 

 

 

 

k =0

 

 

 

 

 

 

 

 

1q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

p

 

=

 

1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1

q)2

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Дважды дифференцируя (2.10) по q и используя формулу (2.6), находим:

DX = q2 . p

37