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

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

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

1 9 0 ПРОХОЖДЕНИЕ СЛУЧАЙНОГО СИГНАЛА ЧЕРЕЗ ДИНАМИЧ. СИСТЕМУ [ГЛ. V

Но

 

г

 

 

 

 

lim

f т (t -f- т — т)) т (t — X) d£ =

 

Т + 00 1 1

 

 

 

 

 

 

 

 

 

Т

 

 

= lim J=r

 

f

m (t-\-T -it- \ — Tl) m ( t) d t = Rm( x ^ - \ — yi).

 

T +oo

y T

 

 

Поэтому

вместо (5.7)

можно написать:

 

 

 

 

СО

с о

 

 

Я , С 0 =

/

d\ f k (7|) k (X) (т + X- 7J) rfTj.

(5.8)

 

— CO

— CO

 

Формула (5.8) представляет собой искомое соотношение между корреляционными функциями входного m(t) и выходного x[t) сиг­ налов.

Найдем теперь соотношение между спектральными плотностями

входного

и выходного

сигналов.

Имеем:

 

 

 

 

 

 

 

СО

 

 

 

 

 

Sx (ш )=

f

Rx (T.)e-i™dx.

(5.9)

Подставляя (5.8) в (5.9),

найдем:

 

 

с о

СО

 

о о

 

 

 

 

Sx iш ) =

f dr J dX J e - ^ f t ( 7 | ) * ( X ) / ? w ( x 4 - X - 7 ]) d 7 , =

 

— CO

—со

— с о

 

 

 

 

 

о о

с о

 

CO

 

 

 

 

=

J" dt

J

dX

J

e~Jm(T+^—1л) g/ш). g-/u)T| X

 

— CO

— CO

— ГХ1

 

 

 

 

 

 

 

 

 

 

 

X Rm(* +

X— 7)) ft (X) ft (7)) dT) =

 

C O

 

 

 

C O

 

 

 

— J k (у) е~/шг>d7)

J

k (X) e^u)) dX X

 

 

— CO

 

 

 

— 0 0

 

 

 

 

 

 

 

 

 

 

CO

 

XJ ” Ят (т -|- X — tj) e-M*+b-i))dT =

o o

CO

= ф ( » Ф ( — M / W e - /W dT = I Ф (уш) I* • Sm (u>). (5.10)

Выражение (5.10) показывает, что спектральная плотность Sx (io) величины на выходе линейной динамической системы равна произ­ ведению квадрата модуля ее передаточной функции Ф(уш) на спек­ тральную плотность 5 т (ш) величины на входе.

3]

ОБЩЕЕ

БЫРАЖЕНИЕ

ДЛЯ

СПЕКТРАЛЬНОЙ ПЛОТНОСТИ ОШИБКИ

191

3.

Общее

выражение

для

спектральной плотности ошибки })

Рассмотрим теперь

более

общий случай, когда линейная

дина"

мическая

система

находится

под

влиянием двух стационарных слу­

чайных воздействий: управляющего

т (() (полезный сигнал) и возму­

щающего

п (t) (помеха),

приложенных к

различным ее точкам.

 

Согласно

(1.93)

имеем:

 

 

СО

 

 

 

 

 

00

 

 

 

 

 

 

 

X (/) =

J* т (f — X) ft (X) d\ 4 -

[ ft (t — X) l (X) rfX.

(5.11)

 

 

 

 

о

 

 

 

 

()

 

Подставляя (5.11) в выражение

(3.12) для корреляционной функ­

ции

Rx (i) переменной х (t), получим:

 

 

 

 

 

j"

Т

 

со

 

 

со

 

Rx {x) =

lim

 

f dt

 

f

 

т — ti) k (tj) dr\ f m (t— X)6(X)dX4~

 

 

 

 

 

T

i

 

 

о

 

 

 

 

 

 

со

 

 

со

 

+/ rf/Jn(/-f-x— i])/(i|)dij Jn(< — X)/(X)dX +

 

 

 

 

- T O

 

0

 

 

 

 

 

 

T

со

 

oo

 

 

 

 

 

4*

f d

t f

m{t+

i — 7|) k (-ф dfj j n (t — X) l (X) dX -|-

 

 

 

 

- r

 

о

 

о

 

 

 

 

 

 

T

 

со

oo

1

 

 

 

 

4 -

f

dt f

m (t X) к (X) d\ j n (t + 1 — ц) l (д) d-ц 1

и,

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

-

Г О

 

0

J

 

 

 

 

 

 

Rx (t) =

 

 

 

 

 

 

 

 

 

CO

CO

 

 

 

 

 

 

 

=

/ л

/{л(Х)/гт (т+ х-71)Л(11) + /(Х)/гя(х + х - 1,)/(1о+

 

--С О

— с о

 

 

 

 

 

 

 

 

+ LO') Rmn(-C

 

X— T|) k (7]) -f - k (X) Rnm(t + X - 1]) / (I))} rfT).

(5.12)

 

Перейдем к определению спектральной плотности Sx (со).

 

 

Умножая обе

части

(5.12) на

e~iwc и интегрируя по т от — оо

до + со, при помощи выкладок,

аналогичных тем, которые были

применены при выводе формулы (5.10) в конце предыдущего

пара­

графа, получим:

 

 

 

 

 

 

 

 

(®) =

I Ф О )

I2

(Ш) ■+ I V (Jio) Р S„ (со)

 

 

 

+

ф * ( М Srnn(со) Y ( » 4 - Ф (/ш) S am(cu) Y* С/т).

(5.13)

где

Y (/’со) — передаточная функция

сигнала помехи.

 

 

4 См. цитированную книгу X. Джеймса и др.

 

192

ПРОХОЖДЕНИЕ

СЛУЧАЙНОГО

СИГНАЛА ЧЕРЕЗ ДИНАМИЧ. СИСТЕМУ

[ГЛ. V

В частном,

 

но весьма часто встречающемся случае, когда

 

 

 

 

 

 

 

Sm»

=

$„„(<“) =

0,

 

(5.14)

т. е.

когда управляющее и возмущающее

воздействие

взаимно

неза­

висимы

(корреляция

между

ними отсутствует), формула (5.13) упро­

щается

и принимает

вид

 

 

 

 

 

 

 

 

 

5 , ((о) = | Ф ( »

|* 5„, (ш) +

1К (у'ш) р S„ (ш).

(5.15)

Найдем теперь выражение для спектральной плотности 55(ш) ве­

личины

ошибки

воспроизведения

е (/). .■

 

 

 

 

Для этого

в выражении

(5.13) необходимо лишь заменить Ф (у'ш)

через

ФЕ(у'ш),

имея

в виду,

что

согласно

(1.190)

 

 

Итак,

 

 

 

 

 

Ф Д » = 1 — Ф О ) .

 

(5.16)

 

 

 

 

 

 

 

 

 

 

 

 

 

(«)) =

| Фе ( »

 

I2 S m(ш) + 1К (у'ш) р S B(ш) +

 

 

 

 

+

 

«Г ( » S mn (ш) Y (у'ш) + Ф, (у'ш) S nm(ш) V* (у'ш).

(5.17)

В частном случае

(5.14) выражение (5.17)

сводится к

виду

 

 

 

5. («О =

I Ф. (

»

I2 S m(ш) 4 - 1Y (у'ш) р S a (ш).

 

(5.18)

П р и м е р .

 

С п е к т р а л ь н а я п л о т н о с т ь о ш и б к и с л е д я ­

щей

с и с т е м ы .

Рассмотрим

следящую

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

антенной радиолокатора *),

в

котором

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

слежение за

целью осуществляется при помощи конической развертки луча, схе­ матически изображенной на рис. 5.2. Передатчик радиолокатора из­ лучает импульсы длительностью порядка долей микросекунды. Ис­ точником излучения служит вибратор, установленный вблизи фокуса параболического рефлектора и вращающийся от мотора со скоро­ стью 1800 об мин. При вращении ось луча отклоняется от оси па­ раболоида на 1,25° и описывает в пространстве коническую поверх­ ность (рис. 5.2).

Когда наблюдаемая цель находится на оси конуса (точка А на рис. 5.2), то условия облучения цели остаются неизменными, не­

смотря на вращение луча. Соответственно

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

этой

цели импульсы имеют постоянную амплитуду. Если

же цель

нахо­

дится не

на оси конуса, например в точке

В, то максимальная ве­

личина

отраженных

импульсов будет

в точке 3,

минимум

будет

в точке

/;

в точках

2 н 4 величина

отраженных импульсов

будет

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

') Приводимое ниже описание принципа действия радиолокатора взято из книги «Наземные американские и английские радиолокационные стан­ ции», Воениздат, Москва, 1947.

3]

ОБЩЕЕ ВЫРАЖЕНИЕ ДЛЯ СПЕКТРАЛЬНОЙ ПЛОТНОСТИ ОШИБКИ

193

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

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

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

ось рефлектора,

причем

антенна вновь

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

цель так,

чтобы последняя

попала

в ось конуса.

Ошибка слежения

за целью

как по углу места, так

и по азимуту не превышает в обычных ус­

ловиях 4 угловых минут.

 

 

 

 

Огибающая отраженных импульсов

выделяется

в приемнике при

помощи детектора и соответствующего фильтра.

Эта огибающая

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

 

 

 

 

Т I1 +e(O sin(u)/-f-tp)],

 

(5.19)

где f — величина, пропорциональная коэффициенту отражения цели, е(£) — угловая ошибка слежения, ojj — угловая частота, с кото­ рой происходит развертка луча, ср — фазовый угол, определяющий31

13 Зак. 1083. В. В. Солодовников 4 4

1 9 4

ПРОХОЖДЕНИЕ СЛУЧАЙНОГО СИГНАЛА ЧЕРЕЗ ДИНАМИЧ. СИСТЕМУ [ГЛ. V

положение цели относительно оси рефлектора

(случай, когда ср — О,

соответствует

наличию только азимутальной ошибки).

 

Для упрощения рассуждений предположим, что цель движется

по

прямой,

проходящей непосредственно

над радиолокатором.

Очевидно, что в этом случае угол азимута цели сохраняется по­ стоянным.

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

исходит лишь в зависимости от

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

оси конуса, описываемого лучом,

т. е. что коэффициент отражения

цели 7 является постоянным. Однако на самом деле это не так. Изменения коэффициента отражения 7 вызываются, как это ука­

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

функцию

времени.

 

 

 

 

 

Можно написать:

 

 

 

 

 

 

 

T = To +

Ti(0 =

T o[l-|-«(01 .

(5-20)

где 7о =

const

и

п (() =

_

 

 

ю

 

 

 

 

 

 

Типичный

вид

функции

n(f),

полученный при помощи

экспери­

ментальной записи сигнала, отраженного от летящего самолета1), приведен на рис. 5.3.

Учитывая

(5.20),

огибающая

(5.19) может быть представлена

с точностью до постоянной 7„ в следующем виде:

 

[1+ п (0][1+ е(081 п (а)^Н -« р )] =

 

 

 

=

1 +

» (0 ■+■ е (0 sin К * +

<Р) 4 - и (0 е (0 sin (ш/ +

ср). (5.21)

Так

как

согласно

эксперименту | е ( 0 1

0,05 J) и

стандартное

значение

величины n(t) (предполагаемой

стационарной

случайной

*) См. цитированную на стр. 189 книгу X. Джеймса и др.

3]

ОБЩЕЕ

ВЫРАЖЕНИЕ

ДЛЯ

СПЕКТРАЛЬНОЙ ПЛОТНОСТИ

ОШИБКИ

 

195

функцией) имеет порядок 0,25, то

последним членом (5.21)

можно

пренебречь

и для

огибающей

написать:

 

 

 

 

 

 

 

 

 

 

l + « ( 0

+

e(0sm (o)/ + cp).

 

(5.22)

Сигнал (5.22) пропускается через разделительный трансформатор,

служащий для отделения постоянной составляющей.

 

 

 

Таким

образом,

напряжение на зажимах вторичной обмотки транс­

форматора

 

пропорционально

величине

 

 

 

 

 

 

 

 

 

l2 (t) — а [п (t) -(- s (t) sin (ws( —j—tp)].

 

 

 

Напряжение l2{t) поступает на коммутатор, который дает на выходе

две

составляющие,

одна

из которых пропорциональна 21г (t) sin u>st,

а другая

пропорциональна

 

 

 

 

 

 

 

2/2 (0 cos о / . Первый сигнал

дщ +^пи )^

хи>

УстройствоупраВ-

поступает

 

на

вход

канала

 

 

l+n(t)+e!t)sinu)sL

пениаантенной

 

управления

а

антенной

по Разделитель

 

 

 

азимуту,

 

второй

 

На

ньш транс

 

 

 

 

 

вход

канала

управления Форматор

n lt)+ £ .(l)sin u js t

 

 

 

антенной

 

по

углу

 

места.

Комму­

 

 

 

 

 

Так как оба канала обычно

татор

2n(t)sin iost+e(t)-e.(t)cos2(itst

 

одинаковы, то в дальней­

 

 

 

шем мы ограничимся рас­

Фильтр

 

 

 

 

 

смотрением лишь

одного

из

 

 

2nIt)sinajst+Eftj

 

 

 

них, например азимуталь­

 

 

 

 

 

 

 

 

 

 

 

 

ного, и для простоты пред­

 

 

 

Рис. 5.4.

 

 

 

положим,

что ® = 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Итак, на входе этого

канала мы будем иметь сигнал, пропорцио­

нальный

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(/) -}- е (t) sin ш/] 2 sin iost =

2ti (t) sin ш /

s (0 — e (t) cos 2wst.

(5.23)

Коммутированный сигнал (5.23) пропускается затем через фильтр

нижних частот,

который

устраняет

пульсации сигнала (5.23)

вида

е (t) cos 2ш/, имеющие

частоту, равную

удвоенной частоте развертки

луча. Сигнал e(t) проходит через

фильтр без изменений, так

как

его

спектр

содержит

лишь

низкие

частоты (обычно

меньше

2

гц).

Составляющая 2n(£)sinu>^ изменяется фильтром в соответствии с ви­ дом его передаточной функции (обычно передаточная функция фильтра выбирается таким образом, что он срезает все частоты выше 10 гц).

Рассмотренный выше ход преобразования сигнала схематически изображен на рис. 5.4.

Схема на рис. 5.4 еще более упрощенно может быть представ­ лена так, как это сделано на рис. 5.5. Здесь предположено, что на

вход

фильтра

действует сигнал

 

 

 

 

2n (t) sin

s (t),

(5.24)

а не

сигнал

(5.23), поскольку, как

это мы видели,.составляющая

13*

1 9 6

ПРОХОЖДЕНИЕ

СЛУЧАЙНОГО

СИГНАЛА ЧЕРЕЗ ДИНАМИЧ. СИСТЕМУ [ГЛ. V

e(t)cos2ust

через

фильтр не

проходит и может не приниматься

во

внимание.

 

 

 

 

Найдем корреляционную функцию /?п(-с) и спектральную плот­

ность

S n (x)

флуктуаций n(t).

 

 

Корреляционная функция Rn (х) может быть найдена из экспери­

ментальной

кривой на рис.

5.3, как это изложено в § 12 гл.

III.

В

результате

обработки

мы получим кривую, показанную

на

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

Rn М

Р« W Ял (0) '

Для определения соответствую-

gffj+r+nftl x!t) Канал

управления

антенной

2n(tlsiniost-He!t)

фильтрI

Рис. 5.5.

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

00

ап(со) = 2 J р„ (т) cos сот dx.

о

Этот интеграл может быть вычислен графически, например методом,

изложенным в § 12 гл. I,

или

аналитически,

если

известен вид

функции р„(т).

 

 

 

 

 

Кривая р„ (т) на рис. 5.6

приближенно может быть

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

в виде

 

 

 

 

 

рп (х) яь е-0,831

[cos 2х —(—0,106sin2|x|],

(5.25)

График функции (5.25) изображен

на рис. 5.6

пунктиром. Под­

ставив (5.25) в выражение для ап (со),

получим:

 

 

(ш)

1,44

—[—11,11

 

(5.26)

о)4 — 6,27 о)2-|-23,64

 

 

 

 

График функции (5.26), являющейся нормированной спектраль­ ной плотностью ол (со) флуктуации n(t), приведен на рис. 5.7.

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

пх (в) = 2n (t) sin со^

сигнала (5.24).

3]

 

ОБЩЕЕ

ВЫРАЖЕНИЕ ДЛЯ СПЕКТРАЛЬНОЙ ПЛОТНОСТИ ОШИБКИ

197

 

По

определению

 

 

 

 

 

 

 

 

 

 

Rni (т) =

4п (t) n{t-\-z ) sin w j sin cos (t + x).

(5-27)

Так

как

n(t)

не

зависит от положения

луча при

развертке,

то сиг­

налы

2«(£)sinco^,

отличающиеся

лишь

фазой, одинаково вероятны.

Поэтому

запишем (5.27) в виде

 

 

 

 

 

 

 

Rn, СО =

4п (t)n(t-\- х) sin (со^ -|-ср) sin [со4

—)—х) —)—ср].

(5.28)

 

Усредняя в (5.28) по всем возможным значениям фазы ср, получим:

 

 

 

 

2тс

 

 

 

 

 

 

 

 

 

 

^ f

sin (wsf +

ср) sin [cos (t -(- x) -f- cp] dcp =

j cos coyt

 

 

 

 

 

0

 

 

 

 

 

 

 

 

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

 

Rn, (0 = 2 cos wsxRn (x).

 

(5.29)

 

 

 

 

 

 

 

 

 

 

Спектральная плотность S„, (со):

 

 

 

 

Sn, (u>) =

 

2 J

R„t (x) cos cox dx = 4 J Rn (x) cos co^x cos cox di =

 

 

 

 

 

о

 

 

 

о

 

 

 

 

 

 

 

 

с о

 

 

 

 

o o

 

 

 

 

=

 

2 J

Rn (x) cos (со — coj x dx -f- 2J

Rn (x) cos (co -(- “ 0 x dx =

 

 

 

 

О

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

— S a (to — “ s) Ч -5 Я(со -f- ®s).

(5.30)

 

Нормированная

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

плотность

коммутированного

сигнала,

 

вычисленная

по

 

 

 

 

формуле (5.30),

изображена

 

 

 

 

на

рис.

5.8.

к

 

вычисле­

 

 

 

 

 

Перейдем

 

 

 

 

 

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

 

 

 

 

ошибки е (/).

 

 

 

 

 

 

 

 

 

Ввиду

отсутствия

кор­

 

 

 

 

реляции

между траекторией

 

 

 

 

движения цели и флуктуа­

 

 

 

 

циями ti(t) взаимные корре­

 

 

 

 

ляционные функции в (5.17)

 

 

 

 

равны

нулю,

и

выражение

 

 

 

jfo

для

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

 

 

 

 

сводится к виду (5.18). Кроме того, необходимо учесть,

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

Ф(/ш) = Y ( » ,

иравенство (5.18) принимает вид

а. И = I ®. О ) I2 s m( « Ж ф I2 Sa(со).

198- ПРОХОЖДЕНИЕ СЛУЧАЙНОГО СИГНАЛА ЧЕРЕЗ ДИНАМИЧ. СИСТЕМУ [ГЛ. V

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

S m(ш) — 28 (ш),

и поскольку функция ФЕ(уш), как это будет ясно из дальнейшего, содержит ш в виде множителя 1), то

I Ф6 (Уш) |2 S m(ш) = 0.

Следовательно, мы будем иметь:

51(ш) = |ФС/«о) |25„, (ш),

(5.31)

и ошибка будет определяться исключительно флуктуациями n(t).

Из (5.31) ясно, что для определения спектральной плотности

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

 

Ф (ус)

К Щ М

(5.32)

1 + KW

(у'ш)

[см. формулу (1.174)1.

 

 

функцию следящей

системы

Для того чтобы найти передаточную

в разомкнутом состоянии KW (у'ш), входящую в (5.32), рассмотрим

упрощенную структурную схему, изображенную на рис. 5.9.

 

Как мы видим, она состоит из блока

1, служащего для

выделе­

ния. сигнала ошибки, фильтра 2, усилителя 3, электромашинного усилителя (управляющая обмотка 4, якорь 5), приводного двигателя 6 , редуктора 7, мостика обратной связи 8 и фильтра обратной связи 9.

*) Это ясно также из того обстоятельства, что рассматриваемая следя* щая система является астатической.

3]

ОБЩЕЕ ВЫРАЖЕНИЕ ДЛЯ СПЕКТРАЛЬНОЙ ПЛОТНОСТИ ОШИБКИ

199

Согласно эксперименту ‘) передаточная функция фильтра 2 может быть представлена в виде

М«а(01

 

kd

(5.33)

L[eAt)}

~

(7> + l)(7Vs + l)a’

 

где е ,(0 — напряжение на

входе фильтра, пропорциональное ошибке;

e2(t) — напряжение на выходе фильтра; k d — постоянный коэффициент; Td, Тт— постоянные времени.

Обозначая через — изменение тока в управляющей обмотке электромашинного усилителя (или сокращенно .ЭМУ), мы можем на­ писать:

L[lc (t)\ = 1 ^ ^ L \ e 2 { t ) \ - T ^ L { e J(t)\,

(5.34) •

где ef (t) — напряжение на выходе цепи обратной связи; kf , kT— по­ стоянные коэффициенты; Тс — постоянная времени управляющей об­ мотки.

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

идругих подвижных частях, то передаточная функция ЭМУ, двигателя

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

&S

L[* (01

N

Гб.35)

£ ['е(0]“

« ( ^ s + lH ^ s +

l) *

*) См. цитированную книгу X. Джеймса и др.

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