книги из ГПНТБ / Живоглядов, В. П. Адаптация в автоматизированных системах управления технологическими процессами
.pdfгде Kui*1), |
(х1)— коэффициенты, |
зависящие |
общем |
||
случае от х 1, или в дискретной |
по |
t и х 1 системе |
|||
|
^[ArJs]=/CUKM[s—к] + /^ к (X[s—к], |
(1. 149> |
|||
|
X |
t |
|
коэффициенты. |
|
где К= Т Х х ’ S== д7 ’ Кик’ |
|
||||
Функцию |
потерь |
примем в виде |
lFs= (g*[s^.T] —KlLU[s\ — |
||
—К^ p[s])2. |
Для |
нахождения |
оптимального |
управления |
a*[s] воспользуемся формулами (1. 39), (1. 41). Можно по казать, что система является нейтральной и оптимальному значению «*[5] ноответствует минимум функции ccls j_. =.
J <7*[з-И1 — Kua[s] —
Q'Vs) s—1 l
- K ^ p[s])aP( (A[s]) Y \ { f ] ЯМк.ЛЗИЬ-к],alj—к] J dQ,
/=1 K= 0
причем для аддитивных гауссовых независимых помех
|
Ply[*j\ I v-[i— K ] M i— K])= |
|
|
1 _ |
f |
(у [ ^ \ —К ^ и - к ] —Кики\1—к]) \ _ |
|
oh\K\V2~ |
Р1 |
2о*А[к] |
J |
После вычислений получаем
|
|
|
(1. 150)- |
где' |
|
S— 1 |
I |
|
^ |
||
а' S + ' ~ |
P(v-s) J |
J I [~}P { y n j IPj—Ki llj —k) J ^ = |
|
■ 1 ' |
Tv"* ^ |
- 1=1 |
K=° |
|
Q(n-i) |
|
|
90
|
|
|
|
1 |
|
s—1 |
|
|
= A J |
... jexp |
А Л |
^ j(lV ~ PHy-i)2A i+ |
|||||
2 a 2 |
||||||||
|
|
|
|
|
|
/=1 |
|
|
s — |
1 |
l |
|
|
|
|
|
|
^ |
^ |
{Укj |
V |
V'j —k |
Kukuj—к)2у]к |
djir ..rfixs, |
||
j —1к= 0 |
|
|
|
|
|
|||
|
|
% = |
0, ^ — » 0<o*A<oo ; |
(1. 151)) |
i4i— постоянный коэффициент; a"s+T
лишь наличием множителя |xs под Введем дополнительно обозначение
отличается от a '^ ^ . знаком интеграла.-
^j=(y\K,j]—Kulcu[j—K]).
Выражение в (1. 151) в фигурных скобках равно: |
|
|
|||
{■}= —{—|*V 0+ tlSiA + .»+PV iJPs-i—гЫЯо+гн)— |
(1. |
152> |
|||
|
5—1 |
1 |
|
|
|
—2\i1(q°1-\-r\i2)—...~2^s- 1q°s- 1+ - ~ Y |
|
|
|
|
|
|
j= |
0 |
л с = 0 |
|
|
где |
Ро =^Л ~(т + Дд); |
|
|
|
|
|
А |
|
|
|
|
A = ^ ^ " ( A 1+P2A i+m ) при |
s—I; |
(1. |
153> |
||
|
P /-i= |
|
|
|
|
( |
I |
|
|
|
|
I
при 5—/> /
/я =
при K — S —1</.
z = 0
9h
Если же |
для |
всех 2==0,1,.../, то |
|
|
т = |( 1+1 |
при |
s—г > / |
|
{ s—i |
при |
k = s — г'<Х |
При нулевых начальных условиях все £К;- равны нулю для ■K>s—j. С учетом этого получим
|
N |
<7i°= |
i j ^ . / = « + ' ПРИ z= 0>b--,s—2, |
|
АГ=0 |
Я S ~ l
1
2а2/г[0]
тде |
|
jV=min{/,s—i—1). |
|
||
"Подставив |
(1. 152) |
в выражение |
для |
a's+T, a"s+T и вы |
|
полнив ряд |
интегрирований, |
находим по формуле (1. 150 |
|||
оптимальное управление m*[s) |
|
|
|||
k*[s]= -jr <7*[s+ t] — |
2а2/гр s_ u |
-1 4" |
|||
|
Лц |
|
|||
|
г |
(^,s—1+••• + |
r |
W . . . ) + |
|
|
P s - 2,1 |
Pol |
|||
P s - 2,1 |
|
1+ p s —3,1 ■(5l,s—2 + •■ Pol |
|||
H- |
|
~+^K,S—1+ |
P s— 1— K,1 -(ZK,s— 2 + |
||
|
|
||||
П^s-l-U |
|
|
|
||
i= l |
|
|
|
|
|
'+ - H — |
~— 5*с,к)...)Ч— |
i-------—--------- (Zl,s—1 + |
П
/=1
Э2
+ |
^ ------в 1в- 2 + - + - £ Н « ) - - ) И . (1- |
154> |
|
|
Ps—2—1,1 |
Poi |
|
где |
Pi.1=Pi— |
(1. |
155> |
|
Pi-i,i |
|
Ро1~Рош
Алгоритм оптимального управления (1. 154) можно запи сать в более компактной форме:
S— 1 |
I |
Ч s+ t |
§Kj(yKj KuiiUj—fc) |
W h P s - i , i ;
j = 0 K—0
(1. 156)'
причем коэффициенты gKS_\ равны
S k .s — 1 |
K 0 , 1 , - . . , / , |
157> |
|
(1. |
ПPs—I—1,1
i= I
аостальные коэффициенты вычисляются по рекуррентной формуле
|
g*j=gK,!+1—^Г------- ■ |
/= s —2,..., 1,0. |
.(1- |
158) |
||
|
|
Pi—к, 1 |
|
|
|
|
При т)к=1 и любых к и / и |
0< р<1 имеет |
место |
не |
|||
равенства: |
|
|
|
|
|
|
|
0<С§ко<С---<СЯку—i <^ - § K j <-ig ic ,s — |
О- |
159)' |
|||
' > |
|
- , ‘г |
«*-■> - i r ^ ^ r ^ o . |
|
||
Доказательство |
этого |
факта опирается на |
результаты,, |
|||
описанные |
выше. |
оптимальный |
алгоритм управления явля |
|||
Таким |
образом, |
ется линейным. Он предполагает использование распределен ного контроля с весовой функцией, убывающей по мере уда ления точки отбора сигнала от точки приложения входных
93
^воздействий, и обработку информации во времени. При этом старая информация берется с уменьшающимся весом. Ин
формация в УУ обновляется. |
режим и |
распространим ре |
|||
Рассмотрим |
стационарный |
||||
зультаты на непрерывную систему |
при |
/( к=1, |
о2Л[^]= аг/,в |
||
к = 0 ,1 .- |
|
|
|
|
|
Естественно |
предположить, |
что |
в УУ можно |
будет выде |
лить блок обработки информации во времени и устройство распределенного контроля с экспоненциальной (или близкой к ней) весовой функцией. При этом можно надеяться на уп рощение и обеспечение технической реализации управляюще го устройства. Посмотрим, в какой мере это предположение оправдается.
из |
Пусть |
интервалы |
квантования |
|
Д t |
и Дли малы и г'-»-со, |
||||
(1- 155) |
находим |
|
|
г2 |
|
|
|
|
||
|
|
Р<*,\ - P i |
, |
|
|
|
||||
|
|
- |
- |
|
|
|
||||
|
|
|
|
|
г 00,1 |
|
|
|
|
|
|
|
РооЛ = -г^-т ( AI+P2 д 1+от)— Г—— , |
(1-160) |
|||||||
|
|
lah |
|
|
|
|
Р00,1 |
|
|
|
так |
как при этом Рп=Раз,\ ~Pi—1.1 |
• |
|
Решая |
уравнение |
|||||
(1. |
160) и учитывая |
условие р 00д > р 0, |
получаем |
|
|
|||||
|
т + Д х+/72Д X+ V (m + Д 1—р2 Д 1)2—4р2Д 21 |
8 |
|
|||||||
Рсо,1 — |
|
|
4°д2 |
|
|
|
|
W |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
(1. |
161) |
Для перехода к непрерывной системе |
введем допущение |
|||||||||
|
|
|
°л2 (2*ySh |
|
|
|
(1- |
162) |
||
|
|
|
|
Д х Д t |
|
|
|
|
|
|
При малых ДгУ+О можно записать |
|
|
|
|
||||||
|
|
Д1=- |
(2* r s h |
|
|
|
|
|
||
|
|
А х { М ) г2аю^ |
• |
|
|
|
||||
Используя |
(1. 161), |
получим |
|
|
|
|
|
|
||
|
|
г |
|
1 |
|
|
|
1 |
|
|
|
РОО,-1 |
т |
■а+1°д / |
|
\ + L ° h t ' |
|
|
|||
|
|
2 Д 7 |
|
|
|
|
|
|
94
где |
. |
. |
Гк A x2aial, 2 |
|
|
|
1К |
|
|
|
||
|
|
|
|
|
|
|
||||||
L ( к А * ) - у |
— 2д)55^- + а*1 ПРИ (s—Q & t< 0 > |
|
|
|||||||||
L°(l&.x) = |
I /~* А -*aigV -u.^ |
при (f— i'A ^ ) > — |
, |
(1. |
163) |
|||||||
|
г |
|
2-2S;i |
|
|
|
у |
|
|
|
|
|
|
|
|
|
к Axt=(s—/) А |
— . |
|
|
|
|
|
||
Аналогично |
вычисляется множитель |
|
|
|
|
|
||||||
|
|
|
|
aia2 |
|
|
|
|
|
|
|
|
|
|
|
М i = 2TCS/j(L(0).+ai) |
|
|
|
(1. |
164) |
||||
Формула |
|
для |
оптимального |
управления |
us* |
прини |
||||||
мает вид |
|
|
|
S— 1 |
I |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
»s*=-^(<7*s + t— M i K p S |
|
|
|
|
|
(1- |
165) |
|||||
|
|
|
|
i = l к = 0 j- j |
[ 1 + £ о ( Х ) д д |
|
|
|
|
|||
|
|
|
|
|
|
X—0 |
|
|
|
|
|
|
Проанализируем и преобразуем выражение для весовых |
||||||||||||
коэффициентов в (1- 165): |
|
|
|
|
|
|
|
|
||||
|
|
|
|
к + i — 1 |
|
1 |
|
/С+1— 1 -L°UW |
|
|
||
Я1!*,*)" |
П l + |
L |
° ( j ) A |
t |
П |
|
|
|
|
|||
|
|
|
|
i=0 |
|
|
|
/= о |
|
|
|
|
|
|
|
|
|
к + i — 1 |
|
|
|
|
|
|
|
|
|
|
|
=ехр{— |
|
Ц / ) Д ^ |, |
|
|
|
|
||
|
|
|
|
|
|
/= о |
|
|
|
|
|
|
|
|
|
|
|
|
|
*'+1>6 |
|
|
|
|
|
|
|
|
|
^ 1(x0)=exp{— J |
L\i\)dt\.\. |
(1. |
1 |
66) |
||||
Перейдем |
к |
непрерывному |
времени, |
устремив |
Ах |
и |
A t |
|||||
к нулю и заменив суммы интегралами: |
|
|
|
|
|
|
|
|
|
|
x + v O |
|
|
|
|
|
|
t la |
|
j |
L°(-q)df] |
|
“ * ( 0 S X |
|
|
|
о о |
|
0 |
( у(*,*-о~ |
|
|
|
|
|
|
|
|
|
|
|
—Ku(x)u(t—0— -£• ) \dxdO. |
|
(1. 167> |
|||||
Весовая |
функция |
УУ |
зависит |
|
от |
пространственной |
||
координаты х и времени 0 и |
представляет |
собой |
некото |
|||||
рую поверхность |
с |
максимумом |
в |
точке (х=0, |
0= 0). |
|||
Пусть |
|
|
fZg — Ol°V |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
2 rSSlt |
|
|
|
|
И Х + ц О > /1Г |
|
|
|
|
|
|
|
|
Тогда |
|
/ц |
|
jt-f-оО |
|
|
||
л+t/O |
|
|
|
|
||||
L°{t])df\= |
^ |
|
dt{ + |
^ |
fl?vj/as/H+ at2 = |
|||
О |
|
о |
|
|
/„ |
|
|
|
= - J r |
T/ M H + a i2)3- |
2 oi3 |
|
V «2/H+ Q'"l+ |
|
|||
3^2 |
|
|
||||||
|
|
|
|
|
|
|
|
-f(* + lO)VGj/jj+a2! = 5 1 + (x+ l0)52,
——B-Lx —B 2vO
|
|
g'(x, 0)= e |
-e |
■e |
. |
(1 . |
168) |
|
Перейдем |
к |
случаю, |
когда |
х + и 0 < /„. |
Этому условию со |
|||
ответствует |
наиболее важная |
часть |
весовой функции. |
|||||
Находим |
|
|
|
|
|
|
|
|
x + v b |
L*(r\)d-q= _2_ |
|
|
|
|
|
|
|
^ |
У(айх+сьВ-Ь а ^)3—Gj3 |
(1. |
169) |
|||||
О |
|
Зс |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Итак, |
при |
лН-Ы)</н |
нельзя |
представить |
g(x, О) |
|||
только от х, |
а другая — только |
от 0. |
Следовательно, |
УУ |
в виде произведения двух функций одна из которых зависит нельзя представить в виде двух отдельных блоков: уст
96
ройства распределенного контроля и устройства временной обработки информации. Управляющее устройство является динамической системой с распределенными параметрами. Этот результат имеет качественный принципиальный ха рактер.
Найдем приближенно оптимальное, более простое уп равление. Напомним, что 0<Z.°(y|)<L°(L) и заменим L°(vj) константой L°(fi*) = D, где 0<т]*</.. Тогда из Ц. 167) легко получаем
|
|
|
__D_ |
||
«*(*) = Ка Ч*У+ Ч ) - М 1К„ \е |
—Ш |
■v |
|||
е |
'[у(х,*?—0)— |
||||
|
|
О |
О |
|
|
Кахи |
t—0- |
л: |
dx 1 dQ |
(1. 170) |
|
v |
|||||
|
|
|
|
Для приближенной реализации закона управления, ока зывается, необходимо иметь устройство распределенного конт роля с экспоненциальной функцией веса и фильтр — аперио дическое звено. При/Гцд.= 1 сигнал u(t) подается на вход фильтра через блок с передаточной функцией
JL . \ |
V |
—р?и |
0 W ~ D + P |
>- |
|
где |
lHD |
|
сс=е |
|
АГ*
Значит, высказанное выше предположение о возможности представления УУ в виде совокупности устройств распреде ленного контроля и фильтра с сосредоточенными параметра ми справедливо лишь для субоптимальной системы.1
1.5. ДУАЛЬНОЕ УПРАВЛЕНИЕ МНОГОМЕРНЫМИ ОБЪЕКТАМИ
СЧИСТЫМ ЗАПАЗДЫВАНИЕМ
1, 5. 1 . Алгоритмы управления
В этом подразделе мы приведем полученные в замкнутой форме общие алгоритмы управления многомерными линейны-
7 |
2247 |
97 |
ми объектами е запаздыванием. Рассмотрим случаи, когда запаздывания по каналам управления одинаковы или различ ны, а возмущение представляет собой векторную случайную величину с коррелированными компонентами или марковский векторный случайный процесс.
Математическая модель объекта и измерительного уст ройства имеет вид
экв
<7[s]r=CP p[s]+/4l/[s—т] = |х [s]~t~A l/[s—т],
1/[s]= F (k[s],s), |
(1, 171) |
y[s] = G?[s]+£A[s], |
8=1,2,... |
где —/-мерные векторы типа
*ф] = || ui [s]...M;[s] || г
l/[s—1 ]= || "yjs— x j,...Vils—X/] || r ,
y[s],/i|sl—тг]— мерные |
векторы, |
размерность |
вектора |х |
равна д; |
|
нелинейный |
оператор; |
F(-)— взаимнооднозначный |
|||
А —(/Х/)-матрица, |
В—(7]Х>]) = матрица, | В | =£0, |
Ср. —/ХЛ)-матрица, G—(гХ/) — матрица.
Матрицы в общем случае зависят от времени s.
Блок-схема объекта управления показана на рис. 1. 11.
Блоки, характеризующиеся матрицами А, |
, G, В и вектор- |
|||
функцией F, обозначены соответствующими буквами. Вектор |
||||
ные сигналы отмечены двойными |
стрелками, п, тг, |
..., п — |
||
звенья с |
чистым запаздыванием Т1< Т г < ...< т /. |
Вначале |
||
считаем, что все запаздывания xt |
одинаковы: т /= т , / = 1 , ..., /, |
|||
a p[s] = p |
— векторная случайная величина. |
|
||
Значения вектора помехи /i[s] |
в различные моменты вре |
мени независимы, но составляющие вектора h[s] могут быть между собой коррелированы. Закон распределения р и h — нормальный- с плотностями вероятности:
|
Рис. 1.11 |
|
|
P(A[s])~A/(0,q 7 ) |
|
|
|
гили |
|
|
|
P{h[s\)=yf |
* 1 e'cpj- y At[s]Q aA[s] j, |
(l. |
172)V |
|
V |
(1. |
173) |
—l —1
•где Qh ,Q (j.— ковариационные матрицы размера (tjX'»]) и <ЛХД); р и Л Is] —статистически независимы.
Критерий качества квадратический общего вида:
rs+T=Hf{^r[C?(9*-^[s+x])(7*-<7[s+T])r] | и [«-1], |
(1. 174) |
|
y|s—1]}=Л4{(?* —^[s+ x])7" |
—ф + т ] ) | / s |
!} |
с симметрической матрицей весовых коэффициентов Сд=СдТ. Буквой Is- 1 обозначена условно информация об объекте, на-
99