УДК 517.958
DOI: 10.14529/mmph240305
СИММЕТРИИНЫИ АНАЛИЗ УРАВНЕНИЯ ТИПА БЕЛЛМАНА
Д.И. Камалетдинова1, В.О. Лукащук2
1 ООО ИК «СИБИНТЕК», г. Уфа, Российская Федерация
2 Уфимский университет науки и технологий, г. Уфа, Российская Федерация E-mail: [email protected]
Аннотация. Исследуется уравнение с квадратичной нелинейностью, которое является частным случаем дифференциального уравнения типа Беллмана. В результате решения задачи групповой классификации установлено, что в случае произвольного вида двух функций, входящих в нелинейные члены, уравнение допускает четырехпараметрическую группу преобразований, которая расширяется до пятипараметрической и шестипара-метрической в случае параметрической и линейной зависимости функций. Построены некоторые инвариантные решения.
Ключевые слова: уравнение типа Беллмана; группа точечных преобразований; инвариантное решение.
Введение
Уравнение Беллмана [1] является одним из основных инструментов теории оптимального управления и динамического программирования. Оно широко используется в технических науках и экономике, является ключевым элементом многих алгоритмов машинного обучения и искусственного интеллекта.
Одной из актуальных задач теории управления является задача оптимальной коррекции траектории движения материальной точки при малых возмущениях. Согласно подходу академика Ф.Л. Черноусько [2, 3], исследование таких процессов может быть сведено к решению краевой задачи для нелинейного уравнения в частных производных второго порядка (уравнения типа Беллмана).
Характерной особенностью уравнения, существенно усложняющей его аналитическое исследование, является его многомерность: даже в простейшей постановке оно содержит три независимые переменные - временную и две пространственные. Одним из эффективных методов анализа нелинейных уравнений является групповой анализ [4]. Данный метод позволяет не только установить ряд качественных свойств уравнения, но и найти его частные решения уравнения. Его преимущество при решении задач оптимального управления лазером, например, подробно описаны в работе [5]. Групповой анализ применяется и в задачах об оптимальном управлении пограничным слоем [6].
В работе рассматривается уравнение с квадратичной нелинейностью, которое является частным случаем дифференциального уравнения Беллмана, вида [7, с. 259]
\2
^t^ x J^x^yy bVJy^y j
где f(t) и g(t) - произвольные функции, f(t) Ф 0, g(t) Ф 0. Решается задача групповой классификации относительно функций f(t) и g(t) по допускаемым этим уравнением однопараметриче-ским группам Ли точечных преобразований. Также строятся некоторые инвариантно-групповые решения уравнения (1).
Следуя работе [8], в дальнейшем используется обозначение d/dt = dt.
1. Групповая классификация
Будем искать для уравнения (1) допускаемую группу преобразований с оператором вида
X =t(t, x, y, u)dt +X(t, x, y, u)d x + £(t, x, y, u)d y +h(t, x, y, u)du. (2)
Действуя вторым продолжением оператора (2) на уравнение (1), получаем определяющее уравнение
\24
utux - f(t)uxuyy - g(t) (uy ) = 0 (1)
X (utux- f(t)uxuyy - g(t) (uy )2 )
= 0.
(1)
Камалетдинова Д.И., Симметрийный анализ уравнения типа Беллмана
Лукащук В. О.
Расщепляя его по независимым переменным их, иу, иху,..., приходим к довольно громоздкой определяющей системе, большая часть уравнений которой легко решается. В результате находим
т = т(0, Х = с\х + С2, С = сз У + С4, Т = С5и + С6, (3)
где С1 - произвольные постоянные, 1 = 1,...,6, а из нерешенных в определяющей системе останутся лишь следующие два уравнения:
¡ВТ'+ 8'т = (2Сз - С1)8,
Д т'+ Д 'т = 2С3Д.
(4)
Таким образом, задача групповой классификации сводится к решению классифицирующей системы на функции ДУ) и 8(У), построенной на основе системы (4)
К(0 + ) = ( 2 А - В) 8(0, Д) + Д '(У) w(t) = 2А), где А, В - произвольные постоянные, w(t) - произвольная функция. Поскольку ) Ф 0,8(У) Ф 0, разделим первое уравнение системы на 8(У), а второе - на ). Переходя к логарифмам функций
вд = 1п|ед|, 8(0=1л|8(0|,
имеем систему
М/1) + 8'(УМу) = 2А - В, ) + £'(У) w(t) = 2А. Вычитая первое уравнение из второго, приходим к соотношению
(Г- = в. (6)
Случай 1. Пусть £' Ф 8', тогда из (6)
В
w = --.
Г- 8'
Вычисляя производную w(t) и подставляя в систему (5), получаем уравнение
в г(г-8')-в (£''- 8'') = 2 а( £'-8')2, (7)
которое связывает функции ),8(У) и постоянные А, В.
Если ) и 8(У) - произвольные функции, то условие (7) выполняется при А = В = 0. Возвращаясь к системе уравнений (4), получаем, что С = 0, С3 = 0, т = 0. Следовательно, уравнение (1) допускает оператор (2) с координатами
т = 0, Х = С2, С = С4, Т = С5и + С6.
Пусть ДУ) и 8(У) удовлетворяют условию (7). Выясним их общий вид, когда А Ф 0, В Ф 0. Обозначим
2 А
1 =—, МО = ДО - 8(0,
В
причем Ъ'(У) Ф 0. Условие (7) перепишется в виде
£ Ъ'- Ъ'' = 1Ъ'2,
откуда
£ = 1Ъ + 1пЪ'+ 1п Ъ, 8 = (1-1) Ъ + 1пЬ'+ 1п Ъ,
где Ъ - произвольная постоянная, Ъ Ф 0. Получаем общий вид функций ДО и 8(0, которые удовлетворяют условию (7) в параметрическом виде:
Д = Ъ Ь'е1Ь, 8 = Ъ Ъ'е(1-1)Ь.
С 1
Подставляя функции ДО и 8(0 в систему (4), находим т = ——, С3 = —С1. Следовательно, в
Ъ'(У) 2
этом случае уравнение (1) допускает оператор (2) с координатами_
Вестник ЮУрГУ. Серия «Математика. Механика. Физика» 33
2024, том 16, № 3, С. 32-37
Математика
С 1
т =——, Х = СлХ + С2, С=—СлУ + С4, п = С5и + С6.
ь 1 ^ ь 2 5 6
Случай 2. Пусть Г' = тогда функции линейно зависимы: g(t) = к:(,), где к Ф 0 - произвольная постоянная. Из системы (4) получаем С1 = 0 и находим функцию т, а координаты допускаемого оператора (2) будут иметь вид 2С3 Г ^)сИ + С7
т=—:-, Х = С,, С = С3 у+С., Л = С ^и+Сб.
Таким образом, доказана справедливость следующего утверждения.
Утверждение. В зависимости от вида функций и g(t) уравнение типа Беллмана (1) допускает четырех-, пяти- или шестипараметрическую группу Ли точечных преобразований. А именно,
1) если ОД, g(t) - произвольного вида, ) Ф 0, g(t) Ф 0, то уравнение (1) допускает четырех-параметрическую группу Ли G4 с операторами
°4: Х1 = Э х, Х 2 = Э у , Х3 = Эи , Х 4 = иЭи ,
2) если : = Ь Ь'e1h,g = Ь Ь'е(1 1Ь, где Ь(,) - произвольная функция, ) Ф со^, Ь, 1 -произвольные постоянные, Ь Ф 0, 1ф 0, то группа Ли G4 расширяется до пятипараметрической оператором
1л л 1
-а. + хд х +—
Щ) ' х 2'
3) если g(t) = к:(.), где к Ф 0 - произвольная постоянная, то группа Ли G4 расширяется до шестипараметрической операторами
^ Э,, Хб = ^; т . 6 т
X5 = ^dt + xd X +1Г yd y,
X 5 = ^Zdt' Хб = «Л dt + УЭ У •
2. Инвариантные решения
Построим некоторые инвариантные решения уравнения (1). Общий инвариант для оператора Х2 есть функция I(,,х,и). Анзац на решение и = ((,, х) приводит к двум возможным инвариантным решениям уравнения (1) вида
и = ((,), и = х).
Для абелевой подалгебры с базисными операторами Х1, Х2 инвариантное решение тоже имеет вид произвольной функции от переменной х: и = ((х). Если взять за основу решения инвариант линейной комбинации операторов аХ2 — ЬХГ, где а,Ь - произвольные постоянные, то, предполагая, что решение имеет вид и = ((,, ах + Ьу), получим линейное редуцированное уравнение
а( = аЬ2^ )( + )(,
где г = ах + Ьу.
Рассмотрим двумерную подалгебру Ь2 ={Х1 + аХ3, Х2 + ЬХ3}, где а, Ь - произвольные постоянные. Из общего инварианта имеем анзац на решение и = ах + Ьу + ((,), подставляя который в уравнение (1), получим редуцированное уравнение
а( = Ь2g(t),
которое легко интегрируется. В результате находится инвариантное решение в виде
и = ах + Ьу + — Г g(t ^ + С,
п
где C - произвольная постоянная.
Камалетдинова Д.И., Симметрийный анализ уравнения типа Беллмана
Лукащук В. О.
Рассмотрим другую двумерную подалгебру {Х1 + аХ4,Х2 + ЪХ4}, из общего инварианта которой имеем анзац на решение вида и = еах+Ъур(У). В результате решения редуцированного уравнения получаем инвариантное решение
ах+Ъу+Ъ2 Г| )+^
и = Се ^ а ; ,
где С - произвольная постоянная.
Для оператора Х5, в случае параметрической зависимости функций Д и 8, анзац на решение
( ЩУ ) ^
а уравнение (1) редуцируется в уравнение с двумя переменны -
имеет вид и = р
хе~Н( >, уе 2
ми
Л
Р( + - + ЪРгРт + ЪР( = 0,
лиц)
где г = хе_И(У),V = уе 2 .
В случае линейной зависимости функций Д и 8 для оператора Х 5 получаем инвариантное решение и = р(х, у), где р удовлетворяет редуцированному уравнению
РууРх + (( = 0.
Для оператора Х6 инвариантное решение ищется в виде и = р( х,^т | ), а редуцирован-
ное уравнение меньшей размерности по пространству имеет вид
где г:
= ^ | ЩсИ.
4кг(2 + (6 г -1) р2 Рх + 4 г 2(хР22 = °
Заключение
В работе исследованы симметрийные свойства уравнения типа Беллмана. Установлено, что в случае произвольного вида функций от времени, стоящих при первых производных по пространству, допускаемая группа преобразований для исследуемого уравнения является четырехпара-метрической. Если функции связаны через параметрическое соотношение, то группа расширяется до пятипараметрической. В случае, когда функции пропорциональны, группа расширяется до шестипараметрической. Для некоторых допускаемых операторов и на их двумерных подалгебрах построены инвариантные решения, которые могут быть использованы в дальнейшем, например, при решении уравнения (1) с помощью численных методов для тестирования вычислительных программ. С помощью ряда операторов удалось понизить порядок уравнения (1) и выписать соответствующие редуцированные уравнения.
Литература
1. Ванько, В.И. Вариационное исчисление и оптимальное управление / В.И. Ванько, О.В. Ермошина, Г.Н. Кувыркин. - М.: МГТУ им. Н.Э. Баумана, 2006. - 488 с.
2. Черноусько, Ф.Л. Автомодельные решения уравнения Беллмана для задач оптимальной коррекции случайных возмущений / Ф.Л. Черноусько // ПММ. - 1971. - Т. 35, № 2. - С. 333-342.
3. Черноусько, Ф.Л Оптимальное управление при случайных возмущениях / Ф.Л. Черноусько, В.Б. Колмановский. - М.: Наука, 1978. - 351 с.
4. Ибрагимов, Н.Х. Практический курс дифференциальных уравнений и математического моделирования / Н.Х. Ибрагимов. - М.: Физматлит, 2012. - 332 с.
5. Полякевич, А.С. Преимущества группового анализа дифференциальных уравнений при решении задач оптимального управления лазерными системами / А.С. Полякевич, Б.Н. Пойзнер // Оптика атмосферы и океана. - 1995. - Т. 8, № 11. - С. 1697-1699.
Вестник ЮУрГУ. Серия «Математика. Механика. Физика» 2024, том 16, № 3, С. 32-37
Математика
6. Гараев, К.Г. Инвариантные краевые задачи оптимально управляемого пограничного слоя / К.Г. Гараев, В.А. Овчинников // Прикладная механика и техническая физика. - 2003. - Т. 44, № 1 (257). - С. 33-38.
7. Полянин, А.Д. Методы решения нелинейных уравнений математической физики и механики / А.Д. Полянин, В.Ф. Зайцев, А.И. Журов. - М.: Физматлит, 2009. - 256 с.
8. Овсянников, Л.В. Групповые свойства дифференциальных уравнений / Л.В. Овсянников -Новосибирск: СО АН СССР, 1962. - 238 с.
Поступила в редакцию 19 декабря 2023 г.
Сведения об авторах
Камалетдинова Дарья Ильшатовна - специалист управления информационных систем Lexema, ООО ИК «СИБИНТЕК», г. Уфа, Российская Федерация.
Лукащук Вероника Олеговна - кандидат физико-математических наук, доцент кафедры высокопроизводительных вычислительных технологий и систем, Уфимский университет науки и технологий, г. Уфа, Российская Федерация, e-mail: [email protected].
Bulletin of the South Ural State University Series "Mathematics. Mechanics. Physics" _2024, vol. 16, no. 3, pp. 32-37
DOI: 10.14529/mmph240305
SYMMETRY ANALYSIS OF THE BELLMAN EQUATION
D.I. Kamaletdinova1, V.O. Lukashchuk2
1 SIBINTEKIC LLC, Ufa, Russian Federation
2 Ufa University of Science and Technology, Ufa, Russian Federation e-mail: [email protected]
Abstract. The optimal correction of a material point's trajectory under small perturbations is an important problem in control theory. The study of such processes can be reduced to solving a boundary value problem for a special nonlinear second-order partial differential equation known as the Bellman equation. This paper deals with a partial differential equation with quadratic nonlinearity, which is a special case of the Bellman equation. The equation contains three independent variables—one temporal and two spatial, and two arbitrary time-dependent functions as multipliers. Multidimensionality is a distinctive feature of this equation, which significantly complicates its analytical study. Therefore, we use Lie group analysis which is an effective technique to analytically study nonlinear partial differential equations. It allows the investigation of not only the symmetry properties of equations, but also to find their particular solutions. The problem of group classification for this Bellman equation is solved with respect to two arbitrary time-dependent functions. It is established that in the case of the arbitrariness of these functions, the equation admits a four-parameter group of point transformations. This group expands to five-parameter and six-parameter groups for the parametric representation and linear dependency of the functions, respectively. Several invariant solutions are also constructed.
Keywords: Bellman type equation; group of point transformations; invariant solution.
References
1. Van'ko V.I., Ermoshina O.V., Kuvyrkin G.N. Variatsionnoe ischislenie i optimal'noe upravlenie (Calculus of Variations and Optimal Control). Moscow, MGTU im. N.E. Baumana Publ., 2006, 488 p. (in Russ.).
2. Chernousko F.L. Self-Similar Solutions of the Bellman Equation for Optimal Correction of Random Disturbances. Prikladnaya matematika i mekhanika, 1971, Vol. 35, no. 2, pp. 333-342. (in Russ.).
3. Chernous'ko F.L., Kolmanovskiy V.B. Optimal'noe upravlenie pri sluchaynykh vozmushcheniyakh (Optimal Control for Random Disturbances). Moscow, Nauka Publ., 1978, 351 p. (in Russ.).
Камалетдинова Д.И., Симметрийный анализ уравнения типа Беллмана
Лукащук В. О.
4. Ibragimov N.Kh. Prakticheskiy kurs differentsial'nykh uravneniy i matematicheskogo modeliro-vaniya (A Practical Course in Differential Equations and Mathematical Modeling). Moscow, Fizmatlit Publ., 2012, 332 c. (in Russ.). [Ibragimov N.H. A Practical Course in Differential Equations and Mathematical Modeling. Higher Education Press Limited Company, 2009, 364 р.]
5. Polyakevich A.S., Poizner B.N. Advantages of Group Analysis of Differential Equations in Solution of Some Problems of Laser System Optimum Control. Optika atmosfery i okeana, 1995, Vol. 8, no. 11, pp. 1697-1699. (in Russ.).
6. Garaev K.G., Ovchinnikov V.A. Invariantnye kraevye zadachi optimal'no upravlyaemogo pogranichnogo sloya (Invariant Boundary Value Problems of an Optimally Controlled Boundary Layer). Prikladnaya mekhanika i tekhnicheskayafizika, 2003, Vol. 44, no. 1 (257), pp. 33-38. (in Russ.).
7. Polyanin A.D., Zaytsev V.F., Zhurov A.I. Metody resheniya nelineynykh uravneniy matematicheskoy fiziki i mekhaniki (Methods for Solving Nonlinear Equations of Mathematical Physics and Mechanics). Moscow, Fizmatlit Publ., 2009, 256 p. (in Russ.).
8. Ovsyannikov L.V. Gruppovye svoystva differentsial'nykh uravneniy (Group Properties of Differential Equations). Novosibirsk, SO AN SSSR Publ., 1962, 238 p. (in Russ.).
Received December 19, 2023
Information about the authors
Kamaletdinova Daria Ilshatovna is Specialist, Department of Information Systems Management Lexema, Sibintek, Ufa, Russian Federation.
Lukashchuk Veronika Olegovna is Cand. Sc. (Physics and Mathematics), Associate Professor, Department of High Performance Computing Technologies and Systems, Ufa University of Science and Technology, Ufa, Russian Federation, e-mail: [email protected].
Вестник ЮУрГУ. Серия «Математика. Механика. Физика» 2024, том 16, № 3, С. 32-37