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

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

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

88

 

Глава 3

 

жет быть счетным или несчетным. Для любого N и любых фикси^

рованных моментов tx,

/2>

tN значения хк((х), xk(t2), ..., xk(tNy

представляют собой N

случайных величин в функции индекса к.

Необходимо, чтобы для каждого N

существовала однозначная

JV-мерная функция распределения. Ансамбль выборочныхфункций,

образующих случайный процесс, изображен на рис. 1.10.

Отдельная выборочная функция xk(t), вообще говоря, не может

описывать весь случайный

процесс

(xfc(0), которому она при­

надлежит. Однако при определенных условиях, описанных нижр, оказывается, что в частном случае эргодических случайных про­ цессов можно получить нужную статистическую информацию о процессе в целом по результатам надлежащего анализа одной^ произвольной выборочной функции. В случае двух случайныхпроцессов (xft(£)} и [ук(1)} соответствующая задача состоит в оценке совместных статистических характеристик случайных процессов по результатам надлежащего анализа одной произволь­

ной пары выборочных функций xk{t) и ук{і)-

Рассмотрим два

произвольных случайных процесса {хА(0}

и !*/*(/)}• Первые

искомые статистические характеристики —

средние значения по ансамблю при произвольных фиксированных

моментах t — определяются в соответствии с формулой

(3.8)

М 0 = М[**(0].

 

М 0 = М Ы 0 1 .

(3.^.7)

Вообще говоря, эти средние значения различны для разных мо­ ментов времени и их следует вычислять отдельно для каждого ЛГ Иначе говоря,

Рх (к) Ф Ѵ"х (к) >если

Ф /2,

 

Ру (іі) Ф Ру (*а). если

к Ф к-

(3-38)

Ковариационные функции при произвольных фиксированных

значениях t x

= t ,

= / + т определяются формулами

 

Cx (t,

t + т) = М f(хк(/) —[X* (f))(*é(t + т) — ( t н- т))],

 

C y ( t ,

t + x ) = M

( y k ( t ) - v. y ( i ) ) ( y k ( t + x ) - v . y ( t + x ) ) \ ,

(3.39)

Сху ік

^ + Т)= М [{хк (I)—(J.J,(t))(yk {t + t) Ру (t + t))J.

(3.40)

В общем случае эти величины различны для разных комбинаций

t i и І г . Заметим, что при т = 0 ( t 1 = t 2

= t )

 

Сх (к

0 = М 1(хк (t) - ц х { т

= Фх(/),

£

сЛк

0 = М| (Ук(і)-Ѵ-ит = Ф у ( і),

(3.41)

Сху (к t) = М Кхк(0 -|Х Ж(Щук ( t ) - {Xy(0)1= СХу (0.

(3.42)

\ y k (t)},
{ y k(t)\

 

 

Математические основы анализа

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

89

ѴГаким

образом,

значения

C x {t,

і)

и

C y (t, t)

ковариационных

функций представляют собой

не что

иное,

как обычные дис­

персии

процессов (дЦ/)}

и \ y k {t)}

при фиксированном

значе­

нии

t,

а C xy(t,

і) — ковариация

процессов {хЛ(/)} и

\ y k (t)}.

Как

и ранее, значения ковариационных функций в общем случае

зависят от t.

По ансамблю можно определить и другие статистические характеристики, зависящие от трех и более фиксированных мо­ ментов времени. Таким образом., вероятностная структура слу­ чайного процесса будет описываться все подробнее. Однако если

при фиксированном значении t

процессы \ x k{t) ) и

\ y k { t ) , подчи­

няю тся двумерному

нормальному распределению,

то и каждый

' из процессов [ xk (t)\

и (//*(/)}

распределен по нормальному за­

кону. В этом случае

определенные выше средние значения и ко­

вариационные функции полностью описывают вероятностную структуру процессов1). Именно по этой причине основное внимание в данной главе уделяется этим двум статистическим характери­ стикам и их связи со спектральной плотностью.

Если средние значения рЦ/) и \by {t), а также ковариационные функции Cx(t, t + т), Cy(t, t + т), Cxß , t + т) оказываются оди­ наковыми для различных моментов времени (т. е. независимыми от начала отсчета), то случайные процессы {xé(/)j и назы­ ваются слабо стационарными2). Если все возможные распределе­ ния вероятностей, описывающие процессы {xft(i)) и не зависят от переносов во времени, то процессы называются строго Iстационарными. Поскольку средние значения и ковариационные функции зависят только от одномерных и двумерных функций распределения, класс строго стационарных процессов представ­ ляет собой подкласс слабо стационарных процессов. Однако для гауссовских случайных процессов существование слабой ста­ ционарности означает одновременно и существование строгой ста­ ционарности, так как все возможные функции распределения полностью определяются средними значениями и ковариацион­ ными функциями. Таким образом, для нормальных случайных процессов оба понятия стационарности совпадают.

х) Это замечание неверно; однако в практических приложениях подтвер­ ждение гипотезы о нормальности одномерных или двумерных распределений

делает правдоподобной гипотезу о нормальности процесса. Вообще же спра-

N

■ледлив следующий факт. Если любые случайные величины ви да^

нор-

1*

ред.

і-і

 

мальны, то и процесс д:(() нормален.— Прим.

/ + т), Cy(t, г/ +

т) от­

2) Независимость от t функций (хД/), Рй10*

Cx(t,

носится к стационарности процессов {.ѵА(0}. {'<’*(0};

независимость-от f фун­

кции CXy{t,i + т) называется свойством стационарной связанности этих про­

цессов (в широком смысле).— Прим. ред.

•90

Глава 3

 

3.2.!.

Корреляционные (ковариационные) функции

£

Для рассматриваемых далее в этой главе стационарных слу­ чайных процессов |хДД) и |(/Д/)} средние значения представляют -собой не зависящие от времени постоянные величины. Для всех t

•справедливы равенства

 

СО

 

ід = м [хк(01=

j Хр (x)dx,

 

___

—со

 

'

со

 

^= М Гі!*(0]=

\ у р Ш у>

(3-43)

где р(х) и р(у) — одномерные плотности распределения, соответст­ вующие случайным величинам хк{1) и ук{(). Ковариационные функции стационарных случайных процессов также не зависят от /.

Определим для произвольных фиксированных / и х функции

Rx (т)=

м [хк(/)хк (I + т)],

 

RVN =

МІ!/л(/)ук(/ + т)],

(3.44)

Rxy(t)=№[xk(t}yk(t + x)I,

(3.45)

где символ R введен вместо С для того, чтобы отличить эти соот-^ ношения от формул (3.39) и (3.40). При средних значениях, шД равных нулю, функции R и С не совпадают. Характеристики'

Rx(x) и R„(x) называются автокорреляционными функциями, а Rxy(x) взаимной корреляционной функцией процессов (хДД)

и|#Д0] ■

Необходимое и достаточное условие того, чтобы ДДт) могла

служить автокорреляционной функцией некоторого слабо стацио­ нарного случайного процесса (хДЛ), состоит в том, что Rx{x) = = Rx(х) и Rx(x) — неотрицательно определенная функция [15]. Можно также доказать, что Rx(x) — непрерывная функция х, если она непрерывна в точке х — 0. Аналогичным образом взаим­ ная корреляционная функция Rxy(x) непрерывна при всех х, если Rx(x) или Ry(x) непрерывна в точке т = 0.

Для

пары стационарных

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

(хД^)}

и

{уДОі совместная

плотность

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

d^

, х2) пары

слу­

чайных

величин

Xi = xk{t),

х2 — xk{t + т)

не

зависит

от

/;

ня*

зависит от^ также и совместная плотность распределения p(ylt у2% соответствующая паре случайных величии ух = yk{t), у* — yk{t 4* + т). Это утверждение справедливо также и для совместной плот­ ности распределения р(хх, г/2), связанной с парой случайных ве-

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

9t

^ичин Xi = xk(i), 1/2 = уkit + т)- Корреляционные функции вы­ ражаются через плотности распределения следующим образом:

оо

Я*(т)= Л * і*2 Р (* і> x2)dxldx2,

 

Ru (т) =

J I У1У1Р Ua

, y<i)dyidy2,

(3.46)

 

 

— CO

 

 

 

 

CO

 

 

I

Rxy (t) =

j j ЧУ-іР

yjdxjdyr

 

 

 

— CO

 

 

При произвольных значениях p* и

ковариационные функции

связаны с корреляционными функциями соотношениями

 

Сх (Ѵ = ЯЛ(t) —

(3.47)

C.ru (Т) = R.XU (Т) — IC'P'y

'Таким образом, при нулевых средних значениях корреляционные функции тождественно равны ковариационным. Заметим, что,

согласно

определению,

два

стационарных случайных

процесса

.некоррелироваиы,

если

их

взаимная ковариационная

функция

^ и(і) =

0 при всех т. Как следует из формул (3.47), это условие

,наполняется, если Rxy{t)

= рд.ру при всех т. Следовательно, когда

Rxy(x) =

О ПРИ всех

случайные процессы будут некоррелирова-

ны только в том случае, если хотя бы одно из средних значений р* или р(/ равно нулю.

Из гипотезы стационарности следует, что автокорреляцион­ ные функции Rx(?) nRy(i‘) представляют собой четные функции т:

Rx(— T)= R x (r)>

 

Ry ( - r ) = R y (x).

(3.48)

Взаимная корреляционная функция не обладает ни свойством четности, ни свойством нечетности, но удовлетворяет равенству

Rxy( - T ) = Ryx(r).

(3-49)

Примеры автокорреляционных и взаимных корреляционных функ-

і О

даны на рис. 1.15 и 1.20 и в табл. 3.2.

’'

Неравенство для взаимной корреляции. Верхняя грани­

ц а взаимной корреляционной функции определяется неравенством

ІЯ « (т) І ’ <Я*(0)Я„(0)>

(3.50)

о о

 

Н е к о т о р ы е а в т о к о р р е л я ц и о н н ы е ф у н к ц и и

 

 

 

Таблица 3.2

 

 

 

 

 

 

 

Т и п

 

 

 

А в т о к о р р е л я ц и о н н а я ф у н к ц и я

 

 

П о с т о я н н о й в е л и ч и н ы

1—

 

Rx (

с 2

 

 

 

 

 

 

 

 

 

 

т ) =

 

 

 

 

 

 

 

 

 

 

О

 

 

 

 

 

 

 

 

 

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

Ѵ / Ѵ

/ ^

 

Я * ( т ) —

о

c o s 2 n f 0T

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

Б е л о г о ш у м а

 

 

,

. '

 

 

 

 

 

 

 

 

 

 

fr

i

R x

U )

=

a ö (

t

)

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

Б е л о г о

ш у м а

н а в ы х о д е

/ -k

 

 

 

 

 

 

( s i n 2xB t \

ф и л ь т р а н и ж н и х ч а с т о т

 

 

л л ( т ) - в в ^

^

j

Б е л о г о

ш у м а

н а в ы х о д е

 

 

 

R a. ( t ) _

 

a ß ^

 

 

 

j c o s 2 j t f „ T

п о л о с о в о г о ф и л ь т р а

 

0

 

 

 

 

 

Э к с п о н е н т а

Э к с п о н е н т а , у м н о ж е н н а я н а к о с и н у с

Э к с п о н е н т а , у м н о ж е н н а я н а с и н у с о и д а л ь н у ю ф у н ­ к ц и ю

-- ^ W - e - a | T l

У 4" * ”

# х ( т ) = е a , T | c o s 2 r t / 0T

■ 04,-

•^*/\] V / N ? " -

R * ( t ) = e— a 1 T 1 ( 6 c o s 2 я / „ т +

?

+ c s i n 2 n / 0 1 x | )

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

93

^ірторое доказывается следующим образом. Для любых действи­ тельных постоянных а и b

М [ (ах (і) + by(l + т))2] > О,

поскольку под знак математического ожидания входит величина, которая всегда неотрицательная. Это неравенство эквивалентно следующему:

a~Rx (0) + 2abRxy (т) + b*Ra (0) > 0.

Следовательно, полагая b Ф 0, можно записать

( ^ ) Ч (0) + 2 ( - f )

w + ^ .v (0) > 0.

Это квадратное уравнение относительно alb не имеет различных действительных корней, поскольку его левая часть неотрицатель­ ная. Поэтому дискриминант уравнения должен быть неположи­ тельным:

Дискриминант = 4R2Xy(т) —4Rx (0)Ry (0) < 0 .

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

.

 

 

R%(i) = \(Rxu('i)\2 ^ R x (0) R y (0),

 

 

 

что и требовалось доказать.

на

x(t) — рх и

y(t

+

т) — py,

Заменяя

x(t)

и

y(t + т)

можно совершенно аналогичнымобразом доказать

неравенство

для взаимной

ковариации

 

 

 

 

 

4

 

 

 

I Сху (т) 12 <

Сх (O)Cj, (0).

 

 

(3.51)

Положив

x(t)

y(t),

придем

к

неравенствам

|# х(т)|

< і?х(0)

и I Cx(t) I

^

СДО)

для

всех т. Из неравенства (3.51)

следует, что

функция коэффициентов корреляции (нормированная взаимная ковариационная функция) может быть определена в виде

Рху(х) = у

Сху И)

(3.52)

Сх (0)С„(0)

Эта функция удовлетворяет при всех х неравенству

 

— 1 Д Рад(т).< и

(3.53)

р£сли средние значения рж и р^ равны нулю, то ржДт)

принимает

вид

R*у (т)

 

Рад(т) = /

(3.54)

R x (0} R y (0)

 

Функция рЖ!/(т)

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

цессов (хД/)) и

(г/Д0) при сдвиге процесса {г/Д0} относительно

4 (0 1 на X.

 

R x (v).

94 Глава 3

Суммируя изложенное, можно сказать, что корреляционный свойства стационарных случайных процессов {хА(/)| и [ук (і/)} могут быть описаны четырьмя функциями: R x {х ), R y (т), R xlДг) и R,jx(f)- Эти функции достаточно определить только для неотри­ цательных т, поскольку их значения при т < 0 определяются свойствами симметрии Іформулы (3.48) и (3.49)1.

Пример 3.6. Автокорреляционная функция гармонического процесса. Пусть (хДЛ) = {Хьіп[2л/0^ + Ѳ(&)|( — гармонический процесс, для которого X и / 0 — постоянные и 0(/е) — случайная величина с равномерной плотностью распределения /?(Ѳ) в интер­ вале (0 ,2л). Найдем автокорреляционную функцию

Очевидно, при любых фиксированных значениях і справедлив

вы равенства

h;'

хк (t)=X sin [2л/У + Ѳ(/г)]=*,(Ѳ),

xk(t + r)=Xsin

|2л/0 (/ -f- т) + 0 (/<•)] = X, (0).

Согласно соотношению (3.44),

Rx (T) = M ІЧ (()xk il + t)1=M fxx (Ѳ)х3 (Ѳ)],

где

Р (0) =

(2л)_1

при 0 < 0 <1 2л,

 

 

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

 

 

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

 

 

 

 

 

 

 

 

 

 

“ЕГ~ J sin (2я/Ѵ + °) sin [2л/о(Н-т) +

0^0=

Г''

 

 

 

 

 

X*

О г

 

(3.55)

 

cos 2л/0т

 

есть автокорреляционная функция гармонического процесса, что и утверждалось в формуле (1.31).

Пример 3.7. Автокорреляционная функция прямоугольного волнового процесса. Рассмотрим физическую ситуацию, когда выборочная функция xk(t) случайного прямоугольного волнового процесса (хД/)} может принимать значение либо с, либо —с, причем число перемен знака в интервале (/, t -}- т) случайно; моменты перемен знака независимы, перемены происходят со средней интенсивностью Я. Допустим также, что поведение про-, цесса внутри интервала (I, t + т) не зависит от поведения про­ цесса вне этого интервала. Определим Ап как событие, состоящее в том, что внутри интервала (/, t + т) будет наблюдаться точно п перемен знака. Такой физический процесс описывается распре,^'-- лением Пуассона [1] и вероятность события А„ равна

Р(А )=JMlL>le~x ИІ.

U

п

п \

 

 

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

95

І^айдем

автокорреляционную функцию

процесса (хА(/)).

сле­

^

Автокорреляционная

функция

может быть вычислена

дующим образом.

Любое данное произведение xk(t)xk[t ф- т) рав­

но

либо

с3, если

знаки

xk(t) и xk(l

+ х)

одинаковы, либо

—с3,

если знаки xk(t) и xk(l + х) различны. Суммарная вероятность

появления

с2 равна Р(Л()) +

Р(А2) + Л(Л4) + ...,

а суммарная

вероятность появления

—с2

есть Р{Аг) +

P(AZ) +

Р(АЬ) +

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

 

 

 

 

 

 

 

 

СО

 

 

 

 

Rx (т) = М \xk{t)xk (t + х)]=с2 2 (

- \УР {Ап) =

 

\

 

со

п=0

 

 

 

= с2е~11т I

 

 

 

 

к

V, (— 1)" Ж і і і 1 = с-е-2ХI' I .

(3.56)

 

 

п—1

nl

 

 

 

 

 

 

 

 

 

Пример 3.8. Автоюрреляционная функция суммы двух стационарных процессов. Пусть случайный процесс (г/*(/)) пред­ ставляет собой сумму двух стационарных случайных процессов (х1>А(0| и 1*2,* (01. так что каждая выборочная функция имеет вид

Ук ( О = 0 Л , *(/) + а2х2к (0,

где % и а2— постоянные. Допустим также, что процессы {хг А(/)} и {*2,*{0} могут быть коррелированы. Найдем автокорреляцион­

ную функцию Ry(r). Согласно формуле (3.44),

К (т)= М \Ук(І)Ук{t + т)]=

= М кaxxhk (/) -I- а2х2>к ( i ) ) ( a (/ + т) + а2хг к(/ + т))] == = аЩ [*м (0*і,k (t + 01 + Oia2M \x1<k(t)x2>k ({+ т)] +

+ fljOjM [x2)Ä (/)jc1>A(t + t)| + a\M [x2)* (t)x2k (t + x)j =

=alRXl (x) -f axa2 (x)+Rx#x(x)j + a\Rx„(x).

Таким образом, для построения автокорреляционной функции суммы стационарных процессов необходимо знать автокорреля­ ционные и взаимные корреляционные функции слагаемых.

'*}.2.2. Спектральная плотность

/

Определения спектральной и взаимной спектральной плот­ ностей приведены в подразд. 1.3.4 и 1.4.3, где рассмотрены прак- т-Аеские аспекты операций фильтрации, возведения в квадрат и осреднения выборочных реализаций. Описанный в этих разделах ^подход приводит к односторонним спектральным плотностям GQ), существующим только в интервале неотрицательных ча­

9fi

Глава 3

стот (0, оо). Более непосредственно спектральные плотности опреу* деляются как преобразования Фурье соответствующих корреля­ ционных функций, рассмотренных в подразд. 3.2.1. Такой подход ведет к двусторонним спектральным плотностям

 

ОО

 

5 (/) =

(x)e~2llilxdx,

(3.57)

— СО

заданным на всей частотной оси о т — оо дооо. Эта функция су­ ществует, если существует функция R(x) и если

СО

 

J I R (т) I dx < ос.

(3.58^

— СО

 

Обратное преобразование Фурье функции S(f)

дает

СО

 

R(т)= JS (/)e 2*/fTd/.

(3.59)

— СО

Предположим теперь, что авто- и взаимные корреляционные функции Rx{x), R„{x) и RxU{т), определенные формулами (3.44 и (3.45), существуют. Тогда двусторонние спектральные и взаим­ ные спектральные плотности имеют вид

 

со

 

 

Sx ( f ) =

J

dx,

 

 

— CO

 

 

 

CO

 

 

Sy{f) =

j*Ry(x)e~2^ lrdx,

(3.60)

 

— CO

 

 

 

0 0

 

 

£ ,,(/) =

f Rxy (т)е-2п^

dx.

(3.61)

 

—оо

 

 

Для решения'практических задач оказывается полезным исследо­ вать функции'Д(т) и 5(/) с помощью дельта-функций. Уравнения (3.60) и (3.61), а также приводимые ниже эквивалентные равен-, ства называются соотношениями Винера — Хинчина.

Из свойств симметрии корреляционных функций стационар­ ного процесса [формулы (3.48) и (3.49)] следует, что

Sx ( - f ) = S x {f),

*

S y { - f) = Su{f),

(3.62),

Sxu{ - f) = S 'Xy(f)=Sux(f).

(3.63)

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

97

»Как видно из этих формул, двусторонние спектральные плотности представляют собой действительные, неотрицательные и четные функции частоты /, а взаимные спектральные плотности — комп­ лексные функции частоты /. Соотношения (3.60) могут быть при­ ведены к виду

 

СО

с о

 

Sx (/) =

J

Rx (т) cos 2я/та'т=2 j' Rx (т) cos 2nfxdx,

 

Sy (f)= j* Ry (t ) c o s 2nfxdx=2^Ru(т) cos 2nfxdx.

(3.64)

s

—CO

0

 

обратно,

 

СО

 

 

 

 

 

 

Я*(т) = 2 j' Sx (/) cos 2njxdf,

 

 

 

Ry{x)=2 Sy ([) cos 2nfxdf.

(3.65)

Односторонние спектральные плотности Gx(f)

и Gy(f) с частотами /

в интервале (0, оо) имеют вид

 

 

 

 

Gx (f)=

2Sx (/)

при 0 <

/ <

оо,

 

0

при других

/,

 

 

 

 

 

_ [ 25у (/)

при 0 <

/ <

о о ,

(3.66)

GJ!)--

0

при других

/.

 

 

 

 

Именно эти характеристики измеряются на практике путем пря­ мой фильтрации. Однако при математических выкладках анализ часто упрощается за счет использования функций Sx(f) и Sy(f), определенных во всей области частот от — со до оо, и экспонент с мнимыми показателями степени. Необходимо уметь правильно обращаться с обоими представлениями спектральных плотностей; здесь будут использоваться они оба. Соотношение между функ­

циями S(l)

и G(/) иллюстрируется рисунком 3.1.

 

Односторонние спектральные плотности Gx(]) и Gy(f) связаны

чс корреляционными

функциями ^ ( т )

и Ru(x) стационарного про­

цесса соотношениями

 

 

 

 

 

СО

 

 

 

 

Gx (/) = 4

j" Rx (т) cos 2nfxdx,

0 <

/ < оо,

 

4i

 

б1

 

 

 

 

 

со

 

 

 

)

Gy (/) = 4

j^ y (т) c o s 2nfxdx,

0 <

Д<;оо.

(3.67)

 

 

о

 

 

 

7^22«

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