книги / Системы экстремального управления
..pdfПусть / (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) |
вида |