Файл: Красовский, А. А. Фазовое пространство и статистическая теория динамических систем.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',