Posit1nov (1)
.pdf
Свойства вероятностей pij : |
|
|
||||
|
|
P |
|
|
|
|
pv |
1. |
pij = pX,i, где pX,i = P (X = xi), |
i = 1, 2, . . . - фиксировано. |
|||
|
|
j |
|
|
|
|
|
|
P |
|
|
|
|
pv |
2. |
pij = pY,j , где pY,j = P (Y |
= yj ), |
j = 1, 2, . . . - фиксировано. |
||
|
|
i |
D - |
P P |
|
|
v |
|
|
|
|
||
pv |
3. Свойство нормировки |
pij = 1. |
|
|||
|
|
|
|
i j |
|
|
p 4. Для любого |
|
борелевского множества на плоскости справедливо соотно- |
||||
|
|
|
|
P |
|
|
шение P (X = (X, Y ) B) = |
pij . |
|
||||
i,j: (xi,yj ) B
Докажем, например, первое свойство. Используя свойства действий над событиями
|
|
|
|
= |
|
= |
|
= |
|
P |
|
|
|
|
P |
|
|
|
|
(см. замечание 2.5 ) будем иметь |
j {X = xi, Y = yj } = |
j ({X = xi} · {Y = yj }) = |
|||||||||||||||||
{ |
P |
|
{X = xi, Y = yj } |
i} · |
|
{ |
|
j несовместны.P |
|
j } |
|
|
|||||||
i} j |
{ |
Y |
|
y |
j } { |
X |
|
Ω = |
|
i} |
, так как |
j { |
Y = y |
= Ω. |
|||||
|
X = x |
|
|
|
|
x |
|
X = x |
|
|
|||||||||
|
События |
|
|
|
|
|
|
при разных |
|
|
Применим аксиому счет- |
||||||||
ной аддитивности |
P |
|
|
|
|
|
|
P |
|
|
|
|
P |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
pX,i = P (X = xi) = P ( j |
{X = xi, Y = yj }) = j |
P {X = xi, Y = yj } |
= j |
pij . |
|||||||||||||||
|
В дискретном случае функция распределеня случайного вектора (X, Y ), вычис- |
||||||||||||||||||
ляется по формуле F (x, y) = |
P |
P |
pij . |
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||
i: xi<x j: yj <y
Пример 5.2 В двух ящиках находится по 6 занумерованных шаров.
В первом ящике: 1 шар - с номером 1; 2 шара - с номером 2; 3 шара - с номером 3. Во втором ящике: 2 - с номером 1; 3 - с номером 2; 1 - с номером 3.
Из каждого ящика случайным образом вынули по шару. Пусть X - номер шара, вынутого из первого ящика и Y - номер шара вынутого из второго ящика. Составить таблицу распределения случайного вектора (X, Y ).
Решение. События {X = i} и {Y = j} независимы при любых i, j в силу незавсимости эксперементов (разные ящики). Следовательно, P {X = 1, Y = 1} = P {X = 1}P {Y = 1} = 1/6 · 2/6 = 1/18. Поступая аналогично для вычисления
других вероятностей, получим таблицу.
X \ Y |
1 |
2 |
3 |
1 |
1/18 |
1/12 |
1/36 |
|
|
|
|
2 |
1/9 |
1/6 |
1/18 |
31/6 1/4 1/12
Рассмотрим теперь случайный вектор X = (X, Y ), где X и Y - непрерывные
случайные величины.
Определение 5.4 Пусть существует такая функция f (x, y) ≥ 0, что для всех вещественных x, y справедливо представление
x y |
|
F (x, y) = Z−∞ Z−∞ f (u, v) du dv, |
(5.1) |
43
тогда случайный вектор X = (X, Y ) называется непрерывным, а функция f (x, y) называется плотностью распределения случайного вектора X = (X, Y ) или совместной плотностью распределения случайных величин X и Y .
Основные свойства плотности распределения случайного вектора:
∞ ∞
R R
fv 1. f (x, y)dxdy = F (+∞, +∞) = 1 - свойство нормировки.
−∞ −∞
fv 2. Положим f (x, y) = ∂2F (x, y) , в тех точках, где эта производная существует.
R R ∂x ∂y
fv 3. P ((x, y) D) = f (x, y)dxdy, для любого D B(R2), где B(R2) - боре-
D
левская σ-алгебра на плоскости. Заметим, что это свойство эвивалентно соотноше-
нию (5.1).
|
|
fv 4. |
∞ |
f (x, y)dy = fX (x), где fX (x) - плотность распределения случайной ве- |
|||||||||
|
|
|
|||||||||||
|
|
|
|
|
−∞ |
|
|
F x |
F |
|
x, |
|
|
личины |
RX |
. Это свойство справедливо в силу того, что |
|
+∞) = |
|||||||||
R |
x |
|
|
R |
∞ |
|
X ( ) = |
|
( |
|
|||
−∞ |
( |
f (u, v) dv) du (см. определение 4.5) и означает, что компоненты непрерыв- |
|||||||||||
|
|
−∞ |
|
|
|
|
|
|
|
|
|||
ного случайного вектора являются непрерывными случайными величинами. |
|||||||||||||
|
|
|
|
|
|
|
∞ |
|
|
|
|
|
|
|
|
fv 5. Аналогично, |
f (x, y)dx = fY (y), где fY (y) плотность распределения слу- |
||||||||||
|
|
|
|
|
|
|
−∞ |
|
|
|
|
|
|
чайной величины Y . R |
|
|
|
|
|
|
|||||||
Доказательство свойства E3 математического ожидания.
Проведем доказательство для непрерывного случайного вектора X = (X, Y ), в
дискретном случае доказательство аналогично. В силу абсолютной сходимости интегралов (по определению математического ожидания) имеем
|
E(X + Y ) = Z∞ Z∞(x + y)f (x, y)dxdy = Z∞ xdx Z∞ f (x, y)dy + |
(5.2) |
||
|
−∞ −∞ |
|
−∞ −∞ |
|
+ |
Z∞ ydy Z∞ f (x, y)dx = |
Z∞ xfX (x)dx + Z∞ yfY (y)dy = E(X) + E(Y ). |
|
|
|
−∞ −∞ |
−∞ |
−∞ |
|
Замечание 5.1 Заметим, что в только что приведенном доказательстве нельзя ограничиться требованием непрерывности случайных величин X и Y , так как из этого требования не следует непрерывность случайного вектора X = (X, Y ). Например, при Y = X значения случайного вектора X = (X, X) сосредоточены на прямой y = x и, если X - непрерывная случайная величина, то его распределение будет сингулярным. То есть вероятность P (X = x, Y = y) = 0, для любой точки плоскости (x, y), множество значений случайного вектора V = {(x, y) = (X(ω), Y (ω)), ω Ω} несчетно, а мера множества V (площадь) равна нулю, следовательно, f (x, y) - плотность распределения, не существует.
С использованием интеграла Лебега-Стилтьеса по мере P = P (dω) приведен-
ное доказательство обобщается на общий случай.
Условные законы распределения вероятностей компонент случайного вектора.
44
Определение 5.5 Условным законом распределения случайной величины X, входящей в систему (X, Y ) называется её закон распределения, вычисленный при
условии, что другая случайная величина Y приняла определенное значение yj0 |
или |
||||||
y0. Обозначается {pi|j0 }, i = 1, 2, ...; F (x|y0); f (x|y0). |
|
|
|||||
По определению условной вероятности имеем при i = 1, 2, . . . |
|
||||||
pi|j0 = P (X = xi|Y = yj0 ) = |
P (X = xi, Y = yj0 ) |
= |
pij0 |
при P (Y = yj0 ) 6= 0, |
(5.3) |
||
P (Y = yj0 ) |
|
pY,j0 |
|
||||
Так как при P (Y = yj0 ) = 0 вероятность P ({X = xi, Y |
= yj0 } {Y = yj0 }) = 0 |
||||||
положим pi|j0 = 0 при P (Y = yj0 ) = 0,
Аналогично, можно показать, что плотности условного распределения имеют
вид |
|
|
|
|
|
|
|
|
f (x|y0) = |
f (x, y0) |
при fY (y0) 6= 0 |
или f (y|x0) = |
f (x0, y) |
при fX (x0) 6= 0, |
(5.4) |
||
|
|
|
|
|||||
fY (y0) |
fX (x0) |
|||||||
и положим f (x|y0) = 0, для y0 : fY (y0) = 0 или f (y|x0) = 0, для x0 : fX (x0) = 0.
Пример 5.3 Пусть задано распределение дискретного случайного вектора (X, Y ):
X |
|
Y |
|
|
-1 |
0 |
1 |
|
|
|
|
-1 |
0.5 |
0.1 |
0 |
1 |
0 |
0.2 |
0.2 |
Найти распределения вероятностей случайных величин X и Y , условное распределение X при условии Y = 1 и условное распределение Y при условии X = −1.
Решение. Используем для нахождения распределений компонент случайного вектора (X, Y ) свойства pv1 и pv 2 вероятностей совместного распределения pij ,
тогда
X
pX1 = P {X = −1} = p1j = P {X = −1, Y = −1} + P {X = −1, Y = 0}+
|
|
|
|
|
|
|
|
|
|
j |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
+P {X = −1, Y = 1} = 0, 5 + 0, 1 + 0 = 0, 6; |
|
|
|
|
|
|
|
|
||||||||||||||||||
аналогично pX2 = P {X = 1} = 0, 4. Такие же вычисления проводятся и для Y. |
|
|||||||||||||||||||||||||||||||
Таким образом, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
или |
|
|
xi |
|
|
-1 |
|
1 |
|
|
|
|
или |
|
yi |
|
-1 |
|
0 |
|
|
1 |
|
|
|
||||
|
|
|
|
|
pX,i |
|
|
0.6 |
|
0.4 |
|
|
|
|
|
pY,i |
0.5 |
|
0.3 |
|
0.2 |
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
Для нахождения условного распределения, используя формулу (5.3) имеем: |
|||||||||||||||||||||||||||||
P |
{ |
X = 1 |
Y = 1 |
} |
= |
|
P {X = −1, Y = 1} |
= 0, |
|
|
P |
{ |
X = 1 |
| |
Y = 1 |
} |
= |
|||||||||||||||
|
|
− | |
|
|
|
|
|
|
P |
{ |
Y = 1 |
} |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
P {X = 1, Y = 1} |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
= 1 или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
P {Y = 1} |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
xi |
|
|
−1 |
|
1 |
. |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
pi|3 = P (X = xi|Y = 1) |
|
1 |
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
||||||||||||
45
Заметим, что в сечении Y = 1 возникло вырожденное распределение, сосредоточенное в точке X = 1.
Аналогично, находится условное распределение случайной величины Y при X = −1
yj |
−1 |
0 |
1 |
(5.5) |
|
pj|1 = P (Y = yj |X = −1) |
1/6 |
0 |
|||
5/6 |
|
Пример 5.4 Система случайных величин (X, Y ) равномерно распределена в тре-
угольнике D с вершинами в точках (0; 0), (0; 1), (1; 0). Найти:
|
|
|
|
|
|
|
|
|
1) |
плотность распределения f (x, y) случайного |
y |
|
6 |
|
|
|
вектора X=(X, Y ); |
||||
|
|
|
|
|||||||
1 |
|
|
|
@y = −x + 1 |
2) |
плотности компонент fX (x) fY (y); |
||||
|
|
|
||||||||
|
|
|
@ |
|
|
|
|
|||
|
|
|
|
D@@ |
|
|
- |
|
3) |
плотности условных распределений f (y|x0). |
0 |
|
|
1 |
x |
||||||
|
||||||||||
Решение. 1) Плотность равномерного распределения постоянна в треугольни-
|
|
|
|
c |
|
const, |
(x, y) |
D, |
||
ке D и равна нулю вне треугольника: f (x, y) = (0,= |
|
|
вне D |
|
||||||
Используя свойство fv 1 совместной плотности f (x, y), найдем c: |
||||||||||
1 |
−x+1 |
1 |
(−x + 1) dx = c |
−2 + 1 = 2 c. |
||||||
1 = Z0 |
Z0 |
c dxdy = c Z0 |
||||||||
|
|
|
|
|
|
1 |
|
1 |
|
|
Таким образом, c = 2, следовательно, f (x, y) = 2 для всех (x, y) D.
2) Для вычисления плотностей компонент используем свойства fv 4, fv 5, тогда:
fX (x) = |
−0 |
2 dy = 2 · (−x + 1 + 0) = 2 − 2x, x (0, 1) ; |
||||
|
|
|
x+1 |
|
|
|
|
|
R |
0, |
x / (0, 1) |
||
fY (y) = |
|
−y+1 |
|
|
||
|
0 |
2 dx = 2 · (−y + 1 + 0) = 2 − 2y, y (0, 1) . |
||||
|
|
R |
|
|
|
|
|
|
|
|
0, |
y / (0, 1) |
|
|
|
|
|
|
|
|
3) Используя 5.4 для вычисления плотностей условных распределений, имеем: |
||||||
| 0 |
|
( |
|
0, |
(x, y) / D |
|
f (x y |
) = |
2/(2 − 2y0) = 1/(1 |
− y0), x (0, 1 − y0), y0 |
(0, 1) ; |
||
| 0 |
|
( |
|
|
|
|
|
|
0, |
(x, y) / D |
|
||
f (y x ) = |
2/(2 − 2x0) = 1/(1 |
− x0), y (0, 1 − x0), x0 |
(0, 1) . (5.6) |
|||
|
|
|
|
|
|
|
В этом примере распределение в любом сечении является равномерным.
46
5.2Независимые случайные величины.
Определение 5.6 Случайные величины X и Y называются независимыми, если для любых боревских множеств A и B выполнено соотношение
P (X A, Y B) = P (X A)P (Y B). |
(5.7) |
Замечание 5.2 Событие (X A, Y B) на языке теории множеств записывается как пересечение множеств: {ω : X(ω) A} T {ω : Y (ω) B)
Замечание 5.3 Свойство независимости на геометрическом языке площадей фигур звучит очень привычно: "площадь"прямоугольника равна произведению "длин"его сторон. Вероятностная мера не всегда обладает этим свойством, если же для вероятностной меры оно справедливо, то породившие эту меру случайные величины называются независимыми, то есть "поведение"одной из них на вторую "не влияет".
Теорема 5.1 Необходимым и достаточным условием независимости случайных величин X и Y является представимость функции их совместного распределения
в виде
FX(x, y) = FX (x)FY (y), при всех x, y R. |
(5.8) |
Доказательство. Необходимость. Если случайные величины X и Y независимы,
то представимость функции их совместного распределения в виде (5.8) следует из определения функции распределения F (x, y) = P (X < x, Y < y) = P (X (−∞, x), Y (−∞, y)) и того, что множества (−∞, x), (−∞, y) B - борелевской σ- алгебре на R1 при любых x, y.
Достаточность. Схема доказательства. Доказательство того, что из (5.8) следует независимость случайных величин X и Y основано на свойстве Fv 6 функции
распределения - из этого свойства следует, что соотношение (5.7) выполнено для любого прямоугольника: пусть D = {(x, y) R2 : a ≤ x < b, c ≤ y < d} -
прямоугольник, тогда
P (X [a, b), Y [c, d)) = P ((X, Y ) D) = F (b, d) − F (b, c) − F (a, d) + + F (a, c)) = FX (b)FY (d) − FX (b)FY (c) − FX (a)FY (d) + FX (a)FY (c) =
=FX (b)(FY (d) − FY (c)) − FX (a)(FY (d) − FY (c)) =
=(FX (b) − FX (a))(FY (d) − FY (c)) = P (X [a, b)) P (Y [c, d)).
Отсюда легко следует, что соотношение (5.7) справедливо для любого множества из алгебры, построенной на базе прямоугольников. Затем доказывается теорема апроксимации о вероятностной связи алгебры и порожденной ею σ-алгебры (то есть пересечения всех σ-алгебр, содержащих данную алгебру), откуда следует свойство (5.7) для всех боревских множеств A и B.
Теорема 5.2 Для дискретных случайных величин необходимым и достаточным условием независимости случайных величин X и Y является представимость
вероятности их совместного распределения в виде
pij = pX,ipY,j , где pX,i = P (X = xi), pY,i = P (Y = yj ) при всех i, j. |
(5.9) |
47
Доказательство. Необходимость. Как и в теореме 5.1 из условия (5.9) независимость случайных величин X и Y следует в силу того, что pij = P (X = xi, Y = yj ), а одноточечные множества {xi} и {yj } являются борелевскими.
Достаточность следует из свойства вероятностей pv 4: для любых борелевских A и B имеем
P (X A, Y B) = P ((X, Y ) {(x, y) : x A, y B}) =
= P ((X, Y ) A × B) = |
X × |
|
iX × |
|
|
|
pij = |
pX,i pY,j = |
|
= |
i,j: (xi,yj ) A |
B |
i,j: (x ,yj ) A |
B |
i: xi A pX,i |
j: yj B pY,j = P (X A)P (Y B). |
|||
|
X |
|
X |
|
Теорема 5.3 Для непрерывных случайных величин справедливы два утверждения.
A)Из независимости непрерывных случайных величин X и Y следует непрерывность случайного вектора X= (X, Y ).
B)Необходимым и достаточным условием независимости непрерывных случайных величин X и Y является представимость плотности их совместного
распределения в виде
f (x, y) = fX (x)fY (y), при всех x, y R. |
(5.10) |
Доказательство пункта А и необходимости пункта B. Пусть X и Y независимые непрерывные случайные величины. Из теоремы 5.1 имеем F (x, y) = FX (x)FY (y), откуда по теореме Фубини для любых вещественных x, y
F (x, y) = FX (x)FY (y) = |
Z∞ fX (u)du |
Z∞ fY (v)dv = |
Z∞ Z∞ fX (u)fX (v)dudv |
|
x |
y |
x y |
По определению 5.4 , при f (x, y) = fX (x)fY (y), отсюда следует непрерывность случайного вектора X= (X, Y ) и необходимость пункта B.
Доказательство достаточности пункта B. Пусть справедливо равенство (5.10), отсюда следует непрерывность случайного вектора X= (X, Y ). Из (5.1), используя
теорему Фубини и определение 4.5 получим
F (x, y) = |
Z∞ Z∞ f (u, v)dudv = |
Z∞ Z∞ fX (u)fX (v)dudv |
(5.11) |
|
|
x |
y |
x y |
|
= |
Z x fX (u)du Z y fY (v)dv = FX (x)FY (y). |
|
||
|
∞ |
∞ |
|
|
По теореме 5.1 из (5.11) следует независимость случайных величин X и Y .
Утверждение 5.1 (Свойство E5 математического ожидания) Пусть случайные величины X и Y независимы и существуют математические ожидания
E(X), E(Y ), тогда E(XY ) = E(X)E(Y ).
48
Доказательство. Докажем это утверждение для дискретных случайных величин, в случае непрерывных случайных величин доказательство проводится аналогично с заменой сумм на интегралы. Используя свойства абсолютно сходящихся рядов
(теорему Фубини для интегралов) будем иметь: |
|
|||||
|
X X |
|
X X |
|||
E(XY ) = |
|
|
xiyj pij |
= |
|
xiyj pX,ipY,j = |
= |
|
i |
j |
j |
i |
j |
i |
xipX,i |
yj pY,j = E(X)E(Y ). |
||||
|
|
X |
X |
|
||
Из этого утверждения и свойства D4 дисперсий вытекает
Следствие 5.1 Пусть случайные величины X и Y независимы и существуют дисперсии D(X), D(Y ), тогда D(X + Y ) = D(X) + D(Y ).
Пример 5.5 Пусть задана плотность распределения случайного вектора:
(
f (x, y) = |
Ae−x−2y , x ≥ 0 и y ≥ 0, |
|
|
0, |
x < 0 или y < 0. |
Найти A,
∞
1). A R e−x
F (x, y), fX (x), fY (y).
∞
dx R e−2y dy = A/2 = 1 - свойство нормировки, откуда A = 2.
00
и |
|
x |
y |
|
|
|
|
|
|
x |
|
y |
|
|
|
|
|
|
|
|
R |
R |
|
|
|
|
|
|
R |
|
R |
|
|
|
|
|
|
|
|
2). F (x, y) = |
|
f (x, y)dxdy = 2 e−udu |
e−2v dv = (1 − e−x)(1 − e−2y ), при x ≥ 0 |
||||||||||||||||
|
−∞ −∞ |
|
|
|
|
|
|
0 |
|
0 |
|
|
|
|
|
|
|
||
|
y ≥ 0; F (x, y) = 0 при x < 0 или y < 0. |
|
|
|
|
|
|
|
|
||||||||||
|
∞ |
|
|
|
|
, y |
∞ |
0 |
|
Y ( ) = 0 |
|
≤ 0 |
|
( |
|
|
) = |
||
Аналогично,R fY (y) = 2e− |
|
|
R |
|
|
|
|
|
|||||||||||
3). fX (x) = 0 |
2e−x−2y dy = 2e−x |
0 |
e−2y dy = e−x, при |
x > 0 и fX (x) = 0, при x |
≤ 0. |
||||||||||||||
|
|
|
|
|
2y |
|
> |
|
и f |
y |
|
, y |
|
. Следовательно, f x, y |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
fX (x)fY (y) и случайные величины X и Y независимы. |
|
|
|
|
|
||||||||||||||
|
Можно |
найти |
функции |
распределения |
по |
определению |
FX (x) |
|
= |
||||||||||
x ∞ |
|
|
x |
|
∞ |
y ≥ 0 и FY (y) = 0, y < 0. С другой стороны, в силу |
|||||||||||||
Аналогично,R R |
FY (y) = 1R− e−R2y , |
||||||||||||||||||
|
f (u, y)dudy = |
du |
|
2 e−u−2y dy |
= |
1 − e−x |
, x |
≥ 0 FX (x) = |
0, x |
< 0. |
|||||||||
−∞ −∞ |
|
|
0 |
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
||
независимости F (x, y) = FX (x)FY (y), и этот результат следует из пункта 2.
5.3Числовые характеристики случайного вектора.
Математические ожидания компонент случайного вектора. Начальные и центральные моменты компонент случайного вектора определяются как и ранее - в подразделе 4.2. Совокупность математических ожиданий mx = E(X) и my = E(Y ) представляет собой характеристику положения системы: (mx, my ) - координаты
центра рассеивания системы. Можно привести формулы для вычисления математических ожиданий и дисперсий компонент случайного вектора X. Для дискретного случайного вектора:
X X |
X X |
|
E(X) = |
xipij ; E(Y ) = |
yipij |
i j |
i |
j |
49
X X |
X X |
||
D(X) = |
(xi − E(X))2pij ; D(Y ) = |
(yi − E(Y ))2pij , |
|
i |
j |
i |
j |
Для непрерывного случайного вектора: |
|
|
|
E(X) = Z∞ Z∞ xf (x, y)dxdy; |
E(Y ) = Z∞ Z∞ yf (x, y)dxdy |
||
|
−∞ −∞ |
−∞ −∞ |
|
D(X) = Z∞ Z∞(x − E(X))2f (x, y)dxdy; |
D(Y ) = Z∞ Z∞(y − E(Y ))2f (x, y)dxdy. |
||
−∞ −∞ |
|
−∞ −∞ |
|
Дисперсии D(X), D(Y ) характерезуют рассеивание случайной точки (X, Y ) соответственно вдоль осей OX и OY
Ковариация, коэффициент корреляции системы случайных величин.
Для многомерного случайного вектора кроме начальных и центральных моментов рассматривают смешанные начальные и центральные моменты.
Определение 5.7 Начальными смешанными моментами порядка k случайного вектора X=(X, Y ) называются числа E(XiY j ), i, j = 1, 2..., i+j = k. Центральными смешанными моментами порядка k случайного вектора X=(X, Y ) называются числа E((X − E(X))i(Y − E(Y ))j ), i, j = 1, 2..., i + j = k.
Рассмотрим более подробно второй центральный смешанный момент.
Определение 5.8 Ковариацией случайных величин X и Y называется матема-
тическое ожидание произведения отклонений этих величин от своих математических ожиданий, обозначается K(X, Y ), cov(X, Y ):
− − ˚˚ ˚ −
K(X, Y ) = E ((X E(X))(Y E(Y ))) или K(X, Y ) = E(XY), где X= X E(X),
˚ −
Y= Y E(Y ) - центрированные случайные величины. То есть, ковариация есть ма-
тематическое ожидание произведения центрированных величин (второй центральный смешанный момент). Ковариация является качественной характеристикой зависимости случайных величин. Если распределение дискретное, то
X X
K(X, Y ) = (xi − E(X))(yj − E(Y ))pij .
ij
В непрерывном случае,
Z∞ Z∞
K(X, Y ) = (x − E(X))(y − E(Y ))f (x, y)dxdy.
−∞ −∞
Свойства ковариации:
K1. K(X, Y ) = K(Y, X) - следует из определения.
50
K2. Для любых вещественных чисел a и b справедливо равество K(aX + b, Y ) = aK(X, Y ) - следует из того, что E(aX + b) = aE(X) + b и определения ковариации.
K3. K(X, X) = D(X), K(Y, Y ) = D(Y ) - следует из определения. K4. K(X, Y ) = E(XY ) − E(X)E(Y ).
˚˚ − − −
Действительно, K(X, Y ) = E(YX) = E((X E(X))(Y E(Y ))) = E(XY
E(X)Y − XE(Y ) + E(X)E(Y )) = E(XY ) − E(X)E(Y ) − E(X)E(Y ) + E(X)E(Y ) =
E(XY ) − E(X)E(Y ).
Следствие 5.2 В силу свойства E5 математического ожидания (см. утверждение 5.1) для независимых X и Y ковариация равна нулю - K(X, Y ) = 0.
Определение 5.9 Если K(X, Y ) = 0, то случайные величины X и Y называются
некоррелированными.
Замечание 5.4 Из следствия 5.2 получаем, что независимые случайные величины являются некоррелированными. Но достаточно просто построить пример (см. пример 5.6 ) некоррелированныx случайных величин, являющихся зависимыми. Для случайных величин X и Y , имеющих нормальное распределение, незави-
симость и некоррелированность эквивалентны.
Если K(X, Y ) = 0, то случайные величины X и Y зависимы.
Утверждение 5.2 Если случайные величины X и Y некоррелированны, то D(X +
Y ) = D(X) + D(Y ).
Доказательство. Равенство E(XY ) = E(X)E(Y ) следует из свойства K4 ковариации, тогда равенство D(X + Y ) = D(X) + D(Y ) следует из свойства D4 дисперсии.
Это свойство можно распространить на любой набор попарно некоррелированных случайных величин. Пусть X1, X2, . . . , Xn набор состоящий из n попарно некор-
релированных случайных величин, тогда D( |
n |
|
n |
|
=1 Xi) = |
i=1 D(Xi). |
|||
Замечание 5.5 В общем случае, используяPковариацию,i Pможно записать |
||||
n |
n |
|
K(Xi, Xj ). |
|
D i=1 Xi = |
i=1 D(Xi) + 2 i>j |
|||
X |
X |
X |
|
|
В частности
D(X + Y ) = D(X) + D(Y ) + 2K(X, Y ).
Коэффициент корреляции.
Определение 5.10 Коэффициентом корреляции случайных величин X и Y назы-
|
|
|
|
K(X, Y ) |
|
|
|
|
|
D(Y ). |
r(X, Y ) = |
, где, как и ранее, σ(X) = pD(X), σ(Y ) = |
|||||
|
σ(X)σ(Y ) |
|||||||
вается величина |
|
|
|
|
|
|
||
p |
|
|
|
|
|
|
|
|
51
Cвойства коэффициента корреляции.
r1. Коэффициент корреляции по модулю не превышает единицу: |r(X, Y )| ≤ 1.
Доказательство. |
Обозначим |
˘ |
|
˘ |
|
X = (X − E(X))/σ(X), Y = (Y − E(Y ))/σ(Y ) |
|||||
- соответствующие случайным величинам X и Y центрированные нормированные |
|||||
случайные величины. Тогда |
|
|
|
|
|
˘ |
˘ |
|
˘ |
˘ |
˘ ˘ |
E(X) = E(Y ) = 0, D(X) = D(Y ) = 1, K(X, Y ) = r(X, Y ). |
|||||
Отсюда |
|
|
|
|
|
˘ |
˘ |
˘ |
|
˘ |
˘ ˘ |
0 ≤ D(X ± Y ) = D(X) + D(Y ) ± 2K(X, Y ) = 2 ± 2r(X, Y ).
Следовательно, r(X, Y ) ≤ 1, r(X, Y ) ≥ −1.
r2. Из следствия 5.2 получаем, что если случайные величины X и Y независимы, то r(X, Y ) = 0.
r3. Равенство |r(X, Y )| = 1 достигается тогда и только тогда, когда величины X и Y связаны линейной зависимостью Y = aX + b, где a 6= 0 и b постоянны. При этом ar(X, Y ) > 0, то есть a и r(X, Y ) одного знака. Доказательство см. например, [1].
Пример 5.6 Рассмотрим систему случайных величин (X, Y ) с плотностью
|
(0, |
|
− |
x2 |
+ y2 |
> a2. |
f (x, y) = |
2(a2 |
− x2 |
|
y2)/(πa4), x2 |
+ y2 |
≤ a2, (a > 0), |
Найти коэффициент корреляции r(X, Y ).
Вычисления будем производить в полярной системе координат. Обозначим A = 2/(πa4), тогда
|
|
|
|
|
|
2π |
a |
|
|
2π |
a |
|
|
|
E(X) = E(Y ) = A Z0 |
dϕ Z0 |
r cos ϕ(a2 − r2)rdr = A Z0 |
cos ϕdϕ Z0 |
r2(a2 − r2)dr = 0. |
||||||||||
Cледовательно, K(X, Y ) = E(XY ) и |
|
|
|
|
||||||||||
r(X, Y ) = 0. |
2π |
|
a |
|
|
|
2π |
a |
|
|
|
|||
R |
R |
|
|
|
R |
R |
r3(a2 − r2)dr = 0, |
|||||||
K(X, Y ) = A |
|
|
dϕ |
|
|
r sin ϕr cos ϕ(a2 − r2)rdr = A sin ϕ cos ϕ |
||||||||
|
0 |
|
0 |
|
0 |
0 |
|
|
|
|||||
В данном случае, случайные величины X и Y являются некоррелированными, |
||||||||||||||
но зависимыми, так как при x (−a, a), y (−a, a) |
|
|
|
|
||||||||||
|
√ |
|
|
|
|
|
|
|
|
|
|
|
||
|
a2−x2 |
|
|
|
|
|
|
|
|
|||||
|
√ Z |
|
|
|
|
8 |
|
|
8 |
|
||||
fX (x) = A |
|
|
|
|
|
(a2 − x2 |
− y2)dy = |
|
(a2 − x2)3/2, fY (y) = |
|
|
(a2 − y2)3/2. |
||
|
|
|
|
|
3πa4 |
3πa4 |
||||||||
|
|
2 |
2 |
|
|
|
|
|
|
|
|
|
|
|
− a −x |
|
|
|
|
|
|
|
|
|
|
|
|||
Вне указанных промежутков обе плотности равны нулю. Следовательно, f (x, y) 6= fX (x)fY (y).
Ковариационная матрица системы случайных величин.
52
