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

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

.pdf
Скачиваний:
11
Добавлен:
30.10.2023
Размер:
19.04 Mб
Скачать

1 2 0

СТАЦИОНАРНЫЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ

[ГЛ. III

4.Значение корреляционной функции R (т) при любом т не мож

превышать ее начального значения R (0), т. е.

 

Я ( 0 ) > |Д ( т )|.

 

(3.19)

Действительно, мы всегда

имеем право написать неравенство

или

[*(/)

+

а-(/ +

т)]2>

0

 

 

 

 

 

 

 

 

**(0 +

* 2(* +

* ) >

+

2* ( /) * (/ +

т).

Беря среднее по

времени

от

обеих

частей

этого неравенства

и учитывая формулу (3.17), получим

формулу (3.19).

Заметим, что неравенство (3.19) непосредственно следует из рас­ смотренного выше смысла функции R (т).

Выражением (3.12) формально можно пользоваться не только применительно к стационарным случайным функциям, но и примени­ тельно к обычным периодическим функциям. Это оказывается прак­ тически необходимым при экспериментальной обработке записей случайных функций.

В виде примера найдем «корреляционную функцию» для случая, когда

х (/) = a sin (со/ -f- ср),

понимая при этом под термином «корреляционная функция» просто результат применения к функции а (/) операции, выражаемой инте­ гралом (3.12).

Имеем:

Я(т) =х(0*(*-И)=

 

 

т

 

 

 

=

-2

f a2sin (со/ -)- ср) sin (со/ -f- сот

ср) dt =

 

 

 

 

 

 

 

1

т

 

 

 

=

Г а2

 

 

 

Yjr

/ -g- [cos сот — cos (2со/ + сот-|-2ср)] dt =

 

 

 

1 Г sin (2“ 7" + «от -)- 2ср) -f- sin (2со7*— сот —2<р) Ц

=

-r

____

{cosu>l —

---------------------------

5=-----------

:-----

(3.20)

При достаточно большом Т с любой степенью точности можно

написать:

а2

R (т) = - у cos сот.

Мы видим, что функция R ( т) имеет тот же период, что и функ­ ция х (/), но в отличие от нее является четной и не зависящей от фазы ср.

41 КОРРЕЛЯЦИОННАЯ ФУНКЦИЙ 121

Точно так же легко показать,

что

для

функции

П

 

 

* ( 0 = во + 2

ak sin ( V

+ <р*)

«корреляционная функция» R (т)

имеет

вид

 

(3.21)

Таким образом, если сигнал х(^) представляет собой, например, стационарную случайную функцию с наложенной на нее постоянной составляющей, то корреляционная функция R (т) будет иметь вид, изображенный на рис. 3.1. При этом ее начальная ордината R ( 0)

Я{Г)

 

 

 

г

 

 

О

 

 

 

Рис. 3.1.

Рис. 3.2.

равна

среднему значению квадрата

х 2 сигнала х (/), а ее конечная

ордината

R (оо) — значению квадрата постоянной составляющей си­

гнала

х (t).

 

Если

же сигнал x(t) представляет собой стационарную случайную

функцию с наложенной на нее периодической составляющей, то кор­ реляционная функция также будет содержать периодическую соста­ вляющую с тем же периодом и, следовательно, будет иметь вид, изображенный на рис. 3.2. Здесь через /^(т) обозначена корреляци­ онная функция, соответствующая случайной составляющей функ­ ции x(t).

Взаимные корреляционные функции в отличие от рассмотренных выше корреляционных функций не являются четными, но, как это

легко видеть

из (3.15),

 

 

Д*уОО = Яу*(— т).

(3.22)

Кроме того,

можно показать, что

 

1 Л ^ ( 0 ) У /г у ( 0 ) > | / г , у (т)|.

( 3 . 2 3 )

122

СТАЦИОНАРНЫЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ

[ГЛ. II!

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

Иногда (обычно в тех случаях, когда х Ф 0) оказывается удоб­ ным вводить в рассмотрение нормированную корреляционную функцию

п м —

[ас(t) — x) [х (г + т) — х]

Н ) ~

Р - Р

Очевидно, что

 

 

Р ( 0 ) = 1 .

б. Спектральная плотность

Спектральная плотность 5 (ш) определяется*) как преобразование Фурье от корреляционной функции R(т), т. е.

00 00

S ( ш )=

J R (т) е- /шт dx =

2

J R (т) cos шт dx

(3.25)

 

—oo

 

0

 

и, следовательно,

 

 

 

 

 

oo

 

oo

 

Я(т) = -^ -

j S (ш) е/шт rfu) =

JL

J S (u>) cos xu) du>.

(3.26)

 

— CO

 

— CO

 

Формулы (3.25) и (3.26) показывают, что корреляционная функ­ ция и спектральная плотность представляют собой преобразования Фурье друг для друга, так же как и переходная и передаточная функции.

Понятие спектральной плотности 5 (ш) может быть связано с по­ нятием текущего спектра Хт(/ш) стационарного случайного про­ цесса x(t).

!) Это определение вытекает из следующей теоремы А. Я. Хинчина> лежащей в основе теории стационарных случайных процессов. Для того, чтобы функция R (т) представляла корреляционную функцию некоторого Непрерывного стационарного стохастического процесса, необходимо и доста­ точно, чтобы ее можно было представить в виде

СО

 

R (т) = J cos сот dF(a).

(3.26а)

— СО

Если предположить, что функция F(a>) является дифференцируемой при всех о, т. е. что d F { ш) = S (со) dco, то мы получим формулу (3.26) как част­ ный случай формулы (3,26а).

5]

 

 

 

 

СПЕКТРАЛЬНАЯ

ПЛОТНОСТЬ

 

123

Текущий

спектр

X j

U

^ ) м о ж н о

 

определить

следующим

образом.

Пусть

 

д:T(t) =

x{t)

при

— T < £ < Т,

 

 

 

 

 

 

 

 

хт(0 — 0

 

при всех остальных значениях t.

 

Тогда

текущим спектром ХтО'ш)

называется преобразование Фурье

для функции

 

т.

е.

 

 

 

 

 

 

 

 

 

 

 

оо

 

 

т

 

 

 

 

Х г (у'ш) =

J

х тЦ)еч ы d t = f xT(t)e~jal dt.

(3.27)

 

 

 

 

- о о

 

 

 

 

 

Введем в рассмотрение функцию

 

 

 

 

 

 

 

 

 

 

т

 

 

 

 

 

 

 

RT(^) = -^f

j x T(t)xT{t-\-i)dt.

(3.28)

 

 

 

 

 

 

-V

 

 

 

 

Преобразование

Фурье

5г(ч>)

для функции /?г(т):

 

 

ОО

 

 

 

 

СО

со

 

 

5 г (а))=

f

RT(i)e~Jmdx = -^r j

e~J™dx j

x T(t) x T(t +

t) dt =

— CO

 

 

 

 

 

“ CO

 

— CO

 

 

 

 

CO

 

CO

 

 

 

 

 

 

 

=

W

I

dz f xtWxtV+

“0 e~JU>V+X) eM dt =

 

 

 

— CO

— CO

 

 

 

 

 

 

 

 

 

CO

 

 

CO

 

 

 

 

=

4 r

f

х т(0е/шГЫ f

* r (H -T )e --'tol/+x)rfT.

(3.30)

 

 

- o o

 

 

 

— oo

 

 

 

 

 

Произведя в

(3.30)

замену

переменной ^—f—-с = X и имея в виду,

что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

СО

 

 

 

 

Хт(у'ш) =

Хт(— yu)) =

J x T(t)e/mtdt,

 

получим:

 

 

 

 

 

 

 

 

— СО

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

со

 

 

 

 

 

 

 

 

Sr (co)=

J

^ ( т ) е ~ ;штйт = -1 г ^ г (ус0) ^ ( » .

(3.31)

 

 

 

— СО

 

 

 

 

 

 

 

Учитывая, что функция /?г (т) четна, найдем:

ОО

ОО

СО

 

f

RT( т ) е~ушт dx = 2 J RT(т) e~,m dx =

2 J RT(т) cos шт dx,

(3.32)

- О О

0

О

124 СТАЦИОНАРНЫЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ [ГЛ. III

и, следовательно, на основании

(3.31) можно

написать:

со

 

со

 

О

 

О

 

В качестве

спектральной плотности 5 (со) можно рассматривать

математическое

ожидание

 

 

 

5 ( с о )=

М [Sr ( с о )]

(3.33)

для совокупности функций x T(i). В случае стационарного случай­ ного процесса совокупность функций x T{t) может быть получена из одной единственной реализации, если взять время наблюдения доста­ точно большим и разрезать реализацию на куски. Отсюда следует, что если заменить в силу свойства эргодичности математическое ожидание в (3.33) средним во времени на достаточно большом ин­ тервале изменения t, то вместо (3.33) можно написать

5 ( с о ) ^ Х т ( |/ с о ) | 2 . (3.34)

Соотношение (3.26) между корреляционной функцией и спектральной плотностью имеет исключительно большое значение для теории стацио­ нарных случайных процессов.

Значение этого соотношения прежде всего заключается в том, что оно позволяет определить функцию спектральной плотности S ( со) по заданной аналитически или в виде графика корреляционной функ­ ции R(~), или обратно, корреляционную функцию /?(т) по заданной функции спектральной плотности 5 (со).

Так, например, если

Я(т) = е-1"Ч

то

если

то

Кроме

того,

значение соотношений (3.25), (3.26) заключается

в том, что

они

позволяют

использовать

при изучении связи между

функциями

R (т)

и 5 (со)

все известные

свойства преобразования

51

 

 

 

СПЕКТРАЛЬНАЯ

ПЛОТНОСТЬ

 

125

Фурье. Так,

например, если R (x)— монотонно убывающая

функция

от

х, то S (ш) — также

монотонно

убывающая функция от

ш. Чем

уже

функция

R(x),

тем

более

пологой и широкой является

функ­

ция

S (ш). Если

/?(т) стремится

к нулю в течение очень короткого

промежутка времени

Д,*то 5'(«О

сохраняет приблизительно

постоян-

ное

значение

до

частоты

 

2к

 

 

порядка -д -. Такого рода спектр 5 (ш) часто

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

®2(*1. t{, хг, t2) = wi (xl, ty) w 2(x2, t2),

и мы имеем так называемый абсолютно случайный процесс.

Если 5 (со) имеет максимум при некоторой частоте шс и обладает

симметрией относительно этой точки, т. е. если

 

 

S(io) =

F(a> — шс)

 

И

F(u) =

F { - и ) ,

 

 

 

ТО

 

 

 

.

 

с о

 

R (т)

cos соус У* F (со) cos сот йш

(3.35)

 

 

— СО

 

при условии, что ширина 5 (со) мала по сравнению с сос. Предполо­ жим, например, что 5 (со) = const в полосе частот 2 , содержащей частоту шс. Тогда

^ (^

i cos ш^ - т 51п;х

(3-36)

при условии, что 2 сой. Таким образом, корреляционная функция будет иметь характер затухающих колебаний с частотой шс.

Пользуясь соотношением (3.25), можно получить выражение для спектральной плотности также и в том случае, когда функция x(t), помимо случайной составляющей, содержит периодические члены или когда среднее значение этой функции отличается от нуля. Очевидно, что если функция х (t) содержит периодическую составляющую частоты (1>с, то при достаточно больших Т

2^ - 1Х т(ju>c) [2= оо

и спектральная плотность 5 (ш) имеет в точке ш = шс разрыв непре­ рывности. Если же х Ф.О, то функция S(w) будет иметь разрыв непрерывности в начале координат, т. е. при ш = 0.

1 2 6 СТАЦИОНАРНЫЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ [ГЛ. III

Математическую трудность, связанную с получением выражения для 5 (со), в указанных случаях можно обойти при помощи введения понятия о дельта-функции 8(ш).

Имея в виду, что

(см. гл. I, §

7)

 

 

 

 

СО

 

 

 

 

 

 

IГ cos сох di — i

8(ш)(

 

 

о

 

 

 

 

(3.37)

/ cos (о0т cos шт di =

 

 

о(со — со0),

 

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

плотности в пред­

положении, что

 

 

 

 

 

 

 

х (t) =

а0+

c o s^ jf +

ср).

(3.38)

Для корреляционной

функции R ( т) в этом

случае получим:

 

R (т) = а* -}-

a'f

 

 

(3.39)

 

cos tOjT.

Подставляя (3.39)

в (3.25) и учитывая (3.37),

найдем:

S (со) =

о»8 (со)+

ai

(со — Ш1).

(3 .4 0 )

- £ 8

Точно так же легко видеть, что если функция x(t) представляет собой некоторый тригонометрический ряд

п

 

(3.41)

*(0 = fl0+ S

aftcos(co^ + <?*),

л-1

 

 

то спектральная плотность будет представлять собой разрывную функцию, состоящую из отдельных линий:

5 (ш) = а\ о (со) 4 - ^

-J о (со — | coft |).

(3.42)

к=1

'

 

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

Покажем, что

х°- =

J S (со) dco.

(3.43)

— СО

si

Имеем:

 

СПЕКТРАЛЬНАЯ ПЛОТНОСТЬ

127

СО

со

 

J | Xj (уш) |2

dm =

^ Xj (— усо)X2 (уto) dm =

 

— со

 

 

— со

 

 

 

 

 

 

со

со

 

 

 

=

J Хт(— усо) dm J

х 7 (/) е~№ dt

 

 

 

 

— оо

— оо

 

 

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

 

 

 

оо

со

 

со

 

со

 

J | Хт (Уш) |2(1ш— J

х т(/) dt J Хт (— У'ш) e-Jwl du>= 2тс J

(t) dt.

— ОО

— ОО

 

— ОО

 

— СО

 

Итак,

 

 

 

 

 

 

со

 

Т

 

со

 

 

f

x \ { t ) d t = j

x \{ t ) d t = ^

f

|X j (у‘ш) |2dm.

(3.44)

—оо

 

T

 

—oo

 

 

Разделив (3.44) на 2T, при достаточно больших Т получим, при­ нимая во внимание (3.34):

Т

оо

оо

** = ■&■/

/ ^

1^ (У ш )|2й ? о )= ^ у >5 (ш)

Т

—СО

—оо

Физический смысл функции S ( со) можно пояснить следующим образом. Если предположить, что x T(t) есть ток, то выражение (3.44)

можно рассматривать как энергию, которую он рассеивает в сопро­ тивлении 1 ом, а выражение

Т оо

р 1 = Т г f 4 ( 0 d t = i

/ 2 7 ! ^ С /® )I2 rf®

(3.45)

—Г

—оо

 

как среднюю мощность этого тока для промежутка времени 2 Т. Опуская в левой части индекс Т, получим согласно (3.44) выра­

жение для средней мощности

СО

Р = - ^ / S(a)da .

(3.46)

— СО

 

Из (3.46) ясно, что спектральная плотность S ( со) в рассматри­ ваемом случае имеет размерность энергии, вследствие чего она иногда называется энергетическим частотный спектром функции x(t).

Заметим, что спектральная плотность 5(ш), как это следует из формулы (3.34), не содержит, так же как. и корреляционная функ­ ция, никаких сведений о сдвигахфаз отдельных составляющих функ­ ции X (t).

128

СТАЦИОНАРНЫЕ

СЛУЧАЙНЫЕ

ПРОЦЕССЫ

[ГЛ. III

Если

функции

x(t) и у (7)

относятся

к двум различным

множе­

ствам, причем

 

т

 

 

 

 

 

 

 

 

 

 

 

ХгО'ш) =

f

x T(t)e-Jw<dt,

 

 

 

 

- г

 

 

 

 

 

 

 

т

 

 

 

 

 

(Уш) =

f

y T(t)e~ial dt,

 

то выражение

 

 

 

 

 

 

 

 

 

 

 

 

 

sxy ( io ) =

2 Х7 *т( j u ) Y Т (

»

(3.47)

называется их взаимной спектральной плотностью. Через

X * (j со)

обозначена функция, комплексно-сопряженная функции X(jw).

Из соотношений

 

 

 

 

 

 

 

X r (jw) =

X r ( — Уш),

 

(3.48)

 

 

У г (у ш ) = К г ( — У ш)

 

 

следует,

что

 

 

 

 

 

 

 

 

(to) = s ; y ( -

ш) = Syx ( - (О).

(3.49)

Заметим, что в том случае,

когда х Ф 0

и функция 5(ш)

может

быть представлена

в виде

 

 

 

 

 

 

 

5(ш) = 00* о (ш)+ 5.(0)),

(3.50)

где 5,(ш) — непрерывная функция, иногда оказывается удобным ввести понятие о нормированной спектральной плотности:

0(ш) = ------\ ^ ------- ,

(3.51)

-jr f S1(<o)d<o

имеющей, очевидно, размерность времени. Легко видеть, что

(х х)г = х 2 — (х)г — J S l (co)du).

(3.52)

о

Действительно,

( х — x f = { х 2 — 2 х х + ( х ) 2) =

2 ( х ) 2 + ( х ) 2 = х * - ( х ) \

61

ГАУССОВ СЛУЧАЙНЫЙ ПРОЦЕСС

129

 

6. Гауссов случайный процесс 1)

 

На практике весьма часто встречается особый вид случайного

процесса,

который обычно называют гауссовым случайным

процес­

сом. Примером могут служить шумы в радиоприемниках, по крайней мере до детекторного каскада, имеющего нелинейную характеристику. Эти процессы характеризуются тем, что все функции распределения вероятности w, о которых мы говорили в § 2, являются гауссовыми, или нормальными.

Предположим, что функция x(t) есть стационарная случайная функция с нормальным распределением. Рассмотрим некоторый доста­

точно большой промежуток

времени

2Т и предположим, что функ­

ция

x{t) является периодической с периодом 2).

Тогда эта функ­

ция

может быть разложена в ряд Фурье:

 

 

 

СО

 

 

 

 

 

х (0 = i

(ак cos иikt -j- b,. sin cukt).

(3.53)

 

ft= i

 

 

 

 

В ряде (3.53) отсутствует постоянный член, так как без ограни­

чения общности можно предположить, что среднее

значение

функ­

ции

x{t) равно нулю. Рассмотрим совокупность

функций

х (.(t)

вида

(3.53).

 

коэффициенты

ак, Ьк, вообще

Каждая из функций x t (t) имеет

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

чениями и с дисперсиями,

зависящими от

числа k,

но одинаковыми

как для ак, так и для Ьк с одним и тем же k.

 

 

Итак,

 

 

 

 

 

 

 

 

 

 

 

Ч

=

h = 0,

]

 

 

 

 

 

ака1— bjpi — сь°ы.

|

 

(3.54)

 

 

 

акЬ,

=

0, .

 

J

 

 

где ок1= 1 при

k — / и ок1 =

0 при k ф I.

 

 

 

Функция распределения

вероятности:

 

 

 

 

 

 

 

 

 

СО

 

2

 

 

та(а„

а2, . . .; Ьи Ь2 . .

) = Т Т — - е

2е* .

(3.55)

 

 

 

 

 

 

М

 

 

 

1)

L a w s o n

and

U h l e n b e c k ,

Threshold

signals.

Radiation

Lab. Se­

ries, t . 24.

от предыдущего случая, когда

предполагалось, что функ­

2)

В отличие

ция х (t) обращается

в нуль вне пределов интервала 27’,

 

 

9

Зак. 1083. В.

В.

Солодовников

 

 

 

 

 

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