Файл: Красовский, А. А. Фазовое пространство и статистическая теория динамических систем.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 23.10.2024

Просмотров: 89

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

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',