Доклады Российской академии наук. Математика, информатика, процессы управления, 2020, T. 495, № 1, стр. 100-106

ОЦЕНКА РОСТА СТЕПЕНИ НЕВЫПУКЛОСТИ МНОЖЕСТВ ДОСТИЖИМОСТИ УПРАВЛЯЕМЫХ СИСТЕМ В ТЕРМИНАХ α-МНОЖЕСТВ

Член-корреспондент РАН В. Н. Ушаков 1*, А. А. Ершов 1**

1 Институт математики и механики им. Н.Н. Красовского Уральского отделения Российской академии наук
Екатеринбург, Россия

* E-mail: ushak@imm.uran.ru
** E-mail: ale10919@yandex.ru

Поступила в редакцию 04.09.2020
После доработки 23.09.2020
Принята к публикации 24.09.2020

Полный текст (PDF)

Аннотация

Изучаются свойства α-множеств, являющихся одним из обобщений выпуклых множеств. Установлена взаимосвязь между α-множествами и слабо выпуклыми множествами по Виалю и по Ефимову–Стечкину. Получена оценка роста с течением времени меры невыпуклости α у множеств достижимости одного класса управляемых систем в двумерном фазовом пространстве.

Ключевые слова: обобщенно выпуклое множество, α-множество, слабо выпуклое множество, множество достижимости, управляемая система

α-множества введены в 2000-х гг. для классификации множеств достижимости управляемых систем по степени их невыпуклости [1]. Они являются одним из видов так называемых обобщенно выпуклых множеств, к которым относятся, например, сильно и слабо выпуклые множества по Виалю и по Ефимову–Стечкину [2], линейно выпуклые множества в пространстве над полем комплексных чисел [3], $\alpha $-паравыпуклые множества Э. Майкла [4, 5], созданные на их основе функционально паравыпуклые множества [6] и многие другие.

Прикладная значимость в теории управления и актуальность изучения α-множеств также основывается на двух предположениях: во-первых, для некоторых управляемых систем степень невыпуклости α у множеств достижимости растет непрерывно с течением времени (аналог теоремы 3.6.2 из [2, гл. 3, §3.6]), и, во-вторых, α-множества обладают свойствами, полезными для решения задач (в качестве аналога отметим работы [79]).

Главная цель настоящего сообщения заключается в доказательстве справедливости первого предположения для линейных по управлению управляемых систем в двумерном фазовом пространстве. Для достижения этой цели установим взаимосвязь между α-множествами и слабо выпуклыми множествами на плоскости, и затем воспользуемся результатом теоремы 3.6.2 из [2, гл. 3, §3.6].

1. СООТНОШЕНИЕ МЕЖДУ α-МНОЖЕСТВАМИ И СЛАБО ВЫПУКЛЫМИ МНОЖЕСТВАМИ ПО ЕФИМОВУ–СТЕЧКИНУ

Будем использовать следующие обозначения.

Через coM обозначим выпуклую оболочку множества M, clM – замыкание множества M, $\left\langle {{{x}_{ * }},x{\text{*}}} \right\rangle $ – скалярное произведение ${{x}_{ * }}$ и $x{\text{*}}$ из ${{\mathbb{R}}^{n}}$, $\left\| {{{x}_{ * }}} \right\| = {{\left\langle {{{x}_{ * }},{{x}_{ * }}} \right\rangle }^{{1/2}}}$ – евклидову норму в ${{\mathbb{R}}^{n}}$, diam(M) = = $\mathop {sup}\limits_{x,y \in M} \left\| {x - y} \right\|$ – диаметр множества M, B(a, r) = = $\{ x\, \in \,{{\mathbb{R}}^{n}}\,:\,\left\| {x\, - \,a} \right\|\, < \,r\} $ – открытый шар радиуса r и с центром в точке a, $\angle ({{x}_{ * }},x*)$ = ${\text{arccos}}\frac{{\langle {{x}_{ * }},x{\text{*}}\rangle }}{{\left\| {{{x}_{ * }}} \right\| \cdot \left\| {x{\text{*}}} \right\|}}\, \in \,[0,\pi ]$ – угол между векторами ${{x}_{ * }}$ и x*.

Под метрической проекцией $p{\text{*}}$ точки x* на множество M мы понимаем ближайшую к x* точку из M. Множество всех проекций точки x* на множество M обозначим через ${{\Omega }_{M}}(x*)$.

В частном случае $\alpha $-множества из двумерного евклидового пространства его можно определить следующим образом [10].

Определение 1. Пусть A – замкнутое множество в двумерном евклидовом пространстве ${{\mathbb{R}}^{2}}$ и $z* \in {{\mathbb{R}}^{2}}{\backslash }A$. Определим функцию

${{\alpha }_{A}}(z*) = \pi \quad {\text{в}}\;{\text{случае}},\quad {\text{если}}\quad z* \in {\text{co}}{{\Omega }_{A}}(z*),$
и

$\begin{gathered} {{\alpha }_{A}}(z*) = \mathop {max}\limits_{x,y \in {{\Omega }_{A}}(z*)} \angle (x - z*,y - z*) \\ {\text{в}}\;{\text{противном}}\;{\text{случае}}. \\ \end{gathered} $

Полагаем ${{\alpha }_{A}} = \mathop {sup}\limits_{z* \in {{\mathbb{R}}^{n}}\backslash A} {{\alpha }_{A}}(z*) \in [0,\pi ]$.

Множество A с мерой невыпуклости $\alpha = {{\alpha }_{A}}$ назовем α-множеством.

Отметим, что в евклидовых пространствах размерности 3 и выше α-множества определены иначе (см. [1, 10]), но в ${{\mathbb{R}}^{2}}$ все определения эквивалентны [10, теорема 1].

Определение 2 [11]. Множество A в линейном пространстве E называется  слабо выпуклым по Ефимову–Стечкину, с постоянной R > 0, если существует непустое множество ${{A}_{1}} \subset E$ такое, что

$A = \bigcap\limits_{a \in {{A}_{1}}} {(E{\backslash }B(a,R))} .$

Определение 3 [2, §1.1]. Пусть в линейном пространстве E с нормой ${{\left\| {\; \cdot \;} \right\|}_{E}}$ заданы две точки x, y и задано число $R \geqslant \frac{{{{{\left\| {x - y} \right\|}}_{E}}}}{2}$. Множество

${{D}_{R}}(x,y) = \bigcap\limits_{a \in E:{\text{\{ }}x,y{\text{\} }} \subset {\text{cl}}B(a,R)} {{\text{cl}}} B(a,R)$
называется сильно выпуклым отрезком.

Определение 4 [2, §1.1]. Множество A в линейном пространстве E с нормой ${{\left\| {\; \cdot \;} \right\|}_{E}}$ называется слабо выпуклым по Виалю с постоянной R > 0, если для любых двух точек $x,y \in A$ таких, что $0 < {{\left\| {x - y} \right\|}_{E}} < 2R$, существует точка $z \in {{D}_{R}}(x,y) \cap A$, не совпадающая с точками x и y.

Под расстоянием от точки x до множества M в линейном пространстве E с нормой ${{\left\| {\; \cdot \;} \right\|}_{E}}$ будем понимать величину $\rho (x,M) = \mathop {inf}\limits_{y \in M} {{\left\| {x - y} \right\|}_{E}}$.

Определение 5 [2, §1.7]. Будем говорить, что множество M в линейном пространстве E с нормой ${{\left\| {\; \cdot \;} \right\|}_{E}}$ обладает  чебышевским слоем толщины R > 0, если для любой точки $x \in E$ такой, что $\rho (x,M) < R$, существует и единственна метрическая проекция точки $x$ на M.

Лемма 1. Пусть заданы две точки p(1), p(2)на плоскости и число $R > \tfrac{{{\text{|}}{{p}^{{(1)}}}{{p}^{{(2)}}}{\text{|}}}}{2}$, где |p(1)p(2)| = ${\text{||}}{{p}^{{(2)}}} - {{p}^{{(1)}}}{\text{||}}$длина отрезка p(1)p(2). Обозначим через $s = \frac{1}{2}$(p(1) + p(2)) середину отрезка p(1)p(2), через $ps$ серединный перпендикуляр к отрезку p(1)p(2), через $B(q,R)$ такой открытый круг с центром в точке q и радиусом R, что ${{p}^{{(1)}}} \notin B(q,R)$, ${{p}^{{(2)}}} \notin B(q,R)$. Пусть некоторая точка $z \in ps \cap B(q,R)$ (рис. 1).

Рис. 1.

Схема расположения отрезка ${{p}^{{(1)}}}{{p}^{{(2)}}}$ и точки $z$.

Тогда, если $\left| {zs} \right| \geqslant R - \sqrt {{{R}^{2}} - \frac{1}{4}{\text{|}}{{p}^{{(1)}}}{{p}^{{(2)}}}{{{\text{|}}}^{2}}} $, то точки z и q лежат в одной полуплоскости относительно прямой $\Pi $, содержащей p(1)и p(2).

Доказательство. Введем декартову систему координат с центром в точке s, ось sy проведем через точку p(2). Без ограничения общности считаем, что точка z находится слева от точки s, как показано на рис. 1.

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

$\begin{gathered} {{p}^{{(1)}}} = (0, - {{y}_{2}}),\quad {{p}^{{(2)}}} = (0,{{y}_{2}}), \\ z = ( - x,0),\quad q = ({{x}_{0}},{{y}_{0}}), \\ \end{gathered} $
где

(1)
${{y}_{2}} = \frac{1}{2}{\text{|}}{{p}^{{(1)}}}{{p}^{{(2)}}}{\text{|}},$
(2)
$x \geqslant R - \sqrt {{{R}^{2}} - \frac{1}{4}{\text{|}}{{p}^{{(1)}}}{{p}^{{(2)}}}{{{\text{|}}}^{2}}} .$

Предположим противное, что точка q находится по другую сторону от прямой Π относительно точки z, т.е. ${{x}_{0}} > 0$. Покажем, что это невозможно.

Действительно, условия ${{p}^{{(1)}}} \notin B$(q, R) и ${{p}^{{(2)}}} \notin B$(q, R) равносильны системе неравенств

$x_{0}^{2} + {{({{y}_{2}} - {{y}_{0}})}^{2}} \geqslant {{R}^{2}},$
$x_{0}^{2} + {{({{y}_{2}} + {{y}_{0}})}^{2}} \geqslant {{R}^{2}},$
из которых (при ${{x}_{0}} > 0$) следует, что

(3)
${{x}_{0}} \geqslant \sqrt {{{R}^{2}} - y_{2}^{2}} .$

С учетом (1), (2) и (3) получаем, что

$\begin{gathered} {{\left| {zq} \right|}^{2}} = {{({{x}_{0}} + x)}^{2}} + y_{0}^{2} \geqslant \\ \geqslant \;\mathop {\left( {\sqrt {{{R}^{2}} - y_{2}^{2}} + R - \sqrt {{{R}^{2}} - \frac{1}{4}{{{\left| {{{p}_{1}}{{p}_{2}}} \right|}}^{2}}} } \right)}\nolimits^2 = {{R}^{2}}, \\ \end{gathered} $
т.е. $z$ не может принадлежать кругу $B(q,R)$.

Таким образом, пришли к противоречию.

Лемма 2. Пусть множество $A \subset {{\mathbb{R}}^{2}}$ слабо выпуклое по Ефимову–Стечкину с постоянной R, связно и его диаметр ${\text{diam}}(A) < 2R$.

Тогда A обладает чебышевским слоем толщины R.

Доказательство проведем от противного. По определению 5 для наличия чебышевского слоя толщины R у множества A все точки, удаленные от A на расстояние, меньшее R, должны иметь единственную метрическую проекцию на A. Предположим, что существует точка $z \in {{\mathbb{R}}^{2}}{\backslash }M$ такая, что $\rho (z,A) < R$, и имеющая по крайней мере две различные метрические проекции p(1) и p(2) на A. Обозначим через $r = {\text{||}}z - {{p}^{{(1)}}}{\text{||}}$. Заметим, что $r = \rho (z,A)$ < R и что в круге $B(z,r)$ не может быть ни одной точки из A.

Обозначим через $s = \tfrac{{{{p}^{{(1)}}} + {{p}^{{(2)}}}}}{2}$ и рассмотрим точку

$\begin{gathered} y = s + \frac{1}{2}\left( {r - \sqrt {{{r}^{2}} - \frac{1}{4}{\text{|}}{{p}^{{(1)}}}{{p}^{{(2)}}}{{{\text{|}}}^{2}}} } \right. + \\ + \;\left. {R - \sqrt {{{R}^{2}} - \frac{1}{4}{\text{|}}{{p}^{{(1)}}}{{p}^{{(2)}}}{{{\text{|}}}^{2}}} } \right)\frac{{{{p}^{{(1)}}} + {{p}^{{(2)}}} - 2z}}{{{\text{||}}{{p}^{{(1)}}} + {{p}^{{(2)}}} - 2z{\text{||}}}}. \\ \end{gathered} $

Заметим (см. рис. 2), что

$\left\| {y - z} \right\| = \left\| {y - s} \right\| + \left\| {s - z} \right\| = $
$\begin{gathered} = \;\frac{1}{2}\left( {r - \sqrt {{{r}^{2}} - \frac{1}{4}{\text{|}}{{p}^{{(1)}}}{{p}^{{(2)}}}{{{\text{|}}}^{2}}} } \right. + \\ + \;\left. {R - \sqrt {{{R}^{2}} - \frac{1}{4}{\text{|}}{{p}^{{(1)}}}{{p}^{{(2)}}}{{{\text{|}}}^{2}}} } \right) + \sqrt {{{r}^{2}} - \frac{1}{4}{\text{|}}{{p}^{{(1)}}}{{p}^{{(2)}}}{{{\text{|}}}^{2}}} < \\ \end{gathered} $
$ < \;r - \sqrt {{{r}^{2}} - \frac{1}{4}{\text{|}}{{p}^{{(1)}}}{{p}^{{(2)}}}{{{\text{|}}}^{2}}} + \sqrt {{{r}^{2}} - \frac{1}{4}{\text{|}}{{p}^{{(1)}}}{{p}^{{(2)}}}{{{\text{|}}}^{2}}} = r$
из-за монотонного убывания функции

f(r) = r ‒ $\sqrt {{{r}^{2}} - \frac{1}{4}{\text{|}}{{p}^{{(1)}}}{{p}^{{(2)}}}{{{\text{|}}}^{2}}} $ при $r > \frac{1}{2}{\text{|}}{{p}^{{(1)}}}{{p}^{{(2)}}}{\text{|}}$.

Следовательно, точка y находится в круге B(z, r), и, значит, $y \in {{\mathbb{R}}^{2}}{\backslash }A$. В силу слабой выпуклости по Ефимову–Стечкину множества A существует открытый круг радиуса R, который содержит точку y, но не пересекается с множеством A.

Рис. 2.

Точка $z$ и ее проекции ${{p}^{{(1)}}}$ и ${{p}^{{(2)}}}$ на множество A.

Однако

$\left\| {y - s} \right\| > R - \sqrt {{{R}^{2}} - \tfrac{1}{4}{{{\left| {{{p}^{{(1)}}}{{p}^{{(2)}}}} \right|}}^{2}}} ,$
и, следовательно, в силу леммы 1, центр любого такого круга, содержащего точку y, но не содержащего точки p(1) и p(2), будет находиться по одну сторону от прямой p(1)p(2) вместе с точкой y (т.е. находиться в одной полуплоскости, образуемой прямой p(1)p(2)). Аналогично, центр любого открытого круга радиуса $R$, не пересекающегося с A, содержащего точку z, но не содержащего точки p(1) и p(2), будет находиться по другую сторону от прямой p(1)p(2) относительно точки y (рис. 2).

Очевидно, что точки p(1) и p(2) невозможно соединить какой-либо кривой из A, не нарушив условия ${\text{diam}}(A) < R$. Следовательно, нарушается условие связности A и наше предположение о существовании нескольких проекций неверно.

Теорема 1. Пусть связное множество $A \subset {{\mathbb{R}}^{2}}$ слабо выпукло по Ефимову–Стечкину с постоянной R и ${\text{diam}}(A) \leqslant 2r$, где $r < R$. Тогда A является α-множеством с числом

$\alpha \leqslant arccos\left( {1 - 2{{{\left( {\frac{r}{R}} \right)}}^{2}}} \right).$

Доказательство. Прежде всего заметим, что, в силу [2, лемма 1.3.4], A является замкнутым, как и любое другое множество, слабо выпуклое по  Ефимову–Стечкину. Тем самым A является α-множеством с числом $\alpha \leqslant \pi $. Оценим степень его невыпуклости $\alpha $ более аккуратно.

Вначале докажем, что ни одна точка $z* \in {{\mathbb{R}}^{2}}{\backslash }A$ не может находиться в выпуклой оболочке своих проекций на множество A. Действительно, предположим, что существует такая точка $z* \in {\text{co}}{{\Omega }_{A}}(z*)$. Тогда множество проекций ${{\Omega }_{A}}(z*)$ должно содержать не менее двух точек. Однако, в силу леммы 2, множество A обладает чебышевским слоем толщины  R, или, иными словами, все точки из $R$-окрестности множества A обладают единственной метрической проекцией на A. Следовательно, расстояние $\rho (z*,{{\Omega }_{A}}(z*)) \geqslant R$ или, иначе говоря, точка z* является центром открытого круга $B(z*,\hat {R})$ с радиусом

(4)
$\hat {R} \geqslant R,$
который не содержит точек из A. В то же время на границе данного круга располагается множество ${{\Omega }_{A}}(z*)$. Если множество ${{\Omega }_{A}}(z*)$ содержит всего две проекции, то в силу соотношения $z* \in {\text{co}}{{\Omega }_{A}}(z*)$ эти проекции есть диаметрально противоположные точки круга $B(z*,\hat {R})$, расстояние между которыми $2\hat {R}$. Это противоречит неравенству ${\text{diam}}(A) \leqslant 2r$. Если множество ${{\Omega }_{A}}(z*)$ содержит большее количество точек, то можно выбрать из них такие точки x(1), x(2), x(3), что $z* \in {\text{co}}\{ {{x}^{{(1)}}},{{x}^{{(2)}}},{{x}^{{(3)}}}\} $. Однако несложными рассуждениями, учитывая линейную связность $A$, приходим к выводу, что данный случай также невозможен из-за неравенства ${\text{diam}}(A) \leqslant 2r$.

Итак, мы получили, что A является α-множеством с числом $\alpha < \pi $. Оценим теперь значение α, ссылаясь на определение 1. Пусть

$\alpha = \angle {{x}_{ * }}z{\text{*}}x* = \mathop {max}\limits_{x,y \in {{\Omega }_{A}}(z*)} \angle xz{\text{*}}y.$

Рассмотрим $\Delta {{x}_{ * }}x{\text{*}}z{\text{*}}$ на рис. 3.

Рис. 3.

Точка $z{\text{*}}$, имеющая более одной проекции на множество A.

По теореме косинусов

${{\left| {{{x}_{ * }}x{\text{*}}} \right|}^{2}} = {{\hat {R}}^{2}} + {{\hat {R}}^{2}} - 2\hat {R} \cdot \hat {R}cos\alpha .$

Отсюда,

(5)
$\alpha = arccos\left( {1 - \frac{{{{{\left| {{{x}_{ * }}x*} \right|}}^{2}}}}{{2{{{\hat {R}}}^{2}}}}} \right).$

Так как ${{\Omega }_{A}}(z*) \subset A$, то

(6)
$\left| {{{p}^{{(1)}}}{{p}^{{(2)}}}} \right| \leqslant {\text{diam}}({{\Omega }_{A}}(z*)) \leqslant {\text{diam}}(A) = 2r.$

Поскольку функция (5) достигает максимума при наибольшем ${\text{|}}{{p}^{{(1)}}}{{p}^{{(2)}}}{\text{|}}$ и наименьшем $\hat {R}$, то из неравенств (6) и (4) следует утверждение теоремы.

2. ОСНОВНОЙ РЕЗУЛЬТАТ

Введем в рассмотрение управляемую систему на промежутке $[{{t}_{0}},\vartheta ]$, описываемую векторным дифференциальным уравнением

(7)
$\dot {x}(t) = f(t,x(t)) + B(t)u(t),\quad x({{t}_{0}}) \in {{X}^{{(0)}}},$
где $t$ – время, $x(t) \in {{\mathbb{R}}^{2}}$ – фазовый вектор, ${{X}^{{(0)}}}$ – ограниченное множество в ${{\mathbb{R}}^{2}}$, $u(t)$ – вектор управления из компакта $P \subset {{\mathbb{R}}^{2}}$, $f:[{{t}_{0}},T] \times {{\mathbb{R}}^{2}} \mapsto {{\mathbb{R}}^{2}}$ – непрерывная функция, $B(t)$ – (2 × 2)-матрица с непрерывными коэффициентами. Под множеством допустимых управлений U системы (7) понимаем множество всех интегрируемых по Лебегу вектор-функций $u(t){\text{:}}\;[{{t}_{0}},T] \mapsto P$.

Предполагаем выполненными следующие условия на систему (7):

1) для любого $t \in [{{t}_{0}},T]$ функция $f(t,x)$ дифференцируема по x,

2) производная $f_{x}^{'}( \cdot ,x)$ ограничена константой M и удовлетворяет условию Липшица с постоянной L,

3) множество X(0) слабо выпукло по Ефимову–Стечкину с постоянной ${{R}_{0}} > \tfrac{1}{2}{\text{diam}}({{X}^{{(0)}}})$,

4) для любого $t \in [{{t}_{0}},T]$ множество B(t) · P = = ${\text{\{ }}B(t)u{\text{:}}\;u \in P{\text{\} }}$ сильно выпукло [2, определение 1.1.4] с постоянной ${{r}_{P}}$,

5) выполнено неравенство ${{r}_{P}} \leqslant M\hat {R}$, где

(8)
$\hat {R} = {{\left( {\frac{L}{{6M}} + \frac{1}{{{{R}_{0}}}}} \right)}^{{ - 1}}}exp( - 6M(T - {{t}_{0}})).$

При этих условиях существуют решения x(t) системы в классе абсолютно непрерывных функций на $[{{t}_{0}},T]$, которые мы называем движениями системы (7).

Обозначим через $X(T) = X(T;{{t}_{0}},{{X}^{{(0)}}})$ множество достижимости системы (7) в конечный момент времени T со стартовым множеством X(0). Относительно $X(T)$ сформулируем и докажем следующие утверждения.

Лемма 3. Справедлива оценка

${\text{diam}}(X(T)) \leqslant {{e}^{{M(T - {{t}_{0}})}}}{\text{diam}}({{X}^{{(0)}}}) + \hat {R}({{e}^{{M(T - {{t}_{0}})}}} - 1).$

Доказательство. Рассмотрим два произвольных движения $x(t)$ и $y(t)$ управляемой системы (7), порожденное допустимыми управлениям $u(t)$ и ${v}(t)$ соответственно, и оценим рассогласование $\left\| {y(T) - x(T)} \right\|$ этих движений в момент T, тем самым мы оценим ${\text{diam}}(X(T))$.

Из уравнения (7) следует, что для этих движений выполняется соотношение

$\dot {y}(t) - \dot {x}(t) = f(t,y(t)) - f(t,x(t)) + {v}(t) - u(t),$
которое после интегрирования по времени $t$ принимает вид

(9)
$\begin{gathered} y(t) - x(t) = y({{t}_{0}}) - x({{t}_{0}}) + \\ + \;\int\limits_{{{t}_{0}}}^t {(f(\tau ,y(\tau )) - f(\tau ,x(\tau )))} d\tau + \int\limits_{{{t}_{0}}}^t {({v}(\tau ) - u(\tau ))} d\tau . \\ \end{gathered} $

Учитывая сильную выпуклость множества $B(t) \cdot P$ и неравенство ${{r}_{P}} \leqslant M\hat {R}$ из условия 5), получаем, что

(10)
$\begin{gathered} \int\limits_{{{t}_{0}}}^t {\left\| {{v}(\tau ) - u(\tau )} \right\|} d\tau \leqslant \int\limits_{{{t}_{0}}}^t {{\text{diam}}} (B(t) \cdot P)d\tau \leqslant \\ \leqslant \;(t - {{t}_{0}}){{r}_{P}}(t - {{t}_{0}})M\hat {R}. \\ \end{gathered} $

Кроме того, применив теорему об оценке конечных приращений (см., например, [12, гл. 1, §3.3, теорема 3.3.2]) для функции $f(t, \cdot )$, получаем для любого момента времени $t \in [{{t}_{0}},T]$ неравенство

(11)
$\left\| {f(t,y(t)) - f(t,x(t))} \right\| \leqslant M\left\| {y(t) - x(t)} \right\|.$

Из (9), (10) и (11) следует оценка

(12)
$\begin{gathered} 0 \leqslant \left\| {y(t) - x(t)} \right\| \leqslant \left\| {y({{t}_{0}}) + x({{t}_{0}})} \right\| + \\ + \;M\int\limits_{{{t}_{0}}}^t {\left\| {y(\tau ) - x(\tau )} \right\|} d\tau + (t - {{t}_{0}})M\hat {R}. \\ \end{gathered} $

Применяя к (12) лемму Гронуолла [13, гл. 1, §2, с. 26], получаем неравенство

$\left\| {y(t) - x(t)} \right\| \leqslant \left\| {y({{t}_{0}}) - x({{t}_{0}})} \right\|{{e}^{{M(t - {{t}_{0}})}}} + \hat {R}({{e}^{{M(t - {{t}_{0}})}}} - 1) \leqslant $
$ \leqslant \;{\text{diam}}({{X}^{{(0)}}}){{e}^{{M(t - {{t}_{0}})}}} + \hat {R}({{e}^{{M(t - {{t}_{0}})}}} - 1),$
которое выполняется, в частности, при t = T. В силу произвольности выбора движений $x(t)$ и $y(t)$ из него следует утверждение леммы.

Замечание 1. Для выполнения утверждения леммы 3 по отношению к системе (7) условие 3) излишне, а условия 4) и 5) можно заменить простым неравенством ${\text{diam}}(B(t) \cdot P) \leqslant M\hat {R}$, где $\hat {R}$ – любая наперед заданная постоянная.

Теорема 2. Если

${{e}^{{M(T - {{t}_{0}})}}}{\text{diam}}({{X}^{{(0)}}}) + \hat {R}({{e}^{{M(T - {{t}_{0}})}}} - 1) < 2\hat {R},$
то

(13)
$\begin{gathered} {{\alpha }_{{X(T)}}} \leqslant \\ \leqslant \;{\text{arccos}}\left( {1\, - \,\frac{{{{{({{e}^{{M(T - {{t}_{0}})}}}{\text{diam}}({{X}^{{(0)}}})\, + \,\hat {R}({{e}^{{M(T - {{t}_{0}})}}}\, - \,1))}}^{2}}}}{{2{{{\hat {R}}}^{2}}}}} \right). \\ \end{gathered} $

Доказательство. Так как ${\text{diam}}({{X}^{{(0)}}}) < 2{{R}_{0}}$ по условию 3), то по лемме 2 множество A0 обладает чебышевским слоем толщины R0, что эквивалентно тому, что множество X(0) является слабо выпуклым по Виалю с постоянной R0 (см. [2, теорема 1.7.4]). Кроме того, в силу [2, лемма 1.3.4] множество X(0) является замкнутым. Но тогда, по [2, теорема 1.8.1], множество X(T) является замкнутым, слабо выпуклым по Виалю множеством с постоянной $\hat {R}$. По [2, теорема 1.3.4] множество $X(T)$ является слабо выпуклым по Ефимову–Стечкину с той же константой $\hat {R}$, а из [2, теорема 1.4.1] следует его связность.

В силу леммы 3 имеем оценку

${\text{diam}}(X(T)) \leqslant {{e}^{{M(T - {{t}_{0}})}}}{\text{diam}}({{X}^{{(0)}}}) + \hat {R}({{e}^{{M(T - {{t}_{0}})}}} - 1).$

В свою очередь, применяя к множеству $X(T)$ теорему 1, получаем оценку (13).

3. ПРИМЕРЫ

3.1. Пример 1

Рассмотрим математический маятник без силы трения, описываемый системой

(14)
$\begin{array}{*{20}{c}} {{{{\dot {x}}}_{1}} = {{x}_{2}},} \\ {{{{\dot {x}}}_{2}} = - sin{{x}_{1}} + {{u}_{2}},} \end{array}$
где ${{x}_{1}} = {{x}_{1}}(t)$ – угол отклонения маятника от вертикального положения, ${{x}_{2}} = {{\dot {x}}_{1}}(t)$ – скорость отклонения маятника от вертикального положения. Управление маятником на промежутке $[{{t}_{0}},T]$ осуществляется при помощи управляющего воздействия $u = u(t)$, стесненного ограничением u2 = = ${{u}_{2}}(t) \in [ - 1,\;1]$.

Определим множество начальных позиций ${{X}^{{(0)}}} = \{ ({{x}_{1}},{{x}_{2}}){\text{:}}\;x_{1}^{2} + x_{2}^{2} \leqslant 1\} $.

Теперь перейдем к рассмотрению системы (14) на предмет выполнения условий 1)–5). Прежде всего, запишем систему (14) в матричной форме

(15)
$\dot {x}(t) = f(x(t)) + u(t),$
где $x(t) = \left( {\begin{array}{*{20}{c}} {{{x}_{1}}(t)} \\ {{{x}_{2}}(t)} \end{array}} \right)$, $f(x)\, = \,\left( {\begin{array}{*{20}{c}} {{{x}_{2}}} \\ { - {\text{sin}}{{x}_{1}}} \end{array}} \right)$, $u(t) = \left( {\begin{array}{*{20}{c}} 0 \\ {{{u}_{2}}(t)} \end{array}} \right) \in P$,

$P = \left\{ {u = ({{u}_{1}},{{u}_{2}}):{{u}_{1}} = 0,{{u}_{2}} \in [ - 1,\;1]} \right\}$.

Поскольку P не является сильно выпуклым множеством, то теорема 2 не применима к системе (14).

3.2. Пример 2

Модифицируем систему (14) так, чтобы она полностью удовлетворяла условиям 1)–5). Для этого заменим ограничение P на управление u на множество $\hat {P}$ = $\{ u = ({{u}_{1}},{{u}_{2}}):u_{1}^{2} + u_{2}^{2} \leqslant 1\} $. Тогда ${{r}_{{\hat {P}}}}$ = 1.

Далее, находим

$f_{x}^{'} = \left( {\begin{array}{*{20}{c}} {\frac{{\partial {{f}_{1}}}}{{\partial {{x}_{1}}}}}&{\frac{{\partial {{f}_{2}}}}{{\partial {{x}_{1}}}}} \\ {\frac{{\partial {{f}_{1}}}}{{\partial {{x}_{2}}}}}&{\frac{{\partial {{f}_{2}}}}{{\partial {{x}_{2}}}}} \end{array}} \right) = \left( {\begin{array}{*{20}{c}} 0&{ - cos({{x}_{1}})} \\ 1&0 \end{array}} \right).$

Обозначим получившуюся матрицу $f_{x}^{'}$ через $\mathcal{A}$. Спектральная норма производной $f_{x}^{'}$, подчиненная евклидовой векторной норме, равна

${\text{||}}f_{x}^{'}{\text{||}} = \left\| \mathcal{A} \right\| = \mathop {sup}\limits_{\left\| x \right\| = 1} \sqrt {(\mathcal{A}x,\mathcal{A}x)} = \sqrt {{{\lambda }_{{max}}}(\mathcal{A}{\text{*}}\mathcal{A})} ,$
где $\mathcal{A}* = {{\bar {\mathcal{A}}}^{t}}$ – матрица, сопряженная к матрице $\mathcal{A}$, ${{\lambda }_{{max}}}(\mathcal{A}{\text{*}}\mathcal{A})$ – максимальное собственное значение произведения матриц $\mathcal{A}{\text{*}}\mathcal{A}$.

Последовательно находим

$\mathcal{A}* = \left( {\begin{array}{*{20}{c}} 0&1 \\ { - cos{{x}_{1}}}&0 \end{array}} \right),\quad \mathcal{A}{\text{*}}\mathcal{A} = \left( {\begin{array}{*{20}{c}} 1&0 \\ 0&{co{{s}^{2}}{{x}_{1}}} \end{array}} \right),$
$\begin{gathered} \left\| {\mathcal{A}{\text{*}}\mathcal{A} - \lambda E} \right\| = \left| {\begin{array}{*{20}{c}} {1 - \lambda }&0 \\ 0&{co{{s}^{2}}{{x}_{1}} - \lambda } \end{array}} \right| = \\ = \;(1 - \lambda )(co{{s}^{2}}{{x}_{1}} - 1),\quad {{\lambda }_{{max}}}(\mathcal{A}{\text{*}}\mathcal{A}) = 1. \\ \end{gathered} $

Отсюда, ${\text{||}}f_{x}^{'}{\text{||}} = 1$ в любой точке из ${{\mathbb{R}}^{2}}$. Следовательно, в качестве ограничивающей норму производной константы можно взять M = 1.

Константа Липшица L должна удовлетворять условию

$||f_{x}^{'}(t,x) - f_{y}^{'}(t,y)|| \leqslant L\left\| {x - y} \right\|,$
где $x = ({{x}_{1}},{{x}_{2}}) \in {{\mathbb{R}}^{2}}$, $y = ({{y}_{1}},{{y}_{2}}) \in {{\mathbb{R}}^{2}}$.

Полагаем

$\mathcal{B} = f_{x}^{'}(t,x) - f_{y}^{'}(t,y) = \left( {\begin{array}{*{20}{c}} 0&{cos({{y}_{1}}) - cos({{x}_{1}})} \\ 0&0 \end{array}} \right).$

Находим

$\begin{gathered} \mathcal{B}* = \left( {\begin{array}{*{20}{c}} 0&0 \\ {cos({{y}_{1}}) - cos({{x}_{1}})}&0 \end{array}} \right), \\ \mathcal{B}{\text{*}}\mathcal{B} = \left( {\begin{array}{*{20}{c}} 0&0 \\ {{{{(cos({{y}_{1}}) - cos({{x}_{1}}))}}^{2}}}&0 \end{array}} \right). \\ \end{gathered} $

Отсюда видно, что ${{\lambda }_{{max}}}(\mathcal{B}{\text{*}}\mathcal{B})$ = (cos(y1) – ‒ ${\text{cos}}({{x}_{1}}){{)}^{2}} \leqslant {{({{y}_{1}} - {{x}_{1}})}^{2}}$, т.е. $\left\| \mathcal{B} \right\| \leqslant \left\| {x - y} \right\|$. Отсюда следует, что $L = 1$.

Поскольку начальное множество X(0) выпукло, то ${{R}_{0}} = \infty $.

Далее находим по формуле (8) величину

$\hat {R} = {{\left( {\frac{L}{{6M}} + \frac{1}{{{{R}_{0}}}}} \right)}^{{ - 1}}}exp( - 6M(T - {{t}_{0}})) = 6\,{{e}^{{ - 6(T - {{t}_{0}})}}}.$

Неравенство ${{r}_{P}} \leqslant M\hat {R}$, необходимое для применения теоремы 2, будет выполняться при Tt0 ≤ ≤ $\tfrac{1}{6}ln6$ = 0.2986...

Наконец, применяя теорему 2, получаем оценку

$\alpha \leqslant arccos\left( {1 - \frac{{{{{({{e}^{{T - {{t}_{0}}}}} + \hat {R}({{e}^{{T - {{t}_{0}}}}} - 1))}}^{2}}}}{{2{{{\hat {R}}}^{2}}}}} \right) = $
$ = \;arccos\left( {1 - \frac{{{{{({{e}^{{7(T - {{t}_{0}})}}} + 6\,({{e}^{{T - {{t}_{0}}}}} - 1))}}^{2}}}}{{72}}} \right).$

Заметим, что величина

$1 - \frac{{{{{({{e}^{{7(T - {{t}_{0}})}}} + 6\,({{e}^{{T - {{t}_{0}}}}} - 1))}}^{2}}}}{{72}} < - 1$
только при $T - {{t}_{0}} > 0.{\text{3245}}\,...$, т.е. действие оценки прекращается из-за нарушения условия 5) раньше, чем аргумент арккосинуса выйдет из области определения арккосинуса.

ЗАКЛЮЧЕНИЕ

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

Список литературы

  1. Успенский А.А., Ушаков В.Н., Фомин А.Н. $\alpha $-множества и их свойства / ИММ УрО РАН. Екатеринбург, 2004. 62 с. Деп. в ВИНИТИ 02.04.2004, № 543-В2004.

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

  3. Зелинский Ю.Б. Выпуклость. Избранные главы. Киев: Институт математики НАН Украины, 2012. 280 с.

  4. Michael E. Paraconvex sets // Mathematica Scandinavica. 1959. V. 7. № 2. P. 312–315.

  5. Ngai H.V., Penot J.-P. Paraconvex functions and paraconvex sets // Studia Mathematica. 2008. V. 184. № 1. P. 1–29.

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

  7. Иванов Г.Е. Слабо выпуклые множества и их свойства // Матем. заметки. 2006. Т. 79. № 1. С. 60–86.

  8. Иванов Г.Е., Половинкин Е.С. Второй порядок сходимости алгоритма вычисления цены линейных дифференциальных игр // ДАН. 1995. Т. 340. № 2. С. 151–154.

  9. Ivanov G.E., Golubev M.O. Strong and weak convexity in nonlinear differential games // IFAC PapersOnline. 2018. V. 51. Iss. 32. P. 13–18.

  10. Ершов А.А., Кувшинов О.А. О свойствах пересечения $\alpha $-множеств // ИМИ УдГУ. 2020. Т. 55. С. 79–92.

  11. Стечкин С.Б., Ефимов Н.В. Опорные свойства множеств в банаховых пространствах и чебышевские множества // ДАН СССР. 1959. Т. 127. № 2. С. 254–257.

  12. Картан А. Дифференциальное исчисление. Дифференциальные формы. М.: Мир, 1971. 392 с.

  13. Лизоркин П.И. Курс дифференциальных и интегральных уравнений с дополнительными главами анализа. М.: Наука, 1981. 384 с.

Дополнительные материалы отсутствуют.

Инструменты

Доклады Российской академии наук. Математика, информатика, процессы управления