Раздел 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