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

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

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

3 8 4

СИНТЕЗ СЛЕДЯЩИХ

СИСТЕМ

ПРИ СЛУЧАЙНЫХ

ВОЗДЕЙСТВИЯХ [ГЛ.

IX

Эти

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

можно

найти также,

построив вначале

по

заданной передаточной функции <D(s) частотные характеристики зам­ кнутой системы, а затем применив соответствующие формулы или номограммы 1), позволяющие найти по этим характеристикам искомые характеристики разомкнутой системы.

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

даточными функциями.

Поэтому следующий шаг должен состоять

в том, чтобы, исходя •

из найденных оптимальных характеристик

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

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

В общем случае необходимо поступать следующим образом. В результате вычитания ординат характеристик объекта из оптималь­ ных характеристик системы найдем оптимальные логарифмические ча­ стотные характеристики последовательного корректирующего устрой­ ства. Затем, имея в виду, что для удовлетворения предъявляемых к системе требований особенно существенную роль играет вид частот­ ных характеристик разомкнутой системы при частотах, меньших частоты среза 2), и что вид частотных характеристик при тех часто­ тах, при которых Z,(cu)^(— 15)-г-(—20) дб, не существенен, заменим найденные оптимальные характеристики желаемыми.

Желаемые характеристики будут отличаться от оптимальных в основном тем, что при частотах, ббльших частоты среза шс, послед­

!) См., например, Основы автоматического регулирования. Теория. Под

редакцией

В. В. Солодовникова, гл. VI, Машгиз, 1954,

или C h e s t n u t ,

M a y e r ,

Servomechanisms and regulating system disign,

т.

I, J. Wiley, 1951.

Имеется

русский перевод: Ч е с т н а т Г. и М а й е р

Р.,

Проектирование

ирасчет следящих систем и систем регулирования, ч. 1, Госэнергоиздат, 1959.

2)Частотой среза называется координата точки пересечения амплитуд­ ной логарифмической характеристики с осью частот, см., например, Основы, автоматического регулирования. Теория. Под редакцией В. В. Солодовникова, Машгиз, 1954.

21 ЛОГАРИФМИЧЕСКИЕ ХАРАКТЕРИСТИКИ СИСТЕМЫ 3 8 5

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

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

вую частотные характеристики корректирующего устройства.

 

Поясним это на примере. Предположим,

что

требуется

найти

передаточную

функцию

разомкнутой следящей

системы с астатизмом

первого порядка в виде

дробно-рациональной

функции

таким

обра­

зом,

чтобы при коэффициентах ошибки Ct = 0,003,

С., =

0,002,

вре­

мени

переходного процесса

Г да 1 сек

и перерегулировании о ^ 4 0 %

обеспечивалась

минимально

возможная

среднеквадратическая ошибка

при наличии помех в виде белого шума с уровнем /V2= 0,6 - 10— Согласно изложенной выше методике импульсную переходную

функцию необходимо искать в следующем виде:

к(1)=^А0- \- А ^ + ^ ,

0 < г < Г .

(9.5)

Для этого примера, пользуясь таблицей 8,1 (случай II), получим Л0= 8,652; Л1 = —34,104; А-, = 28,200; eck = 8 град.

Передаточная функция замкнутой системы будет:

Ф (s) = ^ (1 - в - ST).+ А (1 _ 1 ,в-*г._ Те-ST) +

Импульсной переходной функции (9.5) соответствуют вещественная и мнимая частотные характеристики

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

мая

(пунктир) вещественные частотные характеристики приведены-

на

рис. 9.3. Из рис. 9.3

видно,

что оптимальные и желаемые харак­

теристики разомкнутой

системы

достаточно близки.

 

На рис. 9.4 изображен переходной процесс (сплошная линия),

соответствующий (9.5)

и

(9.6).

Перерегулирование в случае опт и-

25 Зак. I0S3. В. В. Солодоаннков

3 8 6

СИНТЕЗ СЛЕДЯЩИХ

СИСТЕМ

ПРИ СЛУЧАЙНЫХ ВОЗДЕЙСТВИЯХ [ г л . IX

мальной

системы можно найти

путем вычисления корней AQ+- A yt +

—)—у42^2

0, получающихся в

результате

приравнивания выражения

для k (t) нулю и подстановки

меньшего

из двух корней

 

(ь2 =

(— A l ± V А 1 — 4А0А*)/2Л1

в формулу для переходного процесса

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

Желаемая передаточная функция разомкнутой системы имеет вид

KW (s) =

125 (0,37s+l) (0,22s+ 1 )

(9.7)

s (1,21s*+ 0,704s + 1) (0,02s + I)'

Она получена путем замены оптимальных логарифмических частотных характеристик сопрягающимися линейными отрезками.

На рис. 9.4 приведен желаемый переходной процесс (пунктиром).

Для

передаточной

функции (9.7) о = 37%; Г = 1

сек

и С1 — 0,008;

Со =

0,002; есй =

8,5 град.

 

 

 

 

 

 

Сравним частотные характеристики существующей следящей си­

стемы

для управления антенной *) с

оптимальными,

вычисленными

так

же, как

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

выше

примере,

но

для

следующих

данных: Ct =

0,005; С2 = 0,0015;

Т =

1 сек;

(V2

= 0,6-

К Г6.

 

’-)

Основы

автоматического регулирования.

Теория.

Под редакцией

В. В.

С о л о д о в и и к о в а, Машгиз,

1954.

 

 

 

 

 

2]

ЛОГАРИФМИЧЕСКИЕ ХАРАКТЕРИСТИКИ СИСТЕМЫ

387

x(t)

Рис. 9.4.

25*

3 8 8

СИНТЕЗ СЛЕДЯЩИХ СИСТЕМ ПРИ СЛУЧАЙНЫХ ВОЗДЕЙСТВИЯХ [ г л . IX

Рис. 9.6

3]

 

ПЕРЕДАТОЧНАЯ

ФУНКЦИЯ КОРРЕКТИРУЮЩЕГО УСТРОЙСТВА

3 8 9

 

Численные

значения

величии Ct, С2, Т при расчете оптимальных

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

взяты

близкими к значениям этих величин для следя­

щей системы радиолокационной стан­ x(t)

 

ции

SC/?-584.

расчета

 

оптимальных

 

 

Методика

 

 

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

будет,

очевидно, та­

 

кой же, как и в предыдущем при­

 

мере,

поэтому ниже

приводятся

 

только

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

результаты."

 

 

На рис. 9.5, 9.6 и 9.7 приве­

 

дены для

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

 

желаемые (пунктир) характеристики,

 

а

в

таблице

9 .1 — динамические

 

показатели

системы.

 

 

 

 

После того как выбраны же­

 

лаемые

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

корректи­

 

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

 

 

Т а б л и ц а 9.1

 

5Ctf= SS4

Оптимальная система

C,

0,005

0,005

Ct

0,0016

0,0015

T

1,1 сек

1,0 сек

eck

8,5

9,4

a

35 96

37%

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

3. Определение аналитического выражения для передаточной функции корректирующего устройства

и численных значений ее коэффициентов

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

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

одной из заданных

частотных характеристик

(например,

логарифми-

>) L i n v i 11

J.,

The approximation with

rational

functions

of

prescribed

madnitude and

phase characteristics, Proc. of

the

IRE,

v. 40, №

6,

June 1952;

стр. 711—721.

 

 

 

 

 

 

 

390 СИНТЕЗ СЛЕДЯЩИХ СИСТЕМ ПРИ СЛУЧАЙНЫХ ВОЗДЕЙСТВИЯХ [ г л . IX

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

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

эффекта

передвижения этих нулей

и полюсов на

комплексной пло­

скости,

может

оказаться

очень сложным. В этом

случае указанные

поправки

можно

найти не путем подбора, а аналитическим методом.

Нужно,

однако,

заметить,

что при

расчете

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

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

решения

с требуемой точностью

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

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

рифмическим

амплитудной

и фазовой частотным характеристикам,

заданным

графически, и что в качестве первого

грубого приближе­

ния для

этой

передаточной

функции можно выбрать выражение

 

 

К »

(x5 + D ( T adS3 + 2 C rfTd5 +

l )

 

 

K'W(s) =

 

(9.8)

(rs+l)(7ls3+ 2C*r*s + l)'

Задача состоит в том, чтобы уточнить первое приближение соот­ ветствующим выбором величин Т, Тk, i d. С* и Crf.

Преобразуем выражение (9.8) к следующему виду:

K*W (s) =

К, (s — т) (* — 1d) —

(9.9)

 

S (S -

A) (S — XA) ( s — a; )

 

где

 

 

 

 

T — '

 

 

 

 

_

+ y .

v

 

 

x=-

 

 

 

(9.10)

 

T

 

 

**=■ Tk

^ J

T,

'■°-k

 

а через -j*. X* обозначены величины, комплексно-сопряженные с ве­ личинами Хл. Окончательно, пользуясь принятыми обозначениями,

3]

 

ПЕРЕДАТОЧНАЯ

ФУНКЦИЯ КОРРЕКТИРУЮЩЕГО УСТРОЙСТВА

391

нужно

написать:

 

 

 

 

 

 

 

 

K*W (s) =

M

i +

1) ( c

) ( l

 

(9.11)

 

 

 

 

 

 

 

s ( ^ + l ) (pt - I t - ' . K pt - f r + ' ) '

 

Если

исходить

из

выражения

(9.9), то легко

видеть, что

перво*

начально

сформулированная

задача состоит в том,

чтобы путем изме­

нения

положения

нулей f,

-j*

и полюсов X, Хк, X* улучшить точ­

ность приближения функции (9.9) к графически заданным частотным характеристикам.

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

L (со) =

20 lg | K*W ( j со) | = 20 Ig К — 20 lg со 4 - 20 1g

У-м

+1 1

 

20 lg £

+

1

1

420 lg

. tO

О.А

 

^ Pd

Pd ^

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

20 lg

/JL _ ? £ . _ /

 

 

_i_/

,

 

 

 

 

 

 

 

 

J h

P*

J J Pa Pa ^ J

 

6(со) =

Arg {K*W ( » }

=

j

+ arctg^---- arctg^- -+-

 

 

 

 

 

 

 

+

arctg яя» Т,2Tвв

 

i

 

2afecо

 

 

 

 

 

 

2- a r c l g —

j

 

или

 

 

 

 

 

 

 

“d +

Pd —“

“ft

i Pft — “

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L (со) =

20 lg К — 20 lg со +

L.,

( cd) LT(со) +

(со) — A, (со)

 

 

0(ш) = — 4 -. 0t (со) — 0г (со) + 0d (со) — 0* (со).

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ц (со) =

20 lg

>

+

1

 

 

 

 

 

 

 

tO-r

 

 

 

 

 

 

 

Ad(co) =

20 lg

. t o

f lr f

 

j J L — b L A - j

 

 

 

УР 7 - Р 7 - У

 

 

 

y Pd

Pd ^

 

 

 

:(co)=arctg — ,

(9.12)

(9.13)

(9.14)

(9c15)

(9.16)

(9.17)

(9.18)

9d(co) = aretg ■ 2° f ^

(9.19)

“d T Pd — “

3 9 2

 

СИНТЕЗ'СЛЕДЯЩИХ СИСТЕМ ПРИ СЛУЧАЙНЫХ ВОЗДЕЙСТВИЯХ

[ г л . IX

Аналогичные

выражения

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

для Lt (w)> Lk (ш).

G7*(ш), 0* (со). Полученные выражения (9.12), (9.13) или (9.14), (9.15)

удобны

тем, что

они позволяют

легко

анализировать влияние каж­

дого

из

полюсов

и нулей передаточной

функции

 

системы в отдель­

ности

на ее частотные

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

Легко

видеть, что

выраже­

ния (9.12) и (9.13) состоят

из членов двух типов:

соответствующих

комплексно-сопряженным полюсам

и нулям и соответствующих веще­

ственным полюсам и нулям.

 

 

 

 

 

 

 

Найдем теперь влияние

малых

изменений в положении

полюсов

и нулей на вид частотных характеристик (9.14) и (9.15).

 

При

данном

со амплитудная

логарифмическая

частотная

характе­

ристика

(9.14) является

функцией нескольких переменных, а именно,

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

L

( p

* * +

А э : * .

Р , , +

 

д

Р а РА * -

+ сот

- | - Д

и >

т ,

ш г - | - Д

ш г ) =

 

=

L (ad,

ak,

prf.

pftl

toT,

со7.)-[-^ Д ос1У4 - ^ - Д р (1 +

 

 

 

 

+

^

 

+

| ^

? ,

+ ^

+

;

ё

1/» г -

(9 20)

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

Д / .

( с о ) =

dL.

 

dL,

 

дЬл

Даь

6L-,

До)г —-

 

- т —

'

Д с о .

 

do.d

дш]

 

с -

4 '

ощ

 

 

“ А. .

 

 

 

 

 

 

 

 

 

 

дР * - - й г - л«*-

(9-21)

 

 

 

 

 

 

 

 

■ЙР*

 

 

 

Аналогично

 

получим:

 

 

 

 

 

 

 

дО.

 

ва

 

два

 

двт

 

 

 

 

 

д

 

 

 

 

Дб (“) =

Дш*+ ж

Д^ +

1=7 Да* -

5=7 Дш7-

 

 

 

 

 

д?а

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

<9-22>

 

 

 

 

 

 

 

 

dh

 

3] ПЕРЕДАТОЧНАЯ ФУНКЦИЯ КОРРЕКТИРУЮЩЕГО УСТРОЙСТВА • 393

Найдем

выражения для

частных

производных, входящих

в фор­

мулы (9.21)

и (9.22):

 

 

 

 

 

 

 

 

6 L.

 

 

20

&

 

 

 

 

 

 

 

 

 

2,303 Mi)T

 

(9.23)

 

 

 

 

 

 

 

 

 

 

 

 

dLd _

 

ad

 

 

 

 

 

 

a

 

40

 

 

| W d \2

 

(9.21)

"d d*d

 

 

2,303

 

 

 

Pd

/М.,

 

'

40

(М±Ж±(^-^Г)

(9.25)

 

d L d

 

 

 

 

 

 

U^dl2

 

 

<?Pd

 

2,303

 

 

 

 

 

 

, n

 

 

57,3-

 

 

 

 

 

 

 

dOT

 

 

_____ (

 

 

 

 

(9.26)

 

dw-

 

 

1 +

 

 

 

 

 

 

 

 

 

Ш 2 '

 

 

 

 

 

 

 

 

 

 

 

 

Pd did

 

 

2 • 57,3

t H

 

i t r -

m

(9.27)

 

 

 

 

 

 

I W'dl2

 

Pd H r

=

 

- 4 - 5 7 ,3

l

w„ i2

 

(9.28)

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

1 ^ 1

=

 

. (i)

aa

 

 

J Pd

Pd ^ J

(9.29)

 

J V d

- f a - J

Аналогичные выражения можно написать и для остальных шести частных производных, входящих в формулы (9.21) и (9.22).

Кривые

(9.30)

и

<Pi

(9.31)

1 +

иеображени на рие. 9.8 и 9.9.

Пользуясь этими кривыми, легко определить влияние перемеще­ ния Д(вт вещественного нуля и церемещения Дшг .вещественного по-

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