
книги из ГПНТБ / Бендат Дж. Измерение и анализ случайных процессов
.pdf88 |
|
Глава 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) |
|
|
Математические основы анализа |
случайных процессов |
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. |
|
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. |
|
|
||
Следовательно, |
|
|
||
|
|
|
|
|
2Л |
|
|
|
|
“ЕГ~ 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«