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

книги из ГПНТБ / Бендат Дж. Измерение и анализ случайных процессов

.pdf
Скачиваний:
73
Добавлен:
25.10.2023
Размер:
16.48 Mб
Скачать

78 Глава 7

ная функция x(g) есть действительная

п-значная функция g,.t

где п — целое

число. Тогда

 

 

 

 

Р [g<Г S (х (fe)) < g + Agl

 

nP [X <~ X Ik) ^ X-f- Ал-]

 

 

A g

 

A g

 

 

пР [.V<( a (к)

:gi x + A-v]

A.v

(3.12)

 

Ал:

 

 

' Ag

 

 

 

 

Следовательно,

при dg/dx Ф 0

 

 

 

 

p{g) = n p { x )№

up (X)

(3.13)

 

I dg/dx I

 

 

 

 

При использовании этой формулы1) необходимо заменять перемен­ ную X в правой части соответствующим значением g.

Пример3.3. Распределение гармонического колебания. Гар­ моническое колебание постоянной амплитуды X и постоянной частоты /0 можно рассматривать как случайную величину, если начальная фаза 0 = G(£) есть случайная величина. Рассмотрим, в частности, случай фиксированного времени ^0; тогда случайное гармоническое колебание описывается выражением

X (k)=x (Ѳ)=Х sin [2л/0/0 + Ѳ(/?)].

Пусть C(k) имеет равномерную плотность распределения

^ _1(2л)_1 при 0 ^ Ѳ^ 2я,

\ 0 при других Ѳ.

Какова в этом случае плотность распределения р(х) функции

*(*)?

Здесь прямая функция х(Ѳ) однозначна, а обратная функция Ѳ(х) двузначна. Подставив в формулу (3.13) 0 вместо х и х вме­ сто g, получим

2 р (Ѳ)

при

dx

Р (*) = dx/dQ

Ф 0.

где

 

 

^ Г = Х сос (2л/0/0+ Ѳ)= Х-j/1 —s.nz (2л)0/0 + 0) = у/ X*хй.

L) При выгоде данной формулы авторы, по-видимому, предполагают, что величина dg/dx одинакова при всех п значениях х, соответствующих данному g. В общем случае это не так, и правильно записать

Р (Xj)

Р ( е > = Е - ^ - ' dg

I ~г— dxit

где Х[ = xt(g) — функции, обратные функции g(x).— Прим. ред.

Математические основы анализа случайных процессов

79

^Таким образом,

 

р(х) = \ ( я / * " х*)~1 ПРИ \х \ < х >

(3.14а)

\ 0 при 1хI ^ X

 

[см. формулу (1.24)]. Соответствующая функция распределения гармонического колебания имеет вид

0

при X< —X,

X

 

Р(х) = I р(Ю^= = 4 “ ( т _ + агсз'п ' т )

при — Х < х < X, (3.146)

- X

 

1

при X > X.

Пример 3.4. Неравенство Чебышева. Пусть x(k) — произ­ вольная случайная величина со средним значением р.х, средним значением квадрата Ф'2 и дисперсией о2. Пусть одномерная плот­ ность распределения, вообще говоря неизвестная, есть р(х). Тогда

уЛр (x)dx ^

г

х-р (x)dx ^ е2

С

р (x)dx,

Ф* _ I

j

I

—со

 

I X I ^ 8

I л I ^ 8

 

поскольку подынтегральное выражение неотрицательно и х2 ^ е2

6ВО всей области интегрирования [|х(

ej.

Следовательно,

 

Р [ | * ( * ) | > е ] = f

Р{x)dx <

-Hi-

(3.15)

\ x \ ^

e

 

 

 

Заменим теперь x(k) на x{k) — р,^. Тогда Ф2 заменится на а2 и со­ отношение (3.15) примет вид

 

Р ||х ( А ) - а л | ^

е ] < - ^ .

(3.15а)

В частности, положив е — сах, найдем

 

 

Р IIX (k) цх I ^

сах] < - і - ,

(3.156)

что

эквивалентно

 

 

*

Р [ ! * ( * ) —

> 1 ------

(3.15В)

Любая из форм неравенств (3.15) известна под названием неравен­ ства Чебышева.

80

Глава 3

 

3.1.2.

Совокупность двух случайных величин

 

Рассмотрим теперь две случайные величины x(k) и y{k), где k

точки

соответствующего выборочного пространства.

Пусть Р(х)

и Р(у)

— функции распределения, соответствующие

случайным

величинам х{Іг) и у(/г). Совместная функция распределения Р(х, у)

определяется как вероятность, соответствующая подмножеству точек k выборочного пространства, удовлетворяющих одновремен­

но неравенствам x(k) ^ х и y(k) < у.

Это можно записать так:

Р(х, у ) = Р [ к ( к ) ^ х

и y{k)^y] .

(3.16)

Очевидно, совокупность всех точек k удовлетворяет неравенствам

x(k)

^ оо и yik) ^ оо. Отсюда

 

 

 

 

 

 

 

 

Далее, очевидно,

Р(оо, оо) = 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Р ( - о о , у)=Р(х, —оо) = 0.

 

 

 

 

(3.17)

Предполагая, как и ранее, что случайные величины непрерыв­

ны,

можно .

определить

совместную

плотность

распределения

р(х,

у) соотношением1)

 

 

 

 

 

 

 

 

 

 

 

р(х,

 

Р |х < х (ft) <

X -f Дх и у<Су (k)

<

у +

Ду) і

(3.18)

 

у) Пт

 

 

Длі/

 

 

 

 

 

 

 

 

 

Д.ѵ->)

 

 

 

 

 

 

 

 

 

 

 

 

 

Ді/->0

 

 

 

 

 

 

 

 

 

 

Заметим,

что

 

р (х,

у)

^ 0,

 

 

 

 

 

 

(3.19)

 

 

 

 

 

 

 

 

 

 

 

 

 

j

j р (х, y)dxdy= 1,

 

 

 

 

(3.20)

 

 

 

— со

 

 

 

 

 

 

 

 

 

 

 

 

У

 

 

 

 

 

 

 

 

 

 

 

Р(х,

у)—

1 j p ( 6 , ^didrt,

J - \ d-

^ j l \

=

p(x,

у). (3.21)

 

 

 

— СО — СО

 

 

 

 

 

 

 

 

 

 

Одномерные

плотности

распределения

находятся

по

формулам

 

 

 

 

 

СЗО

 

 

 

 

 

 

 

 

 

 

р(х) =

\p(x,y)dy,

 

 

 

 

 

 

 

 

 

 

— СО

 

 

 

 

 

 

 

 

 

 

 

 

 

00

 

 

 

 

 

 

 

 

 

 

 

р (У) =

J

Р (X,

y)dx.

 

 

 

 

(3.22)

_______________________

 

 

— ОО

 

 

 

 

 

 

 

 

*) Из того, что x(k) и y(k) порознь непрерывны, еще не следует существо­ вание р(х, у) в смысле (3.18). Если при любом фиксированном значении x(k) значения y(k) могут сплошь заполнять некоторый интервал, то это на прак­ тике может служить интуитивным основанием предположения о существова­ нии р(х, у).— Прим. ред.

Математические основы анализа случайных процессов

81

ІбрСЛИ

 

р(х, У)=р (х)р(у),

(3.23)

то случайные величины х{іг) и у(/г) называют статистически неза­ висимыми. Отсюда следует, что для статистически независимых случайных величин

Р(х, у)=Р(х)Р(у).

(3.24)

Математическое ожидание действительной однозначной не­ прерывной функции g(x, у) двух случайных величин х(!г) и y(k) определяется выражением1)

№[g(x, у)\ = j j g (х, у)р{х, y)dxdy.

(3.25)

Если, например, g(x, у) = (х(/г) — px)(y(k) — pty), где рх и ру

средние значения случайных величин x(k) и y(k) соответственно, то формула (3.25) определяет ковариацию Сху величин x(k) и y(k)

Сху= М [(х (k)— рх)(у (k) —Ру)\ =

= М [х (k)y (/г)]—М [*(£)] М (£)]=

СО

 

= J J {х— рх){у— ру) Р (х, y)dxdy.

(3.26)

Заметим, что Схх — а |

есть дисперсия случайной величины x{k),

определяемая соотношением (3.11).

 

Между ковариацией и стандартными отклонениями случай­

ных величин x(k), y(k)

существует простое соотношение ^

4

 

\Сху\ < о хоу.)

(3.27)

Таким образом, ковариация случайных величин не превышает произведения их стандартных отклонений. Доказательство при­ водится ниже в подразд. 3.2.1.

Как следует из приведенного неравенства, значение норми­

рованной

величины

 

I

рxy= S « L - ,

(3.28)

 

а.\9ц

 

а) Эта формула справедлива и для разрывных функций g(x, у), во всяком вручав если интеграл в правой части (3.25) имеет смысл и

СО

I f g (х> У) I Р (X, y)dxdy < оо. —Прим. ред.

6 -2 2 4 4

—oo
(3.33)

82

Глава 3

называемой коэффициентом корреляции, находится в предела?^ между --1 и +1. Говорят, что случайные величины х(Іг) и y{k) некоррелированы, если их коэффициент корреляции равен нулю. Понятие некоррелированности следует отличать от приведенного ранее понятия независимости случайных величин. Заметим, что если х(Іс) и y(k) независимы, то, как следует из равенства (3.23),

СО

 

 

м ІХ (k)y (&)] = j* j* xyp (x,

y)dxdy =

 

— CO

 

 

 

CO

 

\ xp (x)dx

I yp (,y)dy = M \x (&)|M (/г)].

(3.29)

— 00

— CO

 

Поэтому Схи и, следовательно, рГІІ равны нулю, так что независи­ мые случайные величины всегда некоррелированы. Обратное утверж­ дение, вообще говоря, неверно: некоррелированные случайные величины не обязательно независимы. Однако в физически важ­ ных ситуациях, когда рассматриваются две или более случай­ ных величии с нормальным (гауссовским) законом распределения, некоррелированность означает также и независимость случайных величин. Доказательство приводится ниже в подразд. 3.1.3.

Пример 3.5. Распределение суммы двух случайных величин. Пусть x(k) и y(k) — случайные величины с совместной плот­ ностью распределения р(х, у). Определим одномерную плотность распределения p{z) суммы случайных величин

z (k)=x (k) + у (li).

(3.30f

Если X +■ у — z, то имеем у = z x.

Отсюда

p(x,' y)=p(x, z— x).

(3.31)

При любом фиксированном значении г величина х может прини­ мать значения в пределах от — оо дооо. Следовательно,'

oo

 

p (z)= \p{x, zx)dx.

(3.32)

— oo

Таким образом, для нахождения одномерной плотности распреде­ ления суммы случайных величин z{k) достаточно знать совмест­ ную плотность распределения слагаемых. Если x(k) и y(k) — независимые случайные величины с одномерными плотностями

распределения pt(x)

и р2(у)

соответственно то,

р(х, у) =

= РЛФгІУ) = Рі (х)Р2

z X)

и

*

 

00

 

 

Математические основы анализа случайных процессов

83

іф.І.З. Гауссовское (нормальное) распределение

Говорят, что случайная величина x(k) подчиняется распреде­ лению Гаусса (нормальному), если ее плотность распределения имеет вид

р (х) = (& / '2 л ) 1е х р w a)-

(3.34)

где а и b — действительные постоянные и Ь > 0. Нетрудно убе­ диться,. что а и b суть соответственно среднее значение и стандарт­ ное отклонение случайной величины x(k). Действительно,

w

М(&)] = I' хр (x)dx = а рл,

*оо

М[(x(k)—а)2] = J а)2р (x)dx=b2= a x2.

Таким образом, плотность распределения случайной величины, подчиняющейся нормальному закону, есть

Р(*)=(Дс V ^

) : ехр

(х-р*)2

 

(3.35a)

2оі

 

цПо определению нормальная функция распределения

имеет вид

 

X

 

 

 

Р (x)=(axf/r2л )_1

j expj^-

2аі

dl.

(3.356)

 

 

 

 

С физической точки зрения важность нормального закона рас­ пределения объясняется, в частности, существованием централь­ ной предельной теоремы [30, 33], которая утверждает, что сумма большого числа совместно действующих независимых случай­ ных величин распределена в довольно общем случае по закону,

близкому к нормальному.

Рассмотрим конкретный пример. Пусть хД/г), x2(k), ..., xN(k) представляют собой N взаимно независимых случайных величин

с любыми законами распределения. Пусть далее рг и of

— сред­

ние значения и дисперсии случайных величин хД/г), t =

1, 2,...

Ц.., N. Рассмотрим случайную величину

 

* ( * ) = 2 аЛ (*)»

 

і=і

 

6*

 

84

 

 

 

 

Глава 3

 

 

 

 

где ог — произвольные фиксированные

постоянные. Как следуе£

из приведенной

ниже

формулы

(3.137),’

среднее

значение р,*

и дисперсия

случайной величины ' x(Ji)

определяются форму­

лами

 

 

 

 

 

 

 

 

 

 

|А* = М [*(£)] =

М

 

 

 

N

 

 

 

N

2

аіхі (/г)

 

М Г**

=

і=і

 

 

_ « = і

_

 

/=і

 

 

 

 

 

 

N

 

 

 

N

о!= М [(*(£)—(ц)2] = М

J ] ai(Xi (k) — щ.)5 = £

щві

 

 

 

 

 

і = 1

 

 

 

 

 

Последнее выражение есть следствие взаимной

независимости

величин xt(k) и

Xj{k)

при і ф у'1).

Центральная

предельная тео­

рема утверждает,

что

при

достаточно общих

условиях и при

N —- оо распределение случайной

величины x(k)

стремится к

нормальному закону

со

средним значением рж и дисперсией ст^.

Рассмотрим

теперь N случайных

величин

хг(к), х2{к), ...,

Xfij(k), вообще говоря коррелированных. Введем следующие обозна­ чения для соответствующих средних значений, дисперсий и кова­ риаций:

\Х; (/?)],

а?=М [(х; (А )-^)*],

 

 

S '

С и — М ІС*І (70 — !J-i) (Xj

],

С „= а ? .

Говорят, что случайные величины

имеют

N-мерный гауссовский

(нормальный) закон распределения, если соответствующая ^-мер­ ная плотность распределения определяется выражением

ѳхр

f

- 2- | С | - ^ 2

I cij I

(xi — щ) (*/ — И/)

р{хъ X2, . . .,Хы)=

- Ч__________

=1__________________=_

V

 

(2n)NU I С 11/2

 

 

 

 

(3.36)

гдеС — ковариационная матрица с элементами

Си , определенная

ниже, I С I— детерминант

матрицы С

и |СІ;-|

— алгебраические;

') Первая

из этих формул справедлива в самом общем случае, лишь бы

М \xi\k)} были конечными,

а вторая 'праведлива для любых попарно некор­

релированных

случайных

величин.— Прим. ред.

Математические основы анализа случайных процессов

85

элемента Сг;-. Точнее

 

Сп С12 . .

C-IN

С= С2і С22 . . C2n

_CfjiCN2 ..

C n n -

а алгебраическое дополнение | Си | элемента Cij определяется как детерминант (N — 1)-го порядка, полученный вычеркиванием і-й строки и /-го столбца из матрицы С и умножением полученного минора на (— \)i+J.

, Наиболее важная особенность /^-мерного нормального закона

распределения состоит в том,

что он

полностью определяется

средними значениями рг и ковариациями Сі}. При N = 1 соотно­

шение (3.36) принимает вид

)

 

 

Р(хі) — (°іѴ

(*і — Ці)2

 

1 е2ха?Р

(3.36а)

 

 

 

т. е. сводится к нормальной плотности распределения, определен­

ной ранее формулой (3.35а).

двумерная нормальная

плотность

При

N =

2

совместная

распределения

имеет вид

 

 

 

Р(х1, *я)=

 

 

 

 

 

 

 

ехр (—[2(1—р?,)]->

,'Хі — Hl \ 2

п ' [Хі — Pi \ (Хъ — (Х2

 

ц=___!___________

 

 

----------------------»

 

 

 

 

2 7 1 0 ^ 2 / 1—Рі2

(3.336)

 

 

 

 

 

 

 

где р12 =

С,2/аі0 2

коэффициент корреляции случайных величин

Xi(k) и х2(А'). Заметим, что

при некоррелированных

величинах

Хі(£) и Xo(k), т.

е.

при р12

=

0 имеем

 

 

 

 

 

р{хъ

хг)= р (х 1)р(хг),

 

откуда следует, что если величины x^k) и х2(£) некоррелнрованы, то они также независимы. Этот результат не остается справедли­ вым для любого закона распределения.

Можно записать аналогичные формулы и для распределений

более высокого порядка при N — 3, 4, 5, ....

Нетрудно показать,

дто при произвольном N и при некоррелированности всех воз­

можных пар

нормально

распределенных

случайных величии

(т. е. при рг;-

= 0,

і Ф /)

эти случайные величины будут и взаим­

но независимыми

в вероятностном смысле:

 

 

р (хг х2, . . . ,

хы)= р (xjp (х2) . . . , p ( x N).

&

 

 

Таблица 3.1

 

Некоторые плотности распределения

 

Виі иаспредслеинп

 

 

Плотность распределения

 

Дискретное

. _

 

р (* ) = Л 6 а) + ( х

Ь) + • • •

 

1 2

%

• • • + N Ö { x - n ) ,

1

 

--------1L

где А - г В + ■■■+ N

Однородное (прямоуголь­ ное)

1

1

( (6 — а ) - 1 при а < х ^

Ь,

а

Ь

Р (А) = 1

 

 

 

1 0 при других X

ц.

Гармонического процесса

 

|Ѵ|7|

р

^ _

{

(~ V X*

при \ х \ < Х ,

 

 

 

-X

X

 

 

ІО при других X

 

Гауссовское (нормальное)

 

 

»

 

 

 

 

 

 

 

 

 

 

 

р (X) =

(а,

/ Ъ - Г

1

е -

 

 

 

 

/ f a :

 

 

 

 

 

 

 

 

Рэлея

| /

\

_____р ix) =

1

^

e

- v2/2f2

при х > 0 ,

^

 

 

О

 

 

 

 

1

0

при других X

 

Максвелла

 

 

 

р

(X) =

 

 

 

 

 

 

 

/

V .

(

^

f / т

-

е-Х 2/2С 2 ПР И

х > 0 ’

 

О

 

 

 

 

ѵ

0 при других X

 

і

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Усеченное

.

Ср ( х )

В предположении,

 

что исходная

плот-

 

 

1

\

ность распределения рх(х) определена

 

 

на интервале ( — с » ,

оо), усеченное рас-

 

а

t ~

пределение

имеет

вид р ( х ) = Срх (х)

 

 

 

 

при а < X < Ь; р

(х) = 0 при другііх х,

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

j р (х) dx = С 1 P i (х) dx = 1

 

 

 

 

 

со

 

 

 

а

 

Математические основы анализа случайных процессов

87

Продолжение табл.

3.1

*

В и д р а сп р е д е л е н и я

П л о т н о с т ь р а с п р е д е л е н и я

 

Срезанное

*

В

лДСгМ'

Л

аі~

В предположении, что исходная плот­ ность распределения рх{х) определена на интервале (—о о , о о ) , срезанное рас­

пределение имеет вид

[ р, (X) при а < X < 6 ,

А5 (ха) при X ~ а,

р (.V) =

Bö (х Ь) при х = Ь,

О при X < а пли X > Ъ,

где У p ( x ) d x ~ j p1(х) dx + А +- B = I

—со

о

Важность /Ѵ-мерного нормального закона распределения в физических проблемах объясняется, в частности, существованием

многомерной центральной предельной теоремы [30, 33J, которая утверждает, что распределение векторной суммы большого числа взаимно независимых N-мерных случайных величин стремится при достаточно общих условиях к іѴ-мерному нормальному заікону.

Ч В табл. 3.1 дана сводка некоторых частных плотностей рас­ пределения, используемых в теории вероятностен.

3.2. Стационарные случайные процессы

Случайный процесс

\xk(l)}, — оо < і < оо (называемый так­

же временным рядом1),

или стохастическим процессом) представ­

ляет собой ансамбль действительных (или комплексно-значных) функций, который задается посредством описания его вероят­ ностной структуры. Из соображений удобства в последующем под. переменной t будет подразумеваться время. Каждая отдельная функция xk (/), где t — переменная, а индекс k фиксирован, назы­ вается выборочной функцией. На практике выборочная функция (или какая-либо конечная реализация выборочной функции вре­ мени) может рассматриваться как результат наблюдения в от-, даіьном опыте. Возможное множество результатов опытов опре­ деляет выборочное пространство индекса k\ это пространство мо-*)

*) Под временным рядом обычно принято понимать случайный процесс с дискретным временем.— Прим, перге.

Соседние файлы в папке книги из ГПНТБ