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

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

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

1 3 0 СТАЦИОНАРНЫЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ t г л . нг

Имеем:

СО

*2(0 =

 

(flft cos2 о)kt

Щ sin2

=

 

 

 

 

 

 

 

 

Й«1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

СО

 

 

 

СО

 

 

 

 

00

 

 

 

 

 

 

=

2

с* =

w 2

5г (ш*} ~ i

—со

^

= с2-

(3-56)

 

 

 

Дг = 1

 

 

к= I

 

 

 

 

 

 

 

Согласно

(3.53)

jc (if)

есть

линейная

функция от

ак, bt при дан­

ном t. Отсюда следует,

что

функция x(t) при любом t имеет нор­

мальное

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

w i (x).

Как мы видим,

время

выпало

из

рас­

смотрения. Этого и следовало ожидать,

так как процесс стационарен.

Рассмотрим

теперь

вторую

функцию

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

w2(x v

х 2, т),

характеризующую

вероятность

нахождения пары

зна­

чений

х

в интервалах x v

x i - \ - d x l и х г,

x 2- \- d x2

в моменты

вре­

мени t x и t2,

причем

t2t l =

’z.

 

 

 

 

 

 

 

Так

как x ( t j и x ( t 2) — линейные функции

от uk,

bL, то мы полу­

чим двумерное

нормальное

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

л-2 ((,) и л:2(^2) по-прежнему

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

формулой

(3.56) и

 

 

х ( t j a:

(t2) =

^

(в» cos

 

cos шА/3 -f- b2 sin cofc'/t sin сол,7.) —

 

 

 

 

fc = i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

GO

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

^

с\ cos шкх

 

 

J

S (ш) cos тсо rfco = R (х) =

р (т) с2;

(3.57)

 

к= 1

 

 

 

- с о

 

 

 

 

 

 

 

 

 

здесь р(т) — нормированная корреляционная функция.

 

 

Как

мы видим, корреляция зависит только от разности т =

f21{,

как это и должно быть в силу стационарности процесса.

 

 

Согласно

(2.50) функция

распределения вероятности w2(x {, х 2, т)

определяется

формулой

 

 

 

 

 

 

 

 

 

 

 

 

 

w2(x1, х 2, Т):

 

 

 

 

 

a c a

( i - p«' J(^) + 4 -

2p.vr v3)

(3.58)

 

 

Ъ к : *

/ 1

р 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Нормированной корреляционной функции р(т) соответствует нор мированная спектральная плотность

 

 

 

— ОО

 

 

 

Легко

видеть,

что

третья функция

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

вероятности

w 3(xl, t{,

х 2, /2;

x s. t3)

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

трехмерное

нормальное

распределение, зависящее только от t?

 

t312.

 

Точно так же ясно, что все функции

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

вероятности

являются

нормальными

и зависят только

от

с2 и р(х).

 

7]

ТЕЛЕГРАФНЫЙ

СИГНАЛ

 

131

7. Корреляционная функция

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

 

для телеграфного сигнала

 

 

Вычислим

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

корреляционную

функцию

и спек­

тральную плотность функции x(t),

представленной на рис. 2.10. Пред­

положим, что

вероятность числа

перемен знака р.

в течение

проме­

жутка времени Т определяется распределением Пуассона, причем функция х (t) сохраняет постоянное абсолютное значение а.

Для вычисления корреляционной функции необходимо рассмотреть

среднее значение произведения х (t) х (t -j- т).

Это произведение

равно

-f-a2

или — а2 в зависимости от того, имеют ли

выбранные

значе­

ния х

один и тот же или противоположный

знак.

В первом

случае

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

знака, равное нулю)

в промежутке (t , ^ + т), а во втором — нечетное.

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

 

х (t) х (/-|-т ) = а2 X

[вероятность четного числа перемен знака в про­

межутке (t , / — т)] — а2 X [вероятность нечетного числа перемен

знака в промежутке (t , £-[-т)]. (3.59)

Так как по предположению вероятность перемены знака в тече­ ние интервала т не зависит от того, что происходит'вне пределов этого интервала, то отсюда следует, что это предположение будет оставаться справедливым для любого интервала вне зависимости от начала его отсчета. Поэтому вероятности в (3.59) не зависят от t и могут быть получены из (2.36), если положить 7' = | т | . Итак, на основании (2.36) найдем:.

х (0 х ( t ■ z ) == а2 [Р (0) -|- Я (2) -j- Я (4) -j- .

. . ] -

 

_ a 2[P(l)_ j_ P(3)-j_ P(5)-(_

.

(3.60)

Таким образом, для корреляционной функции на основании (3.13) получим:

У?(т) = А - аЯ'1.

(3.61)

Функция спектральной плотности S ( со) согласно (3.25) будет равна

СО

 

5 (ш) = 2й2 f е - ^ с о з ш =

. (3.62)

о

 

Корреляционные функции и спектральная плотность вида (3.61), (3.62) часто встречаются на практике. Заметим, например, что спек­ тральная плотность тока, протекающего через последовательно сое­ диненные индуктивность и сопротивление, к которым приложено напряжение с «белым» спектром, имеет вид (3.62).

!)

См. Rice

S. О., Mathematical

analysis of Random noise. BSTJ, t. 23,

№ 3,

1944; t . 24,

№ 1, 1945...................

'

9*

132

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

[ГЛ. III

8. Спектральная плотность типичного сигнала на входе следящей системы

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

УЩ

изменения угловой скорости могут соответствовать маневрам цели. Кривая изменений угловой координаты цели относительно радиоло­ катора y{t) в этом случае будет иметь вид, изображенный на рис. 3.3.

Легко видеть, что в рассматриваемом нами случае функция y(t) не представляет собой стационарной случайной функции. Поэтому, вместо того чтобы искать спектральную плотность функции у (t), мы

найдем спектральную плотность

*

dy

(ш) для ее производной у (t) —

— ,

относительно которой (рис. 3.4) в отличие от y(t) можно предпо-

!) Д ж е й м с X., Н и к о л ь с Н., Ф и л л и п с Р., Теория следящих систем,

ИЛ, 1951.

8]

СИГНАЛ НА ВХОДЕ СЛЕДЯЩЕЙ СИСТЕМЫ

133

дожить, что она является стационарной случайной функцией. Функ­

ция y{t) в первом приближении сохраняет постоянное значение в тече­ ние каждого из последующих интервалов. Более точное определение

функции y (t) может быть дано следующим образом:

У (0 = ап. * „ < * < * « +!. « = 0 , ± 1 . ± 2 .........

(3.63)

где через ап обозначены независимые случайные переменные, имею­ щие одно и то же распределение. Интервалы

ХЛ~ ^Я +1 ^II

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

Определим, прежде всего, вероятность w(x) того, что значения функций y{t) и _у(^ + т) лежат в одном и том же интервале. Из определения функции y(t) ясно, что эта вероятность будет равна вероятности отсутствия перемены знака в случае телеграфного сигнала, рассмотренного в § 7 гл. III. Ит£к,

w(x) = е~^ Hi,

(3 .64)

где через р, обозначена величина, обратная среднему значению х интервалов хп, т. е.

1

г

Найдем теперь корреляционную функцию для

y(f).

 

В зависимости от того, находятся

ли моменты

t и t- \ - x в одном

и том же или в двух различных интервалах хп, мы будем

иметь:

или

У (0 y (t - \ - x) = а„-,

 

 

 

 

 

 

 

У У (t + т) =

апап+к.

 

(3.65)

Задавшись определенным распределением моментов времени tn,

найдем вначале

среднее по совокупности для у (t) у (t -t-т).

Так как ап

не зависит от

ап+к, то

 

 

 

y ( t ) y ( t + х) = а2,

если моменты t и t -\ -x относятся к одному и тому же интервалу, и

> (0 3 '(^ Н -т) = (а)2,

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

134

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

СЛУЧАЙНЫЕ

ПРОЦЕССЫ

[ГЛ.

III

Но

вероятность того, что

t и £ -|-т

лежат в

одном и том

же

интервале, равна да(т). Поэтому, усредняя по всем возможным рас­

пределениям tn и имея в виду, что вероятность

нахождения t и t -|-т

в разных интервалах

равна

1 — -иу(т),

получим:

 

 

 

R СО = У У V +

х) =

a2w (х) 4 - ( Я ) * (1 - да (х)).

(3.66)

Пользуясь

теперь

формулой

(3.37)

для

спектральной

плот­

ности S ( ш), получим:

 

СО

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

(ш) = 2 (а2 — (a)2) Jw C O costo-rf- -\-(а)2о(ш).

(3.67)

 

 

 

 

о

 

 

 

 

 

В том случае, когда точки

tn имеют распределение Пуассона,

для

которого

а — 0

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

выражением (3.64),

фор­

мула

(3.66) сводится

к виду

 

 

 

 

 

 

 

 

S ■' (со) =

а2

 

 

(3.68)

 

 

 

- j - jj.3

 

 

 

 

 

 

 

0)2

 

 

 

9. Спектральная плотность последовательности равноотстоящих импульсов, имеющих одну и ту же форму,

но различную амплитуду

Предположим, что случайная функция y(t) с плотностью распре­ деления вероятности w(y) представляет собой последовательность

равноотстоящих

импульсов

с

периодом повторения Т, имеющих одну

 

 

 

 

 

 

п

и ту же форму, но различ­

л __[1

 

 

 

 

 

ную амплитуду (рис. 3.5).

n

 

 

 

 

Функция y(t),

очевидно,

 

 

 

 

может

быть

представлена

т

 

 

 

 

 

в виде суммы двух состав­

 

 

Рис. 3.5.

 

 

 

 

 

 

 

 

ляющих, одна

из

которых,

 

 

 

 

 

 

 

у ! ((), — периодическая и со­

стоит из последовательности

равноотстоящих импульсов с постоянной

амплитудой

у,

равной

среднему значению y{t),

т. е.

 

 

 

 

 

, y =

f

y > w ( y ) d y ,

 

 

(3.69)

а другая,

у 2(0- — флуктуирующая,

состоит из последовательности

равноотстоящих импульс'ов со

случайной амплитудой

и со средним

Значением, равным нулю.

 

 

 

 

 

 

 

Представим

периодическую

составляющую

y ^ t )

в

виде ряда

Фурье

 

 

«

 

 

со

ЗпЫ

 

 

 

 

 

 

 

 

 

 

 

 

 

• ■

Уг(0 =

2

c j ~ ,

 

 

(3.70)

 

 

_.........

 

- c q .............

 

.......................

9]

ПОСЛЕДОВАТЕЛЬНОСТЬ ИМПУЛЬСОВ РАЗЛИЧНОЙ АМПЛИТУДЫ

135

где

о

последовательно, F (у‘и>) можно рассматривать как непрерывный спектр отдельного импульса с единичной амплитудой. Спектральная плотность 5, (ш), соответствующая периодической функции y l {t):

5 1 (a.) = 2 u ^ f 2

8 ( “ - ^ г ) | / ? С/'“ )1г -

Х3 ' 72)

к = - с о

 

 

Вычислим корреляционную функцию Rz (т) для флуктирующей составляющей y 2(t). Прежде всего, заметим, что если t лежит в пре­ делах одного периода повторения, а t-\-x — в пределах другого, то

R-2СО — У2Уг + "О— О-

так как отдельные импульсы флуктуируют независимо. Таким обра­ зом, допустимые значения для т лежат в пределах (— Г, Т) и ясно, что спектр, соответствующий отдельному импульсу, является непре­ рывным. Флуктуирующий импульс удобно рассматривать состоящим из импульса у 2(Ь> умноженного на «постоянную» составляющую, колеблющуюся относительно нуля с вероятностью, определяемой w (у). Корреляционная функция этой «постоянной» равна с2, где с2 =

= { у У)2- Спектральная плотность, соответствующая этой «постоян­ ной» составляющей:

СО

S'2(m) = J /?2(T)cosmTdT==0»— 302&(a>)==c28((o). : (3.73)

— СО

Спектральная плотность самого импульса равна

Но спектральная плотность S2(a>) произведения двух корреля­ ционных функций, которым соответствуют спектральные плотности

5'(ш), 5" (ш), .равна

 

СО

 

 

5 ,(ш )=

J

5"(ю — u)S2(u)du.

(3.74)

Итак,

— ОО

 

 

 

 

 

СО

 

 

 

S2 (ш) = с2 /

~

}и) р 8 (ц) d u = ^ p (уш) |2,

(3.75J

—СО

Действительно,
части к «мощности» линейной части приблизительно равно

136

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

ПРОЦЕССЫ

[ГЛ. ill

Таким

образом, полная спектральная

плотность

 

 

5 («>) =

 

(3.76)

состоит из двух частей (рис. 3.6): 1) из непрерывного спектра

у 1 ^ ( » | 2.

имеющего такую же форму, как и спектр F (j со) отдельного импульса,

и 2) из линейного спектра, состоящего из линий, имеющих «мощность»

ТТ

расположенных друг от друга на расстоянии - у .

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

Интересно отметить, что отношение «мощности» непрерывной с2

V

СО

(3.77)

— о У

О ) у -

СО

10.Спектральная плотность последовательности

импульсов,

имеющих одну и ту же форму

и амплитуду,

но различную частоту повторения

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

10]

ПОСЛЕДОВАТЕЛЬНОСТЬ ИМПУЛЬСОВ РАЗЛИЧНОЙ ЧАСТОТЫ

137

Пусть Тср — средний период повторения импульсов, а — откло­ нение k-vo промежутка между двумя последующими импульсами

от Гер, причем

т = 0. Пусть

 

 

 

 

 

®(т) — функция

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

_

_

 

 

 

вероятности для величины т и

н

 

 

 

 

пусть

 

I II

II

I 1

И II II

..

со

 

 

 

 

 

i

W ( j ш )= J w { t ) e - ivizdx.

 

 

Рис. 3.7.

 

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

S (ш) =

— 1 ( у ш) р]

1w и») р

тср

ср

 

(3.78)

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

1) См. Mac F а г 1a n е G. G., Ргос. IRE, т. 37, № 10, 1949.

ГЛ А В А IV

ОПРЕДЕЛЕНИЕ СТАТИСТИЧЕСКИХ ХАРАКТЕРИСТИК

ИЗ ЭКСПЕРИМЕНТАЛЬНЫХ ДАННЫХ

1.Введение

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

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

Внастоящей главе описываются специальные вычислительные

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

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

Ниже излагаются некоторые возможные методы решения этой задачи,

21

ФОРМУЛЫ ДЛЯ ВЫЧИСЛЕНИЯ КОРРЕЛЯЦИОННОЙ ФУНКЦИИ

139

2. Расчетные формулы для вычисления корреляционной функции по экспериментальным данным

Для того чтобы пояснить способ определения корреляционной функции /?(т), соответствующей заданной экспериментальной кри­ вой x(t), рассмотрим формулу:

т

(4Ч1)

-V

Разделим промежуток времени Т на N весьма малых интервалов Д

так,

чтобы функция

х (/)

мало изменялась

на протяжении

интер­

вала Д (рис. 4.1), т.

е. положим, что

 

 

 

 

 

 

T = NA.

 

(4.2)

Будем придавать

t и т дискретные

значения,

кратные Д:

 

 

 

 

t = уД,

v =

1, 2, . .'.

 

(4.3)

 

 

 

т = р.Д,

[а = 0 , 1,

 

 

 

 

 

 

 

Очевидно, что при сделанных допущениях интеграл в формуле (4.1)

можно заменить

знаком суммы

и написать:

 

 

 

 

 

 

 

 

N

 

 

 

R СО = R (Ид) « 2АГ1-)- 1 S * № х + Р) Д1

(4-4)

 

 

 

 

 

 

v= —N

 

 

при

r = [iД, (1 =

0,

1. . . .

Если ввести обозначения:

 

 

 

 

 

/?([АД ) = Я ([А ),

 

 

 

 

 

 

х (уД) = jc„,

 

(4.5)

 

 

 

* 10»+

Р)Д] = *, +*.

 

 

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