книги из ГПНТБ / Красовский, А. А. Фазовое пространство и статистическая теория динамических систем
.pdf170 |
РА В Н О В ЕС Н Ы Е Р А С П Р Е Д Е Л Е Н И Я В СИСТЕМАХ |
[ГЛ . |
IV |
|
следящей системы, описываемой уравнениями |
|
|
||
+ a u zi + «12^2 = |
+ « i + I i , х 2 — Хх = 0, |
(4.32) |
||
где знак |
минус — при хй |
0 и знак плюс — при х г < |
0. |
|
Роль шума | х здесь может играть случайный момент на валу двигателя следящей системы.
Если 4 = а г = 0, 4° 4 2) ф 0, т. е. релейное
управление присутствует в уравнении, в котором нет шу ма, то согласно (4.31) условие существования кусочно нормального распределения имеет вид
«1 = — х — • <1ц + а22
2. Система с переменной структурой второго порядка с одной линией переключения описывается уравнениями
|
Х\ + |
йиХ1+ ап х2 = ёт, |
|
Х2 |
4l*^l ' 1,"Л22Х2 — 0 |
при ахг |
х% 0, |
(4.33) |
|
|
|
|
Х1+ аПХ1 + а<12 — ^2, |
|
|
&2 |
®21^1 "Ь 4ю — 0 |
при axi + |
х2<[ 0. |
|
Так же как в предыдущем случае, определяем коэффициенты с двумя индексами:
Л(V |
„(v) -L М |
i(4vi)+ 4 v2,)42) |
|
О И ^ й22 ^(v) __ |
|
||
Л11 |
|
M V |
|
|
|
(4.34) |
|
|
KV4V + (4V)2- 4V4V] (4V+4V) |
||
4 V = |
|
||
|
■sn («£«)* |
2. |
|
|
v = 1, |
|
|
Последние две группы условий существования (4.23) для данного случая имеют вид
— 2а1Л(112) + ^ 2 2 = 4 i — 2а 1А п |
+ А 22(2) 9 |
|
a i [4V —ai4V + — ^ii(4i —ai-4i2)j+ |
(4.35) |
|
+ ( 4 V — a 14 V ) = a j [ 4 21 — « 1 4*2 |
||
+ |
||
, 2) |
„(2К |
+ — S n (^n —®1-4а) I+ (4V —®4i b)'
§ 4.2] Ф П К -У Р А В Н Е Н И Е И С И Н ТЕ З СИСТЕМ У П Р А В Л Е Н И Я 171
Подставляя сюда выражения (4.34), после преобразований получаем
2d! [a^a&ViV-+ Оаа)(ваУ— |
+4а) ( ай) г]= |
|
= (4 1))2 (ail4- Ли ) [ЛиЛм—afMV—(ajg)8—(4i )2] — |
|
|
— (оЙ))*(аЙ)+ в « ))[виви — л'иЛи — (Ли )2 — (л^)2], |
(4.36) |
|
д ( ! ) |
а(2) |
|
31 |
-‘г! |
|
= afi — аЧ1
В отличие от рассмотренной релейной системы, условия (4.36) существования кусочно-нормального равновесного распределения в СПС включают два соотношения. Усло виям существования (4.36) можно удовлетворить, подби рая коэффициент наклона ах прямой переключений и обеспечивая дополнительно одно соотношение между ко эффициентами уравнений в обоих полупространствах.
Если = fl4i\ а выражение в квадратных скобках второго из условий (4.36) отлично от нуля, то условия существования упрощаются: oti=0,
(«и + «и) [л^ли — а$а!$ — (д&У — (ац)2] =
—(Ли 4" а<2ъ) [а П а<21 — ЛцЛад — (Л м )2 — (Ли )2]-(4-37)
Вданном случае линия переключений совпадает с осью %. При других предположениях о коэффициентах уравнений в обоих полупространствах из (4.3б) можно получить большое число частных форм условий существования рав новесного кусочно-нормального распределения в рассма триваемой системе.
§4.2. ФПК-уравнение и синтез систем управления
ифильтрации
Втеории аналитического конструирования систем управления доказаны следующие положения [4.6] — [4.8]. Для объекта, описываемого уравнениями! |
*1 + f t К» •••» *п» 0 = ««. *= 1, 2, |
(4.38) |
172 |
РА В Н О В ЕС Н Ы Е РА С П Р Е Д Е Л Е Н И Я В СИСТЕМАХ [ГЛ . IV |
где ut — синтезируемые управления, могущие быть функ циями всех фазовых координат xL, . . ., хп (полная сте
пень наблюдаемости), оптимальными в смысле минимума функционала
I = V3[Xj (t%), •. ., хп(^2)] -Ь
tг |
t, п |
г |
|
+ j e ( * x , . . . , * n,f ) * |
+ 4 - j 2 |
( т - ) dt, |
(4.39) |
где Va — заданная функция значений фазовых координат в момент времени t — t2^> tx, Q — заданная функция фазовых координат и времени, kt —- заданные коэффици
енты, являются управления
|
|
|
= |
|
|
|
|
|
(4.40) |
где V (а^, |
. . ., хп, t) — решение |
уравнения |
в |
частных |
|||||
производных |
|
|
|
|
|
|
|
|
|
dt |
f |
dV |
1 |
v |
/ь |
|
- Q |
|
(4.41) |
Zj |
— |
2 |
Zj \^ д х , |
|
|||||
|
i=l |
* |
|
i=l \ |
1 / |
|
|
||
при граничном условии |
|
|
|
|
|
|
|
||
|
V (хц . . ., xn, £a) |
= |
V3 (x*, |
. . ., xn). |
|
(4.42) |
|||
Далее, для того же объекта (4.38) |
оптимальными в смысле |
||||||||
минимума |
функционала |
|
|
|
|
|
|
|
|
I = V3[хх (г2), . . . , х п (*г>] + |
j |
Q (хь . . . , x n,t)dt + |
|
||||||
+ х |
j (**■) ^ + |
|
|
f( k i Щ ) |
d t |
(4,43) |
|||
являются управления (4.40), где V, однако, определяется
линейным уравнением в частных производных
п
(4.44)
при том же граничном условии (4.42).
§ 4.2J Ф П К -У РА В Н Е Н И Е И С И Н ТЕ З СИСТЕМ У П РА В Л Е Н И Я 173
Синтез систем управления. Сопоставим уравнение (4.44) с ФПК-уравнением (1.22) для свободного движения (uj = 0) того же объекта:
д Inр |
71 |
71 |
|
V 4 д In Р |
VI j |
(4.45) |
|
|
|
|
Пусть для момента t = t2 задана логарифмическая плот
ность распределения
In /> (* !,.. ., хп, g = In р3 (Хи . . . , Х п ) -
= - М * 1 . • • М*п) (4-46)
в виде отрицательно определенной функции, могущей служить оценкой (нормой) приближения к невозмущенно му состоянию и удовлетворяющей условию нормировки. Пусть определено решение ФПК-уравнения (4.45) при условии (4.46). Тогда, как видно из (4.42)—(4.46), спра ведливо следующее положение.
Если для неуправляемого объекта (м,- = 0) без шумов известно текущее распределение плотности вероятности в фазовом пространстве р (хг, . . ., хп, t), обращающееся в
рг (а*, . . ., ж„) = |
exp [—F3 (®г, |
. . ., хп)) |
|
|
в момент времени t = t2, |
то оптимальными в смысле ми |
|||
нимума функционала |
U п |
|
|
|
I — Ув lxi (h)i ■• • >Х„(<2)] + |
~дх~^ |
|
|
|
J 2 |
|
|
||
|
+ |
? № |
) ■ * |
<«•«> |
являются управления |
|
|
|
|
Ui==A ? i i ^ . |
|
(4.48) |
||
Функционал (4.47) является полуопределенным в том смысле, что назначаемыми по желанию конструктора яв ляются составляющие Уя и
174 |
РА В Н О В ЕС Н Ы Е РА С П Р Е Д Е Л Е Н И Я В СИСТЕМАХ [ГЛ . IV |
Составляющая
определяется |
характеристиками |
/ г исходного |
объекта, а |
||||
составляющая |
|
|
|
|
|
|
|
71 |
tt |
|
lt n |
d In |
y |
|
|
12_ 2J(*, |
d In p \ 2 |
Ш * . |
d t |
||||
dx{ |
) |
||||||
i=l (i N |
/ |
u i= l \ |
|
|
|
||
раскрывается после решения задачи об определении лога рифмической плотности вероятности в фазовом простран стве неуправляемого объекта. Однако сумма
Л
заведомо является неотрицательной функцией фазовых
координат, а сумма 2 |
для объектов с устойчивым |
i=i |
xi |
невозмущенным состоянием чаще всего также неотрица
тельна. Об этом свидетельствует, |
в частности, |
выражение |
|||||
(1.19) для производной общей энтропии, |
которое для сво |
||||||
бодного движения |
неуправляемого |
объекта |
имеет |
вид |
|||
|
|
. п |
|
|
|
|
|
|
|
н = —м| |
|
|
|
|
|
|
|
|
|
XI |
а/{ |
= 0 и мини- |
|
Для обобщенно консервативных систем 2 j |
i |
||||||
|
|
|
|
t=i |
|
|
|
мизируемый функционал принимает форму |
|
|
|||||
I = |
V3[xx (f2) , .. .,xn (t2)]+ |
|
|
|
|
|
|
|
+ - г i |
!(£■ )*+ |
1=1 f, ' |
* |
/ |
|
|
|
i=l i, |
' * / |
|
||||
Для |
линейных систем величина |
2 |
-зг- |
постоянна |
или |
||
|
|
|
7=1 |
дх. |
|
|
|
является заданной функцией времени и наличие соответ
5 4.2] Ф П К -У РА В Н Е Н И Е И С И Н ТЕ З СИСТЕМ У П РА В Л Е Н И Я 175
ствующего члена в минимизируемом функционале никакой роли не играет, так что этот функционал также сводится
к(4.49).
Вработах [4.7]—[4.9] практическая процедура анали тического конструирования рассматривается как итера ционная. При таком подходе сначала на основе техниче ских соображений ориентировочно задается назначаемая часть функционала и определяются оптимальные для дан
ного функционала управления. Далее тем или иным путем, но обычно посредством интегрирования уравнений зам кнутой системы на ЭВМ определяются процессы в синте зированной системе. Если эти процессы оказываются в каком-либо отношении неудовлетворительными, произ водится более или менее целеустремленная трансформа ция назначаемой части функционала и вновь осуществля ется синтез. Процедура повторяется до получения прием лемых результатов. При такой трактовке задачи аналити ческого конструирования синтез на основе минимизации функционала (4.47) оказывается целесообразным.
Для получения оптимальных управлений в явной фор ме необходимо решить уравнение (4.45) при граничном условии (4.46). Если неуправляемый объект (ut = 0) от
носится к классу систем без прямых нелинейных связей и известна полная система первых интегралов для этого объекта, то согласно (1.25) общее решение уравнения (4.45) имеет вид
U п
lnp = Y (фх, • • •, Ф„) — f 2
( i=i i
где произвольная функция определяется из условия
¥ (Фи • • • , Фп)м. = — Уг (*1>• • •, ®п). |
(4.50) |
Таким образом, в данном случае оптимальные управления равны
, _ |
ь» ЗУ |
/ ЗУ |
Зт]д |
ЗУ di|;n |
) . (4.51) |
1 ~ |
1 дхА |
\ 3^i |
dxi ‘ |
Зфп dxi |
Если для объекта того же класса известно р <^п линейно
независимых первых интегралов, но условие
¥ (фх, . . ., фр)(=(1 = - 7з(Жх, . . ., хп) |
(4.52) |
176 |
РА В Н О В ЕС Н Ы Е РА С П Р Е Д Е Л Е Н И Я В СИСТЕМАХ 1ГЛ. IV |
может быть удовлетворено подбором функции Y, то опти мальные управления имеют вид
\ |
дх.г + |
■ д'У дфр \ |
|
+ 5г|>р дх{ ) |
' |
Если первые интегралы не зависят от времени, что бывает для обобщенно консервативных систем, то при удовлетво рении граничного условия (4.50) или (4.52) равенство Y = — V3 имеет место при любом t и управления равны
В этом случае весь синтез управлений сводится к назна чению функции Ляпунова V3.
Для случая линейного исходного объекта согласно (1.53) решение уравнения (4.45), удовлетворяющее усло вию (4.50), имеет вид
In р =
n n п
= |
— Vs[ 2 u>u(h,t)xk, •. •, 2 ^nk(h,t)xk] — j |
2 akkdt- |
|||||||
|
|
k= l |
|
|
|
fc=l |
|
f |
1 |
Таким |
образом, |
оптимальные |
в указанном |
смысле |
|||||
управления для линейного объекта равны |
|
||||||||
|
|
и |
|
|
|
|
•• |
|
|
|
h |
F3[ 2 |
wlk {tit t) xkt • • • i 2 |
(tit t) ^-Jt] • (4.53) |
|||||
|
|
k—l |
|
|
|
k=l |
|
|
|
Если |
V3 |
задана |
в |
виде |
канонической |
квадратичной |
|||
формы |
|
|
|
п |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
У з |
— 2 |
Р р р х р |
> |
|
|
то |
|
|
|
|
р=1 |
|
|
|
|
|
|
п |
|
п |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
щ — — 2 & i |
2 |
х к |
2 wPi {t^ t) Р ppWpk {tit t). |
( 4 . 5 4 ) |
||||
|
|
|
k=X |
p=l |
|
|
|
|
|
Если рассматривать всю совокупность нелинейных динамических систем, то знание первых интегралов яв
§ 4.2] Ф П К -У РА В Н Е Н И Е И С И Н Т Е З СИСТЕМ У П Р А В Л Е Н И Я |
177 |
ляется редким исключением. В общем случае, когда эти интегралы неизвестны, для отыскания решения уравне ния (4.45) можно использовать изложенный выше метод рядов. Уравнения для коэффициентов логарифмической плотности вероятности и выражения этих коэффициентов
через квадратуры весовых |
функций линейного |
прибли |
|||||
жения представлены выше (см. (2.19), (2.23)). |
|
|
|||||
Особенность |
рассматриваемой |
задачи заключается в |
|||||
том, что граничное условие (4.46) |
записано для |
конечно |
|||||
го момента |
времени t = t2. |
Будем полагать, что фун |
|||||
кция |
V3 задана в виде некоторого полинома или сходя |
||||||
щегося |
во всей |
рассматриваемой области фазового |
прос |
||||
транства ряда |
|
|
|
|
|
||
|
|
П |
|
П |
|
|
|
V3 — |
2 |
PiifX&k Н |
2 |
РiklmxixHxlxm "Ь • |
• • |
||
По условию функция V3 должна быть положительно опре
деленной, поэтому в разложении присутствуют лишь члены четных степеней. Коэффициенты pih, рШ т, . . ., симмет
ричные относительно всех индексов, считаются задан ными величинами. Тогда, в соответствии с (2.8), (2.19), (4.46), коэффициенты ряда логарифмической плотности вероятности получаются как решение бесконечной систе мы линейных уравнений (2.19) при граничных условиях следующего вида:
(4.55)
Интегрирование уравнений при подобных граничных усло виях удобно производить в обратном времени т = —t. При этом момент t2 считается за начальный. Для устой
чивых по линейному приближению объектов это выгодно с точки зрения сходимости метода. Соответствующие со ображения приведены в главе II.
Для выражения коэффициентов оптимальных управ лений через квадратуры весовых функций линейного приближения можно использовать выражения (2.23), за менив в них момент времени 0 на t2 и начальные значения
РА В Н О В ЕС Н Ы Е РА С П Р Е Д Е Л Е Н И Я В СИСТЕМАХ [ГЛ . IV
178
Л0у(л ...х на Получаем
|
|
|
п |
U |
п |
|
|
|
|
|
|
1 |
|
At (f) = |
— 2 2 |
|
J 2 |
арр'> |
|
>t)dt , |
|
||||||
|
|
|
|
||||||||||
|
|
|
v = l |
t |
p = l |
|
|
|
|
|
|
||
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
AV(0 = |
— 2 |
|
Pv^vi (^. t) |
(<2, t) — |
|
|
|||||||
|
|
|
v ,p = l |
n |
|
It |
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
— 2 |
2 |
|
j* [з |
2 |
app*v-(O ~b |
|
|
|||
|
|
|
n |
v ,p = l |
t |
P = 1 |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+ |
2 |
aPW(O |
(*')] |
(*'. t) w\4 (*'. t) dt'> |
|||||||
|
|
|
p=1 |
|
|
|
|
|
|
|
|
|
|
Aijk...t (t) = |
|
|
|
|
|
|
|
|
|
|
|
||
'~N~' |
|
n |
|
|
|
|
|
|
|
|
|
|
|
— — |
|
Pvp.-.x^vi (^2> t) Wy.j{ti, t). . . |
t) |
|
|||||||||
|
2 |
|
|||||||||||
|
V - .....x—l |
I T |
|
|
|
|
|
|
|
|
|||
— |
|
2 |
|
J |
|
|
Ч- 1) 2 |
арр*Н"--.х ( о ”Ь |
(4.56) |
||||
|
V, (1 |
X=1 |
t |
|
|
|
P= 1 |
IV+2 |
|
|
|||
|
|
|
+ iv |
2 |
apvp—x ( o -^p( o |
+ |
|
|
|||||
|
|
|
|
P=1 |
N + l |
|
|
|
|
|
|||
+ 2 |
i a p j ^ x ( О л р ^ ( о + |
• • • + |
|
v ) A p x ( o h - |
|
||||||||
p = l |
|
N |
|
|
|
|
|
|
|
N |
|
|
|
+ |
^ |
|
•‘4 |
[a pt.--X ( O |
^ p v p ( O |
+ |
• • • |
|
|
||||
TV— 1 |
|
s-'T' |
|
|
|
|
|
|
|||||
|
|
|
P=1 |
N —1 |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
• • • Ч- flpv...ti |
) -^ppx (01 “b ■• |
• |
|||||
• • . + |
yy 2 |
i~ |
2 |
|
l flppx ( ^ ) |
^ PVp-jJ ( O "b |
• • • |
|
|||||
|
|
|
|
|
P = 1 |
|
|
|
|
|
|
|
|
|
|
. . . -j- |
Hpvp (£ ) |
-4реф...х |
)]} X |
|
|
||||||
X W4i (t',t) Wy.j (t\t) . . . U>XP (<', <)df',
§ 4.2] Ф П К -У Р А В Н Е Н И Е И С И Н ТЕ З СИСТЕМ У П Р А В Л Е Н И Я |
179 |
где Ь изменяется от <; t2 до t2- Сами оптимальные в
смысле минимума функционала (4.47) управления имеют вид
П |
|
щ —— (А\ -(- 2 AiicXft |
2 Aftix^Xi - ) - ...) (4.57) |
Jt=1 |
К1=1 |
Этот ряд при t = t2 по условию сходится.
Для объектов без прямых нелинейных связей и обоб щенно консервативных объектов выражения для коэффи циентов оптимальных управлений принимают вид
A i |
( t ) =О, |
|
|
|
|
|
|
|
|
|
П |
|
|
|
|
|
|
АН(0 = |
— 2 |
Р**»* {h, t) |
(<г, t) , |
|
||||
|
|
v,H=l |
|
|
|
|
|
|
|
|
n |
(| |
n |
|
) A p v ( t r) + |
|
|
А ц к |
(0 = — 2 |
j*{ 2 |
|
|
||||
|
|
v,9,t=l |
t |
p=1 |
|
|
||
|
+ |
flpvt (f) Ap\x (t') + |
flpvji (t') Ape (<')]} X |
|
||||
|
|
|
X Ulvi (t\ t) |
(t\ t) Wek (t\ t) dt', |
||||
Aijkl (t) — |
|
|
|
|
|
|
||
= |
— |
2 |
|
|
(t2, t) W p . j ( t 2 , t) Wtk(t2, t) X |
@.58) |
||
|
v, 9.,*, 11= 1 |
n |
|
f, |
n |
|
|
|
x Wni ( t2, t ) |
|
oo A p t ( t )-(-... |
||||||
2 |
|
J { 2 |
||||||
|
|
«.И1.* .4 = 1 |
t |
P=1 |
|
|
||
|
|
...- ( - |
apvp.t (t') App (<')] + |
|
||||
|
n |
|
|
|
|
|
|
|
+ |
2 |
[apen(0 -4pvn(0 +•••+ ЯрмЦ (^0 ^4рец(<')]) X |
||||||
|
p=i |
(t\ t) ww (t’, t) W t k |
(t’, t) wnl (f, t) dt’, |
|||||
|
|
X W t i |
||||||
Если невозмущенное состояние исходного |
объекта при |
|||||||
ut = |
0 |
асимптотически |
устойчиво по Ляпунову, т. е. |
|||||
wik (t2, t) -*■ 0 |
при t -> |
00, то для ограниченных по моду |
||||||
лю |
коэффициентов |
|
|
|
|
|||
&ijkо • • •
