Необходимые условия оптимальности второго порядка в одной задаче управления с переменной структурой Текст научной статьи по специальности «Математика»

ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА 2018 Управление, вычислительная техника и информатика № 43

УПРАВЛЕНИЕ ДИНАМИЧЕСКИМИ СИСТЕМАМИ

УДК 517.977

DOI: 10.17223/19988605/43/1

К.Б. Мансимов, А.А. Аликберов

НЕОБХОДИМЫЕ УСЛОВИЯ ОПТИМАЛЬНОСТИ ВТОРОГО ПОРЯДКА В ОДНОЙ ЗАДАЧЕ УПРАВЛЕНИЯ С ПЕРЕМЕННОЙ СТРУКТУРОЙ

Рассматривается одна задача оптимального управления, описываемая совокупностью дифференциальных и интегральных уравнений. Установлены необходимые условия оптимальности первого и второго порядков. Ключевые слова: дифференциальное уравнение; интегральное уравнение типа Вольтерра; необходимое условие оптимальности; аналог уравнения Эйлера; аналог условия Лежандра-Клебша; необходимое условие оптимальности второго порядка.

На практике очень важным вопросом является адекватное описание изучаемого процесса. Во многих случаях это приводит к необходимости изучения многоэтапных процессов (составных или ступенчатых), которые очень широко распространены на практике. Такие процессы возникают в космонавтике, теории локомоционнных процессов, химической технологии и др. [1-4].

В работах [1-10] и других изучен ряд задач оптимального управления многоэтапными процессами, описываемых на различных отрезках времени (или же в различных областях) разными дифференциальными (разностными) уравнениями. Подобные задачи оптимального управления называются также задачами оптимального управления с переменной структурой. В предлагаемой работе исследуется одна многоэтапная задача оптимального управления, описываемая совокупностью дифференциальных и интегральных уравнений. При предположении открытости области управления установлен аналог уравнения Эйлера [11-13]. Выведены необходимые условия оптимальности второго порядка.

1. Постановка задачи

Допустим, что на заданном отрезке времени T = T ^ T (T = [t0, t1 ], T = [t,12 ]) управляемый

процесс описывается системой уравнений

t

x(t) = j f (t, s, x(s),u (s))ds, t e T, (1)

t0

y = g(t, y, v), t e T2, (2)

У (ti ) = G(x(t)). (3)

Здесь f (t, s,x,u) (g (t,y,v)) - заданная n (га)-мерная вектор-функция, непрерывная по совокупности переменных вместе с частными производными по (x, u) ((y, v)) до второго порядка включительно, G (x) - заданная дважды непрерывно дифференцируемая га-мерная вектор-функция, t0, ^, t2 заданы, причем t0 < tj < t2, u (t) (v (t)) - г(д)-мерный кусочно-непрерывный (с конечным числом точек разрыва первого рода) вектор управляющих воздействий со значениями из заданного непустого, ограниченного и открытого множества U (V), т.е.

u ( t)eU с Rr, t g T = [t0, 4 ],

v (t )gV с Rq, t gT2 , t2 ]. Пару (u° (t), v° (t)) с вышеприведенными свойствами назовем допустимым управлением. Предполагается, что каждому допустимому управлению (u° (t ), v° (t )) соответствует единственное непрерывное (кусочно-гладкое) решение x° (t) (y° (t)) уравнения (1) (задачи Коши (2)-(3)).

На решениях системы (1)-(3), порожденных всевозможными допустимыми управлениями, определим терминального типа функционал

I (u,v)=ф (x (t1 )) + ф2 (y (t2 )) . (5)

Здесь ф1 ( x ), ф2 ( y ) - заданные дважды непрерывно дифференцируемые скалярные функции.

Допустимое управление (u° (t),v° (t)) , доставляющее минимум функционалу (5) при ограничениях (1)-(4), назовем оптимальным управлением, а соответствующий процесс (u° (t),v° (t),x° (t),y° (t)) - оптимальным процессом.

Нашей целью является вывод необходимых условий оптимальности при предположении открытости области управления.

2. Вариации функционала и неявные необходимые условия оптимальности

Пусть (u°(t ), v°(t ), x°(t ), y°(t )) - фиксированный допустимый процесс. Через

(u (t) = u° (t) + Au (t ), v (t) = v° (t) + Àv(t ), x (t ) = x° (t) + Àx (t ), y (t) = y° (t ) + Ay (t )) обозначим произвольный допустимый процесс и, используя формулу Тейлора, запишем приращение критерия качества

(5):

AI (u 0, v°) = I (u, v )-1 (u 0, v°) = ^ ( x ( ti ))-ф1 ( x°(ti ))] + [ф2 ( y (t2 )) Ф2 ( y °(t2 ))] = ^Ф1( x°(ti )) 1 02ф1 ( x° (ti )) 0ф2 (y°(t2 ))

=-^^ Àx ( ti ) + ^ Àx'(ti )-Àx ( ti ) +-Ày (t2 )+ (6)

+1 Ày'(t2 Ày (t2 ) + Oi (|Àx (ti )||2 ) + 02 (|Ày (t2 )||2 ).

Здесь и в дальнейшем ||а|| есть норма вектора а = (а, а2,..., а„) , в Rn определяемая формулой

n

liai = SlaI, штрих для матриц - операция транспонирования, а для векторов - знак скалярно произве-

¿=1

дения, а о(р2) есть величина более высокого порядка малости, чем р2, т.е. о(р2)/р2 ^ 0 при р ^ 0 .

Далее ясно, что (Àx(t ), Ày(t )) будет решением задачи

t

Àx (t ) = J[ f (t, x ( 5 ), й ( 5 ))- f (t, x0( 5 ), u 0( 5 ))] ds , t G T, , (7)

*0

Aj(0 = g(t,y(t),v(t)) -g(t,y° (t),v° (t)), (8)

Ày (ti ) = G ( x (ti ))-G ( x° (ti )) . (9)

Предположим, что y0 (t) , p0(t) -пока неизвестные n и га-мерные вектор-функции.

Умножая обе части соотношения (7) ((8)) слева скалярно на (V) (р° (()), а затем интегрируя обе части полученного тождества по Т (Т ), будем иметь

}у°С)Ах(V)Л = }у°С) }[/(V,5,X(5),и(5))-/(V,х°(я),и°(5))]

¿0 ^ _

Ж,

(10)

(11)

Полагая N(р°,х) = р° (^ ) О(х) и используя формулу интегрирования по частям, в определенном интеграле получим

= (12)

Ч к

Учитывая тождества (10)—(12), из формулы приращения (6) получим

Ы(и°,у°) = аф'()) Ах(^ ) + 5ф2 (У ()) Ау (г2) + }у°'(^) Ах(V)Ж -

дх ду

-} }(5)_/(5,t,х (г),и (Г)) - /(5,Г,х° (Г),и° (Г))]Ж

¡0 _ í

дМ'(р°,х°к))

Ж + р°' (^2) Ау (^2)-

дх

- Ах(^) -(Г) -(Г)) -г^З'0 (0)]^ +

1 д2ф^ х °(г1 )) 1 д2 N (р°, х°( ^ ))

_1 ( *1)-Ь- Ах ( 1 )--- Ах (^ )"

+1 АУ (*2 )д 2ф2 (у °('2)) АУ (^ ) + °1 (||Ах (^ )||2 ) + °2 (|| Ау ( .

2 ),, ,-°3 (||АхС,

Из (7) ясно, что

ч

Ах(^) = }[/(*1,*,х(0,и (0)-/(х° (0,и°(0)] .

'o

Принимая во внимание (13), из (12) будем иметь

А1 (и °, ^^^^^^ [ / ((1, V, х (V), и (V))-/(Г1, V, х°(г), и ))] Ж +дф2 ( ^^2)) Ау (, 2 )-

(12) (13)

+КС)Ах(V)Ж-} }у°(5)[/(5,V,х(V),и (V))-/(5,V,х),и))]Ж

Ж + р° ) Ау (^ )-

Ч -

(14)

Ч

V (V))-£ (V, у 0('), у0^))] Ж + - Лх'(г1)

(х°(1))

дх2

Ах (^ )-

-1 Ах'( 1 )дМр°хМ Ах (^) +1 Ау'( ^ Ау (^ ) + °1 (||Ах (^ )||2) +

г2 ) )-°з (||Ах(V

2

2

Введем аналоги функции Гамильтона-Понтрягина:

Н (ч,Х,И, у0) =---/ г, X, и ) + (5 ) / (5, г, X, и ) Л +--^ / (Ч, г, X, и ) ,

М (г, у, у, р°) = р0' g (г, у, V). С учетом введенных обозначений формула приращения (14) примет вид:

М (и°, V0) = -|[н (г, X (г), И (г), (г)) - Н (г, х° (г) ,и ° (г), (г ))] -

0 ^ (15)

-}[М (г, у (г), V (г), р°(г))-М (г, у °(г), ^(г), р°(г))] Л +1 Ах'(ч Ах (ч) +

Ч

+1 Ду'(г2 Ф2 )) Ау (г 2)-1 А*'« ^(р^,2х°(г1)) Ах (г) +(г )Ау (г) Лч + р°'(г2 )Ау (г2)-

у <0

-}р°( г) Ау (г) Л + 5ф2 ( У0(<2 )) Ау (ч2 ) + °1 (|Ах (г )||2 ) + °2 (||Ау (ч2 )||2 )-°з (||Ах (г,4"2'

Из (15), используя формулу Тейлора, будем иметь

Ы(и°у)=—^——Ау(^) + р° (г2)ДУ^2)-)р° (04+ (г)Дх(г)Л-

ду к ч

Н'и (г, х° (г),и° (г), (г)) Аи (г)Л -1Н'х (г, х° (г),и° (г), (г)) Ах (г)Лг -

1 г1 г

- |[Ах'(г) Нхх (г, х° (г), и °(г), у°(г ))Ах (г) + 2 Аи'(г) Н„ (г, х°(г), и °(г), у°(г ))Ах (г)

2 ч

ч

+Ах'( г) НИИ (г, х° (г), и ° (г), (г)) Аи (г)] Л -1 М[ (г, у ° (г), V0 (г), р° (г)) Аv (г) Лг -

к

-{му (г,у° (г), V0 (г), р° (г)) Ау (г)Л -1 {[Ау'(г)Муу (г,у° (г),V0 (г),р° (г)) Ау (г) +

2

+2 Аv'(г)Му (г,у°(г),^(г),р°(г))Ау (г) + Аv'(г)(г,у°(г),^(г),р°(г))Аv(г)]Л +

+°1 (||Ах ( Ч1 )||2 ) + °2 (|| Ау ( Ч2 2 )-°з (|Ах (Ч1 2) + \ Ахи( Ч1 Ах (,1) +

1 д2 Ф2 (у °(г2)) 1 д2 N (р°, х°(<1 ))

+ 2 АУ (Ч2)-- АУ(Ч2)-1 Ах (г1)-^--Ах(г1 )-

-1 °4 ([|Ах (г )| + ||Аи (г )|]2) Лг-} °5 ([||Ау (г )|| + |^ (г )|]2) Л.

Предположим, что (у°(<), р° (г)) удовлетворяет соотношениям

у°(г) = Нх (г, х°(г), и° (г), у(г)), (17)

р°({) = -Му({,у°({)У({),р°({)), (18)

р -(,2 . (,9)

ду

Соотношение (17) есть линейное неоднородное интегральное уравнение типа Вольтерра относительно (г), а соотношение (28) - линейное однородное дифференциальное уравнение относительно

р° (г) . Систему (17)-(19) назовем сопряженной системой для задачи (1)-(5).

При выполнении соотношений (17)-(19) формула приращения (16) примет вид:

А/ (

(u°,v°) = -j[HU (t,x° (t),u° (t),(t)) Au (t)dt - _[M'v (t,y° (t),v° (t),p° (t)) Av (t)dt -

to tl

/ \\ „ / \\ „ i

i / ча2Ф1 (xo(tj )) , ч i , ча2Фг (yo(t2)) , ч i , y2N(p°, xo(t1)) / ч +J Ax'(tj) Ф1 ^2(1)) Ax (tj) + J Ay'(t2) Ф2 (y 2 (2)) Ay (t2)-J Ax'( tj)-(^ (l)) Ax (tj)-

1 4

- f[Ax'(t) Hx (t, x °(t), u°(t), (t)) Ax (t) + 2 Au'(t) Нш (t, x °(t), u °(t), y°(t ))Ax (t) +

2 to

+Ax'(t) Huu (t, x °(t), u°(t), V°(t)) Au (t)] dt - - j[Ay'( t) My (t, y °(t), v °(t), p°(t ))Ay (t) +

2 t

+2 Av'( t) Mvy (t, y °(t), v°(t), p°(t)) Ay (t ) + Av'(t) Mw (t, y °(t), v °(t), p°(t ))Av (t)] dt + ^ (Au, Av ),

(20)

где по определению

(Au, Av ) = °i (||Ax (tj )||2 ) + °2 (||Ay (t212 )-°з (||Ax (t- )||2 )-jo4 ([||Ar (t )|| + ||Au (t )||]2) dt -

to

-1°5 ([l|Ay (t )| + |Av (t )||]2) dt.

(21)

_,5|l НЛК/УЛИЛ-НЛК/УЛИ"12

t1

Из (7)-(9), используя условию Липшица, получаем, что

t

IAx(t| < L 1[||Ax(x)|| + ||Au (x)||]dx , (22)

to

t

|Ay (0|< L21 [|| Ay (x)| + ||Av (x)||] dx + L2\Ax (tj )||, (23)

to

где L — const > o, i — 1,3 - некоторые постоянные.

Применяя к неравенству (22) аналог леммы Гронуолла-Беллмана будем иметь

>1

I Ax (t )||< La A Au (x)|| d x . (24)

to

А из (23) следует, что

ч

||Ay (t )|| < Ls 1||Av (x)|| d x + L611 Ax (tj )||, (25)

h

где L4, L5, L6 - некоторые положительные постоянные. С учетом (24) из (25) получим

||Ay (t)|| < L5 )\|Av(x)||dx + L6 J||Au (x)||dx, (26)

t t где L — const > o - некоторое постоянное.

По предположению, U и V открытые множества. Поэтому специальное приращение управления (u° (t), v° (t)) можно определить как

AuE (t) = s5u(t), t e T, Avs(t) = s5v(t), t e T2,

где s - достаточно малое по абсолютной величине число, а 5u (t)e Rr, t e T (5v(t)e Rr, t eT2) -произвольная кусочно-непрерывная (с конечным числом точек разрыва первого рода) r (д)-мерная вектор-функция - допустимая вариация управляющей функции u° (t) (v° (t)) .

Через (Дхе (t), Aye (t)) обозначим специальное приращение траектории (x° (t),y° (t)), отвечающее приращению (27) управления (u° (t), v° (t)) . Из оценок (24), (26) следует, что

||A*s(i)||<L s, t eT, ||Ays( 1)||< h s, t e T2, где h, L = const > 0 - некоторые постоянные.

Используя оценки (18) и учитывая (27), при помощи (7)-(9) по схеме, например, из [14] доказы-

(28)

вается

Лемма 1. Для (ДхЕ ('),АуЕ (')) справедливы разложения

ДхЕ (') = в8х(') + °(в; '), Ау„(г) = в5у(г) + °(е; г), где (5х ('), 5у (')) - вариация траектории, являющаяся решением уравнения в вариациях

5х(') =}_/х (М^ (5),и° (^(5) + / (', 5, х° (5),и° (5))би (5)] , ' £ Т ,

'o

5у (О = 8у (*,/ (О,V0 + ^ (*,/ (()у (ОНО > т2,

8у(tl) = Ох (х*^))8у(tl). Учитывая (27) и разложения (29), (30) в формуле приращения (20) получаем, что

А/е (и°, v°)= I (и°+Аие, V0 + А^ )-1 (и °, v°) =

}Ни (V,х° (V),и° (V),(V))8и (V)Ж + }М'„ (V,у° (V),V° (V),р° (V))8v(V)Л

к

(29)

(30)

(31)

(32)

(33)

= -s

s2

+—

2

5x (tl )-- 5x(tl )-5X (tl )-^72- 5x(tl ) + 5У (t2 )---5У (t2 )-

8x

8x

8y2

(34)

-}[8х' (V) НхХ (V, х°( V), и °(' ), у° (V ))8х (V) + 2 8и' (V) Нш (V, х° (V), и °(V), у°(' ))8х (V) +

t0

+8х ' (') Нии (' , х° (' ), и °(' ), у°(' ))8и (')] Ж - }_8у ' (' ) Му (' , у °(' ), V °(' ), р °(' ))8у (' ) +

+28v'(')Му (г,у0 ('),V0 ('),р° ('))8у (') + 8v'(')Мт (',у° ('),V° ('),р° ('))8v(')]Ж] + °(в2).

Из разложения (34) следует, что первая и вторая вариации (в классическом смысле) функционала качества (5) имеют соответственно вид:

57(u°, v°; 5u,5v) = - |М (t,x° (t),u° (t), (t))Su (t)dt - |m; (t,y° (t), v° (t),p° (t))Sv(t)dt, (35)

д2 ф( х°(Ч)) д2 N (р°, х°(Ч)) 521(и°,V0; 5и,5v) = 5х'(Ч) 2 (1 П 5х(4)-5х'(ч)-(1''5х(<)-

-}[5х'(<)Нхх (<,х° (<),и° (<),(<))5х(<) +

<0

+25и'(<)Них (<,х° (<),и° (<),(<))5х(<) + 5х'(<)Нии (<,х° (<),и° (<),(<))5и (<)]Л + (36)

+5у'(<2)д2ф2 (у°(<2))5у(<2)-}[5у'(<)Муу (<,у°(<),),р°(<))5у(<) +

у <0

+25v'(<)Му (<,у° (<),V0 (<),р° (<))5у (<) + )(<,у° (<),V0 (<),р° (<))5v(<)]Л. Из результатов классического вариационного исчисления (см. напр.: [11-14]) с учетом (35), (36) следует, что вдоль оптимального процесса (и° (<), v0(<),х°(<),у°(<)) для всех 5и (<)е Яг, < е Т1,

5v(<) е Кг, г е Т выполняются соотношения

<1 <2

¡Ни и,х° (О,и° (*),(*))5и(ОЛ + |м; (*,у° (*),V0 (*),р° ^))5v(*)Л = 0, (37)

<1

5х (<1)-8x2-5х(<1 )-5х (<1)--5х(<1 ) + 5у (<2)--5у(<2)-

<1

-|[5х'(<)Нхх (<,х° (<),и° (<), (<))5х(<) + 25и'(<)Них (<,х° (<),и° (<), (<))5х(<) +

<0

<1

+5х'(<)Нии (<,х°(<),и° (<),у°(<))5и(<)]Л-|[5у'(<)Муу (<,у°(<),v°(<),р° (<))5у(<) +

<0

+25v'(<)М„у (<,у°(<),v0(<),р°(<))5у(<) + 5v'(<)М„ (<,у°(<),v°(<),р° (<(<)]Л > 0.

(38)

3. Необходимые условия оптимальности

Тождество (37) является неявным необходимым условием оптимальности первого порядка, а неравенство (38) есть неявное необходимое условие оптимальности второго порядка.

Используя произвольность и независимость вариаций 5и (<), 5v (<) управляющих функций

и° (<) и V0 (<), при помощи (37) доказывается

Теорема 1. Для оптимальности допустимого управления (и°(< ), v0(<)) в задаче (1)-(5) необходимо, чтобы выполнялись соотношения:

Ни (е, х°(е), м° (е), у°(е)) = 0, (39)

для всех точек непрерывности ее[<0, <1) управления и°(<);

Му (е, у° (е), V0 (е), р°(е)) = 0, (40)

для всех точек непрерывности ее[^,<2) управления V°(?).

Система соотношений (39), (40) есть необходимые условия оптимальности первого порядка и представляют собой аналог уравнения Эйлера.

Каждое допустимое управление (и° (<), V0 (<)) , удовлетворяющее необходимым условиям оптимальности (39), (40), назовем классической экстремалью. Ясно, что оптимальное управление (если оно существует) находится среди классических экстремалей.

Для сужения множества классических экстремалей, подозрительных на оптимальность, надо иметь необходимые условия оптимальности второго порядка, выраженные непосредственно через параметры задачи (1)-(5). С этой целью будем использовать неявное необходимое условие оптимальности второго порядка (38).

Пусть Я(ч,х) (п х п) матричная функция, удовлетворяющая матричным интегральным уравнениям типа Вольтерра [15-17]:

г

Я (', т) = } Я (', 5 ) / (х, х°(т), и °(т)) Лз + / (', х, х°(х), и °(х)),

х

г

Я (', х) = } /х (X, х° (х), и° (х))Я (х)Ж + /х (X, х° (х), и° (х)).

X

Через ¥(', х) обозначим (т х т) матричную функцию, являющуюся решением матричного дифференциального уравнения

¥ (', х) = -¥ (', х) Яу (х, у°(х), v0(х))

с начальным условием

¥ (г, г ) = Е,

где Е - (т х т) единичная матрица.

Решение 8х(') интегрального уравнения (31) допускает представление [15-17]

8х (' ) = |я (', х)

}/и (5,х,х° (х),и° (х))8и(5)Ж + /и и,х,х° (х),и° (х)) 8и(х)

Лх. (41)

Преобразуя правую часть представления (40), используя формулу Дирихле (см. напр.: [11]), получим

Ч

8х (г ) = } б (', х)8и (х) Лх, (42)

'0

где по определению

г

б (', х) = /и (г, х, х°(х), и°(х)) + } Я (', 5 ) /и( 5, х, х°(х), и°(х)) Ж .

х

А решение задачи Коши (32)-(33) допускает представление [11, 18]

г

8у(г) = ¥(',¡1 )8у(¡1) + }¥(',х)^ (х,у° (х),V0 (х))8v(х)Лх .

'1

Следовательно,

ч

8у (') = ¥(','1)Ох (х°('1 ))8х('1) + }¥(',х)^ (х,у°(х), v°(х))8v(х)Лх . (43)

'1

С учетом (42) представление (43) записывается в виде ч ч

8у(') = }¥(г,')Ох(х0(ч))б('1,х)8и(х)Лх + }¥(г,х) Vу°(х),^(х)^(х)Лх . (44)

г ч

Учитывая произвольность вариаций 8и (') и 8v ('), предположим, что 8v (' ) = 0. Тогда из представления (44) будем иметь

8у (') = }Ь (', х)8и (х)Лх , (45)

ч0

где по определению

Ь(',х) = ¥(','1)Ох(х°('1))б('1,х) .

А неравенство (36) примет вид

,^д2ф1 ( х°( '1))

8х'(ч )

4

8х (' ) -8х'(ч )

д 2 N ( р°, х°('1)),

8х ('1 ) + 8у'('2)

д2ф2 ( у °('2 ))

8у ('2)-

дх дх ду

}_8х'(г) Нхх (', х° (г), и ° ('), (г ))8х (') + 2 8и'( г) Н„ (', х° (г), и ° (г), (г ))8х (г)

+

'1

+8х'(г) Нии (', х° (г), и° (г), (г ))8и (г)] Л - }8у'(г) Муу (г, у ° (г), V0 (г), р ° (г ))8у (г) > 0.

'0

Используя представление (42), доказывается, что

д2ф, (х°)) ^ д2ф(х°))

8х'('1) ф1 а/ 1 8х (Ч ) = } } 8и' (х) б ( Ч,х) ф1 1)) б ()8и (5 ) Л х,

ч0 ч0

4

}8и'( г) Н„ (г, х° (г), и ° (г), (г ))8х (г) Л =

'0

4 ['1

= } |8и' (х) Нх (х, х°(х), и°(х), у°(х)) б (х, г) Лх 8и (') Л,

д22 (р°,х°(1 )) ч1ч1 д22 (р°,х°(1 )) 8х'('1)-((| ))8х(0 = }}8и'(х)б('1,х)-^ (| ))б('1,)8и(5)Лх ,

ч0 ч0

4

}8х'( г) Н^ (г, х° (г), и ° (г), (г ))8х (г) Л =

1 1

-}}8и'(х)

} б'(',х)Нхх (',х° ('),и° (г),(г))б(',5)Ж

X (х

8и ( 5 ) х,

Далее при помощи представления (45) доказывается, что

-2( у°)) '¡Л -2 у°°г.))

8у'('2) ф2(у2(2))8у('2) = }}8и(х)Ь('2,х) ф2 ^2 ( 2))Ь('2,5)8и(5)Лх,

У 4) 4) У

ч

}8у'(г) Му (г, у ° (г), V0 (г), р° (г ))8у (г) Ж =

4 4

'0 '0

}}8и'(х) }Ь'(г,х)Му (г,у° (г),V0 (г),р° (г))Ь(г,5)Ж

8и (5)х.

Введя обозначение

д2ф(х°(г. )) д2у° ))

К (х, 5 ) = -б'('1, х) ф1(х2(1 )) б ( '1,5 )-Ь'( '2, х)-^2^^ Ь ( '2,5 ) +

+ } б'(',х)Нхх (',х°(г),и° (г),у(г))б('^Лг -

X(х,5)

-б'( '1, х)

д2 N ( р°, х°( '1))

дх2

ч

б ('1,5) + }Ь'(г, х)Муу (г,у° (г), V0 (г), р° (г))Ь (г, 5)Л

и учитывая тождества (47)-(50), из неравенства (46) получим, что

}}8и'(х)К(х,5)8и (5)х + 2} }8и'(х)Нт (х,х°(х),и°(х),у°(х))Ь(х,г)Лх

8и (г) Ж +

+}8и'(г) Нш (г, х° (г), и° (г), (г ))8и (г) Ж < 0.

(47)

(48)

(49)

(50)

(51)

(52)

0 '0

Теперь предположим, что 5и (< ) = 0, 5v (< ) ^ 0. Тогда из представлений (42), (44) получим, что

5х(<) = 0, < еТ ,

г

5у (<) = |^(<, т)^ (т,у°(т),V0 (т)^(т)Лт .

При этом неравенство (36) примет вид д2Ф2 (у°(<2 ))

Sy'( t2 )-

dy2

Sy (t2) - 1[5y'(t)My (t,y° (t),v° (t),p° (t))Sy (t) +

+25v'(<)Му (<,у°(<),v0(<),р°(<))5у(<) + 5v'(<)М„ (<,у°(<),v°(<),р°(<))5v(<)]Л > 0. Используя представление (53), доказывается, что

5у (<2 )-Г1-5У (<2 ) =

ду

\ \ / ч д 2ф( у °(<2)) . .

= ¡¡5и'(т) ^ (т, у° (т), V0 (т)) ^'(<2, т) ^2 Р (<2,5)gv (5,у° (5), V0 (5))5и (5)Л5 ЛТ,

4 4

^ 2

1Sv'(t)My (t,y° (t),v° (t),p°(t))Sy (t):

: 1 JSv' (x)My (x, y° (x), v° (x), p° (x)) F(x, t)dx

,(t, y°(t), v°(t ))Su (t) dt,

h h

15y'(t)My (t,y° (t),v° (t),p° (t))Sy (t) — HSu'(x)g'v (x,y° (x),v° (x)) x

X

тах (т, 5

Введя обозначение

1 F'(t,x)My (t,y° (t),v° (t),p° (t))F(t,s)dt

52Ф2 ( y °(t2))

, (s,y° (s),v° (s))5u (s)dsdx.

M (x, s ) — F ' (t2, x)

F (t2, s ) + 1 F' (t, x)Myy (t, y°(t), v° (t), p° (t)) F (t, s ) dt

fy -< л

■s max (x, s)

и учитывая тождества (55)-(57), неравенство (54) записывается в виде

h h

11«v'(x) g; (x, y° (x), v° (x)) M (x, s ) gv (s, y 0 ( s ), v° (s ))Su (s ) dsd x +

tl tl

+21

(53)

(54)

(55)

(56)

(57)

1«v'(x) gv (x, y° (x),v° (x))My (x,y° (x), v° (x), p° (x))F (x, t)dx

+15v ; (t) Mv (t, y° (t), v° (t), p° (t ))Sv (t) < o.

,(t, y °(t), v°(t ))Sv (t) dt + (58)

Сформулируем полученный результат.

Теорема 2. Для оптимальности классической экстремали (и° (<),V0 (<)) в задаче (1)-(5) необходимо,

чтобы неравенства (52), (58) выполнялись для всех 5и (< )е Яг, < е Т, 5v (<) е Яг, < е Т2 соответственно.

Заметим, что полученный результат является довольно общим. Из него, используя произвольность вариаций 5и (<), 5v(<) управляющих функций и° (<) , v0(<) , можно получить ряд более легко

проверяемых условий оптимальности и, в частности, исследовать особые в классическом смысле [14, 19, 20] управления.

Заключение

В статье изучается одна задача оптимального управления, описываемая совокупностью дифференциальных и интегральных уравнений. При предположении открытости области управления установлен аналог уравнения Эйлера. Выведено конструктивно проверяемое необходимое условие оптимальности второго порядка.

ЛИТЕРАТУРА

1. Габелко К.Н. Последовательное улучшение многоэтапных процессов // Автоматика и телемеханика. 1974. № 11. C. 72-80.

2. Кириченко С.Б. Оптимальное управление системами с промежуточными фазовыми ограничениями // Кибернетика и си-

стемный анализ. 1994. № 4. C. 104-111.

3. Величенко В.В. Оптимальное управление составными системами // Доклады АН СССР. 1976. T. 176, № 4. C. 754-756.

4. Никольский М.С. Об одной вариационной задаче с переменной структурой // Вестник Московского университета.

Сер. Вычислительная математика и кибернетика. 1987. № 1. C. 36-41.

5. Ащепков Л.Т. Оптимальное управление разрывными системами. Новосибирск : Наука, 1987. 226 с.

6. Тадумадзе Т.А., Авалишвили Н.М. Регулярные возмущения в оптимальных задачах с переменной структурой // Оптималь-

ные задачи в системах с переменной структурой : сб. Тбилиси : Изд-во Тбилисского гос. ун-та, 1985. C. 100-154.

7. Харатишвили Г.Л. Принцип максимума в оптимальных задачах с переключением // Труды ИСУ АН ГССР. 1980. T. 19,

№ 1. C. 5-17.

8. Магеррамов Ш.Ф., Мансимов К.Б. Оптимизация одного класса дискретных ступенчатых систем управления // Журнал вы-

числительной математики и математияческой физики. 2001. № 3. C. 360-366.

9. Багирова С.А., Мансимов К.Б. Особые управления в одной ступенчатой задаче управления // Автоматика и вычислительная

техника. 2007. № 3. C. 74-81.

10. Мансимов К.Б., Насияти М.М. Необходимые условия оптимальности в одной многоэтапной дискретной задаче управления // Математическое и компьютерное моделирование. Сер. физико-математических наук. 2011. № 5. C. 162-179.

11. Алексеев В.М., Тихомиров В.М., Фомин С.В. Оптимальное управление. М. : Наука, 1979. 429 с.

12. Демьянов В.Ф. Условия экстремума и вариационное исчисление. М. : Высшая школа, 2005. 335 с.

13. Габасов Р., Кириллова Ф.М. и др. Методы оптимизации. Минск : Четыре четверти, 2011. 472 с.

14. Габасов Р., Кириллова Ф.М. Особые оптимальные управления. М. : Наука, 1973. 256 с.

15. Абдуллаев А.А., Мансимов К.Б. Необходимые условия оптимальности в процессах, описываемых системой интегральных уравнений типа Вольтерра. Баку : Изд-во ЭЛМ, 2013. 224 с.

16. Васильева А.Б., Тихонов А.Н. Интегральные уравнения. М. : Изд-во МГУ им. М.В. Ломоносова, 1989. 156 с.

17. Цалюк З.Б. Интегральные уравнения Вольтерра // Итоги науки и техники. Сер. математический анализ. 1977. T. 15. C. 131-138.

18. Габасов Р., Кириллова Ф.М. Оптимизация линейных систем. Минск : Изд-во Белорус. ун-та, 1973. 256 с.

19. Мансимов К.Б. Необходимые условия оптимальности особых процессов в задаче оптимального управления : автореф. дис. ... д-ра физ.-мат. наук. Баку, 1994. 42 с.

20. Мансимов К.Б. Особые управления в системах с запаздыванием. Баку : Изд-во ЭЛМ, 1999. 174 с.

Поступила в редакцию 2 октября 2017 г.

Mansimov K.B., Alekberov A.A. (2018) THE SECOND ORDER NECESSARY OPTIMALITY CONDITIONS IN A CONTROL PROBLEM WITH VARIABLE STRUCTURE. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie vychislitelnaja tehnika i informatika [Tomsk State University Journal of Control and Computer Science]. 43. pp. 4-15

DOI: 10.17223/19988605/43/1

The paper deals with the problem of optimal control by minimizing functional of a terminal type

I (u,у) = ф1 (X (t, )) + ф2 (y (t2 )) (1)

with constraints

u (t)eU с Rr, t eTj =[t0, t, ],

v(t)eV сRq, t eT2 = [t,t2],

t

X(t) = J f (t,s,x(s),u(s))ds, t e T, (3)

y = g(t,y,v),teT2 (4)

y (A) = G(x(t,)) . (5)

Here f (t, s, x, u), (g (t, y, v)) are give n (ra)-dimensional vector-functions, respectively, continuous with respect to all the variables together with partial derivatives with respect to (x,u) ((y,v)) up to the second order inclusive, G(x) is ra-dimensional twice continuously differentiate vector-function, t0, ^, t2 are given and t0 < ^ < t2, u(t) (v(t)) are r (g)-dimensional piecewise-continu-ous (with a finite number of points of discontinuity of the first kind) vector-functions of control actions, values from a given nonempty, bounded, and open sets U (V) , ^ (x) , (y) are given twice continuously differentiable in Rn (Rm) scalar functions. Here

^ (x) and (y) are given m-dimensional vector functions, respectively, continuous with respect to all the variables together with partial derivatives up to the second order inclusive.

Keywords: differential equations; Volterra type integral equations; necessary optimality conditions; analog Euler equation; analog Lejandr-Klebch conditions; second order necessary optimality conditions.

MANSIMOV Kamil Bairamali ogly (Doktor of Physics and Mathematics, Professor, Baku State University, Institute of Control Systems of Azerbaijan National Academy of Sciences, Baku, Azerbaijan). E-mail: [email protected]

ALEKBEROVAydyn Abdulla ogly. (Institute of Control Systems of Azerbaijan National Academy of Sciences, Baku, Azerbaijan). E-mail: [email protected]

REFERENCES

1. Gabelko, K.N. (1974) Consistent improvement of multi-stage processes. Automation and Remote Control. 11. pp. 72-80. (In Russian).

2. Kirichenko, S.B. (1994) Optimal'noye upravleniye sistemami s promezhutochnymi fazovymi ogranicheniyami [Optimum control of

systems with intermediate phase constraints]. Kibernetika i sistemnyy analiz - Cybernetics and Systems Analysis. 4. pp. 104-111.

3. Velichenko, V.V. (1976) Optimal'noye upravleniye sostavnymi sistemami [Optimal control of composite systems]. Doklady AN SSSR.

176(4). pp. 754-756.

4. Nikolskiy, M.S. (1987) Ob odnoy variatsionnoy zadache s peremennoy strukturoy [On a variational problem with a variable struc-

ture]. VestnikMoskovskogo universiteta. Ser. Vychislitel'naya matematika i kibernetika. 1. pp. 36-41.

5. Ashchepkov, L.T. (1987) Optimal'noye upravleniye razryvnymi sistemami [Optimum control of discontinuous systems]. Novosi-

birsk: Nauka.

6. Tadumadze, T.A. & Avalishvili, N.M. (1985) Regulyarnyye vozmushcheniya v optimal'nykh zadachakh s peremennoy strukturoy

[Optimum problems in systems with variable structure]. In: Optimal'nyye zadachi v sistemakh s peremennoy strukturoy [Optimum problems in systems with variable structure]. Tbilisi: Tbilisi State University. pp. 100-154.

7. Kharatishvili, G.L. (1980) Printsip maksimuma v optimal'nykh zadachakh s pereklyucheniyem [The maximum principle in optimal

switching problems]. Trudy ISU AN GSSR. 19(1). pp. 5-17.

8. Magerramov, Sh.F. & Mansimov, K.B. (2001) Optimization of a class of discrete step control systems. Zhurnal vychislitel'noy

matematiki i matematiyacheskoy fiziki - Computational Mathematics and Mathematical Physics. 3. pp. 360-366. (In Russian).

9. Bagirova, S.A. & Mansimov, K.B. (2007) Singular control in a single step control problem. Automation and Remote Control. 3.

pp. 74-81.

10. Mansimov, K.B. & Nasiyati, M.M. (2011) Neobkhodimyye usloviya optimal'nosti v odnoy mnogoetapnoy diskretnoy zadache upravleniya [Necessary optimality conditions in one multi-stage discrete control problem]. Matematicheskoye i komp'yuternoye modelirovaniye. Ser. fiziko-matematicheskikh nauk. 5. pp. 162-179.

11. Alekseyev, V.M., Tikhomirov, V.M. & Fomin, S.V. (1979) Optimal'noye upravleniye [Optimal control]. Moscow: Nauka.

12. Demyanov, V.F. (2005) Usloviya ekstremuma i variatsionnoye ischisleniye [Extremum conditions and calculus of variations]. Moscow: Vysshaya shkola.

13. Gabasov, R., Kirillova, F.M. et al. (2011)Metody optimizatsii [Methods of optimization]. Minsk: Chetyre chetverti.

14. Gabasov, R. & Kirillova, F.M. (1973) Osobyye optimal'nyye upravleniya [Singular optimal controls]. Moscow: Nauka.

15. Abdullayev, A.A. & Mansimov, K.B. (2013) Neobkhodimyye usloviya optimal'nosti vprotsessakh, opisyvayemykh sistemoy inte-gral'nykh uravneniy tipa Vol'terra [Necessary optimality conditions in processes described by a system of integral equations of Volterra type]. Baku: ELM.

16. Vasilyeva, A.B. & Tikhonov, A.N. (1989) Integral'nyye uravneniya [Integral equations]. Moscow: Moscow State University.

17. Tsalyuk, Z.B. (1977) Integral'nyye uravneniya Vol'terra [Integral equations of Volterra]. Itoginauki i tekhniki. Ser. matematicheskiy analiz. 15. pp. 131-138.

18. Gabasov, R. & Kirillova, F.M. (1973) Optimizatsiya lineynykh system [Optimization of linear systems]. Minsk: Belarusian State University.

19. Mansimov, K.B. (1994) Neobkhodimyye usloviya optimal'nosti osobykhprotsessov v zadache optimal'nogo upravleniya [Necessary conditions for the optimality of singular processes in the optimal control problem]. Abstract of Physics and Mathematics Dr. Diss. Baku.

20. Mansimov, K.B. (1999) Osobyye upravleniya v sistemakh s zapazdyvaniyem [Singular controls in systems with time delay]. Baku: ELM.