CC BY
31
15. Ладыженская О.А., Солонников В.А., Уральцева Н.Н. Линейные и квазилинейные уравнения параболического типа. М.: Наука, 1967.
16. Casas E., Raymond J.-P., Zidani H. Pontryagin's Principle for Local Solutions of Control Problems with Mixed Control-State Constraints // SIAM J. Control Optim. 2000. V. 39. № 4. P. 1182-1203.
БЛАГОДАРНОСТИ: Работа поддержана грантом в рамках соглашения от 27 августа 2013 г. №02.В.49.21.0003 между Министерством образования и науки РФ и Нижегородским государственным университетом им. Н.И. Лобачевского.
Поступила в редакцию 9 июня 2015 г.
Gorshkov A.A. REGULARIZED PONTRYAGIN MAXIMUM PRINCIPLE IN OPTIMAL CONTROL FOR A PARABOLIC EQUATION WITH PHASE CONSTRAINTS IN LEBESGUE SPACES
The stable with respect to the errors in the initial data sequential Lagrange principle and Pontryagin maximum principle in a optimal control problem are considered. The target functional for this problem isstrictly uniformly convex, the control is distributed, the phase constraints are pointwised for a parabolic equation. The control is set from Lebesgue space of summable functions with p e (2, degree. The restriction operators images are put to the Lebesgue space of summable functions with s e (1,2) degree.
Key words: optimal control; parabolic equation; dual regularization; stability; point-wise phase constraint; Lebesgue space; Lagrange's principle; Pontryagin's maximum principle.
Горшков Андрей Александрович, Нижегородский государственный университет им. Н.И. Лобачевского, г. Нижний Новгород, Российская Федерация, аспирант, e-mail: tiger-nn@mail.ru
Gorshkov Andrey Aleksandrovich, Lobachevsky State University of Nizhni Novgorod, Nizhni Novgorod, the Russian Federation, Post-graduate Student, e-mail: tiger-nn@mail.ru
УДК 517.977
ПОСТРОЕНИЕ ОПТИМАЛЬНОГО УПРАВЛЕНИЯ СИНГУЛЯРНО ВОЗМУЩЕННОЙ СИСТЕМОЙ С ЗАПАЗДЫВАНИЕМ ПРИ ИНТЕГРАЛЬНЫХ ОГРАНИЧЕНИЯХ
© И.В. Гребенникова, А.Г. Кремлев
Ключевые слова: сингулярно возмущенная система с запаздыванием; оптимальное управление; фундаментальная матрица.
Рассматривается задача управления по минимаксному критерию для сингулярно возмущенной системы с запаздыванием при интегральных квадратичных ограничениях на ресурсы управления. Предлагается процедура построения управляющего воздействия, аппроксимирующего оптимальное решение с заданной степенью точности относительно малого положительного параметра.
Рассматривается управляемая сингулярно возмущенная система с запаздыванием Н > 0 (по состоянию):
(х(г)/(г = Ап(г)х(г) + (г)у(г) + С1(г)х(г — Н) + в^г, ^)и(г), ^)
^(у(г)/(г = А21 (г)х(г) + А^Шг) + ^(г)х(г - Н) + ад,^)и(г), ( )
где г € Т = [¿о,г1]; х € Еп, у € Ет; и € Ег — управление. Начальное состояние системы (1) х(г) = ф(г), г0 — Н ^ г <г0, х(г0) = х0, у(г0) = у0 точно неизвестно и заданы лишь
ограничения xo € Xo, yo € Yo , где Xo, Yo — выпуклые компакты в соответствующих пространствах, ф(t) € ^(t), t0 — h ^ t <t0, ^(t) — заданное многозначное отображение со значениями в виде выпуклых компактов, непрерывное по t в метрике Хаусдорфа. Управление u(t) , t € T — измеримые по Лебегу функции, удовлетворяющие условию u() € P, P — слабо компактное выпуклое множество в L2(T) . В данном случае
P = |u(-) | У u'(t)R(t)u(t)dt < Л21 , Л = const > 0, (2)
R(t) — симметричная, положительно определенная матрица с непрерывными элементами; штрих — знак транспонирования.
Рассматривается минимаксная задача управления [1]: среди управлений u() € P найти оптимальное u0 = u0(-), доставляющее
e°(ti,ß) = J(u0) = min J(u(-)), (3)
«(•)eP
J(u(-)) = max max <p(z(ti; u(^),zo,^(^)U, z0ez0 ^(-)еФ(-) V ' j
где <(■) заданная выпуклая функция (с конечными значениями); Z = (x',y'), z(t, u(), z0, ф( )) , t € T — решение системы (1), исходящее из Z0 = X0 х Y0 при некотором ф() € Ф(-) и фиксированном u() € P . Выполнено условие экспоненциальной устойчивости для подсистемы быстрых переменных. Запишем систему (1) в виде:
dz(t)/dt = A(t, ß)z(t) + G(t, ß)z(t — h)+ B(t, ß)u(t),
где матрицы A(t, ß), B(t, ß), G(t, ß) имеют следующий блочный вид:
Aii(t) Ai2(t) \ , _ ( Bi(t,ß) \ nu , _ ( Gi(t) 0
Л(ь,»] Ч А21 (ь)/» ЛиШ») , В(ь,»] ч , с(ь,ц) = V С2(1)/^ о
Пусть 2[Ь, т] есть фундаментальная матрица решений системы (1) (при и = О ), причем 2[т, т] = Е, 2[Ь, т] = О при т > Ь . Матрицу 2[Ь, т] представим в следующем блочном виде:
2[Ьт] = ( 2ц[1,т] 212[Ь,т] 2 [ ,т] V 221[Ь,т] 222% т]
здесь 211[Ь,т] , 212[Ь,т\ , 221[Ь,т\ , 222[Ь,т] — матрицы с размерами соответственно п х п, п х т, т х п, т х т.
Для реализации итерационной процедуры [2] построения оптимального решения задачи (3) важно правильно выбрать начальную асимптотику [3]. Аппроксимация оптимального решения при ограничениях (2) существенно зависит [2] от вида разложения матрицы В2(Ь,») по параметру » (О < » ^ .
Проведем исследование для случая В1(Ь,») = В1(Ь), В2(Ь,») = Л/»В2(Ь) .
В основе предлагаемого метода лежат идеи выделения асимптотики ансамбля траекторий сингулярно возмущенной системы [4] и представления фундаментальной матрицы решений, разбитой на блоки в соответствии с размерностями быстрых и медленных пре-менных, в виде равномерно сходящейся последовательности [3].
В [3] приведены рекуррентные формулы для вычисления блоков 2(1к\1,т] (г,] = 1, 2),к = = 0,1, 2,... , определяющих асимптотику матрицы 2[Ь,т] относительно малого параметра »> 0:
^11+1)[г,Т] = X[г,т] — / ((^М/^А-1^А21(5)^11)[5,Г] + ^(^Ц^ — Н, ТО(5;
т
4к+1)[г,т ] = г [г,т ] + / ^ма^г [5,т ](5;
т
[г, Т] = / ^11)[г,5]А12(5)Г[5,Т](5;
т
^ [г, Т ] = (1/^) / Г [г, 5] (А21(5)^11) [5, Т ] + С2(5)^11) [5 — Н, Т ])(в;
т
причем я^т] = х[г,Т],^2°2)[г,Т] = г[г,т],
где X[г, т] — фундаментальная матрица решений вырожденной системы (система (1) при ^ = 0 ), X[т, т] = Е , X[г, т] = 0 при Т > г,
Г [г, т ] — фундаментальная матрица решений системы = А22(г)у , Г [т, т ] = Е.
Вычисляя в соответствии с [2], при 0 < ^ ^ ^0 , № достаточно мало, имеем:
е0(г1 ) = ^(¿1) + о(/+1),
= тах{х(к)(р, 9) | р € Лп, 9 € Лт} = Х(к)(р(к), 9(к)), (4)
Х(к)(Р, 9) = —Н^к)(Р, 9) — (р, 9))1/2, к = 0,1, 2,...; Н(к) (р, 9) = (р, 9) — р [р,^(1) [г1, г0] + [г1, ^ I Xо] —
—р[р/^(2)[г1,г0]+ ^[¿ЪМ I Г>] — / р(г9к)(т,гьр,9)|Ф(т — н))(т,
г0
^ (р, 9) = / )(т,г1,р,9; ^)Я-1(т)г( ) (т,гьр,9; ^)(т+
¿0
0 2 2
(к)
где а0(^) > 0 : а0 ^ 0, ^ при ^ ^ +0; функции г (т,г1,р,9;г = 1, 2 ,
определяются следующим образом: при г0 ^ т ^ г ^ г1 — а0(^),
г1к)(т, ¿1,р, 9; = (р'^[¿1,т] + 9/^21)[г1,т])Б0(т,^) + (1/^А^)9/Г[¿ьт^т) — £ [р/^(2-1)[г1,т] + (1/^) /9'Г[г1,а]А21(а)^(2-1)[а,т](^А-21(т)в2(т),
т
при 0 ^ 5 < а0
г2к)(5, ¿1,р, 9; = — ^5,^1,р,9;
д,(т,^) = в1(т) — ^ДА12(Т )А-21 (Т )В2(Т),
г!к)(т,г,р,9) = (р'^(1)[г,т ] + 9/^21) [г, т ]) ад— 1Гр'^(2-1)[г,т] + (1/^) Т9'Г[М^да^^т](5)а2-21(т)с2(т),
йт \
с0(т ) = С1(Т) — А12(Т )а-21(т )С2 (т ); _
<^*(р, 9) — функция, сопряженная к <^(р, 9) ; Н**(р, 9) = (со Н)(р, 9) — замыкание выпуклой оболочки функции Н(р, 9) ; р(9^) — опорная функция множества X на элементе 9 . Рассмотрим управляющее воздействие и() (•) :
и(к) (Т) = / и(к) (Т), г0 < Т < г1 — а0 М,
I (г1 — т)/и), г1 — одМ <Т <гь
и(к)(-), -и(к)(-) определяются условиями:
при т € [Ьо,Ь1 - ао(»)},
и(к)(т) = -ХЕ-1(т)т1к) (т,Ь1,р(к),д(к)){ак (р(к) ,д(к)))-1/2; при § € [0, а0(»)/»],
у(к)(в) = -ХЯ-1(Ь1 - »в)г{^(8,Ь1,р(к),д(к)){ак(р(к) ,д(к)))-1/2.
Теорема. Пусть выполнены следующие условия:
1. Вырожденная система — система (1) при » = 0, относительно управляема [5] на Т;
2. Для любого Ь € Т гапк{В2(Ь1,»), Л22(Ь1)В2(Ь1,»),..., Лт,12-1(Ь1 )В2(Ь1,»)} = т;
3. Максимум в (4) достигается на векторе (№)' = (р(к ,д(к ) таком, что г1 (т, Ьър(к),д(к)) = 0, д(к) =0 .
Тогда задача (3) разрешима, причем при 0 < » ^ »0 , »0 достаточно мало, управляющее воздействие и^ (■) доставляет оценку
е0(11) = 3 (и0(■)) = 3 (и^ (■)) + 0(»к+1).
ЛИТЕРАТУРА
1. Куржанский А.Б. Управление и наблюдение в условиях неопределенности. М.: Наука, 1977.
2. Гребенникова И.В., Кремлев А.Г. Об итерационном методе построения оптимального управления сингулярно возмущенными системами с запаздыванием при квадратичных ограничениях // Известия Саратовского университета. Новая серия. Серия Математика. Механика. Информатика. Саратов, 2011. Т. 11. Вып. 3. Ч. 1. C. 8-15.
3. Гребенникова И.В. Аппроксимация решения в минимаксной задаче управления сингулярно возмущенной системой с запаздыванием // Известия вузов. Математика. Казань, 2011. № 10. С. 28-39.
4. Кремлев А.Г. Асимптотические свойства ансамбля траекторий сингулярно возмущенной системы в задаче оптимального управления // Автоматика и Телемеханика. Москва, 1993. № 9. С. 61-78.
5. Кириллова Ф.М., Чуракова С.В. Относительная управляемость линейных динамических систем с запаздыванием // Докл. АН СССР. Москва, 1967. Т. 174. № 6. С. 1260-1263.
Поступила в редакцию 10 июня 2015 г.
Grebennikova I.V., Kremlev A.G. CONSTRUCTION OF OPTIMAL CONTROL FOR SINGULARLY PERTURBED SYSTEM WITH DELAY WITH INTEGRAL CONSTRAINTS
The problem of control for the singularly perturbed delay system with integral quadratic constraints on the control resources with the minimax criterion is considered. Procedure of constructing control response that approximates the optimal solution with given accuracy with respect to a small positive parameter is proposed.
Key words: singularly perturbed system with delay; optimal control; fundamental matrix.
Гребенникова Ирина Владимировна, Уральский федеральный университет, г. Екатеринбург, Российская Федерация, старший преподаватель кафедры информационных систем технологий, e-mail: giv001@mail.ru
Grebennikova Irina Vladimirovna, Ural Federal University, Ekaterinburg, the Russian Federation, Senior Lecturer of the Information Systems Technology Department, e-mail: giv001@mail.ru
Кремлев Александр Гурьевич, Уральский федеральный университет, г. Екатеринбург, Российская Федерация, доктор физико-математических наук, профессор кафедры мультимедиа технологий, e-mail: kremlev001@mail.ru
Kremlev Aleksandr Gurevich, Ural Federal University, Ekaterinburg, the Russian Federation, Doctor of Physics and Mathematics, Professor of the Multimedia Technology Department, e-mail: kremlev001@mail.ru
