Загребаев Лектсии по теории вероятностеы и математическоы статистике 2015
.pdfdϕ |
2 |
|
|
|
D[Y ] = D[ϕ(mx )] + |
dx |
|x=mx × D[(X − mx )] = |
|
|
|
|
|
|
|
dϕ |
|
2 |
]. |
(9.31) |
= |
|x=mx D[X |
|||
dx |
|
|
|
|
Этот подход легко обобщается, когда Y является функцией многих |
||||
переменных |
|
|
|
|
Y = ϕ(X1,..., Xn ). |
|
(9.32) |
||
Тогда функцию Y = ϕ(X1,..., Xn ) также раскладывают в ряд и ограничиваются линейным членом
|
Y ≈ ϕ(m |
|
,..., m |
|
|
) + |
n |
ϕ′ |
(m |
|
,..., m |
|
)(X |
|
− m |
|
|
|
|||||
|
x1 |
xn |
|
x1 |
xn |
i |
xi |
). |
(9.33) |
||||||||||||||
|
|
|
|
|
|
xi |
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
i=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Применяя аппарат числовых характеристик, получим |
|
||||||||||||||||||||||
|
|
|
|
M[Y ] = ϕ(mx ,..., mx |
); |
|
|
|
|
|
|
|
(9.34) |
||||||||||
|
|
|
|
|
|
|
|
|
|
1 |
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
dϕ 2 |
|
|
|
|
dϕ |
dϕ |
|
|
|
||||||||||
|
D[Y ] = |
|
|
D[X ] + |
2 |
|
|
|
|
|
|
|
Kij |
, |
(9.35) |
||||||||
|
dx |
|
dx |
|
|
||||||||||||||||||
|
|
i=1 |
|
|
|
|
i< j |
|
|
|
dx |
|
|
|
|
|
|||||||
|
|
|
i |
|
|
|
i |
|
j |
|
|
|
|||||||||||
ˆ |
K ... K |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
11 1n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
где K = .............. |
|
|
− корреляционная матрица аргументов. |
|
|||||||||||||||||||
|
Kn1... Knn |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Если при разложении функции в ряд оставить еще одно слагаемое
dϕ |
|
|
1 d 2ϕ |
|
( X − mx )2 , |
|
Y ≈ ϕ(mx ) + dx |
|x=mx |
( X − mx ) + |
2 dx2 |
|x=mx |
(9.36) |
то результаты линеаризации можно уточнить. Действительно, применяя к выражению (9.34) операцию нахождения математического ожидания, получим
M[Y ] = ϕ(mx ) + |
1 d 2ϕ |
|x=mx |
D[X ]. |
(9.37) |
2 dx2 |
Поправка к дисперсии просто выглядит, когда случайная величина распределена по нормальному закону
81
D[Y ] = |
dϕ |
| |
2 |
D[ X ] + |
1 |
d 2ϕ |
| |
2 |
[X ]. |
||||
|
|
x=mx |
|
|
|
2 |
|
D2 |
|||||
|
dx |
|
|
2 |
|
dx |
|
x=mx |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|||
Возможно, конечно, обобщение данных результатов на многомерный случай, однако формулы в этом случае выглядят очень громоздкими.
Лекция 10. ЗАКОН РАСПРЕДЕЛЕНИЯ ФУНКЦИЙ СЛУЧАЙНЫХ АРГУМЕНТОВ
Рассмотренные выше способы получения статистических характеристик функции от случайных аргументов имеют недостатки. Вопервых, они справедливы лишь для линейных функций. Если же функция нелинейная и используется метод линеаризации, то трудно оценить ошибку определения числовых характеристик. Во-вторых, иногда необходимо знать не только числовые характеристики, но и закон распределения функции, знание которого решает все эти проблемы. Действительно, если задана функция Y = ϕ(X ) и известна
плотность распределения случайной величины Y, которую мы обозначим g( y), то числовые характеристики – функции Y = ϕ(X ) −
могут быть найдены из следующих соотношений:
|
+∞ |
|
|
|
|
M[Y ] = my = |
y g( y)dy; |
(10.1) |
|
|
−∞ |
|
|
|
|
+∞ |
|
|
|
|
D[Y ] = Dy = ( y − my )2 g( y)dy. |
(10.2) |
||
|
−∞ |
|
|
|
Для |
определения функции g( y) рассмотрим участок |
на |
оси |
|
(a, b), |
на котором лежат все |
возможные значения |
X, |
т.е. |
P(a < X < b) = 1. На рис. 10.1 показан участок монотонно возрастающей функции y = ϕ(x).
Если случайная величина X принимает значение x, то случайная величина Y принимает значение y = ϕ(x). Иными словами, если на
82
плоскости изобразить случайную точку (X, Y), то она может лежать только на кривой y = ϕ(x).
а |
б |
Рис. 10.1. К выводу закона распределения: а − для монотонно возрастающей функции; б – для монотонно убывающей функции
Обозначим функцию распределения случайной величины Y как G( y) = P(Y < y). Из рис. 10.1,а видно, что случайная величина Y
примет значение меньшее, чем y в том случае, если случайная величина X попадет в интервал (a, x). Это означает, что справедливо соотношение
x |
|
G( y) = P(Y < y) = P(a < X < x) = f (x)dx. |
(10.3) |
a
При этом верхний предел интегрирования является функцией от y. Действительно, имея монотонную исходную функцию y = ϕ(x),
можно найти ей обратную, выражающую зависимость x от y, т.е. функцию x = ψ( y), тогда
ψ( y) |
|
G( y) = f (x)dx. |
(10.4) |
a
От функции распределения G( y) можно перейти к функции плотности распределения, в соответствии с (4.7) получим
g( y) = dG |
= f (ψ( y)) |
dψ( y) |
. |
(10.5) |
|
||||
dy |
|
dy |
|
|
83
Если функция y =ϕ(x) является монотонно убывающей, как на
рис. 10.1,б, то случайная величина Y примет значение меньшее, чем y в том случае, если случайная величина X попадет в интервал (x, b). Это означает, что справедливо соотношение
|
|
|
|
|
|
|
|
|
|
|
b |
|
|
|
|
b |
|
|
|
|
|
|
G(y) = P(Y < y) = P(x < X < b) = f (x)dx = |
|
f (x)dx. |
(10.6) |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
ψ( y) |
|
|
|
Взяв производную по переменному нижнему пределу, получим |
|||||||||||||||||
|
|
|
|
|
|
|
g( y) = dG = − f (ψ( y)) |
dψ( y) |
. |
|
(10.7) |
|||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
dy |
|
dy |
|
|
|
|||
|
Формулы (10.5) и (10.7) можно объединить в одну: |
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
dψ( y) |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
g( y) = f (ψ( y)) |
|
. |
|
|
(10.8) |
||||
|
|
|
|
|
|
|
dy |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
Действительно, |
для |
монотонно |
возрастающей |
функции |
|||||||||||||
|
dψ( y) > 0, |
|
dψ |
|
= dψ и справедливо соотношение (10.5), для |
|||||||||||||
тогда |
|
|
||||||||||||||||
|
dy |
|
|
|
dy |
|
|
dy |
|
|
|
|
|
|
|
|
||
монотонно |
убывающей |
функции |
dψ( y) < 0 |
в |
этом |
случае |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
dy |
|
|
|
||||
|
dψ |
|
|
= − dψ |
и справедливо соотношение (10.7). |
|
|
|
||||||||||
|
|
|
|
|
||||||||||||||
|
dy |
|
|
dy |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Приведем пример определения закона распределения функции от случайных аргументов.
Пример. Найти закон распределения линейной функции Y = a X + b от аргумента, распределенного по нормальному закону.
Решение данной задачи приведено в табл. 10.1, где в левом столбце находятся теоретические зависимости в общем виде, а в
правом − конкретно для данного случая.
Выражение для g( y) можно также записать в следующем виде:
|
1 |
|
|
( y − (amx + b)) |
2 |
|
|
|||
g( y) = |
|
|
|
|
(10.9) |
|||||
|
|
|
exp − |
2 |
|
2 |
|
. |
||
| a | σ |
|
2π |
σ |
|
||||||
|
x |
|
2 | a | |
x |
|
|
|
|||
|
|
|
|
|
|
|
|
|
||
84
Определение закона распределения функции |
|
Таблица 10.1 |
|||||||||||||||||
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
f (x) |
|
|
|
1 |
|
exp{− |
(x − mx ) |
2 |
|
|
|
|
|||||||
|
|
|
|
|
} |
|
|
|
|||||||||||
|
|
|
σx |
2π |
|
2σ2x |
|
|
|
|
|||||||||
y = f (x) |
|
|
|
|
y = a x + b |
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x = ψ( y) |
|
|
|
|
|
x = |
y − b |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
a |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
ψ′( y) |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
a |
|
|
|
|
|
|
|
|
|
|
| ψ′( y) | |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
| a | |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
g( y) = f (ψ( y)) | ψ′( y) | |
|
|
|
|
|
|
|
|
|
|
|
y − b |
|
|
2 |
|
|||
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
− mx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a |
|
|
|
|||||
|
g( y) = |
|
|
|
|
exp |
{− |
|
|
|
|
} |
|||||||
| a | |
σx 2π |
|
|
2σ2x |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Из выражения (10.9) видно, что |
g( y) представляет |
собой |
|||||||||||||||||
плотность распределения случайной величины с нормальным законом распределения с параметрами
my = a mx + b; σ y =| a | σx .
Закон распределения функции двух случайных величин
Пусть имеется система двух случайных величин (X , Y ) с плотностью распределения f (x, y) и случайная величина Z, связанная с
(X , Y ) функциональной зависимостью Z = ϕ(X , Y ). |
Требуется |
определить закон распределения случайной величины Z. |
|
Как следует из рис. 10.2, |
|
G(z) = P((X ,Y ) D) = f (x, y) dxdy. |
(10.10) |
D( z) |
|
85
Рис. 10.2. К выводу закона распределения функции случайных величин
Величина z в это выражение входит через предел интегрирования. Дифференцируя G(z) по z, получим плотность распределения
случайной величины Z.
Закон распределения суммы двух случайных величин. Композиция законов распределения
Пусть имеется система двух случайных величин (X , Y ) с плотностью распределения f (x, y). Требуется найти закон распределения случайной величины Z = X + Y. Поверхность z = x + y − плос-
кость, проходящая через начало координат.
В заштрихованной области лежат те значения аргументов при которых X + Y < z. Тогда из рис. 10.3 следует
|
|
|
|
+∞ z− x |
|
|
|
|
|
|
|
|
|
|
|
|
|
G(z) = |
f (x, y)dxdy = |
|
|
|
(10.11) |
|||
|
|
|
f (x, y)dy dx. |
|||||
|
D |
|
|
−∞ −∞ |
|
|
|
|
Плотность распределения случайной величины Z есть |
|
|||||||
|
g(z) = dG |
+∞ |
|
|
|
|
|
|
|
= |
f (x, z − x)dx. |
|
(10.12) |
||||
dz −∞
Исходя из соображений симметрии, можно получить
86
g(z) = dG |
|
+∞ |
|
= |
f (z − y, y)dy. |
(10.13) |
|
dz |
|
−∞ |
|
|
|
|
Рис. 10.3. К выводу закона композиции
Таким образом, если известен закон распределения системы f (x, y), то закон распределения суммы легко получается.
Если требуется найти закон распределения суммы независимых случайных величин, то говорят о композиции законов распределе-
ния.
Рассмотрим случай композиции двух законов распределения. Пусть имеются две случайные величины X и Y, соответственно под-
чиненные f1(x) и f2 ( y). Требуется найти плотность распределения
суммы этих двух независимых случайных величин. В силу независимости X и Y плотность распределения системы имеет вид
|
|
f (x, y) = f1(x) f2 ( y), |
(10.14) |
|
+∞ |
+∞ |
|
тогда |
g(z) = |
f (x, z − x)dx = f1(x) f2 (z − x)dx. |
(10.15) |
|
−∞ |
−∞ |
|
Операцию нахождения композиции двух законно распределения сокращенно записывают так: g(z) = f1(x) f2 ( y).
87
Композиция нормальных законов
Пусть каждая из двух независимых случайных величин X и Y подчинена нормальному закону распределения с плотностью рас-
|
|
|
|
|
|
|
|
|
|
1 |
|
− (x−mx )2 |
|
|
|
|
пределения, |
соответственно, |
f (x) = |
|
|
|
2σ2x |
и f |
2 |
(x) = |
|||||||
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
1 |
σx |
2π |
|
|
|
|
||
|
|
|
|
|
( y−my )2 |
|
|
|
|
|
|
|
||||
|
|
1 |
|
− |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
2σ2 |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
= |
|
|
|
|
|
y . Требуется произвести композицию законов, т.е. |
||||||||||
σ y |
|
2π |
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
найти плотность распределения g(z) = f1(x) |
f2 ( y), |
где Z = X + Y. |
||||||||||||||
|
Применим формулу (10.15): |
|
|
|
|
|
|
|
|
|
||||||
+∞
g(z) = f (x, z −
−∞
1 +∞
2πσyσx −∞
+∞
x)dx = f1(x) f2 (z − x)dx =
|
−∞ |
|
|
|
|
|
|
(x−m |
)2 |
|
( y−my )2 |
|
|
− |
x |
|
− |
|
|
|
|
2σ2y dx = |
|
||||
e |
2σ2x |
|
|
(10.16) |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(x−m )2 |
|
(z−x−my )2 |
|
|
|
|
||||||
|
|
|
|
1 |
|
|
|
|
+∞ |
− |
|
|
x |
− |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
2σ2x |
|
|
|
2σ2y |
|
|
|
|
|||||||||||
|
|
|
= |
|
|
|
|
|
|
e |
|
|
|
|
|
|
|
dx. |
|
|
|
||||||||
|
|
|
2πσ |
|
σ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
y |
x −∞ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Если в формуле (10.16) сделать замену переменных: |
|||||||||||||||||||||||||||||
|
1 |
σ2x + σ2y |
|
|
|
|
|
|
|
m |
|
|
|
|
|
z − my |
|
|
|
|
|
|
m2 |
|
(z − my )2 |
||||
A = |
|
|
; B |
= |
|
|
|
|
x |
|
|
+ |
|
|
|
|
; |
|
|
C == |
x |
+ |
|
, |
|||||
|
2 |
σ2x σ2y |
|
|
|
|
|
|
2σ2x |
|
|
|
|
|
2σ2y |
|
|
|
|
|
|
2σ2x |
2σ2y |
||||||
то формулу (10.16) можно привести к виду |
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
+∞ |
|
|
2 |
|
|
|
|
|
|
|
||
|
|
|
g(z) = |
|
|
|
|
|
|
|
|
|
|
e− Ax |
+2Bx−C dx. |
|
(10.17) |
||||||||||||
|
|
|
|
2πσ |
|
|
σ |
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
y |
x −∞ |
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Из курса математического анализа известно, что |
|
|
|||||||||||||||||||||||||||
|
|
|
|
+∞ |
− Ax |
2 |
+ 2Bx−C |
|
|
|
π |
|
− AC− B2 |
|
|
|
|||||||||||||
|
|
|
|
e |
dx = |
|
|
e |
A |
. |
|
(10.18) |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
A |
|
|
|||||||
|
|
|
|
−∞ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
88
Тогда, проводя необходимые преобразования, получим
|
1 |
|
− |
(z−mz )2 |
|
|
|
|
|
2(σ2 |
+σ2 ) |
|
|
||
g(z) = |
|
|
e |
x |
y |
. |
(10.19) |
σ2x + σ2y |
|
||||||
|
2π |
|
|
|
|
||
Как видно из выражения (10.19), случайная величина Z = X + Y также подчиняется нормальному закону с параметрами:
mz = mx + my |
и Dz = Dx + Dy . Если имеется n случайных величин |
||||||
|
|
n |
|
|
|
|
|
X1,..., Xn , то |
случайная величина |
Z = Xi , |
|
также |
подчиняется |
||
|
|
i=1 |
|
|
|
|
|
|
|
n |
|
и Dz |
n |
|
|
нормальному закону с параметрами |
mz = mx |
= Dx . |
От- |
||||
|
|
i=1 |
i |
|
i=1 |
i |
|
|
|
|
|
|
|
||
метим тот факт, что даже если нормально распределенные случайные величины коррелированны, то их сумма всё равно распределена по нормальному закону. Только выражение для дисперсии будет
|
n |
|
+ 2 Kx y |
|
|
|
иметь вид |
Dz = Dx |
j |
. Таким образом, при композиции |
|||
|
i=1 |
i |
i< j |
i |
|
|
|
|
|
|
|
||
нормальных законов получается нормальный закон.
Нетрудно понять, что закон распределения линейной функции от
n
нормально распределенных аргументов Z = ai Xi + b также будет
i=1
нормальным с числовыми характеристиками
n |
|
n |
|
+ 2 aia j Kx y |
|
|
mz = aimx |
+ b и Dz = ai2 Dx |
. |
||||
i=1 |
i |
i=1 |
i |
i< j |
i j |
|
|
|
|
|
|||
Предельные теоремы теории вероятностей, о которых мы будем говорить ниже доказывают, что при композиции большого числа любых законов, закон распределения суммы близок к нормальному закону.
Предельные теоремы теории вероятностей
Предельные теоремы теории вероятностей можно условно разделить на две группы. Первая группа объединяется под названием
«закон больших чисел», а вторая − под названием «центральная предельная теорема».
89
Предельные теоремы описывают свойство устойчивости массовых случайных явлений, суть которой в том, что конкретные особенности каждого отдельного случайного явления почти не сказываются на среднем результате массы таких явлений. В узком смысле под законом больших чисел понимается ряд теорем, в каждой из которых для тех или иных условий установлен факт приближения средних характеристик большого числа опытов к некоторым определенным постоянным. Например, теорема Бернулли утверждает, что при большом числе опытов, если вероятность события не изменяется от опыта к опыту, частота события сходится по вероятности к вероятности появления этого события. Мы уже упоминали об этом, когда говорили о статистическом подходе к определению вероятности. Приведем еще один пример с бросанием монеты. Если производится n бросков и «герб» появляется в m случаях, то отно-
шение w(n) = mn является частотой появления «герба». Понятно,
что для «правильной» монеты вероятность выпадения «герба» равна вероятности выпадения «решки», равна p = 0,5 и при каждом
броске одинакова. Пусть, например, произведено 12 серий испытаний, число бросков в каждой серии показано в табл. 10.1
Таблица 10.1
Исходные данные для испытаний
Номер |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
12 |
серии |
|
|
|
|
|
|
|
|
|
|
|
|
Число опытов |
10 |
20 |
30 |
40 |
50 |
60 |
70 |
80 |
90 |
100 |
1000 |
10000 |
в серии, n |
|
|
|
|
|
|
|
|
|
|
|
|
Результаты численного эксперимента показаны на рис. 10.4. Частота выпадения «герба» mn получена с использованием функции
генерации случайного числа из программы Excel (Сервис → Ана-
лиз данных → Генерация случайных чисел → Распределение
Биномиальное). Зависимость P m > n |
− 1 |
есть вероятность того, |
|
|
2 |
|
|
90
