Материал: Загребаев Лектсии по теории вероятностеы и математическоы статистике 2015

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

dϕ

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

Источник: https://studfile.net/preview/16708854/