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

книги из ГПНТБ / Красовский, А. А. Фазовое пространство и статистическая теория динамических систем

.pdf
Скачиваний:
35
Добавлен:
19.10.2023
Размер:
7.02 Mб
Скачать

170

РА В Н О В ЕС Н Ы Е Р А С П Р Е Д Е Л Е Н И Я В СИСТЕМАХ

[ГЛ .

IV

следящей системы, описываемой уравнениями

 

 

+ a u zi + «12^2 =

+ « i + I i , х 2 — Хх = 0,

(4.32)

где знак

минус — при хй

0 и знак плюс — при х г <

0.

Роль шума | х здесь может играть случайный момент на валу двигателя следящей системы.

Если 4 = а г = 0, 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, п

г

 

+ 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о • • •

Соседние файлы в папке книги из ГПНТБ