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

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

Раздел 10. Преобразования случайных величин

10.1Преобразование одной случайной величины

Мы будем рассматривать только преобразования случайных величин с абсолютно непрерывными распределениями. Пусть с. в. имеет функцию распределения F (x) и плотность распределения f (x). Построим с помощью функции g : R ! R случайную величину = g( ). Требуется найти функцию распределения и, если существует, плотность распределения .

Замечание 18. Плотность распределения случайной величины = g( ) существует далеко не при любых функциях g. Так, если функция g кусочно-постоянна, то с. в. имеет дискретное распределение, и плотность ее распределения не существует.

Упражнение 16. Привести пример случайной величины с абсолютно непрерывным распределением и непрерывной функции g такой, что = g( ) имеет

а) дискретное распределение; б) невырожденное дискретное распределение.

Плотность распределения g( ) заведомо существует, если, например, функция g монотонна («строго монотонна»). Вспомним, что означает «найти плотность распределения , если она существует».

По определению, если мы представим (для любого x) функцию распределения

x

R

в виде F (x) = h(y) dy, где подынтегральная функция h(y) неотрицательна,

1

то плотность распределения с. в. существует и в точности равна подынтегральной функции: f (x) = h(x).

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

Теорема 23. Пусть имеет функцию распределения F (x) и плотность распределения f (x), и постоянная a отлична от нуля. Тогда случайная величина = a + b имеет

плотность распределения f (x) =

1

 

 

f

 

 

x b

.

jaj

 

a

 

 

 

 

Доказательство. Пусть сначала a > 0.

 

 

a+b

 

 

 

 

 

 

 

 

 

a

 

 

 

 

a

 

 

 

 

(x b)=a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z

 

F

(x) = F

 

 

(x) = P(a + b < x) = P <

x b

= F

 

x b

=

1

f (t) dt =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Замена: t =

y b

,

dt =

dy

 

 

 

3

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y =

 

 

 

Z

a

 

 

 

a

 

 

 

2 t = x b

y = x, t =

 

 

 

 

 

 

 

 

=

4

 

 

 

 

 

a

 

 

a

 

 

5

=

 

1

 

 

f

 

 

y

b

 

 

dy;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7!

 

 

1 7!

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

то есть при a > 0 случайная величина = a + b имеет абсолютно непрерывное распределение с плотностью

 

a

 

 

a

jaj

 

a

f (x) =

1

 

f

 

x b

 

=

1

 

f

x b

:

 

 

 

 

 

 

50

Пусть теперь a < 0.

 

 

 

 

a+b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z

 

 

F

(x) = F

 

 

(x) = P(a + b < x) = P

>

x b

=

(x b)=a

f

(t) dt =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

Замена: t =

y b

,

dt =

dy

,

t =

x b

y = x,

t =

 

y =

1 i

=

a

 

 

 

h

 

 

 

 

 

 

 

a

 

 

a

7!

 

 

 

 

1 7!

 

 

 

 

 

 

 

 

1

 

 

 

 

a

 

x

jaj

 

 

a

 

 

 

 

 

 

 

 

 

Zx a

 

 

Z

 

 

 

 

 

 

 

=

 

1

 

 

f

 

 

y b

dy =

1

 

 

f

y b

 

dy;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

то есть при a < 0 случайная величина = a + b имеет абсолютно непрерывное распре-

 

 

jaj

 

 

a

деление с плотностью f

(x) =

1

 

f

 

x b

.

 

 

 

 

Для произвольной монотонной функции g (то есть такой, что для любых x1 < x2 ëèáî g(x1) < g(x2) (монотонно возрастающая функция), либо g(x1) > g(x2) (монотонно убывающая функция)) справедливо аналогичное теореме 23 утверждение.

Теорема 24. Пусть имеет функцию распределения F (x) и плотность распределения f (x), и функция g : R ! R монотонна. Тогда случайная величина = g( ) имеет плотность распределения

f (x) =

 

1

 

f g 1(x) :

jg0

g 1

(x) j

 

 

 

 

 

Здесь g 1 — функция, обратная к g, è

 

1

 

= g 1(x) 0 — производная функции

g0

g 1(x)

g 1(x).

Упражнение 17. Доказать теорему 24.

Следующие утверждения сразу следуют из теоремы 23. Первое из них мы уже доказывали непосредственно. Доказать остальные.

Следствие 8. Åñëè = N0;1, òî = + a = Na;2 .

 

 

 

 

 

1

 

x

 

a

1

 

 

(x a)2

 

 

 

 

 

 

 

 

 

p2

 

 

Доказательство. Действительно, f

(x) =

 

f

 

 

 

=

 

 

e

2 2 .

 

 

 

 

 

 

 

Следствие 9.

Åñëè = N

 

2 , òî =

a

= N

 

.

 

 

 

 

 

 

 

 

a;

 

 

 

 

0;1

 

 

 

 

 

 

 

Следствие 10.

Åñëè = E , òî = = E1.

 

 

 

 

 

 

 

 

10.2Функции от двух случайных величин

Пусть 1; 2 — случайные величины с плотностью совместного распределения f 1;2 (x1; x2), и задана функция g : R2 ! R. Требуется найти функцию (а если существует, то и плотность) распределения случайной величины = g( 1; 2).

51

Пользуясь тем, что вероятность случайному вектору попасть в область можно вы- числить как объем под графиком плотности распределения вектора над этой областью, сформулируем утверждение.

Теорема 25. Пусть x 2 R, и область Dx R2 состоит из точек (x1; x2) таких, что g(x1; x2) < x. Тогда случайная величина = g( 1; 2) имеет функцию распределения

F (x) = P g( 1

; 2) < x = P ( 1; 2) 2 Dx

 

= ZZ f 1;2 (x1; x2) dx2 dx1:

 

 

Dx

Всюду далее в этой главе предполагается, что случайные величины 1 è 2 независимы, то есть f 1;2 (x1; x2) = f 1 (x1)f 2 (x2).

Следствие 11 (Формула свертки). Åñëè ñ. â. 1 è 2 независимы и имеют абсолютно непрерывное распределение с плотностями f 1 (x1) è f 2 (x2), то плотность распределения суммы 1 + 2 равна «свертке» плотностей f 1 è f 2 :

 

1

 

1

 

 

f 1+ 2 (t) =

Z

f 1 (u) f 2 (t u) du =

Z

f 2 (u) f 1 (t u) du

(13)

 

1

 

1

 

 

Доказательство. Воспользуемся утверждением теоремы 25 для функции g(x1; x2) = x1 +x2. Итегрирование по области Dx = f(x1; x2) : x1 +x2 < xg можно заменить последовательным вычислением двух интегралов: наружного – по переменной x1, меняющейся от 1 до 1, и внутреннего – по переменной x2, которая при каждом x1 должна быть меньше, чем x x1. Òî åñòü Dx = f(x1; x2) : x1 2 (1; 1); x2 2 (1; x x1)g. Поэтому

F 1+ 2 (x) = ZZ f 1;2 (x1; x2) dx2 dx1

=

ZZ

f 1 (x1) f 2 (x2) dx2 dx1 =

 

 

независ-ть

 

 

1 0 x x1f 1 (x1) f 2 (x2) dx21 dx1

Dx

 

Dx

=

 

 

 

 

Z

@

Z

A

 

 

 

 

1

1

Сделаем в последнем интеграле замену переменной x2 íà t òàê: x2 = t x1. Ïðè ýòîì x2 2 (1; x x1) перейдет в t 2 (1; x), dx2 = dt. В полученном интеграле меняем, наконец, порядок интегрирования:

F 1+ 2 (x) =

1 0

x f 1 (x1) f 2 (t

 

x1) dt1 dx1

=

x

0

1 f 1 (x1) f 2 (t

 

 

x1) dx1

1dt:

 

Z

Z

 

 

 

A

 

Z

 

 

Z

 

 

 

 

A

 

1 @ 1

 

 

 

1 @ 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|

 

x

{z

 

 

 

 

 

}

 

 

 

 

 

 

 

 

 

 

 

f 1

+ 2 (t)

 

 

 

 

 

 

 

 

1

 

 

 

1

1

(x)

R1

 

(t) dt, ãäå

Итак, мы представили функцию распределения

F + 2

â âèäå

f 1+ 2

 

f 1+ 2 (t) =

Z

f 1 (x1) f 2 (t x1) dx1 =

Z

f 1 (u) f 2 (t u) du:

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

Второе равенство получается либо из первого заменой переменных, либо если всюду в доказательстве поменять местами индексы 1 è 2.

Следствие 11 не только предлагает формулу для вычисления плотности распределения суммы, но и утверждает (заметьте!), что сумма двух независимых случайных величин с абсолютно непрерывными распределениями также имеет абсолютно непрерывное распределение. Для тех, кто уже ничему не удивляется, упражнение: привести пример двух случайных величин с абсолютно непрерывными распределениями, таких что их сумма имеет вырожденное распределение.

52

Если даже одна из двух независимых случайных величин имеет дискретное, а вторая

– абсолютно непрерывное распределение, то их сумма тоже имеет абсолютно непрерывное распределение, как показывает следующее упражнение.

Упражнение 18. Пусть с. в. имеет таблицу распределения P( = ai) = pi, с. в. имеет плотность распределения f (x), и эти величины независимы. Доказать, что +

P

имеет плотность распределения f + (x) = pif (x ai).

i

10.3Примеры использования формулы свертки

Пример 27. Пусть независимые случайные величины и имеют стандартное нормальное распределение. Докажем, что их сумма имеет нормальное распределение с параметрами 0 и 2.

Доказательство. По формуле свертки, плотность суммы равна

 

 

1

 

 

1

 

 

1

 

2

 

x2

 

f

(x) =

Z

1

e 21 u2 e 21 (x u)2 du =

Z

1

e 21 (2u2+x2 2xu)du =

Z

1

e u

+

2

xu du:

2

2

2

 

+

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

1

 

 

 

 

 

Выделим полный квадрат по u в показателе экспоненты: u2 + x22 Тогда

f + (x) = e

x2

1

2 e (u

2 )

2

du = e

x2

1

2 e v

dv =

2p e

x2 1

4

Z

 

4

Z

Z

 

 

1

x

 

 

 

1

2

 

1

 

 

4

xu = u x2 2 + x42 .

1 2

 

1

 

 

 

x2

 

 

 

4

p

 

e v

dv =

2p

 

e

 

:

 

 

 

1 1 1

Последнее равенство верно поскольку под интегралом стоит плотность нормального

1

распределения с параметрами 0 и p2 , так что интеграл по всей прямой равен 1. Итак, мы получили, что плотность суммы есть плотность нормального распределения с параметрами 0 и 2.

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

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

Лемма 4.

Пусть

случайные

величины

=

è

=

независимы.

Тогда

+ = + .

 

 

 

 

 

 

 

 

Лемма 5.

Пусть

случайные

величины

= Bn;p

è

= Bm;p

независимы.

Тогда

+ = Bn+m;p.

 

 

 

 

 

 

 

Лемма 6.

Пусть случайные величины = Na1; 12

è = Na2; 22

независимы.

Тогда

+ = Na1+a2; 12+ 22 .

 

 

 

 

 

 

 

Лемму 6 мы докажем позднее, используя аппарат характеристических функций, хотя при некотором терпении можно попробовать доказать ее с помощью формулы свертки, как в примере 27.

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

53

Определение 38. Случайная величина имеет гамма-распределение ; с параметрами > 0; > 0, если она имеет плотность распределения

 

 

 

f (x) =

(c x 1e x;

ïðè x > 0;

 

 

 

 

 

0;

ïðè x 6 0;

 

 

 

 

1

1

 

 

 

R

R

где постоянная

c вычисляется из условия

f (x) dx = c x 1e x dx = 1, òî åñòü

 

 

 

 

 

1

0

c =

 

.

 

 

 

 

( )

 

 

 

 

 

 

 

1

 

 

 

 

 

 

Здесь ( ) = R0

x 1e x dx = ( 1) ( 1) — гамма-функция Эйлера;

при k целых (k) = (k 1)! и (1) = 1.

Заметим, что показательное распределение E есть гамма-распределение ;1.

Лемма 7. Пусть независимые случайные величины 1; : : : ; n имеют показательное распределение E = ;1. Тогда 1 + + n = ;n.

Доказательство. Докажем утверждение по индукции. При n = 1 оно верно в силу равенства E = ;1. Пусть утверждение леммы справедливо для n = k 1. Докажем, что оно верно и для n = k. По предположению индукции Sk 1 = 1 + + k 1 = ;k 1, то есть имеет плотность распределения

 

 

fSk 1 (x) =

80; k 1

k

2

 

x

 

ïðè x 6 0;

 

 

 

<

 

x

e

 

;

ïðè x > 0:

 

 

 

(k 2)!

1

 

 

 

:

 

 

k

 

 

 

k равна

Тогда по формуле свертки плотность суммы S

 

=

+

 

+

 

 

 

1

 

 

 

 

1

k

1

 

 

fSk (x) =

Z

fSk 1 (u)f k (x u) du =

Z

 

 

uk 2e uf k (x u) du:

(k 2)!

 

1

 

 

 

 

0

 

 

 

 

 

 

Òàê êàê f k (x u) = 0 при x u < 0, то есть при u > x, то плотность под интегралом отлична от нуля, если переменная интегрирования изменяется в пределах 0 6 u 6 x. Поэтому

x

(k 2)!

x

(k 2)! uk 2 du =

(k 1)! xk 2e x:

fSk (x) = Z0

uk 2e u e (x u) du = e x Z0

 

k 1

 

k

k

Òî åñòü Sk = ;k, что и требовалось доказать.

54