Файл: Живоглядов, В. П. Адаптация в автоматизированных системах управления технологическими процессами.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 23.10.2024
Просмотров: 77
Скачиваний: 0
жает работать с постоянными ненастраиваемыми параметра ми. Структуры блоков УБ и ВУ-1 считаем известными. Возмо жен и другой путь, когда все устройства УБ, ВУ-1, ВУ-2 реали зуются программно на УВМ.
Критерий оптимальности работы системы R определяется значениями выходного показателя до, согласно (3. 1).
Функция потерь W — выпуклая дифференцируемая. Пред полагаем, что в допустимой области изменения параметров а функция /?(а) имеет единственный экстремум. Условия су ществования и единственности экстремума должны быть рас смотрены особо, применительно к отдельным классам систем. Задача синтеза состоит в следующем. Требуется найти алго ритм функционирования ВУ-2, обеспечивающий в условиях случайных помех сходимость параметров к оптимальным зна чениям а*, при которых
/?(а*) = т!п/?(а). |
(3. 11) |
а |
|
Воспользуемся методами стохастической аппроксимации и теории чувствительности. Рассмотрим вначале простой одно мерный случай, а затем проведем обобщение результатов на многомерные линейные системы [3. 151.
3. 3. 1. Настройка коэффициентов распределенного контроля в замкнутой схеме
Пусть объект с распределенными параметрами устойчив и описывается следующими уравнениями в операционной
форме: |
|
|
д(х,р) = Ф(х,р){и(р)+р°<р)]+А1111(х,р),(г(х,р)}; |
(3. |
12) |
(0 < х < /н, 0<х</„) |
|
|
и>(р)=Фи(р)и(Р)+Фр (Р)\>-0(Р)+1{ Р )+ М Н *’Р)Л'{Р)\- |
(3. |
13) |
Здесь Фи (р ), Ф,А(р), Ф(х, р) — передаточные функции, при
чем Ф{х, р) зависит не только от оператора р, но и от прост ранственной координаты х. Возмущающие воздействия ц°(/), K(t) считаем стационарными случайными функциями. Распре деленное возмущение ц(х, /) — случайное поле, стационарное по аргументу /; А\ и А2 — некоторые линейные операторы, ха рактеризующие связь между ц(х, /) и выходными сигналами
16»
q(x.. t) и w(t)\ контроль q (x, i) осуществляется в дискретном ряде точек х i (i = О, 1, /). Управление определяется сле дующим выражением:
/ |
|
« ( » = —'Фр№)ЛсчФ;,р)+Фр(р)в(р), |
(3. 14) |
г=0 |
|
гдеФ ^д) — передаточная функция регулятора. Все входные и выходные переменные считаем скалярами. Предполагаем также, что инерционность и запаздывание в канале «управ ление и — выход w» значительно больше, чем в контуре ре гулирования по косвенным показателям. Это естественное предположение обычно оправдывается на практике. Одновре менно это одно из условий эффективного применения косвен ных показателей для целей управления.
Определим градиент функции потерь W по а:
Va W =W ,(e)S. |
(3. 15) |
где 5 — вектор функций чувствительности. |
|
5 = УЗ Д .....Si..... 5/ || 5/== |
. |
(/=0,1,...,/)
В соответствии с [3. 15] запишем алгоритм настройки пара метров распределенного контроля внепрерывной форме:
“57T(0W'(e)S. |
(3. 16) |
где у (/) — матрица переменных коэффициентов, удовлетво ряющих условиям Роббинса—Монро. Обычный путь получения функций чувствительности S/ [3. 16] состоит в построении модели чувствительности, на вход которой подается сигнал ге. Учитывая спорадический характер поступления сигнала w и возможность случайного перерыва, откажемся от такого пути ипопытаемся сформировать функциичувствительности 5 /, используя лишь косвенные показатели qi. При этом можно также надеяться на ослабление требований к точности реали зации модели чувствительности и упрощение системы.
169
Из уравнений (3. 12) — (3. 14) находим
I
^ ) = - Ф ц(р)Фэ(^)|5]аД ф (д ;/)р)1Ао(/7)+Л{^(х1/;),(7(^.,1р ) ) ] -
/=0
(3. 17)
— В{р) 1+ Ф^ (р)|1о(/»)+Х(/0+Л2{|*(х,/?),«;(/>)],
/
^•(Р) = 9 (^ у ./0 = -Ф (л:у.Р )Ф э(^{2 а/1Ф(^./7)|А0(Р)+ /=*0
( 3. 18)
+Mv-(*>P)q(xi,P)}l-f>(p)) + (&(xj.p)\>.0(p)+Al{\>.(-A,p),q{Xj,p)},
где
ФР(р) |
(3. 19) |
Фэ(Р) = |
|
1 + Ф ^ ) 2 |
«/Ф(XI,р) |
Дифференцируя w по а/ и учитывая [3. 18], получаем функ ции чувствительности:
Slip) = |
= — Фи(Р)Фэ(Р)<1(хьР) = -®o(P)<Jlix iP)> |
|
|
5/(/) = - 7 - 1|?(х/,р)Ф0(/;)), |
(3. 20) |
где Z,- 1 [ - ] — |
обратное преобразование Лапласа. |
Алгоритм |
настройки определяется формулами (3. 16) и (3. 20). Струк турная схема системы управления приведена на рис. 3. 6, где указаны передаточные функции отдельных блоков. Приняты обозначения: J— интегратор; ~i(t) — усилители с переменным
коэффициентом т (0; X — множительное звено; W' (•) — безынерционный преобразователь с соответствующей характе
ристикой; А — блок, определяющий |
характер воздействия |
||
распределенного возмущения |
ц (х, t) |
на сигналы |
q(x, i) и |
w(t). Случайную помеху h(t) |
с нулевым средним |
в данной |
|
схеме можно исключить из рассмотрения, изменив |
соответ |
170
ствующим образом характеристики возмущения X(t). Выпи~ шем для сравнения алгоритмы настройки коэффициентов рас пределенного контроля при квадратичной функции потерь:
а) в замкнутой схеме управления
- ^ - = ^ t ) ( w —w*)L-4q(xi,p)<t>Q(p)), i = 0 |
, \ , (3. 21> |
б) в разомкнутой схеме |
|
- ^ = 4 { t ) { w —w)L-4q{x,p)<$>qJw Ср)}. |
(3. 22> |
' .г |
. |
i . v |
ill |
l
л |
(3. 23) |
w = |
|
i= 0 |
|
Здесь <t>g.w (p) — передаточная функция между |
сигналами |
л
<l{xi, р) и w; w — предсказанное значение выхода w. Из срав нения (3. 21) и (3. 22) видно существенное различие алгорит мов, несмотря на их внешнее сходство. Разомкнутая схема распределенного контроля настраивается на максимум тес ноты связи между выходом w и косвенными показателями q{x , t). Настройка замкнутой схемы управления осуществля
ется на минимум тесноты связи между отклонениями |
выхода |
от заданного значения (ш—ш*) и сигналами q(X{, t). |
|
Выше предполагалось, что задание 0 регулятору |
(р)— |
постоянная величина. Управление происходит по |
сигналам |
■q(X[, i), выходной сигнал ш используется только для адапта ции системы. Однако информацию, содержащуюся в w, и ал горитм (3. 16) можно применять и для коррекции задания 0. Достаточно лишь включить 0 в вектор настраиваемых пара метров а, соответственно увеличив размерность а.
Сформулируем условия сходимости скалярного параметра а к точке (в область) оптимума [3. 17]. Справедлива сле дующая
Теорема
Пусть при любом фиксированном aeQ(a) процесс !/(£,/)= у а W,/>0,EeQ является сепарабельным стохастически
непрерывным случайным процессом, задан на вероятностном
пространстве |
(Й£, о, Р) с вероятностной мерой |
Р на сигма- |
|
алгебре множеств а случайного пространства |
и справед |
||
ливо условие |
|
|
|
|
t |
|
|
|
|
(3. |
24) |
|
о |
|
|
где |
Г(а=а*)=0; |
|
|
|
г(а)>0 При а>а*, |
|
|
|
г(а)<0 при «<а*, |
|
|
|
'■(ai)>'"(aa) при <4 > а 2, |
(3. |
25) |
172
I r(cc) I <A /=const при — CO < « < o o .
Тогда с вероятностью единица существует решение уравне-.
ния (3. 16):
do.(l,t)
—т(ОК(?,о
dt
и имеет место равенство |
|
ЖЕ: 11щ»(«,0=а*} = 1. |
(3. 23) |
/ —» о о |
|
Доказательство теоремы здесь не приводим. Для алгорит ма (3. 21) многопараметрической настройки условия сходи мости подобны (3. 24), (3 25). Для линейных стационарных систем условия сходимости, отказавшись от строгой записи, выразим через спектральные плотности случайных сигналов. gzS (со). Потребуем, чтобы с вероятностью единица выполня
лось равенство
СО
^ g esHflfo>=ar(a), |
(3. 27). |
--СО
гДе 5 es(w) — взаимная спектральная плотность е и со, кото рая легко вычисляется, если известны спектральные плотности входных возмущений и амплитудно-фазовые характеристики системы.
Итак, полученный алгоритм при выполнении определен ных условий обеспечивает асимптотическую сходимость а к точке оптимума даже при наличии большого чистого запазды вания в объекте [3. 17]. Однако качество процесса настройки, скорость сходимости существенно зависят от динамики си стемы по каналу «а—Va^»- Для улучшения качества настрой
ки можно модифицировать алгоритм (3. 16), используя не которые результаты раздела 1: ввести, например, в алгоритм обратную связь по аналогии с (1. 84).
3. 3. 2. Автоматическая оптимизация многомерных распределенных систем
Обобщим метод синтеза алгоритмов параметрической оп тимизации на многомерные линейные системы. Уравнения объекта и регулятора в операционной форме имеют вид
173.
Ч р ) = ^ и { р )ф ) + % |
(р)1^)+ М »)+ Л 2{[*(*,/?),ш); |
(3. |
28) |
|||||
^)= O i(o )« (/’) + |
c^'1J.(p)p-0(p)+/l1{ii(A-,p),9(^)}; |
(3. |
29 |
|||||
|
u(o) = - % p(pjRq(p), |
|
|
(3. |
30) |
|||
где w(p), q(p), и(р), [Л |
X, ji—векторы |
типа |
|
|
|
|||
q = |
|д0,....<7; II T,w= |
II wu...,wr |г , u= |ul,...,um |т-, |
|
|||||
|
|
|
т— знак транспонирования; |
|
|
|||
Фи(Р).ф1 (Р) и Ф|х(Р)>ф11х (Р) — матрицы |
передаточных |
функ |
||||||
|
|
|
ций объекта для управляющих |
|||||
|
|
_ |
и |
возмущающих |
воздействий |
|||
|
|
соответственно; |
|
|
|
|||
|
|
Фр(д)—матрица |
передаточных |
функ |
||||
|
|
|
ций регулятора |
размера |
||||
|
|
|
_ /л(/+1); |
|
|
|
|
|
|
|
|
/<— матрица |
весовых |
коэффици |
|||
|
|
|
|
ентов размера ш(1-И ). |
|
30) |
||
После некоторых преобразований из (3. |
28) — (3. |
|||||||
получим |
|
|
|
|
|
|
|
|
w(p) = - Фы[ I + Фр кФ гП Ф рК ( < Г ^ ( р ) + А 1М**).Ч}) + |
||||||||
|
+Ф(л (А0(р)+^(/’)_(г^2|н-(л:>/?)>ау); |
|
(3. |
31) |
||||
я(р) = - Ш I + % W i ПФр f - lI ) ( % ^ G p ) + |
|
|
||||||
|
|
+Mp(x,p),q(p)}> |
|
|
|
|
||
где I — единичная матрица. |
|
лишь |
отклонением |
|||||
Критерий качества |
R определяется |
|||||||
вектора |
w от вектора |
|
заданных значений w* |
|
|
|
||
|
R=M[W* |
(w,w*)}. |
|
|
|
|||
Введем |
вектор настраиваемых |
параметров |
|
|
|
|||
|
**= |
II а1|...,ак... ап II г - |
|
|
|
174