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

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

..pdf
Скачиваний:
8
Добавлен:
19.11.2023
Размер:
33.28 Mб
Скачать

Пусть / (2) = g sin 0)2, т. e. имеет место гармоническая модуляция. Тогда выход первого инерционного звена будет иметь вид

у (2) = Wx (р) (х -\- g sin 0)2) = х 4- gGx (o))sin[o)2 — <px(o))], (7.4.3)

где Gx (со) — амплитудно-частотная характеристика пер­ вого звена, а фг (о)) — фазово-частотная характеристика этого же звена. Эти функции определяются из следующих соображений.

Пусть передаточная функция Wx (р) имеет вид

W 1(p)

Ътр т +

•••+

bip + 1

(7.4.4)

 

апРП +

• • • +

Р + 1

 

где, как обычно, р = ^ — оператор дифференцирования.

Подставляя v = g sin 0)2 в уравнение и = Wx (р)и, ищем решение в виде и = gG (со) sin (о)2 — Фх (ш))- После необхо­ димых преобразований получаем выражения:

 

 

Gx (0 ) =

[ Р »

+

S2 (<D)]V«,

 

(7.4.5)

 

 

 

 

44(") = ?£$ >

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

Р С ) - 5

| у .

 

 

 

 

 

 

 

 

о / \

Ьс4ad

 

 

 

 

 

 

 

 

5 (“ ) =

7 ф ? .

 

 

 

 

 

 

 

 

а (о)) =

1 — Ь2оо2

 

Ь4©4 —

 

 

(7.4.6)

 

 

b (©) =

bx(ù 6303+ Ь5©5 —

 

 

 

 

 

с (о)) =

1 — а2(о2 +

<z4o)4 —

 

 

 

 

 

d (о)) =

ax(ù — а3©3 + а5о)6 —.

 

 

Так,

для

апериодического

звена

(7.1.4),

для

которого

ах =

Тх,

at = 0 (2

=

2,

.,

п),

bj =

0 (/

= 1,

., ni),

получаем

а (©) =

с (со) =

1,

b (©) =

0,

d (©) = Тха> и

 

 

Gx(со) =

 

1

;

tg ф! (w) =

TX(Ù.

(7.4.7)

У 1 + т2©2

Далее, разрешая (7.4.1), имеем

q (t) — к [у (t) — у*]2 + д*

или

q (t) = к {х у* + gGj (to) sin [coi— <px (©)]}2 -f q*. (7.4.8)

Переходим к выходу объекта:

Q (t) = W2 (p) q(t) = к (x y*)2 + q* +

+ 2kg (x y*)Gi (<Ù)G2(o))*sin Ы — ф^ю) — <pa (to)] +

+ % kg2G\ (со) — Y kg2G* )G\ (2©) • cos [2ooi — 2q)x (со) —

- Фг (2co)], (7.4.9)

где G%(со) и ф2 (со) — амплитудная и фазовая характерис­ тики второго инерционного звена W2 (р), которые вычис­ ляются аналогично по формулам (7.4.5) и (7.4.6).

Осредняя выход синхронного детектора и восполь­

зовавшись соотношением sin ( a t — ф) sin at ~ 4- fcos ф—

и

— cos (2coi— ф)1, получаем

В = kg2 (x — y*)Gx Ù)G2 ( со) - C O S [фх (со) + ф2 (со)]. (7.4.10)

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

Так, при Фх(ш) + ф2 (м) <С Y обратная связь должна

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

ТГ = -

(7.4.11)

где р, — величина обратной связи.

Если же я фх (со) + Фа (<а) > - у - , то сигнал В перед интегрированием уже не Следует инвертировать, т. е.

£ = ИЯ.

(7.4.12)

Данное обстоятельство порождает определенные трудности при экстремальном управлении инерционными объек­ тами с неизвестными динамическими свойствами.

Рис. 7.4.2. Блок-схема экстремального управления методом син­ хронного детектирования с введением инерционного авена.

Ф1

Эти

трудности

можно

преодолеть,

если

величина

+ Фг в процессе поиска

изменяется

не

более чем на

я,

т. е.

(фх +

ф2)тах

— (фх — Фг)пИп < Л.

 

(7.4.13)

 

 

 

 

 

Пусть

ф* — среднее

значение

этой

суммы

 

 

 

 

 

т

 

 

 

 

 

 

 

 

 

 

ф* = П т

(ф! +

ф2)

 

[(ф! +

ф2)тах +

(ф! +

ф2)ттЬ

 

 

 

 

 

 

 

 

 

 

 

 

 

(7.4.14)

=

Будем

на синхронный

детектор

подавать

не

/ (t) =

g sin

cat, а его задержанное на ф* значение. Для этого

достаточно ввести инерционный блок W (р), фазовый

сдвиг

которого

на

выбранной

частоте

о>

совпадает

с

(7.4.14). Тогда на выходе фильтраY,.

 

получаем

(рис. 7.4.2) сигнал,

пропорциональный величине

 

 

 

 

cos (фх +

ф2 — ф*),

 

 

 

 

которая при условии (7.4.14) сохраняет положительный знак. Этот факт позволяет фиксировать знак обратной связи.

Теперь перейдем к случайному возбуждению. Пусть / (t) — стационарный случайный процесс с

нулевым математическим ожиданием, автокорреляционная функция которого равна R (т). Из (7.4.1) получаем для

выхода

фильтра

Ф (см. рис.

7.4.2)

 

т

 

 

 

В = lim 4rr\f(t)W г(р) {к[х +

Wi (р)/ (t) — у*]2+ g*}dt,

Т -ю о

■* J

 

 

 

 

 

 

 

(7.4.15)

Преобразуем это

выражение

к

виду

 

т

В = lim|[/с(я -

уУ + g*] 4"S / ( 0 ^ + 2fc(* - у*) X

Т->СЮ

Q

т

x - r \ f ® W i ( r ) w , ( p ) № & +

О

Т

+ -£ -$ № W t (p) [ ^ ( р Ш О Р * } .

О

Как легко заметить, первое и последнее слагаемые этого выражения равны нулю. В результате получаем

т

В = 2к(х — у•) lim-i-Ç f (t) W1(р) W2 (р) f (t) dt. (7.4.16) г-.» J

Рассмотрим подробнее этот интеграл. Пусть h (т) — им­ пульсная переходная функция инерционного звена с передаточной функцией Wx (p)W2 (р)-

оо

h (т) = ~Тл \ Wi(i(ù)W2(Uû)eiax d(ù, где i = Y — 1.

Тогда имеем

оо

Wl (p)W2(p)f(t) = $ h{ t )f (t - x)dx .

о

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

Тоо

- j r jj f(t) ^ h (x)f (t T) dx dt

0

T

= [ A M [ 4 - $ H 0 / ( < - X ) d t \ i x .

0 0

Как видно, в квадратных скобках заключена оценка

автокорреляционной функции

R (т) возбуждения

/ (t)

на базе

Т.

Эта

функция предполагается

известной.

В результате

получаем следующее выражение для

В:

 

 

 

 

оо

 

 

 

 

 

В =

2к(х -

у*)^ h (т) R (т)dx,

(7.4.17)

 

 

 

 

0

 

 

 

где автокорреляционная

функция R (т) имеет вид R (т) =

1

т

 

 

 

 

 

 

I*

x)dt. Так,

в случае апериодических

= lim-^7

\

Т-too 1

О

 

 

 

 

 

 

звеньев

о

 

 

 

 

 

 

(7.1.4)

 

 

 

 

 

h( т) =

1

 

 

 

1'г

 

и для экспоненциальной автокорреляционной функции R (т) = с2<гат,

где а- — дисперсия возмущения / (t), получаем следую­ щее выражение для В :

В = 2к (х у ) а2

(74 _ ï\) (cCTi—1) +

+ _____ 3 _____1

^ (Г1-Т2)(аГа—1) J*

Из этого выражения хорошо видно, что сигнал В умень­ шается с приближением возбуждения / (t) к белому шуму (напомним, что случай а —>оо соответствует белому шуму). Это хорошо согласуется с известным положением о том, что инерционные объекты лучше откликаются на медленные процесссы ( а —►0), чем на быстрые (çt—> 00).

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

В = 2к(х — у*)<?[-^--

ctTiTz — Ti — 'Гг

\

(aî'i — 1) (аТг — 1)

J‘

Несложный анализ показывает, что для любых значений

параметров а,

и Т2 выражение в квадратных скобках

всегда

больше

нуля.

Поэтому в данном случае можно не опасаться появ­

ления

неустойчивости за счет изменения параметров

объекта Тх и Тг.

Любопытно, что при гармоническом возбуждении ана­ логичного объекта возможно появление неустойчивости при больших значениях постоянных времени Тг и Т2 (фазовый сдвиг в объекте при этом превысит я/2). Сле­ довательно, использование случайного возбуждения не предъявляет строгих требований к динамике экстремаль­ ного объекта, что особенно важно при изменении его динамических свойств.

§ 7.5. Улучшение процессов экстремального управления инерционными объектами (объекты первого рода)

Как показано выше, инерционность вносит значи­ тельное запаздывание в процесс экстремального управле­ ния. И чем больше инерционность, тем медленнее процесс поиска экстремума. Можно ли преодолеть это препятствие при управлении процессами с большой инерционностью?

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

у = W t (р) х.

(7.5.1)

Здесь х — вход объекта, у — вход нелинейного

звена

Q — Q (y). Q выход этого звена — показатель качества объекта.

Задача ускорения процесса поиска сводится к быстрей­

шему переводу величины у из

исходного

состояния у0

в требуемое ух с заданными

нулевыми

значениями

производных у —ÿ = .= 0 в этом состоянии. Этоусловие гарантирует неподвижность системы в состоянии yv Кроме этого, в процессе перехода системы из состояния у0 в состояние ух управляющее воздействие х (t) не должно выходить за заданные пределы, которые определяются

физическими особенностями задачи (так,

ускорение

х (t)

не может быть равным бесконечности,

 

скорость

х (t)

также конечна и т. д.). Таким образом,

решение

этой

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

Для конкретности рассмотрим сначала апериодическое

звено первого порядка (7.1.4):

 

ИМ?) = 7 т Г + Т '

<7-5-2>

которое реализуется следующим дифференциальным урав­ нением:

ТгУ Ч- У = х (t).

(7.5.3)

Задача заключается в том, чтобы перевести у из состояния у0 в за минимальное время, причем управление х (t) не должно выходить за заданные пределы:

Æmin ^ X (t) #тах« (7.5.4)

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

печить управляющему

воздействию х (2)

произвольное

значение.

 

 

Например, при управлении печыо у — температура

в рабочей зоне печи,

а ж — температура

управляемого

нагревателя, которая не может превосходить температуру

плавления

нагревателя,

т. е.

х (t) <[ 2ПЛ.

 

В исходном состоянии, которое предполагается стаци­

онарным,

имеем

 

 

 

 

Уо = У (0) =

х0 и

у (0) = 0.

(7.5.5)

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

на рисунке). В момент, когда значение у (t) совпадает с требуемым и заданным уг, следует управление х (t) сделать равным ух. Это стабилизирует процесс, так как в

этом

случае

ÿ (t) =

0

и

получим у (t)

=

у1 для

всего

последующего времени*.

 

 

 

 

 

 

 

Таким образом, алго­

 

 

 

 

 

 

 

ритм

ускоренного

сме­

 

 

 

 

 

 

щения у0 в Ух сводится

 

 

У{ >

 

 

 

 

к следующему:

(0 ^

 

 

 

 

 

 

Первый

этап

 

£р=Уо

 

 

 

 

< t< tx )

 

 

 

 

 

 

 

 

 

 

x(t) =

\ ^тахПриу1>1/о,

 

 

о

 

 

 

 

 

1жш тприу1<1/0.

 

Рис. 7

 

 

 

 

 

 

(7.5.6)

 

 

£

 

 

 

 

Этот этап продолжается до момента совпадения у

(t) и

У1 , т. е. у (tj =

 

*х)

 

 

 

 

 

 

 

Второй этап (t >

 

 

 

 

 

 

 

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

 

х (t) = yv

 

 

(7.5.7)

 

 

 

 

 

 

 

 

 

 

 

 

 

Г

* m a x

П р и

У (* ) <

Ух,

(7.5.8)

 

 

Х^

“"1

 

Vi

при y{t) = yx.

Аналогично

 

 

 

..

 

Гхт1п при y(i)> yi,

 

 

 

 

 

(7.5.9)

 

 

*(*)="{

 

Ух при

у (*) =

ух.

 

 

 

 

 

 

 

 

Блок-схема устройства, решающего эту задачу аппара­

ту р а

показана на

рис.

7.5.2.

Здесь

уставка ух сопо­

ставляется с текущим

значением у (t)

и

разность

а —

— ух у (t)

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

О,

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

 

 

|

*тах при а > 0 ,

(7.5.10)

 

 

1 ^min при а < 0 .

 

 

 

При а <С 0

получаем х — xmin-

Если

задача выполнена

и

у (t) = Ух, то переключатель

П

устанавливается

в

нравом

положении

(а = 0) и в системе реализуется

стационарное управление д; = ух, которое стабилизирует си­ стему до тех пор, пока у = ух, т. е. а — 0. При новом значении уставки у2, отличающемся от предыдущего, процесс перевода ух в у2 произойдет также с минимальной затратой времени, так как используются при этом экстре­ мальные значения управления жтах или a:min.

Рассмотренная схема может быть использована лишь при возможности наблюдения величины выхода апериоди­ ческого звена у (t). Однако часто в экстремальных объектах

Рис. 7.5.2. Блок-схема оптимального уп­ равления апериодического звена.

эту величину невозможно непосредственно измерить. В этом случае можно воспользоваться одним из двух путей.

Уравнение экспоненты у (t) при х = хтах (см. пунктир на рис. 7.5.1) имеет вид

- J L

(7.5.11)

У (t) = Уо + (Zmax — Уо) (1 — « т‘)-

Определим момент времени £lt когда происходит переклад­ ка управления с хШах на yv Получаем

t1 = — T1ln f l - - - - - Vl~ y? ) .

(7.5.12)

I

хтах~Уо )

V

'

Теперь, вводя блок, реализующий это преобразование, получаем управление без обратной связи, т. е. без изме­ рения выхода у (г) инерционного звена. Для этого доста­ точно экстремальное значение æmaxподдерживать в течение времени tv после чего следует перейти на х = ух. Такое программное управление оказалось возможным за счет знания зависимости (7.5.12). Однако практическая реа-

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

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

Рис. 7.5.3. Блок-схема оптимального управления с моделью инер­ ционной части объекта.

выход модели ум (t). Образованное таким образом управ­ ление х (t) используется для управления инерционным объектом.

Блок-схема такой системы показана на рис. 7.5.3. Здесь управляющее устройство УУ реализует алгоритм синтеза оптимального управления х (t) на основе устав­ ки уг и состояния модели у№(t). Такое УУ работает опи­ санным выше образом (на рис. 7.5.2 оно очерчено пунк­ тиром).

Рассмотрим теперь звено второго порядка

W i { p ) = ^ r ,

(7.5.13)

которое эквивалентно уравнению второго порядка

Tjÿ = х (t)

 

(7.5.14)

с начальными условиями

 

 

У(0) = Уо; У(0) =

0

(7.5.15)

и ограничениями на управление х

(t)

вида