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

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

.pdf
Скачиваний:
95
Добавлен:
20.10.2023
Размер:
8.65 Mб
Скачать

Г л а в а 7

СЛУЧАЙНЫЕ ПРОЦЕССЫ В АВТОМАТИЧЕСКИХ СИСТЕМАХ

§ 7.1. ПОСТАНОВКА ЗАДАЧИ

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

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

Воздействие может принимать с

течением времени

самые раз­

нообразные

случайные

значения.

К таким

воздействиям отно­

сят

различные помехи,

шумы, произвольное

изменение нагруз­

ки,

дрейфы в электронных усилителях и пр.

Эти

случайные ве­

личины можно оценить

только с точки зрения

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

появления в тот или

иной

момент

времени, в той

или

иной

форме.

 

системы

с

воздействиями

случайного

 

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

характера

может служить

автопилот (см.

рис

1.3), который

реагирует на различные изменения

атмосферных условий

(на­

пример, ветра), тяги, напряжения питания усилителей и

руле­

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

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

136

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

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

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

чем пренебречь этими случайными

величинами невозможно.

В этом случае требуется определить

величину дополнительных

отклонений системы от установившихся значений, возникающих

от помех' и шумов,

рассматриваемых как

случайные возму­

щения.

основа

статистических

вероятностных ме­

Математическая

тодов расчета разработана

А. Н. Колмогоровым, Н. Винером,

В. С. Пугачевым и др.

 

 

§ 7.2. СТАТИСТИЧЕСКИЕ ХАРАКТЕРИСТИКИ СЛУЧАЙНЫХ ПРОЦЕССОВ

Случайный процесс — есть случайная функция x[t) от неза­ висимой переменной t. Каждое испытание дает определенную функцию хДО. которая называется реализацией процесса.

Случайный процесс можно рассматривать либо как сово­ купность реализаций xt{t), либо как совокупность случайных величин X[j(t), т. е. выборочных значений случайной функции в дискретные моменты времени. В последнем случае достаточно полное описание случайного процесса дает описание соответ­ ствующего множества случайных величин.

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

137

случайного процесса. Подробное изложение понятии, связанных

с числовыми

характеристиками и

их свойствами, дается

в кур­

сах математики. В настоящем

параграфе

приводится

лишь

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

и функций, исполь­

зуемых для

определения

установившихся

ошибок

при воздей­

ствии на АС случайных

сигналов, а также

краткое

пояснение

их математического и физического смысла.

При этом за основу

приняты эргодическне стационарные случайные процессы.

1. Математическое ожидание случайной величины МЫ, ил статистическое среднее ее значение, является простейшей и в то

же время важнейшей числовой характеристикой

случайной

величины:

 

М (х ]= х = J xf(x)dx,

(7.1)

где f ( х ) — плотность распределения непрерывной случайной величины.

Величина х определяется как среднее значение для всей со­ вокупности случайных величин данного случайного процесса-

водин и тот же момент времени.

2.Среднее значение случайной величины но времени х:

1

(7.2)

л' = П т_2 7 г ) x(t)di.

У—оо

у.

 

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

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

x=F.

(7.3)

3. Среднеквадратическое значение случайной величины хск определяется как корень квадратный из среднего квадрата этой

случа й иой величины:

 

 

*ск= У х2.

(7.4)

4.

Дисперсия D (от латинского

слова — рассеяние) ест

средний квадрат отклонения случайной величины от ее среднего значения:

?

*

1

Г

(7.5)

D = J

(x—x)2f{x)d x= [im 1jjr^

(x—x)2d(.

138

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

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

a=rVD. (7.6)

Часто случайные процессы приводят к центрированным, т. е. к таким, для которых математическое ожидание равно нулю

(х = 0). В таких процессах колебания случайной величины про­ исходят как бы относительно ее нулевого значения. Для центри­ рованных процессов дисперсия равна среднему квадрату слу­ чайной величины.

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

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

которая характеризует' степень связи, т. е. корреляцию между значениями двух функций x\(t) и д-2 (/), отстоящими друг от друга по времени на величину т:

(7.8)

Величина Я\д (0) определяет эту связь в один и тот же мо­ мент времени.

7. Спектральная плотность S (со) по своему физическом смыслу есть величина, которая пропорциональна средней мощ­ ности процесса в интервале частот от со до co-f dm.

Хотя среднее значение случайной . величины, дисперсия и корреляционная функция достаточно полно характеризуют слу-

139

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

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

Математическое определение спектральной плотности слу­ чайной функции хцД) может быть представлено в следующем виде:

 

—г

 

Для случайной безграничной во времени функции лц

(/) не

применимо

преобразование Фурье, поэтому вместо

функ­

ции Х\ (t)

при определении спектральной плотности использует­

ся функция ли r(t), которая удовлетворяет условиям xir{t)=x1(t) при t , лежащем в пределах ±7\ и х, (t) =0 при |/|>7\

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

четной функцией часто­

ты и) (7.9). Спектральные плотности

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

одного какого-либо стационарного

эргодического случайного

процесса, подобно их корреляционным функциям, одинаковы.

8.

Используется также понятие взаимной спектральной пло

ности S|,2 (/ со)

для двух

коррелированных процессов, которая

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

 

 

 

 

 

 

S i,2 (/со) =

1im Yjr X iT{—/ш) X2 т(/<“)•

(7.10)

Взаимная спектральная плотность не является четной функ­

цией и

поэтому

имеет

вид не

вещественной, а

комплексной

функции.

/со на —/со и наоборот выражение для взаимной

При замене

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

меняется

на сопряженное,

т. е.

 

 

Si,2 (/<«)= S-2,1 ( /И).

(7.11)

Аналитическое определение спектральной плотности для эргодического стационарного процесса возможно через корре­ ляционную функцию этого процесса. Для такого процесса пре-

МО

образование Фурье от корреляционной функции есть спектраль­ ная плотность.

В соответствии

с этим

 

 

 

 

 

+ со

 

 

 

 

 

S(U>)=

j

R (-)e

; ,,IXс/т=2 J R (т) cos <» x dt\

(7.12)

 

 

 

 

о

 

 

 

5|.г(/“ )=

j

/?i,2 (t) e JraTс/т,

(7.13)

т. e. взаимная спектральная

плотность

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

Фурье взаимной корреляционной функции.

 

Если выполняются все условия существования

интеграла

Фурье, то имеют место и формулы обращения:

 

/?(T) = J _ +j

S(co)cE“Tc/co=^jS(co)coscoTc/u>,

(7.14)

 

—- со

 

 

U

 

 

/?,.2( т ) = ^ г j

S1,2(y<U) / ' "

rf<„.

(7.15)

В соответствии с формулами (7.12) и (7.14) в приложении 3 даны некоторые функции R(x) и их изображения Фурье S(co).

Важным свойством спектральной плотности является то, что интегрирование ее по всем частотам от — оо до +оо дает сред­ ний квадрат исходной функции времени:

(7Л6)

§ 7.3. СЛУЧАЙНЫЕ ПРОЦЕССЫ В ЛИНЕЙНОЙ АВТОМАТИЧЕСКОЙ СИСТЕМЕ

Рассмотрим линейную АС (рис. 7.1) с частотной характери­ стикой Ф (/со), находящуюся под воздействием на входе слу­ чайной функции x{t). Если это входное воздействие х (t) являет­ ся стационарным случайным процессом, то на выходе линей­ ной АС также будет стационарный случайный процесс y{t).

x(t)

W(p)

V ( t )

 

 

 

Рис. 7.1

 

141

Предположим, что рассматриваемая система устойчива. Частотную характеристику системы в этом случае можно опре­

делить как предел

отношения

функций

выхода YT(j<д) и вхо­

да Х-Л/о.)

п[1и tТ- >

, вычисляемых но формулам

 

 

 

 

 

 

 

 

Хт(/со) =

+'{'

-j.»i

 

 

 

(7.17)

 

 

 

 

\x{t)e

dt,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Yr(j(o)=

jisit

 

 

 

(7.18)

 

 

 

 

[ y(t)e

dt,

 

 

 

т. е.

 

 

 

 

 

 

(/“ )

 

 

 

 

 

 

 

 

 

ф 0‘») = 1i m X r (H '

 

 

 

(7.19)

Будем

считать,

что

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

входного

воз­

действия

Sj(qi) определена

одним из

методов

(например,

по

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

функции).

выходного

воздействия

находится

Спектральная

плотность

следующим образом:

 

 

 

 

 

 

 

 

 

 

 

 

 

И =

l i m

^

(-/со) YT(/«>).

 

 

(7.20)

Умножив и

разделив

выражение,

находящееся

под знаком

предела,

на величину Хт ( — /1|>) X т(/со)

и

сгруппировав

члены,

получим:

 

YT(—j^)Yr (/щ)

lim

 

 

 

 

 

 

S., (со)- lim

 

 

 

 

(7.21)

Хт{—/<и) Хг(/ш)

Г - <х> -ур Х т(—/со) Xr(j<»).

 

Т - с с

 

 

 

 

 

 

 

Первый сомножитель этого

выражения

может быть

пред­

ставлен как Ф{— /со) Ф (/со) = |Ф (/со) |2,

а второй как

S^co).

 

Отсюда окончательно:

 

 

 

 

 

 

 

 

 

 

 

 

S2 N =

Iф (/“ )|2S iH .

 

 

 

(7.22)

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

выходного

воздействия АС

равна

произведению спектральной плотности входного воздействия на квадрат модуля частотной характеристики АС.

Зная величину

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

выходного

воздей­

ствия АС,

в соответствии с выражением (7.16) можно

опреде­

лить среднее значение

квадрата выходного

воздействия АС по

формуле:

 

_

I

+~

 

 

 

 

 

 

Из двух последних

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

среднее

значение

квадрата

выходного воздействия

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

плотности

входного

воздействия:

 

 

 

7 ( 0 = - ^ - J

^ (H P S .H c /c o .

(7.24)

142

Формула (7.24) дает возможность определить характер изменения выходного воздействия АС, при случайном входном воздействии, если известны частотная характеристика системы н спектральная плотность входного воздействия.

Взаимная спектральная плотность входного и выходного воздействий линейной АС по определению равна

 

Si,2

(/со) = lim

L _ A r (—/ю)УУ(/<а).

 

 

 

2Т

 

 

 

Умножая п деля

выражение,

стоящее под

знаком

предела,

па Аг(/<°) и группируя члены, получим:

 

 

■S1.2 ( Н -

 

lim

Х г ( - /ш ) Х 7-(/ш).

(7.25)

 

Г- ООА т1

/

X.1

 

 

Заменяя

сомножители в

формуле (7.25)

их значениями,

окончательно

имеем

 

 

 

 

 

 

*Si,2 (уш)= ф (/to) S, (со).

 

(7.26)

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

Выражение

(7.26)

позволяет

по взаимной

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

плотности

входного

и

выходного

воздействий

и

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

плотности

входного

воздействия

найти

частотную

характери­

стику АС.

 

 

 

позволяет

найти

связь

между соответ­

Это же выражение

ствующими корреляционными функциями.

 

 

 

Если переходная характеристика системы представлена

функцией h(t),

причем

Л(0)=0,

а импульсная

 

переходная ха-

 

 

I

 

..

,,, dh(t)

и известна

корреля­

рактеристика — функцией

g(t)=

^

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

R\»W = [ g ( t ) R i b - t ) d t ‘

(7.27)

о

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

143

§ 7.4. УСТАНОВИВШИЕСЯ ОШИБКИ СИСТЕМЫ ПРИ ВОЗДЕЙСТВИИ СЛУЧАЙНЫХ СИГНАЛОВ

Рассмотрим замкнутую АС (рис. 7.2). к входу которой при­ ложены случайный полезный сигнал (планируемая составлн- I! тая входного воздействия) х п (7) и случайная помеха (непла-

нпруеман составляющая входного воздействия) д'н (().

Необходимо, чтобы выходное воздействие у (/)

наиболее точ­

но соответствовало планируемой составляющей

входного воз­

действия хп (t).

В этом случае ошибка системы

(7.28)

 

 

s (t)=xn (t)— y{t),

а входное воздействие

 

При воздействии

*(0 = *п (*)+*н(0-

(7.29)

на АС случайных сигналов дс,,(^) и хн (t)

ошибка системы

е (0

все время меняется. Установившееся зна­

чение ошибки отсутствует. Поэтому определяют не саму ошибку

системы или ее установившееся

значение, а величину

среднего

значения

квадрата

ошибки е2((). Для определения

е2 (t) тре­

буется знать величину

спектральной плотности сигнала ошиб­

ки

(о>)

(7.16),

которая может быть определена как

 

 

 

 

 

(ш) =

Нш

Ег ( — / со) Е 7.(/со),

(7.30)

где

 

 

 

 

 

dt,

 

 

 

е 7.(-(«>)= У

 

}и>Т

 

 

 

в (/) е

dt.

 

 

 

 

 

 

 

 

 

 

В соответствии с определением ошибки системы (7.28) имеем

 

 

 

 

 

Ег (/ш)=Хпт (/со) — YT(/со),

(7.31)

где

А''пт (/со) =

+ Т

 

— ]ш Т

 

 

j

хп (/) с?

dt,

 

 

 

 

 

—Т

 

 

 

 

 

 

 

_ . шГ

 

 

 

Yп (/св)=

j

y(t)e

1Ш dt.

 

 

144

Аналогичным образом преобразуется Ё?.(— /<■>). Учитывая уравнения (7.30) и (7.31), можно записать:

5= (ш) = П т -^р= г |Х пт (— /со) — У г ( — / ш)] [ У п т (/со)— К г (/ю )] —

сс£ I

=Ит-^г ( |X пт (/со) 2+|Кг(/со)|2—А’пт(—/«>) У Г ( П ~

- Х Пт(/о))Ут-(/со)1.

. (7.32)

Преобразуем далее отдельные слагаемые выражения (7.32), находящиеся под знаком предела:

а) в соответствии с уравнением (7.9) первое и второе слага­ емые будут иметь вид

Нш

2^г|Хпт(/‘о)Г =5п(со) и

Пш - ^ г |K7 (/co)|2= S 4(co),

(7.33)

причем второе слагаемое на основании выражения (7.32)

может

быть

приведено

к виду

 

 

 

 

 

 

 

 

 

 

 

 

S 2 H = S 1(со) |ф (/со) 2.

 

 

 

 

 

 

 

 

 

 

 

 

У т(/<*>)

 

 

 

б)

умножая и деля третье слагаемое на Хт(/«) ’

группируя

члены и преобразуя

в

соответствии

с

уравнениями

(7.19) и

(7.9),

получим

 

 

 

 

 

 

 

 

 

 

Игл -Д р гХпт (—И У т ( П = ; l i m

 

lim

p L

Xm (—/со) Хг(/со)=

\]№)

Нш

 

 

 

 

У’-r

 

 

 

 

/'-* оо А 7

I -► со

^ i

 

 

 

 

 

 

 

 

 

= 0(/io)Sm(/co);

 

 

 

 

(7.34)

в)

 

после

аналогичных преобразований

четвертое

 

слагаемо

будет иметь вид

 

 

 

 

 

 

____

 

Игл

А'пт (Н У г (—/со)=ф(—/со) Sm (/со).

 

 

(7.35)

Подставив

значения

отдельных

слагаемых,

полученные

после преобразований, в выражение (7.32),

найдем

использу­

емую

при

расчетах формулу

 

 

 

 

 

 

 

 

 

 

Sc (co)=Sn И + |Ф (/<») |2

(со)

 

 

 

 

 

 

 

Ф (/со) Sm (/u>)—Ф (—/со) Sin (/<“)•

 

 

(7.36)

Имея

значение

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

плотности

сигнала

ошибки,

можно получить среднее значение квадрата случайной

ошибки

в соответствии с выражением (7.16) по формуле

 

 

 

 

 

 

 

Г2( ^ ) = ~ +|

S«(co) cfco.

 

 

 

(7.37)

10 Учебник

145

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