Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Учебное пособие 1603

.pdf
Скачиваний:
17
Добавлен:
30.04.2022
Размер:
1.46 Mб
Скачать

Если же факт принятия величиной ξ конкретного значения x изменяет вероятность принятия СВ η значения y, то ξ и η являются зависимыми. Для работы с подобными величинами часто оказывается удобно использовать условные ФРВ и плотности вероятности. Знакомство с условными законами распределения разумно начать с анализа следующего предела

lim P{η < y | x ξ < x +

x}= lim

P{x ξ < x + x,

η < y}

=

 

 

 

 

x0

 

 

 

 

 

 

 

x0

P{x ξ < x + x}

 

 

 

 

x0

xx

y

Wξη (x, y)dydx

 

xx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= lim

 

 

 

 

 

 

Wξη (x, y)dydx =

 

 

 

 

 

 

 

x

−∞

 

 

 

 

x

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

Wξη (x, y)dy

 

 

 

y

 

 

 

 

 

 

 

−∞

 

 

 

 

Wξη (x, y)dy

 

Wξη (x, y)dy =

 

 

 

 

= lim

 

x

x

 

 

 

.

 

 

 

 

x0

 

 

 

 

 

 

 

 

 

 

 

Wξ (x)

 

 

 

 

 

−∞

 

 

 

 

−∞

 

 

 

 

 

 

 

Проанализированная выше вероятность, определяющая ФРВ СВ η, вычисленную при условии, что принятое в опыте величиной ξ значение фактически известно, называется условной функцией распределения величины η

F (y | x) = lim P η < y | x ξ < x +

x =

y

W

x, y

)

dy W

x

) (6.22)

η

x0 {

}

−∞

ξη (

 

ξ (

 

Аналогично, условная плотность вероятности величины η, зависящей от величины ξ, определяется выражением

Wη (y|x) = lim

P{yη<yy|xξ<xx}

=

Fη (y|x)

=

Wξη (x, y)

(6.23)

y

y

 

Wξ (x)

x0

 

 

 

 

y0

 

 

 

 

 

 

 

и характеризует возможность наблюдения значений СВ η из бесконечно малой окрестности аргумента y в тех опытах, где величина ξ находилась в бесконечно малой окрестности x.

Для условных плотностей вероятности тоже выполняется условие нормировки

150

 

1

Wξ (x)

 

 

Wη (y | x)dy =

 

Wξη (x, y)dy =

 

=1 ,

(6.24)

Wξ (x)

Wξ (x)

−∞

 

−∞

 

 

а приведенные ниже выражения являются аналогами формулы полной вероятности

 

Wη (y)= +∞Wξη (x, y)dx =

+∞Wη (y | x) Wξ (x)dx ,

(6.25)

 

 

 

 

 

−∞

 

 

−∞

 

 

 

 

 

 

 

(y

)=

+∞

 

y0

 

 

+∞

 

(y | x) W

(x)dx .

 

F

W

(x)

W (y | x)dy dx =

F

(6.26)

η

0

 

ξ

 

η

 

η

0

ξ

 

 

 

 

−∞

 

 

−∞

 

 

−∞

 

 

 

 

 

В итоге, аналогом формулы Байеса для непрерывных случайных величин ξ и η может служить правило

Wξ (x | y)=

 

Wη (y | x) Wξ (x)

.

 

 

Wη (y | x) Wξ (x)dx

(6.27)

 

 

 

−∞

 

 

6.6.Числовые характеристики системы двух СВ

Кнаиболее часто используемым числовым характерис-

тикам системы СВ относят смешанные начальные и централь-

ные моменты распределения. Различающиеся порядками моменты могут иметь заметно отличающийся «физический смысл»; в связи с этим, начнем с формального определения числовых характеристик, а их свойства и практическое назначение будем исследовать для каждого конкретного случая в отдельности.

151

6.6.1. Определения и общие свойства моментов распределения системы СВ

Смешанным начальным моментом распределения порядка k, r для системы случайных величин {ξ, η} называют константу

mkr {ξ,η}=

xk yr Wξη (x, y)dxdy , k 0, r 0 .

(6.28)

 

−∞ −∞

 

В этот набор констант входят как величины, характеризующие особые свойства системы, так и «классические» моменты распределений, характеризующие свойства отдельных компонент. Действительно, с учетом (6.19) для произвольной системы СВ при k = 0 получим

m ξ,η

}

=

y

r

0

 

dy = m η

.

 

x W (x, y)dx

 

(6.29)

0r {

 

 

ξη

r { }

 

 

 

 

−∞

 

−∞

 

 

 

 

 

Аналогично,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mk 0 {ξ,η} = mk {ξ}

 

 

 

(6.30)

и, в частности,

 

 

 

 

 

 

 

 

m10 {ξ,η} = m1 {ξ} = Mξ ; m01 {ξ,η} = m1 {η} = Mη .

(6.31)

Если вместо значений самих СВ использовать при усреднении отклонения величин от их математических ожиданий, то получаемые числовые характеристики называют центральными моментами распределения системы СВ

μkr {ξ,η}= (x Mξ )k (y Mη )r Wξη (x, y)dxdy .

(6.32)

−∞ −∞

 

Нетрудно заметить, что для независимых случайных величин, у которых плотность вероятности системы распадается на произведение законов распределения составляющих систему СВ, как начальные, так и центральные смешанные моменты пред-

152

ставляют собой просто произведение числовых характеристик отдельных СВ

kr {

}

 

ξ

η

k { }

r { }

 

=

;

m

ξ,η

 

 

xkW (x)dx

 

yrW ( y)dy = m ξ

m η

 

 

 

 

−∞

 

−∞

 

 

 

(6.33)

μkr {ξ,η} = μk {ξ} μr {η}.

Врезультате, информацию о системе (т.е. характеризующую особенности зависимости между величинами) несут, прежде всего, смешанные моменты порядков k 1, r 1.

6.6.2.Корреляционные характеристики случайных

величин

Моменты распределения системы СВ порядка (1,1) важны для многих практических приложений теории вероятностей и имеют собственные названия.

Смешанный начальный момент порядка (1,1) называют корреляцией величин ξ и η. Корреляцию Bξη рассчитывают

по формулам

 

 

 

 

n

m

 

 

 

 

∑∑xi y j pij

(для ДСВ)

B

= m

{ξ,η}=

i=1 j=1

(6.34)

ξη

11

 

 

 

 

 

 

xy Wξη (x, y)dydx

(для НСВ)

 

 

 

−∞ −∞

 

Смешанный центральный момент порядка (1,1) называют ковариацией величин ξ и η. Для системы непрерывных СВ ковариацию Kξη рассчитывают по правилу

Kξη = μ11 {ξ,η}= (x Mξ ) ( y Mη ) Wξη (x, y)dydx . (6.35)

−∞ −∞

153

Раскрыв в (6.35) скобки и проанализировав четыре получающихся интеграла, нетрудно убедиться, что ковариация Kξη

и корреляция Bξη связаны простым соотношением

Kξη = Bξη Mξ Mη ,

(6.36)

а значит для величин ξ и η, хотя бы одна из которых обладает нулевым математическим ожиданием, они просто совпадают.

Наконец, безразмерное отношение

rξη =

Kξη

=

Bξη Mξ Mη

(6.37)

σξ ση

Dξ Dη

 

 

 

называют коэффициентом корреляции величин ξ и η.

Примечание: Зафиксированное после (6.36) совпадение значений величин Kξη и Bξη привело к неоднозначности в терминологии;

часть авторов в своих книгах две эти величины меняют местами. В их книгах название коэффициент корреляции для rξη выглядит бо-

лее логичным, чем при предложенных выше обозначениях, однако при более широком рассмотрении это приводит к иным нелогичностям и хуже согласуется с трудами иностранных авторов.

Рассмотрим физический смысл понятия корреляции случайных величин. Для этого с учетом (6.19) и (6.35) проанализируем приведенный ниже интеграл:

 

 

 

∞ ∞ (x Mξ )

 

( y Mη ) 2

(x, y)dydx

 

 

 

 

 

 

 

 

 

 

J =

 

 

 

 

±

 

 

 

=

 

 

 

 

 

 

 

 

 

σ

 

σ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Wξη

 

 

 

 

 

 

 

 

 

 

 

−∞ −∞

 

ξ

 

 

η

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

=

 

(x Mξ )2Wξ (x)dx +

( y Mη )2Wη (y)dy ±

 

 

 

 

 

 

σ 2

σ 2

 

 

 

 

 

 

 

 

 

 

ξ

 

−∞

 

 

 

 

 

 

η −∞

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

Dξ

 

Dη

 

 

Kξη

±

 

 

 

 

(x Mξ ) ( y Mη ) Wξη (x, y)dydx =

 

+

± 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

σ

ξ

σ

η

 

 

σ2

σ 2

σ

σ

 

 

 

 

 

 

−∞ −∞

 

 

 

 

 

 

 

 

 

ξ

 

η

 

 

ξ

 

η

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

154

 

 

 

 

 

 

 

 

 

 

 

Итак,

J = 2 (1± rξη ) ,

но т.к. исходный интеграл J, очевидно, не может быть отрицательным, то коэффициент корреляции всегда удовлетворяет соотношению

rξη

1.

(6.38)

Кроме того, случаю rξη =1, очевидно, соответствует J = 0 ,

но такое возможно лишь если для всех наблюдаемых пар значений величин ξ и η соблюдается равенство

x Mξ = y Mη ,

σξ ση

т.е. если значения СВ функционально связаны между собой линейной зависимостью.

Рассмотрим также ряд представленных на рис. 34 случаев, характеризующих разные варианты взаимозависимостей между величинами системы {ξ,η}. Значения xi и yi, принятые в i-м опыте величинами ξ и η, определяют расположение на плоскости очередной i-й точки, а наиболее вероятным сочетаниям значений соответствуют области плоскости с наибольшей густотой расположения точек.

Левый рисунок соответствует случаю объединения в систему независимых случайных величин. В результате, число точек, принадлежащих каждой четверти плоскости (относительно линий x = Mξ = a и y = Mη = b ), оказывается вполне сопос-

тавимым, а потому результаты усреднения произведений (xi Mξ ) ( yi Mη ) с учетом их отмеченных на рис. 34 знаков,

дают нулевое значение ковариации Kξη , а значит и нулевой коэффициент корреляции rξη .

155

«–» y

«+»

b

(xi, yi)

 

«+»

a «–» x

rξη

0

y

 

«+»

y

(xi, yi)

 

b

 

(xi, yi)

b

 

 

 

 

 

 

«+»

a

x

 

a

x

 

 

 

 

 

rξη ≈ −0,5…−0,8

rξη ≈ +0,8…+0,9

 

Рис. 34. Варианты взаимозависимости величин ξ и η

На центральном рисунке наблюдается тенденция, при которой увеличение значений xi СВ ξ сопровождается, как пра-

вило, уменьшением значения yi парной СВ η. Как следствие,

число усредняемых слагаемых, имеющих отрицательный знак, существенно превосходит число положительных слагаемых, а потому как ковариация Kξη , так и коэффициент корреляции rξη ,

оказываются отрицательными.

Наконец, на правом рисунке представлен случай «жёсткой», функциональной зависимости между величинами. При показанном, близком к линейно нарастающему характере зависимости отрицательных слагаемых не наблюдается совсем, а коэффициент корреляции rξη принимает значения близкие к 1,0.

Итак, коэффициент корреляции rξη характеризует степень линейной зависимости между двумя случайными величинами: при жесткой, функциональной зависимости rξη =1; ослабевание связи влечет уменьшение абсолютной величины rξη . Случайные величины, для которых rξη = 0 ( Kξη = 0 ), называются

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

156

Согласно (6.33), для независимых СВ

K

ξη

= μ

{

ξ,η

}

= μ

ξ

}

μ η

0

,

(6.39)

 

1,1

 

1 {

 

1 { }

 

 

т.е. статистически независимые случайные величины являются некоррелированными. Обратное в общем случае неверно, т.е. некоррелированные СВ вполне могут оказаться зависимыми.

6.7. Двумерный нормальный закон распределения

Случайные величины ξ и η считаются распределенными по двумерному нормальному (гауссовскому) закону, если плотность вероятности системы {ξ, η} может быть записана в виде

Wξη (x, y) =

1

×

2π σξση 1r2

 

1

 

(x aξ )

2

(x aξ )(y aη )

 

(y aη )

2

 

(6.40)

 

 

 

 

 

 

,

×exp

 

 

 

2r

 

+

 

 

 

2(1r2 )

σξ2

 

ση2

 

 

 

 

σξση

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где параметры aξ , aη могут принимать произвольные значения, а остальные параметры удовлетворяют условиям σξ > 0, ση > 0 и r 1.

Найдем условную плотность вероятности величины η, входящей в двумерное нормальное распределение

Wη (y | x) =Wξ1 (x) Wξη (x, y)=

2πσξ exp +

(x aξ )2

 

 

1

×

2σξ2

2π σξση 1r2

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

(xaξ )

 

 

 

(xaξ )(yaη ) (yaη )

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

×exp

 

 

 

 

 

σξ2

 

2r

 

 

σξση

 

+

ση2

 

=

 

 

 

 

×

 

 

 

 

 

 

 

 

 

 

 

 

2(1r2 )

 

 

 

 

 

 

 

 

 

 

 

2π ση 1r2

 

 

 

 

 

(

ξ )

 

 

 

 

(

ξ )

 

 

(

 

 

 

η ) (

η )

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

ξ )(

 

2

 

 

 

1

 

 

 

 

 

 

 

(1r

2

)

 

 

 

 

 

 

 

 

 

 

 

 

 

=

×exp

 

 

 

 

 

σξ2

 

 

 

 

σξ2

 

2r

 

σξση

+

ση2

 

 

 

 

 

 

 

 

 

 

 

 

2(1r2 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

157

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

y aη

 

x aξ

 

2

 

 

 

 

 

 

 

=

 

 

exp

 

 

 

r

 

 

 

.

2π ση 1

r2

 

ση

σξ

 

 

 

2(1r2 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из сопоставления полученного выражения с одномерной нормальной плотностью вероятности (3.24) видно, что при r = 0 условная плотность вероятности СВ η приобретает вид

 

1

 

 

 

 

( y aη )2

 

Wη ( y | x) =

 

 

 

exp

 

 

.

(6.41)

π

σ

 

2

σ 2

 

2

 

η

 

 

 

η

 

 

Она уже не зависит от x и имеет классическое нормальное распределение (3.24). Как следствие, при r = 0 можно записать

Wξη (x, y) =

1

 

 

 

exp

 

(x aξ )2

 

×

 

 

 

 

 

 

 

 

 

 

 

 

 

π

σ

 

2 σ 2

 

 

 

2

 

 

ξ

 

 

 

 

 

 

ξ

 

 

 

(6.42)

 

 

 

 

1

 

 

 

 

 

 

 

( y aη )2

 

 

 

 

 

 

 

 

 

 

 

 

×

 

 

 

 

 

 

 

exp

 

 

 

 

=Wξ (x) Wη ( y)

 

 

 

π

σ

 

 

2

σ 2

 

 

2

 

 

η

 

 

 

 

 

 

 

η

 

 

и, в соответствии с (6.21), величины ξ и η оказываются независимыми.

Отсюда следует важное свойство: для нормальных случай-

ных величин некоррелированность означает также и их статистическую независимость.

6.8. Функциональное преобразование системы СВ

Пусть система СВ {ξ1 ,ξ2 } преобразуется в новую систему {η1 ,η2 } в соответствии с некоторым известным правилом

η = f

(ξ ,ξ

2

)

 

 

 

1

 

1

1

 

.

(6.43)

 

 

=

f2

(ξ1,ξ2 )

η2

 

 

 

 

 

 

 

 

 

 

 

Случаи с бесконечнозначными обратными преобразованиями, подобными обсуждавшимся в п. 5.2.2, для систем СВ хотя и возможны, но не являются типичными для радиотехнической

158

практики. Поэтому в дальнейшем будем полагать, что обратные функции, позволяющие вычислить значения x1, x2 входных СВ, необходимые для получения на выходе преобразователя значений y1, y2, являются конечнозначными (т.е. состоят из одной или конечного числа ветвей). В подобном случае можно доказать, что взаимосвязь между плотностями вероятности систем СВ на входе и выходе функционального преобразователя может быть записана в виде

Wη1η2 ( y1, y2 ) = Wξ1ξ2 (ϕ1i ( y1, y2 ),ϕ2i ( y1, y2 )) | Ji | ,

(6.44)

i

 

где функции ϕ1i (...),ϕ2i (...) - определяют одну из ветвей обратного преобразования

ξ

=ϕ

(η ,η

2

)

 

 

 

1i

1i

1

 

,

(6.45)

 

 

 

 

 

 

ξ2i =ϕ2i (η1,η2 )

 

 

 

 

 

 

 

 

 

 

а | Ji | – модуль якобиана обратного преобразования – поправоч-

ный коэффициент, характеризующий растяжение/сжатие элементарных областей при переходе от одного набора аргументов к другому и численно равный значению определителя, состоящего из частных производных обратных функций

 

 

 

 

 

 

 

 

ϕ1i ϕ1i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ji =

 

(x

, x

)

 

=

 

y1

 

y2

 

.

(6.46)

 

 

 

 

 

1i

2i

 

 

 

ϕ2i

 

ϕ2i

 

 

 

( y1, y2 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y1

 

y2

 

 

 

Пример 3: Найти функцию распределения полярных координат случайно расположенной на плоскости точки, декартовы координаты которой независимы и подчиняются нормальному распределению с нулевым математическим ожиданием и одинаковым среднеквадратическим отклонением σ .

159