Научная статья на тему 'Альфа-множества в конечномерных евклидовых пространствах и их приложения в теории управления'

Альфа-множества в конечномерных евклидовых пространствах и их приложения в теории управления Текст научной статьи по специальности «Математика»

CC BY
162
15
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
α-МНОЖЕСТВО / ВЫПУКЛАЯ ОБОЛОЧКА / CONVEX HULL / ХАУСДОРФОВО РАССТОЯНИЕ / HAUSDORFF DISTANCE / УПРАВЛЕНИЕ / CONTROL / БЫСТРОДЕЙСТВИЕ / PERFORMANCE / УРАВНЕНИЕ ГАМИЛЬТОНА-ЯКОБИ / HAMILTON-JACOBI EQUATION / α-SET

Аннотация научной статьи по математике, автор научной работы — Ушаков Владимир Николаевич, Успенский Александр Александрович, Ершов Александр Анатольевич

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Ушаков Владимир Николаевич, Успенский Александр Александрович, Ершов Александр Анатольевич

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Alpha-sets in finite-dimensional Euclidean spaces and their applications in control theory

In this paper, a technique for investigating nonconvex sets that occur when describing the evolution of wave fronts, in the construction of generalized solutions of boundary value problems for equations of Hamilton-Jacobi type, in the formation of resolving structures in the problems of dynamic control is developed. An estimate is obtained for the Hausdorff distance between such sets and their convex hulls. The estimate is based on the concept of a measure of nonconvexity α. It is shown that for small α, nonconvex α-sets are close to convex. An example of a solution of the optimal control problem on the basis of α-sets is give.

Текст научной работы на тему «Альфа-множества в конечномерных евклидовых пространствах и их приложения в теории управления»

2018 ВЕСТНИК САНКТ-ПЕТЕРБУРГСКОГО УНИВЕРСИТЕТА Т. 14. Вып. 3

ПРИКЛАДНАЯ МАТЕМАТИКА. ИНФОРМАТИКА. ПРОЦЕССЫ УПРАВЛЕНИЯ

ПРОЦЕССЫ УПРАВЛЕНИЯ

УДК 514.74

МБС 52А27, 49Л5, 49J52

Альфа-множества в конечномерных евклидовых пространствах и их приложения в теории управления

В. Н. Ушаков, А. А. Успенский, А. А. Ершов

Институт математики и механики им. Н. Н. Красовского Уральского отделения РАН, Российская Федерация, 620990, Екатеринбург, ул. Софьи Ковалевской, 16

Для цитирования: Ушаков В. Н., Успенский А. А., Ершов А. А. Альфа-множества в конечномерных евклидовых пространствах и их приложения в теории управления // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления. 2018. Т. 14. Вып. 3. С. 261-272. https://doi.org/10.21638/11702/spbu10.2018.307

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

Ключевые слова: а-множество, выпуклая оболочка, хаусдорфово расстояние, управление, быстродействие, уравнение Гамильтона—Якоби.

Введение. Понятие а-множества было введено в начале 2000-х годов в работе [1], оно сформировалось при изучении множеств достижимости управляемых систем [2]. Множества достижимости, как правило, невыпуклы (см., например, [3]). Для одних систем они отличаются от выпуклых мало, для других — весьма ощутимо. В связи с этим возникла естественная потребность в наведении определенной классификации таких множеств по степени невыпуклости. Отметим, что есть другие классы обобщенно-выпуклых множеств, например а-паравыпуклые [4], функционально па-равыпуклые [5], а также семейство слабовыпуклых множеств ([6] и др.) Очевидно существование некоторой монотонной функции : [0, п] ^ [0,1] такой, что любое а-множество является ^(а)-паравыпуклым. Нахождение указанной функции — отдельная задача, для решения которой может быть полезна полученная здесь оценка.

© Санкт-Петербургский государственный университет, 2018

Данная работа посвящена изучению зависимости степени невыпуклости плоского a-множества M от его параметра а. Для частного случая, когда все «прогибы» границы представимы в виде графиков функций, показывается, что если а мало, то соответствующие a-множества близки к выпуклым множествам, а именно,

а

d(M, соМ) < А(М) • tg—, где d(■, ■) — хаусдорфово расстояние; X(M) = sup \\z — z*\\ — диаметр множества M.

z,zteM

Развиваемая техника исследования невыпуклых множеств, опирающаяся на понятие меры невыпуклости а, практически применяется при описании эволюции волновых фронтов [7], построении обобщенных решений краевых задач для уравнений Гамильтона—Якоби [8], конструировании функций оптимального результата в задачах управления [9, 10].

1. Оценка расстояния между a-множеством и его выпуклой оболочкой в частном случае.

1.1. Формулировка результата. Приведем определение a-множества из [11].

Пусть A — замкнутое множество в n-мерном евклидовом пространстве R" и z* €

R"\A. Под проекцией p(z*) точки z* на A понимаем ближайшую к z* точку в A.

Полагаем, что Qa(z*) = {p(z*)} — множество всех проекций p(z*) точки z* на A; coQa(z*) — выпуклая оболочка множества Qa(z*); con(co^A(z*) — z*) = {h = X(z — z*) : Л ^ 0,z € co^A(z*)} — конус в R", натянутый на множество coQa(z*) — z* = {z — z* : z € co^a(z*)}; (h*,h*) — скалярное произведение h* и h* из R",

1/2 (h h*\ ||/i*|| = (h*,h*) ; Z(K,h*) = arceos ' G [0,7г] — угол между векторами h*

\\h*\\\\h \\

и h*; aA(z*) = max Z(h*,h*) € [0,эт]; aA = sup aA(z*) € [0,эт].

Определение 1. Пусть a a = a. Тогда множество A назовем a-множеством.

Отметим, что 0-множество является обычным выпуклым множеством. К эт-мно-жествам относится кольцо или множество, состоящее из двух различных точек.

Определение 2. Пусть M — множество в R2. Назовем лакуной участок границы y С dM, гомеоморфный отрезку прямой и такой, что его крайние точки дy содержатся в д co M, а его внутренние точки Y\dY — в intco M.

Целью п. 1 является доказательство следующей теоремы.

Теорема. Пусть M — a-множество (a < -к) в R2 такое, что все лакуны границы dM представимы в виде графиков непрерывных функций, для которых ось абцисс параллельна отрезку, соединяющему крайние точки лакун.

Тогда

a

d(M,со М) < А(М) -tg-. (1)

1.2. Вспомогательные утверждения.

Лемма 1. (О пересечении альфа-графика с вертикальной полосой.) Пусть a ^ c < d < b, f € C[a,b], Г = {(x,y) : y = f (x),a < x < b}, y = {(x,y) : У = f (x), c ^ x ^ d}.

Тогда

ar = sup ar(z*) ^ aY = sup aY(z*). (2)

z*el2\r z*£R2\7

Доказательство. Прежде всего заметим, что так как Г — график однозначной функции, то аг (г*) всегда совпадает с величиной максимально возможного угла между двумя проекциями на Г, выпущенными из точки г * .То же верно и для 7.

Для доказательства неравенства (2) для каждой точки О Е М2\7 найдем такую точку О" Е М2\Г, что а1 (О) < аг (О").

Итак, возьмем произвольную точку О Е М2\7. Если а7(О) = 0, то в качестве О" подойдет любая точка из М2\Г.

Пусть оказалось, что а1 (О) > 0. Тогда изучим следующие случаи.

Первый случай. Проекции точки О на 7 также являются проекциями точки О на Г. В таком случае в качестве точки О" можно взять саму точку О.

Второй случай. Предположим, что это не так (рис. 1).

н-1-1-1—^

ас ъ

Рис. 1. Проекции точки О на 7 и точки О" на Г

Рассмотрим график 7 = {(х,у) : у = /(х),х Е [с, ¿]}. Пусть в некоторой точке О существует более одной проекции на 7. Обозначим величину максимального угла

ОА2 между проекциями через а. Это означает, что существует круг О с центром в точке О и радиусом \ОА1\ такой, что с1О не содержит никаких точек графика 7, кроме А1, Л2 и, может быть, некоторых точек на дуге А1Л2 (рис. 1).

Теперь покажем, что существует такая точка О" Е М2\Г, имеющая проекции О"В1 и О''В2 на множество Г (т. е. В1 ,В2 Е Г) такие, что АВ1О"В2 > А1ОЛ2.

Действительно, множество ГПс1 О не может быть расположено только на границе дО круга О, так как тогда мы находимся в рамках первого случая. Следовательно, множество Г П с1О размещено не только на дО, т. е. р(О, Г) < \ОА1 \.

Без ограничения общности будем считать, что абциссы точек А1 и А2 удовлетворяют соотношению хА1 < хА2. Обозначим через 71 = Г П с1О П{(х,у) : х ^ хА1}, 72 = Г П с1О П {(х,у) : х ^ хА2}. Заметим, что 71 ,у2 = 0, так как А1 Е 71, А2 Е 72.

Поскольку р(О, Г) < \ОА1 \, то хотя бы одно из расстояний р(О,/у1), р(О,^2) меньше радиуса круга О. Без ограничения общности будем считать, что расстояние р(О,^2) < \ОА1 \. Возьмем точку О' = О и будем сдвигать ее по направлению к точке А1. Заметим, что круг О' с центром в точке О' и радиусом \О'А1 \ всегда будет лежать внутри круга О. Кроме того, в некоторый момент р(О', у2) = \О'А1 \. При этом должна существовать хотя бы одна точка В2 Е у2 такая, что р(О',у2) = \О'В2\. Зафиксируем какую-нибудь одну такую точку В2. Докажем, что /А1 О'В2 ^ /А1ОА2 = а.

Во-первых, заметим, что точки Ах, А2 лежат либо не выше точки О, либо не ниже, иначе график перестает быть изображением однозначной функции. Без ограничения общности будем считать, что обе точки Ах, А2 размещены не выше точки О.

Во-вторых, точка В2 не может находиться внутри двух сегментов окружности О, заштрихованных на рис. 2, где точка А'2 является отражением точки А2 относительно прямой Ах О, поскольку точки А!2 и В2 всегда расположены по разные стороны от вертикальной прямой А2 А'2.

Возможны две ситуации относительно расположения точек Ах, А2 и вертикальной прямой V, проходящей через центр О круга О (рис. 2).

Рис. 2. Геометрическое место точек, в котором не может находиться точка В2 а — ситуация 1; б — ситуация 2.

Ситуация 1. Точки Ах и А2 расположены по разные стороны от вертикали V (т. е. абциссы точек Ах, О и А2 удовлетворяют соотношению хА1 < хо < хА2 ). В этой ситуации луч О'В2 пересекает отрезок ОА2, и очевидно, что /Ах О'В2 > /Ах ОА2.

Ситуация 2. Точки Ах и А2 размещены по одну сторону от V. Докажем от противного, что в данной ситуации также /А\О'В2 ^ /А\ОА2. Предположим, что в = /А\О'В2 < /А\ОА2 = а. Как ранее было отмечено, точка Вх не может находиться вне угла, образованного лучами Ах А2 и Ах А[ (заштрихованная область на рис. 2). Так как треугольник АА\ОА2 равнобедренный, то ¿.ОА\А2 = АА\А20 = Треуголь-

ник А Ах О' В2 также равнобедренный, поэтому АО'АхВ2 = ААхВ20' = Значит, АО'Ах В2 > /ОАхВ2. Но тогда точка В2 будет вне угла, образованного лучами Ах А2 и Ах А[. Следовательно, высказанное предположение неверно и /А\О'В2 ^ /А\ ОА2.

Рассмотрим теперь круг О' е центром в точке О' и радиусом \О'Ах \ = \О'В21. По построению О' С О, 72 и тШ' = 0. Если также 7х и тШ' = 0, то из неравенства /АхО'В2 > /АхОА2 следует, что а1 (О') > а^(О).

Если оказалось, что 7х и тШ' = 0, то тем же образом будем сдвигать точку О' по направлению к точке В2, пока не найдем новый центр круга О'', который касается множества 7х в некоторой точке Вх, с радиусом \О''Вх \ = \О''В2 \. При этом будет выполняться неравенство

ZBiO"B2 > ZAiO'B2 >

Таким образом, установили, что для каждой точки O G R2\y можно найти такую точку O" G М2\Г, что aY(O) ^ аг(O"), и тем самым доказали неравенство (2). □ Определим над- и подграфик функции f следующим образом:

Г+ = epi f = {(x,y) : y > f (x), a < x < b},

Г- = hypo f = {(x, y) : y < f (x), a < x < b}. Тогда будет верно следующее утверждение.

Следствие. Пусть f G C[a, b], a ^ c < d ^ b, Г = {(x, y) : y = f (x), a ^ x ^ b}, Y = {(x, y) : y = f (x), c < x < d}. Тогда ar- ^ aY- .

Доказательство. Возьмем произвольную точку O G R2\r-. Заметим, что если ar- (O) > 0, то Qr- (O) С Г. Тем самым попадаем в ситуацию леммы 1. □

Определение 3. Пусть f G C[a, b]. Сечением надграфика Г+ функции f, отвечающему значению y0 G R, назовем множество S(yo) = {(x, y) : y = yo}n Г+; сечением внутренности надграфика Г+ — множество Sint(y0) = {(x, y) : y = y0}nintr+, его диаметром —

HSint(yo)) = sup \x2 — xi\,

x1,X2:{(x,y):x1<x<X2,y=yo}CSint(y0)

т. е. длину наибольшего интервала в сечении внутренности надграфика. Заметим, что

HSint(y)) = max \x2 — xi\.

Xl,X2:f (xi) = f (X2) = y, X\<X<X2^f (X)<y

Кроме того, несложно показать, что если \(Sint(y)) = 0, то y ^ min f (x). Лемма 2. Пусть f G C[a,b], f (a) = f (b), Г = {(x,y) : y = f (x),a ^ x ^ b}, ar- = a < п.

a

Тогда f(a) — min f(x) ^ \b — a|tg—.

[a,b] 2

Доказательство. Вначале докажем оценку

a + £

/(a) - min f(x) < Ib- a|tg——, (3)

[a,b] 2

где число £ таково, что 0 < £ < п — a.

Построим монотонно убывающую числовую последовательность {yk}. Возьмем

yi = f (a). Если \(Sint(yi)) = 0, то min f (x) = f (a) и утверждение леммы выполняет-

[aM

ся. Если X(Sint(yi)) > 0, то тогда выберем точки Ai = (ai,yl), Bi = (bi,yl) из Г такие, что ai < bi и f (x) < yi для всех x G (ai, bi). (Тем самым отрезок [Ai, Bi] С cl Sint(yi).)

( ai + bi bi — ai (a £\\ Пусть точка Оi = —-—H---— • ctg + -JJ. Тогда угол ZAiOiBi =

a+£. Через B(Oi, OAi\) = {(x, y) : \ \(x,y) — Oi\ \ < OAi\} обозначим открытый круг радиуса \OiAi \. В силу следствия (для участка между точками Ai и Bi) и ar- < a+£ вытекает, что Г n B(Oi, \OiAi\) = 0.

Выберем точку C2 = (x2,y2) из множества Г П clB(Oi, \OiAi\) c минимальной ординатой y2 (рис. 3). Так как f (x) < yi для всех x G (ai, bi), то y2 < yi.

Рис. 3. Построение последовательности точек C2, C3, ...

Если Sint(y2) П {ai ^ y ^ bi} = 0, то выберем некоторые точки A2 = (a2,y2),

B2 = (b2,y2) из Г такие, что f (x) < y2 для всех x G (a2,b2). Возьмем точку 02 =

(a2 + b2 b2 - a2 (а £\\

2/2 H--^--ctg - + - и точку C3 = (ж3, уз) = argmm у.

2 \ 2 2// (х y) ¿=rncl в(Oo I

V 2

(х, y)ernc\ в (O2,\O2 A2\)

Продолжая действовать тем же образом, получим последовательность точек {Сп = (хп,уп)}. По построению ух > у2 > ...

Предположим без ограничения общности, что \АхС2\ < \C2В2\. Пусть АхК — касательная к окружности с центром в точке Ох и радиусом \ОхАх\. Через точку С2 проведем перпендикуляр НК2 к отрезку АхВх. Здесь через К обозначим точку пересечения этого перпендикуляра и касательной.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Из прямоугольного треугольника ААхНК получаем, что

\Ах Н \ а + е

=

\HK\

2

Но \AiBi\ - \A2B2\ > \AiH\, \HK\ > yi - y2. Таким образом,

а + £

Аналогично,

\AiBi \-\A2B2\ > (yi - y2) • ctg \A2 B2 \-\A3 B3 \ > (y2 - У3 ) • ctg

2

а + £ 2 '

(4)

(5)

Сложив неравенства (4), (5), ..., находим, что

\AiBi \ - \AnBn \ > (yi - yn) • ctg

а + £ 2

Таким образом, для любого n выполняется выражение

Уп > yi - \ b - a \tg

а + £ 2

(6)

Убывающая числовая последовательность {уп} ограничена снизу, следовательно, она сходится к некоторому значению у0, причем, в силу (6), выполняется оценка

yi - Уо < \ b - a \tg

а + £ 2

Кроме того, заметим, что монотонно возрастающая последовательность {an} ограничена сверху, а монотонно убывающая последовательность {bn} — снизу. То есть An — Ao = (ao, yo) и Bn — Во = (bo, yo) при n — ж. Докажем от противного, что Sint (yo) П [ao, bo] = 0.

Действительно, предположим, что Sint(yo) П [ao, bo] = 0. Но тогда можно выбрать

точки A* = (a*,yo), В* = (b*,yo) из Г такие, что b( — a* > 0, f (x) < yo для всех

л и П (a* + b* . b* — a* а + £\ x G (a*,o*). По точкам А * и построим точку О* = ^---,Уо~\-----ctg—-—J

(так, чтобы ZA*O*В* = а + £). В силу того, что участок Г между точками A* и В* является а-множеством, существует некоторая точка C* = (x* ,y*) из множества Г П В(O*, O* A* |) и такая, что y* < yo. Рассмотрим последовательности yn — yo, An — Ao, Вп — Во. Поскольку ao < a*, bo > b*, C* E int В (Oo, IOo Ao |), то существует

AT Л л w ^ (aN + bN , bN — aN

некоторое N такое, что С* E I П£>(Одг, |ОдгЛдг|), где Одг = ^---, Vn~\--^-х

ctg —Значит, числовая последовательность {уп} не может сходиться к у о > г/*.

Вследствие полученного противоречия Sint (yo) П [ao ,bo] = 0.

Теперь заметим, что на каждом шаге выбор точек Ak и Вк, принадлежащих границе Sint(yk), был достаточно произвольным. Если рассмотрим всевозможные такие последовательности, то придем к заключению, что все точки (xn, yn) из epiT удовлетворяют неравенству (6). Тем самым доказана оценка (3). Переходя в неравенстве (3) к пределу при £ — 0, получаем утверждение леммы 2. □

1.3. Доказательство теоремы. Рассмотрим множества М и co М. Заметим, что M — односвязное множество, так как ам < п. Следовательно, дМ состоит из совокупности участков, на которых она совпадает с дco М, и так называемых лакун. Именно на лакунах возможно отдаление множеств М и co М в хаусдорфовой метрике. Пусть участок границы между точками P\ и P2 — одна из таких лакун (рис. 4).

Рис. 4- Множество M и его выпуклая оболочка co M

Заметим, что множество М лежит по одну сторону от прямой PiP2, в противном случае точки Pi, P2 не были бы крайними точками лакуны.

Рассмотрим прямоугольную систему координат Oxy такую, чтобы ось Ox была параллельной Pi P2, а ось Oy была направлена от множества M. Пусть в данной системе координат точки Pi = (xi ,yi), P2 = (x2 ,yi). По условию теоремы участок границы дМ между точками PiP2 представим в виде графика Г = {(x,y) : y ^ f (x), xi < x < x2} некоторой непрерывной функции f (x).

Вследствие того, что все точки множества М лежат ниже прямой Pi P2, следует, что Г_ = hypo f является a-множеством, в противном случае множество М — не ^-множество.

При этом длина |Р1Р2| ^ Л(М) для любых Р1,Р2 € дМ, где Л(М) — диаметр множества М.

В силу леммы 2, выполняется следующая оценка для прогиба границы на участке между точками Р\ и Р2:

а

yi - min /(ж) < \х2 -xi|tg-,

[xi,Х2] 2

Здесь \x2 — x\I ^ Л(М). Заметим также, что хаусдорфово расстояние d(M, co М) не превосходит максимального прогиба границы дМ.

Таким образом, теорема доказана. □

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

x =Vl' (7)

У = V2,

v = (v\,v2), IIvII ^ 1, с целевым замкнутым множеством М С R2, имеющим нетривиальную геометрию.

Отметим, что u = u(x,y), (x,y) € R2\М, является минимаксным решением [12] краевой задачи Дирихле для уравнения Гамильтона—Якоби—Беллмана

min (у, Vu(x,y)) + 1 = 0, u|r = 0. (8)

v: | |v|

Здесь Уи(х/у) = ——, —— — градиент функции и(х,у), (х/у) G М2\М. Краевое \dx ду )

условие определено на границе Г = дМ замкнутого множества М С R2.

Структура функции оптимального результата для задачи управления c динамикой (7), как, впрочем, и структура минимаксного решения задачи (8), известна [8]:

u = p((x, у), М), (x,y) € Ж2\М,

где p((x,y), М) — евклидово расстояние от точки (x,y) до замкнутого множества М С R2. Функция u = р(^,у),М) является супердифференцируемой функцией на множестве Ж2\М (см. [13]).

Примем в качестве целевого множества в задаче управления двухкомпонентное множество

М = М1 U М2,

в котором первая компонента М1 — двухзвенная ломаная с узлами в точках Ai = (0,2), A2 = (0,0), A3 = (3,0), звеньями ломаной М1 являются отрезки [A1,A2], [A2, A3], вторая компонента М2 = {A4} — одноточечное множество, A4 = (1, 3).

Синтезируем функцию оптимального результата u = u(x,y), (x,y) € К2\М, осуществив надлежащим образом разбиение фазового пространства. Для этого рассмотрим задачу параметрического программирования

inf p(z, m) = p(z,A4), z = (x,y) € Ш2\М. (9)

Решением задачи (9) является кривая Г = grf (рис. 5), равноудаленная от множеств М1 и М2. Здесь

/ (х) =

' —х + 3,

3 - л/2х - 1,

х2 -2х + 10

б ' 2 1

х ^ 1,

1 < х < 4 - л/6, х > 3.

4

3.5 3

2.5 2

1.5 1

0.5 0

-0.5

X

-1

0

1

Рис. 5. а-гиперплоскость, сильно разделяющая компоненты Ы\ и М2 целевого множества

Рис. 6. Карта линий уровня функции оптимального результата

График Г = gr/ является а-прямой (а-гиперплоскостью), сильно отделяющей (см. [11]) множества М1 и М2. Мера невыпуклости подграфика выпуклой функции определяется линейными функциями, участвующими в склейке. Угловая величина а = агссо8(1/\/2б) ~ 1.3734. Кривая Г = gr/ гомеоморфна прямой и делит плоскость на две а-полуплоскости: Ф- = {(х,у) £ М2 : х £ М,у < /(х)} и Ф+ = {(х,у) £ М2 : х £ М,у > /(х)}. В рассматриваемой задаче функция и(х,у) = р((х,у),М), и поскольку а-прямая Г = gr/ равноудалена от компонент целевого множества, то

и(х у) = / р((х,у)М), и(х,у р((х,у ),М2 ),

(х,у) £ Ф1, (х,у) £ Ф2 и Г.

Основную проблему при конструировании решения составляет построение сингулярного множества — биссектрисы (см. [14]) целевого множества М. Для преодоления этой проблемы применим технику решения задач быстродействия, опирающуюся на свойства множеств симметрии [9, 15]. Невыпуклая компонента М2 целевого множества М имеет одну псевдовершину А2 = (0, 0) — точку негладкой склейки звеньев ломаной. Псевдовершина порождает ветвь сингулярного множества функции оптимального результата, совпадающую с открытым отрезком Ь\ биссектрисы первого координатного угла. Крайними точками Ь1 являются А2 = (0, 0) и Аъ = {4-л/6,4-л/6). Таким образом, в рассматриваемой задаче сингулярное множество Ь = Ь1 и Г. При этом точка А5 сочленения кривых Ь1 и Г — точка бифуркации. Это единственная точка сингулярного множества, имеющая не две, а три ближайшие точки на целевом множестве М.

На рис. 6 представлена геометрия сингулярного множества и сечений множества управляемости. Отметим, что с ростом расстояния до целевого множества мера невыпуклости а соответствующих линий уровня функции оптимального результата не увеличивается (см. [11, теорема 2]).

Заключение. Ценность полученной в п. 1 оценки (1) состоит в том, что при малых а можно заменять невыпуклые а-множества выпуклыми с точностью до небольшой погрешности. Аналогичным образом часто применяется известная теорема Шепли—Фолкмана [16].

Литература

1. Ушаков В. Н., Успенский А. А., Фомин А. Н. «-множества и их свойства. Екатеринбург: Ин-т математики и механики Урал. отд. РАН, 2004. 62 с.

2. Ли Э. Б., Маркус Л. Основы теории оптимального управления / пер. с англ. Л. Л. Леонтьевой; под ред. Я. Н. Ройтенберга. М.: Наука, 1972. 574 с. (Li E. B., Markus L. Foundation on optimal control theory.)

3. Пацко В. С., Пятко С. Г., Федотов А. А. Трехмерное множество достижимости нелинейной управляемой системы // Изв. РАН. Теория и системы управления. 2003. Т. 42, № 3. С. 8—16. DOI: 10.1134/S1064230706030075

4. Michael E. Paraconvex sets // Math. Scand. 1959. Vol. 7, N 2. P. 312-315. DOI: 10.7146/ math.scand.a-10583

5. Семенов П. В. Функционально паравыпуклые множества // Мат. заметки. 1993. Т. 54. Вып. 6. С. 74-81.

6. Иванов Г. Е. Слабо выпуклые множества и функции: теория и приложения. М.: Физматлит, 2006. 352 с.

7. Успенский А. А., Лебедев П. Д. Геометрия и асимптотика волновых фронтов // Изв. вузов. Математика. 2008. Т. 52, № 3. С. 27-37. DOI: 10.3103/S1066369X08030031

8. Лебедев П. Д., Успенский А. А., Ушаков В. Н. Построение минимаксного решения уравнения типа эйконала // Труды Ин-та математики и механики Урал. отд. РАН. 2008. Т. 14, № 2. С. 182-191. DOI: 10.1134/ S0081543808060175

9. Успенский А. А., Лебедев П. Д. Построение функции оптимального результата в задаче быстродействия на основе множества симметрии // Автоматика и телемеханика. 2009. Т. 79, № 7. С. 50-57. DOI: 10.1134/S0005117909070054

10. Ушаков В. Н., Успенский А. А., Лебедев П. Д. Геометрия сингулярных кривых для одного класса задач быстродействия // Вестн. С.-Петерб. ун-та. Сер. 10. Прикладная математика. Информатика. Процессы управления. 2013. Вып. 3. С. 157-167.

11. Ушаков В. Н., Успенский А. А. «-множества в конечномерных евклидовых пространствах и их свойства // Вестн. Удмуртск. ун-та. Математика. Механика. Компьютерные науки. 2016. Т. 26, № 1. С. 95-120. DOI: 10.20537/vm160109

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

12. Субботин А. И. Минимаксные неравенства и уравнения Гамильтона—Якоби. М.: Наука, 1991. 214 с.

13. Демьянов В. Ф., Васильев Л. В. Недифференцируемая оптимизация. М.: Наука, 1981. 384 с.

14. Успенский А. А., Лебедев П. Д. Аналитическое и численное конструирование функции оптимального результата для одного класса задач быстродействия // Прикладная математика и информатика: Труды факультета высшей математики и кибернетики Моск. гос. ун-та им. М. В. Ломоносова. 2007. № 27. С. 65-79.

15. Успенский А. А. Необходимые условия существования псевдовершин краевого множества в задаче Дирихле для уравнения эйконала // Труды Ин-та математики и механики Урал. отд. РАН. 2015. Т. 21, № 1. С. 250-263.

16. Starr R. M. Quasi-equilibria in markets with non-convex preferences // Econometrica. 1969. Vol. 37. Iss. 1. P. 25-38. DOI: 10.2307/1909201

Статья поступила в редакцию 21 января 2018 г.; принята к печати 14 июня 2018 г. Контактная информация:

Ушаков Владимир Николаевич —д-p физ.-мат. наук, проф., гл. науч. сотр.; [email protected]

Успенский Александр Александрович — д-p физ.-мат. наук, ведущий науч. сотр.; [email protected]

Ершов Александр Анатольевич — канд. физ.-мат. наук, науч. сотр.; [email protected]

Alpha-sets in finite-dimensional Euclidean spaces and their applications in control theory

V. N. Ushakov, A. A. Uspenskii, A. A. Ershov

Krasovskii Institute of Mathematics and Mechanics Ural Branch of Russian Academy Sciences, 16, S. Kovalevskaya ul., Yekaterinburg, 620990, Russian Federation

For citation: Ushakov V. N., Uspenskii A. A., Ershov A. A. Alpha-sets in finite-dimensional Euclidean spaces and their applications in control theory. Vestnik of Saint Petersburg University. Applied Mathematics. Computer Science. Control Processes, 2018, vol. 14, iss. 3, pp. 261-272. https://doi.org/10.21638/11702/spbu10.2018.307

In this paper, a technique for investigating nonconvex sets that occur when describing the evolution of wave fronts, in the construction of generalized solutions of boundary value problems for equations of Hamilton—Jacobi type, in the formation of resolving structures in the problems of dynamic control is developed. An estimate is obtained for the Hausdorff distance between such sets and their convex hulls. The estimate is based on the concept of a measure of nonconvexity a. It is shown that for small a, nonconvex a-sets are close to convex. An example of a solution of the optimal control problem on the basis of a-sets is give.

Keywords: a-set, convex hull, Hausdorff distance, control, performance, Hamilton—Jacobi equation.

References

1. Ushakov V. N., Uspenskii A. A., Fomin A. N. a-mnojestva i ih svoistva [a-sets and their properties]. Yekaterinburg, Institute of Mathematics and Mechanics, Ural Branch of the Russian Academy of Sciences Publ., 2004, 62 p. (In Russian).

2. Lee E. B., Markus L. Foundations of optimal control theory. New York etc., Wiley Press, 1967, 576 p. (Russ. end.: Lee E. B., Markus L. Osnovy teorii optimal'nogo upravlenija. Moscow, Nauka Publ., 1972, 574 p.)

3. Patsko V. S., Pyatko S. G., Fedotov A. A. Trehmernoe mnozhestvo dostizhimosti nelinejnoj upravljaemoj sistemy [Three-dimensional reachability set for a nonlinear control system]. Proceedings of the Russian Academy of Sciences. Theory and Systems Control, 2003, vol. 42, no. 3, pp. 8—16. (In Russian) DOI: 10.1134/S1064230706030075

4. Michael E. Paraconvex sets. Math. Scand., 1959, vol. 7, no. 2, pp. 312-315. DOI: 10.7146/ math.scand.a-10583

5. Semenov P. V. Funkcional'no paravypuklye mnozhestva [Functionally paraconvex sets]. Mathematical Notes, 1993, vol. 54, iss. 6, pp. 1236-1240. (In Russian)

6. Ivanov G. E. Slabo vypuklye mnozhestva i funkcii: teorija i prilozhenija [Weak convex sets and functions: theory and applications]. Moscow, Fizmatlit Publ., 2006, 352 p. (In Russian)

7. Uspenskii A. A., Lebedev P. D. Geometrija i asimptotika volnovyh frontov [Geometry and asymptotics of wavefronts]. Proceedings of Higher educational Institutions. Mathematics, 2008, vol. 52, no. 3, pp. 27-37. (In Russian) DOI: 10.3103/S1066369X08030031

8. Lebedev P. D., Uspenskii A. A., Ushakov V. N. Postroenie minimaksnogo reshenija uravnenija tipa jejkonala [Construction of a minimax solution for an eikonal-type equation]. Proceedings of the Steklov Institute of Mathematics and Mechanics Ural Branch of the Russian Academy of Sciences, 2008, vol. 14, no. 2, pp. 182-191. (In Russian) DOI: 10.1134/S0081543808060175

9. Uspenskii A. A., Lebedev P. D. Postroenie funkcii optimal'nogo rezul'tata v zadache bystrodejstvija na osnove mnozhestva simmetrii [Construction of the optimal outcome function for a time-optimal problem on the basis of a symmetry set]. Automation and Remote Control, 2009, vol. 70, iss. 7, pp. 50-57. (In Russian) DOI: 10.1134/S0005117909070054

10. Ushakov V. N., Uspenskii A. A., Lebedev P. D. Geometrija singuljarnyh krivyh dlja odnogo klassa zadach bystrodejstvija [Geometry of singular curves of a class of time-optimal problems]. Vestnik of Saint Petersburg University. Series 10. Applied Mathematics. Computer Sciences. Control Processes, 2013, iss. 3, pp. 157-167. (In Russian)

11. Ushakov V. N., Uspenskii A. A. a-mnozhestva v konechnomernyh evklidovyh prostranstvah i ih svojstva [a-sets in finite dimensional Euclidean spaces and their properties]. Vestnik of Udmurtsk

University. Mathematics. Mekhanics. Computer Sciences, 2016, vol. 26, iss. 1, pp. 95—120. (In Russian) DOI: 10.20537/vm160109

12. Subbotin A. I. Minimaksnyye neravenstva i uravneniya Gamil'tona—Yakobi [Minimax inequalities and the Hamilton—Jacobi equations]. Moscow, Nauka Publ., 1991, 214 p. (In Russian)

13. Dem'yanov V. F., Vasil'yev L. V. Nedifferentsiruyemaya optimizatsiya [Nondifferentiable optimization]. Moscow, Nauka Publ., 1981, 384 p. (In Russian)

14. Uspenskii A. A., Lebedev P. D. Analiticheskoe i chislennoe konstruirovanie funkcii optimal'nogo rezul'tata dlja odnogo klassa zadach bystrodejstvija [Analytical and numerical design of the optimal result function for one class of performance problems]. Applied Mathematics and Informatics. Proceedings of the Faculty of Computational Mathematics and Cybernetics M. V. Lomonosov Moscow State University, 2007, no. 27, pp. 65-79. (In Russian)

15. Uspenskii A. A. Neobhodimye uslovija sushhestvovanija psevdovershin kraevogo mnozhestva v zadache Dirihle dlja uravnenija jejkonala [Necessary conditions for the existence of pseudo-vertices of a boundary set in the Dirichlet problem for the eikonal equation]. Proceedings of the Institute of Mathematics and Mechanics Ural Branch of the Russian Academy of Sciences, 2015, vol. 21, no. 1, pp. 250-263. (In Russian)

16. Starr R. M. Quasi-equilibria in markets with non-convex preferences. Econometrica, 1969, vol. 37, iss. 1, pp. 25-38. DOI: 10.2307/1909201

Author's information:

Vladimir N. Ushakov — Dr. Sci. in Physics and Mathematics, Professor; [email protected]

Aleksander A. Uspenskii — Dr. Sci. in Physics and Mathematics; [email protected]

Aleksander A. Ershov — PhD Sci. in Physics and Mathematics; [email protected]

i Надоели баннеры? Вы всегда можете отключить рекламу.