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

НЕОБХОДИМЫЕ И ДОСТАТОЧНЫЕ УСЛОВИЯ ВНУТРЕННЕЙ УСТОЙЧИВОСТИ ЛИНЕЙНЫХ ФОРМАЦИЙ

А. В. Лакеев *

Институт динамики систем и теории управления им. В.М. Матросова Сибирского отделения Российской академии наук
Иркутск, Россия

* E-mail: lakeyev@icc.ru

Поступила в редакцию 21.11.2019
После доработки 21.11.2019
Принята к публикации 22.11.2019

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

Аннотация

Получены необходимые и достаточные условия внутренней устойчивости формаций, динамика которых определяется линейными дифференциальными уравнениями. При этом в качестве классов допустимых управлений для лидеров выбраны программные управления, а для объектов, имеющих ведущих – аффинные обратные связи, зависящие от состояния самого объекта и состояний его ведущих. Полученные условия легко проверяемы и состоят из требований стабилизируемости пары матриц для уравнений ведомых объектов, гурвицевости и совпадения матриц для лидеров в случае многолидерности, разрешимости некоторых линейных уравнений и ограничений типа равенств на вектора, задающие требуемое взаимное расположение между ведомым и ведущим. Кроме того, описан весь класс управлений, обеспечивающих выполнение свойства линейной внутренней устойчивости. Опираясь на полученные условия, удалось показать, что внутренней устойчивостью могут обладать практически только однолидерные формации. В классе формаций с одним лидером выделен подкласс ( формации граф которых является входящим деревом), в котором не возникает ограничений типа равенств, являющихся основным препятствием для внутренней устойчивости многолидерных формаций.

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

1. Введение. Понятие внутренней устойчивости формаций нескольких управляемых движущихся объектов возникло при изучении задачи управления группой автономных роботов [1]. Под формацией понимается конечное множество взаимосвязанных движущихся объектов, связи между которыми (типа “ведущий–ведомый”) задаются некоторым ориентированным бесконтурным графом. Каждой паре связанных объектов сопоставляется некоторый вектор, определяющий их требуемое (в идеале) взаимное расположение, которое реально может быть выдержано только с некоторой точностью. Изучается задача существования управлений, обеспечивающих малые отклонения формации от желаемого идеального состояния при малом начальном отклонении и малых возмущающих внешних воздействиях. Это свойство, названное внутренней устойчивостью формаций, видимо, впервые появилось в работах Г.Дж. Таннера, Дж.Дж. Паппаса и В. Кумара (H.G. Tanner, G.J. Pappas, V. Kumar) (2002) [25] (см. также [6, 7]).

В этих работах за оценку малости отклонения была выбрана концепция “устойчивости от входа к состоянию” (input to state stability). Данное понятие впервые появилось в работах Э.Д. Cонтага (E.D. Sontag) (1989) [8, 9] и активно исследовалось в работах многих авторов, при этом в основном использовался метод функций Ляпунова. Некоторые итоги этих более чем 20-летних исследований содержатся в обзоре [10].

Отметим, что имеются и различные другие подходы к исследованию поведения формаций. Например, в работах [1113] также на основе метода функций Ляпунова изучается свойство типа диссипативности.

Однако практически во всех работах при исследовании внутренней устойчивости формаций вопрос существования управлений, обеспечивающих свойство внутренней устойчивости формации, в явном виде не ставился.

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

2. Основные понятия и определения. Перейдем к точным формулировкам. Следующие определения взяты из работ [25].

Определение 1. Будем говорить, что задана формация из l взаимосвязанных объектов, если: 

1) задан ориентированный бесконтурный граф $G = (V,E)$ без кратных ребер и петель с множеством вершин $V = {\text{\{ }}1,2,\; \ldots ,\;l{\text{\} }}$ и множеством ребер $E \subseteq V \times V$, называемый в дальнейшем графом формаций. При этом вершины графа отождествляются с движущимися объектами формации, а ребра указывают на наличие взаимосвязей между некоторыми из них. Если $(i,j) \in E$, то i-й объект называется ведомым, а j-й объект – ведущим (ребро $(i,j) \in E$ считаем ориентированным от i к j);

2) динамика движения i-го объекта задается дифференциальным уравнением с управлением

(1)
$\begin{gathered} {{{\dot {x}}}_{i}} = {{f}_{i}}(t,{{{\text{x}}}_{i}},{{{\text{u}}}_{i}}),\quad {{{\text{x}}}_{i}} \in {{R}^{n}}, \hfill \\ {{{\text{u}}}_{i}} \in {{R}^{m}},\quad t \in R_{0}^{1} = [0, + \infty ); \hfill \\ \end{gathered} $

3) каждому ребру $(i,j) \in E$ сопоставлен вектор ${{d}_{{ij}}} \in {{R}^{n}}$, задающий требуемое взаимное расположение (в идеале) i-го и j-го объектов относительно друг друга, при этом в идеале хотелось бы иметь соотношение ${{{\text{x}}}_{i}} + {{d}_{{ij}}} = {{{\text{x}}}_{j}}$;

4) если динамика движения всех объектов задается линейным дифференциальным уравнением с управлением, т.е. ${{\dot {x}}_{i}} = {{A}_{i}}{{{\text{x}}}_{i}} + {{B}_{i}}{{u}_{i}}$, где ${{A}_{i}} \in {{R}^{{n \times n}}}$ и ${{B}_{i}} \in {{R}^{{n \times m}}}$ – матрицы соответствующих размерностей, то формацию будем называть  линейной.

Будем обозначать Li = ${\text{\{ }}j \in V|(i,j) \in E{\text{\} }}$ = = ${\text{\{ }}{{j}_{{i1}}}, \ldots ,{{j}_{{i{{s}_{i}}}}}{\text{\} }}$ – множество ведущих для i-го объекта. Обозначим также ${{V}_{0}} = {\text{\{ }}i \in V|{{L}_{i}} = \phi {\text{\} }}$ – множе-ство  лидеров формации.

Если $(i,j) \in E$, то обозначим ${{{\text{z}}}_{{i,j}}}\, = \,{{{\text{x}}}_{i}} + {{d}_{{ij}}} - {{{\text{x}}}_{j}}$ – отклонение (ошибку) от идеального состояния между i-м и j-м объектами и ${\text{z}} = ({{{\text{z}}}_{{i,j}}})$ – ($l \times l$)-матрица ошибок,

${{{\text{z}}}_{{i,j}}} = \left\{ \begin{gathered} {{{\text{x}}}_{i}} + {{d}_{{ij}}} - {{{\text{x}}}_{j}},\quad {\text{если}}\quad (i,j) \in E, \hfill \\ 0,\quad {\text{если}}\quad (i,j) \notin E. \hfill \\ \end{gathered} \right.$

Для оценки величины отклонений используем стандартные классы функций [14]:

1) $\mathcal{K}$ – класс функций Хана: $\gamma \in \mathcal{K}$, если $\gamma :[0, + \infty ) \to [0, + \infty )$ непрерывная, строго возрастающая, $\gamma (0) = 0$;

2) функции класса $\mathcal{K}\mathcal{L}$: $\beta \in \mathcal{K}\mathcal{L}$, если $\beta :[0, + \infty ) \times [0, + \infty ) \to [0, + \infty )$, $\beta ( \cdot ,t) \in \mathcal{K}$ при фиксированном t, $\mathop {lim}\limits_{t \to \infty } \beta (a,t) = 0$, при фиксированном $a \in [0, + \infty )$, $\beta (a,t)$ непрерывна по (a, t).

Пусть для каждого $i \in V$ задан класс допустимых управлений ${{\mathcal{U}}_{i}}$ для i-го объекта.

Определение 2. Формация называется внутренне устойчивой, если существуют функции ${{\gamma }_{i}} \in \mathcal{K}$, $i \in {{V}_{0}}$ и функция $\beta \in \mathcal{K}\mathcal{L}$ такие, что для любых допустимых управлений ui${{\mathcal{U}}_{i}}$, $i \in {{V}_{0}}$, найдутся управления ${{u}_{j}} \in {{\mathcal{U}}_{j}}$, $j \in V{\backslash }{{V}_{0}}$ такие, что при любых начальных значениях ${{{\text{x}}}_{{0i}}}\, \in \,{{R}^{n}}$, $i = 1,2, \ldots ,l$ и любых $t \geqslant 0$ выполняется оценка

(2)
$\left\| {z(t)} \right\| \leqslant \beta (\left\| {z(0)} \right\|,t) + \sum\limits_{i \in {{V}_{0}}} {{{\gamma }_{i}}} (\mathop {sup}\limits_{\tau \in [0,t]} \left\| {{{u}_{i}}} \right\|),$
где z(t) – матрица ошибок, полученная на решениях уравнений (1) при заданных управлениях ${{u}_{i}} \in {{\mathcal{U}}_{i}}$ и заданных начальных значениях xi(0) = = ${{{\text{x}}}_{{0i}}} \in {{R}^{n}}$, $i = 1,2, \ldots ,l$.

Здесь и всюду в дальнейшем нормы векторов и матриц считаются евклидовыми.

Для линейных формаций будем рассматривать следующие классы допустимых управлений:

1) если $i \in {{V}_{0}}$, то ${{\mathcal{U}}_{i}} = \mathcal{U}_{i}^{{pr}} = {\text{\{ }}{{u}_{i}}|{{u}_{i}}{\text{:}}\,R_{0}^{1} \to {{R}^{m}}$, ui – непрерывно};

2) если $i \in V{\backslash }{{V}_{0}}$ и ${{L}_{i}} = {\text{\{ }}{{j}_{{i1}}},\; \ldots ,\;{{j}_{{i{{s}_{i}}}}}{\text{\} }}$, то ${{\mathcal{U}}_{i}} = \mathcal{U}_{i}^{{aff}}$ = = $\{ {{u}_{i}}|{{u}_{i}}{\text{:}}\,{{R}^{{n({{s}_{i}} + 1)}}} \to {{R}^{m}}{\text{\} }}$,

$\begin{gathered} {{u}_{i}}({{{\text{x}}}_{i}},{{{\text{x}}}_{{{{j}_{{i1}}}}}},\; \ldots ,\;{{{\text{x}}}_{{{{j}_{{i{{s}_{i}}}}}}}}) = \\ = \,\,{{S}_{i}}{{{\text{x}}}_{i}} + \;\sum\limits_{r = 1}^{{{s}_{i}}} {{{K}_{{i{{j}_{{ir}}}}}}} {{{\text{x}}}_{{{{j}_{{ir}}}}}} + {{k}_{i}} = {{S}_{i}}{{{\text{x}}}_{i}} + \sum\limits_{j \in {{L}_{i}}} {{{K}_{{ij}}}{{{\text{x}}}_{j}}} + {{k}_{i}}, \\ \end{gathered} $
где ${{S}_{i}},{{K}_{{ij}}} \in {{R}^{{m \times n}}}$, ${{k}_{i}} \in {{R}^{m}}$.

Определение 3. Линейную формацию будем называть линейно внутренне устойчивой (ЛВУ), если свойство из определения 2 выполняется при ${{\mathcal{U}}_{i}} = \mathcal{U}_{i}^{{pr}}$ для $i \in {{V}_{0}}$, и при ${{\mathcal{U}}_{i}} = \mathcal{U}_{i}^{{aff}}$ для $i \in V{\backslash }{{V}_{0}}$.

В кванторном виде это свойство более кратко можно записать следующим образом:

$\begin{array}{*{20}{c}} \begin{gathered} ({\text{ЛВУ}}) \equiv (\exists {{\gamma }_{i}} \in \mathcal{K}\,{\text{:}}\,\,\,i \in {{V}_{0}})(\exists \beta \in \mathcal{K}\mathcal{L}) \times \\ \times \;(\forall {{u}_{i}} \in \mathcal{U}_{i}^{{pr}}\,{\text{:}}\,\,\,i \in {{V}_{0}})(\exists {{u}_{j}} \in \mathcal{U}_{j}^{{aff}}{\text{:}}\,j \in V{\backslash }{{V}_{0}}) \\ \end{gathered} \\ \begin{gathered} (\forall {{{\text{x}}}_{{0i}}} \in {{R}^{n}}{\kern 1pt} {\text{:}}\;\,i \in V)(\forall t \geqslant 0)(\left\| {z(t)} \right\| \leqslant \beta (\left\| {z(0)} \right\|,t) + \\ + \;\sum\limits_{i \in {{L}_{0}}} {{{\gamma }_{i}}} (\mathop {sup}\limits_{\tau \in [0,t]} \left\| {{{u}_{i}}(\tau )} \right\|)). \\ \end{gathered} \end{array}$

Будем рассматривать еще два свойства линейных формаций:

${{({\text{ЛВУ}})}_{1}}$ – получается из $({\text{ЛВУ}})$ перестановкой кванторов $(\forall {{u}_{i}} \in \mathcal{U}_{i}^{{pr}}$: $i \in {{V}_{0}})$ и $(\exists {{u}_{j}} \in \mathcal{U}_{j}^{{aff}}$: $j \in V{\backslash }{{V}_{0}})$ или в кванторном виде:

$\begin{array}{*{20}{c}} \begin{gathered} {{({\text{ЛВУ}})}_{1}} \equiv (\exists {{\gamma }_{i}} \in \mathcal{K}{\kern 1pt} :\;i \in {{V}_{0}})(\exists \beta \in \mathcal{K}\mathcal{L}) \times \\ \times \;(\exists {{u}_{j}} \in \mathcal{U}_{j}^{{aff}}{\kern 1pt} :\;j \in V{\backslash }{{V}_{0}})(\forall {{u}_{i}} \in \mathcal{U}_{i}^{{pr}}{\kern 1pt} :\;i \in {{V}_{0}}) \\ \end{gathered} \\ \begin{gathered} (\forall {{{\text{x}}}_{{0i}}} \in {{R}^{n}}{\kern 1pt} :\;i \in V)(\forall t \geqslant 0)\left\| {z(t)} \right\| \leqslant \beta (\left\| {z(0)} \right\|,t) + \\ + \;\sum\limits_{i \in {{V}_{0}}} {{{\gamma }_{i}}} (\mathop {sup}\limits_{\tau \in [0,t]} \left\| {{{u}_{i}}(\tau )} \right\|); \\ \end{gathered} \end{array}$
${{({\text{ЛВУ}})}_{0}}$ – получается из $({\text{ЛВУ}})$ при нулевых управлениях лидеров, т.е. ${{u}_{i}}(t) \equiv 0$ для всех $i \in {{V}_{0}}$ (при этом в неравенстве (2) второе слагаемое тождественно равно нулю), или в кванторном виде:
$\begin{array}{*{20}{c}} \begin{gathered} {{({\text{ЛВУ}})}_{0}} \equiv (\exists \beta \in \mathcal{K}\mathcal{L})(\exists {{u}_{j}} \in \mathcal{U}_{j}^{{aff}}{\kern 1pt} :\;j \in V{\backslash }{{V}_{0}}) \times \\ \times \,(\forall {{u}_{i}} \equiv 0{\kern 1pt} :\;i \in {{V}_{0}}) \\ \end{gathered} \\ {(\forall {{{\text{x}}}_{{0i}}} \in {{R}^{n}}{\kern 1pt} :\;i \in V)(\forall t \geqslant 0){\text{||}}{{z}^{0}}(t){\text{||}} \leqslant \beta ({\text{||}}{{z}^{0}}(0){\text{||}},t),} \end{array}$
где z0(t) – матрица ошибок, полученная на решениях уравнений (1) при заданных управлениях ${{u}_{j}} \in \mathcal{U}_{j}^{{aff}}$ для $j \in V{\backslash }{{V}_{0}}$, ${{u}_{i}}(t) \equiv 0$ для $i \in {{V}_{0}}$ и заданных начальных значениях ${{x}_{i}}(0) = {{{\text{x}}}_{{0i}}} \in {{R}^{n}},$ i = 1, 2, ..., l.

Очевидно, что

${{({\text{ЛВУ}})}_{1}} \Rightarrow ({\text{ЛВУ}}) \Rightarrow {{({\text{ЛВУ}})}_{0}}.$

Рассмотрим более подробно граф формации, при этом будем придерживаться терминологии из [15, глава 16].

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

Во-вторых, хорошо известно, что ориентированный граф $G = (V,E)$ является бесконтурным тогда и только тогда, когда он допускает перенумерацию вершин с нижней строго треугольной матрицей смежности [15, с. 235, теорема 16.3]. Поэтому с точностью до перенумерации вершин можно считать, что ${{L}_{i}} \subseteq {\text{\{ }}1,2,\; \ldots ,\;i - 1{\text{\} }}$. В частности всегда ${{L}_{1}} = \phi $, т.е. $1 \in {{V}_{0}}$. Также будем считать, что ${{V}_{0}} = {\text{\{ }}1,2,\; \ldots ,\;{{l}_{0}}{\text{\} }}$.

В дальнейшем нам также понадобится функция $p:V \to V$, выбирающая для каждого $i > {{l}_{0}}$ одного из его ведущих. Для определенности будем считать, что это ведущий с наибольшим номером, а если $i \leqslant {{l}_{0}}$, то $p(i) = i$, т.е.

$p(i) = \left\{ \begin{gathered} i,\quad {\text{если}}\quad i \leqslant {{l}_{0}}, \hfill \\ max{\text{\{ }}j|j \in {{L}_{i}}{\text{\} }},\quad {\text{если}}\quad i > {{l}_{0}}. \hfill \\ \end{gathered} \right.$

Используя эту функцию, определим рекуррентным образом следующие вектора ${{D}_{i}} \in {{R}^{n}}$, которые также понадобятся нам в дальнейшем:

${{D}_{i}} = \left\{ \begin{gathered} 0,\quad {\text{если}}\quad i \leqslant {{l}_{0}}, \hfill \\ {{d}_{{ip(i)}}} + {{D}_{{p(i)}}},\quad {\text{если}}\quad i > {{l}_{0}}. \hfill \\ \end{gathered} \right.$

3. Необходимые условия для свойства ${{({\text{ЛВУ}})}_{0}}$.

Путь для формации выполняется свойство ${{({\text{ЛВУ}})}_{0}}$. Выберем и зафиксируем до конца данного раздела функцию $\beta \in \mathcal{K}\mathcal{L}$ и управления для объектов с номерами из $V{\backslash }{{V}_{0}}$, обеспечивающих это свойство, т.е. матрицы ${{S}_{i}} \in {{R}^{{m \times n}}}$, наборы матриц ${{{\text{\{ }}{{K}_{{ij}}}{\text{\} }}}_{{j \in {{L}_{i}}}}}$, ${{K}_{{ij}}} \in {{R}^{{m \times n}}}$, и векторы ${{k}_{i}} \in {{R}^{m}}$ такие, что если обозначить $x_{i}^{0}(t,{{{\text{x}}}_{{10}}},\; \ldots ,\;{{{\text{x}}}_{{i0}}})$ – решения задачи Коши для следующей системы дифференциальных уравнений:

$\begin{gathered} {{{\dot {x}}}_{i}} = {{A}_{i}}{{{\text{x}}}_{i}},\quad {\text{если}}\quad i = 1,2,\; \ldots ,\;{{l}_{0}}, \\ {{{\dot {x}}}_{i}} = ({{A}_{i}} + {{B}_{i}}{{S}_{i}}){{{\text{x}}}_{i}} + \\ + \;{{B}_{i}}\left( {\sum\limits_{j \in {{L}_{i}}} {{{K}_{{ij}}}{{{\text{x}}}_{j}}} + {{k}_{i}}} \right),\quad {\text{если}}\quad i > {{l}_{0}}, \\ \end{gathered} $

с начальными условиями $x_{i}^{0}(0,{{{\text{x}}}_{{10}}},\; \ldots ,\;{{{\text{x}}}_{{i0}}}) = {{{\text{x}}}_{{i0}}}$, $i = 1,\,\,2,\,\, \ldots ,\,\,l$, и, соответственно, ошибки $z_{{i,j}}^{0}(t)$ = = $z_{{i,j}}^{0}(t,{{{\text{x}}}_{{10}}}, \ldots ,{{{\text{x}}}_{{i0}}})$ = $x_{i}^{0}(t,{{{\text{x}}}_{{10}}}, \ldots ,{{{\text{x}}}_{{i0}}})$$x_{j}^{0}(t,{{{\text{x}}}_{{10}}}$, ... ..., xj0) + dij, при $(i,j) \in E$, то для всех $t \geqslant 0$ выполняется оценка

(4)
${\text{||}}{{z}^{0}}(t){\text{||}} \leqslant \beta ({\text{||}}{{z}^{0}}(0){\text{||}},t).$

Из неравенства (4) получаем следующее

Необходимое условие 1. $\mathop {lim}\limits_{t \to \infty } z_{{i,j}}^{0}(t) = 0$ при любых ${{x}_{i}}(0) = {{{\text{x}}}_{{i0}}} \in {{R}^{n}}$, $i = 1,2,\; \ldots ,\;l$.

Рассмотрим теперь поведение решений системы (3) при нулевых начальных условиях. Обозначим эти решения $\tilde {x}_{i}^{0}(t) = x_{i}^{0}(t,0,\; \ldots ,\;0)$. Так как очевидно, что $\tilde {x}_{i}^{0}(t) \equiv 0$ при $i = 1,2,\; \ldots ,\;{{l}_{0}}$, то индукцией по $i$ из необходимого условия 1 получаем следующее.

Предложение 1. Для всех $i = 1,2,\; \ldots ,\;l$ существует конечный предел $\mathop {\lim }\limits_{t \to \infty } \tilde {x}_{I}^{0}(t) = - {{D}_{i}}$, а если $j \in {{L}_{i}}$, то $\mathop {lim}\limits_{t \to \infty } \tilde {x}_{i}^{0}(t) = - ({{d}_{{ij}}} + {{D}_{j}})$ и, следовательно, ${{d}_{{ij}}} = {{D}_{i}} - {{D}_{j}}$.

Замечание 1. Отметим, что из этого предложения получаем жесткое ограничение на ${{d}_{{ij}}}$ в случае $\left| {{{L}_{i}}} \right| > 1$. Если же $\left| {{{L}_{i}}} \right| = 1$, то ${{L}_{i}} = {\text{\{ }}p(i){\text{\} }}$ и равенство ${{d}_{{ip(i)}}} = {{D}_{i}} - {{D}_{{p(i)}}}$ выполняется по определению Di.

Получим теперь еще одно необходимое условие для того, чтобы выполнялось свойство ${{({\text{ЛВУ}})}_{0}}$. Для этого при любом ${{{\text{x}}}_{0}} \in {{R}^{n}}$ обозначим ${\text{x}}_{{i0}}^{ * } = {{{\text{x}}}_{0}}$Di, $x_{i}^{ * }(t,{{{\text{x}}}_{0}}) = x_{i}^{0}(t,{\text{x}}_{{10}}^{ * },\; \ldots ,\;{\text{x}}_{{i0}}^{ * })$ и $z_{{i,j}}^{ * }(t,{{{\text{x}}}_{0}})$ = $x_{i}^{ * }(t,{{{\text{x}}}_{0}})$ – ‒ $x_{j}^{ * }(t,{{{\text{x}}}_{0}})$ + dij, $i = 1,2,\; \ldots ,\;l$, $(i,j) \in E$. Тогда из равенства dij = ${{D}_{i}} - {{D}_{j}}$ следует, что $z_{{i,j}}^{ * }(0,{{{\text{x}}}_{0}})$ = = ${{{\text{x}}}_{0}} - {{D}_{i}} - ({{{\text{x}}}_{0}} - {{D}_{j}})$ + dij = 0 и поэтому из неравенства (4) получаем следующее

Необходимое условие 2. $z_{{i,j}}^{ * }(t,{{{\text{x}}}_{0}}) \equiv 0$ при любых ${{{\text{x}}}_{0}} \in {{R}^{n}}$, $t \geqslant 0$, $i,j = 1,2,\; \ldots ,\;l$.

Дифференцируя это тождество по t и после этого подставляя t = 0, получаем следующее утверждение.

Предложение 2. Если $(i,j)$ребро графа $G$ ($(i,j) \in E$), то матрицы ${{S}_{i}},{{S}_{j}} \in {{R}^{{m \times n}}}$, векторы ${{k}_{i}},{{k}_{j}} \in {{R}^{m}}$ и наборы матриц ${{{\text{\{ }}{{K}_{{is}}}{\text{\} }}}_{{s \in {{L}_{i}}}}},$ ${{{\text{\{ }}{{K}_{{j\nu }}}{\text{\} }}}_{{\nu \in {{L}_{j}}}}}$ Kis, ${{K}_{{j\nu }}} \in {{R}^{{m \times n}}}$, удовлетворяют следующей системе линейных уравнений:

${{A}_{i}} + {{B}_{i}}\left( {{{S}_{i}} + \sum\limits_{s \in {{L}_{i}}} {{{K}_{{is}}}} } \right) = {{A}_{j}} + {{B}_{j}}\left( {{{S}_{j}} + \sum\limits_{\nu \in {{L}_{j}}} {{{K}_{{j\nu }}}} } \right),$
$\begin{gathered} {{B}_{i}}\left( {{{k}_{i}} - {{S}_{i}}{{D}_{i}} - \sum\limits_{s \in {{L}_{i}}} {{{K}_{{is}}}} {{D}_{s}}} \right) - {{A}_{i}}{{D}_{i}} = \\ = \;{{B}_{j}}\left( {{{k}_{j}} - {{S}_{j}}{{D}_{j}} - \sum\limits_{\nu \in {{L}_{j}}} {{{K}_{{j\nu }}}} {{D}_{\nu }}} \right) - {{A}_{j}}{{D}_{j}}, \\ \end{gathered} $
где для $j \in {{V}_{0}}$ считаем, что ${{S}_{j}} = 0$, $\sum\limits_{\nu \in {{L}_{j}}} {{{K}_{{j\nu }}}} = 0$, ${{k}_{j}}$ = 0.

Из слабой связности графа $G$ и предложения 2 получаем следующие

Предложение 3. Для всех $i = 1,2,\; \ldots ,\;{{l}_{0}}$ все матрицы Ai одинаковы, т.е. ${{A}_{i}} = {{A}_{1}}$, а для $i > {{l}_{0}}$ матрицы ${{S}_{i}} \in {{R}^{{m \times n}}}$, векторы ${{k}_{i}} \in {{R}^{m}}$ и наборы матриц ${{{\text{\{ }}{{K}_{{is}}}{\text{\} }}}_{{s \in {{L}_{i}}}}}{{K}_{{is}}} \in {{R}^{{m \times n}}}$ удовлетворяют следующей системе линейных уравнений:

(5)
${{B}_{i}}\left( {{{S}_{i}} + \sum\limits_{s \in {{L}_{i}}} {{{K}_{{is}}}} } \right) = {{A}_{1}} - {{A}_{i}},$
(6)
${{B}_{i}}\left( {{{k}_{i}} - {{S}_{i}}{{D}_{i}} - \sum\limits_{s \in {{L}_{i}}} {{{K}_{{is}}}} {{D}_{s}}} \right) = {{A}_{i}}{{D}_{i}}.$

Из необходимых условий 1 и 2 получаем также следующее утверждение о свойствах матриц ${{\tilde {A}}_{i}} = {{A}_{i}}$ + BiSi и ${{A}_{1}}$.

Предложение 4. Для всех $i > \;{{l}_{0}},\,\,...,\,\,l$ матрицы ${{\tilde {A}}_{i}} = {{A}_{i}} + {{B}_{i}}{{S}_{i}}$гурвицевы и, следовательно, пары матриц $({{A}_{i}},{{B}_{i}})$стабилизируемы, а если ${{l}_{0}} > 1$, то и матрица A1– гурвицева.

Замечание 2. Сделаем несколько замечаний о разрешимости уравнений (5), (6). Обозначая ${{N}_{i}} = {{S}_{i}} + \sum\limits_{s \in {{L}_{i}}} {{{K}_{{is}}}} $ и ${{\tilde {k}}_{i}} = {{k}_{i}} - {{S}_{i}}{{D}_{i}} - \sum\limits_{s \in {{L}_{i}}} {{{K}_{{is}}}} {{D}_{s}}$, получаем, что ${{N}_{i}} \in {{R}^{{m \times n}}}$ и ${{\tilde {k}}_{i}} \in {{R}^{m}}$ являются решениями следующих уравнений:

${{B}_{i}}{{N}_{i}} = {{A}_{1}} - {{A}_{i}},$
${{B}_{i}}{{\tilde {k}}_{i}} = {{A}_{i}}{{D}_{i}}$
для всех $i = 1,2,\; \ldots ,\;l$.

Обратно, если ${{N}_{i}} \in {{R}^{{m \times n}}}$ и ${{\tilde {k}}_{i}} \in {{R}^{m}}$ являются решениями предыдущих уравнений, ${{S}_{i}} \in {{R}^{{m \times n}}}$ – любая матрица, а набор матриц ${{{\text{\{ }}{{K}_{{is}}}{\text{\} }}}_{{s \in {{L}_{i}}}}}$, ${{K}_{{is}}} \in {{R}^{{m \times n}}}$, удовлетворяет соотношению $\sum\limits_{s \in {{L}_{i}}} {{{K}_{{is}}}} $ = ${{N}_{i}} - {{S}_{i}}$, то, полагая ${{k}_{i}} = {{\tilde {k}}_{i}} + {{S}_{i}}{{D}_{i}} + \sum\limits_{s \in {{L}_{i}}} {{{K}_{{is}}}} {{D}_{s}}$, получаем решение уравнений (5), (6).

4. Необходимые и достаточные условия линейной внутренней устойчивости.

Набор условий для наличия свойства ${{({\text{ЛВУ}})}_{0}}$ из предложений 1–4 оказался не только необходимым, но и достаточным. Более точно, верно следующее утверждение.

Теорема. Для линейной формации выполняется свойство линейной внутренней устойчивости тогда и только тогда, когда выполнены следующие условия:

1) для всех $i > {{l}_{0}}$ пары матриц $({{A}_{i}},{{B}_{i}})$стабилизируемы, т.е. существуют $(m \times n)$-матрицы Si такие, что матрицы ${{A}_{i}} + {{B}_{i}}{{S}_{i}}$гурвицевы;

2) для всех $i > {{l}_{0}}$ существует $(m \times n)$-матрица Ni и векторы ${{\tilde {k}}_{i}} \in {{R}^{m}}$, удовлетворяющие следующим линейным уравнениям:

(7)
${{B}_{i}}{{N}_{i}} = {{A}_{1}} - {{A}_{i}},$
(8)
${{B}_{i}}{{\tilde {k}}_{i}} = {{A}_{i}}{{D}_{i}},$

3) если $(i,j) \in E$, то ${{d}_{{ij}}} = {{D}_{i}} - {{D}_{j}}$;

4) если ${{l}_{0}} > 1$, то ${{A}_{i}} = {{A}_{1}}$ для всех $i = 1,2,\; \ldots ,\;{{l}_{0}}$ и матрица A1– гурвицева.

При этом если Si, Ni и ${{\tilde {k}}_{i}}$ удовлетворяют условиям 1), 2) теоремы, то управления ведомых, обеспечивающие внутреннюю устойчивость, можно выбрать в виде

${{u}_{i}}({{{\text{x}}}_{1}},\; \ldots ,\;{{{\text{x}}}_{i}}) = {{S}_{i}}{{{\text{x}}}_{i}} + \sum\limits_{s \in {{L}_{i}}} {{{K}_{{is}}}} {{{\text{x}}}_{s}} + {{k}_{i}},\quad i > {{l}_{0}},$
где набор матриц ${{{\text{\{ }}{{K}_{{is}}}{\text{\} }}}_{{s \in {{L}_{i}}}}}$, ${{K}_{{is}}} \in {{R}^{{m \times n}}}$, удовлетворяет соотношению $\sum\limits_{s \in {{L}_{i}}} {{{K}_{{is}}}} = {{N}_{i}} - {{S}_{i}}$, а ki = ${{\tilde {k}}_{i}} + {{S}_{i}}{{D}_{i}}$ + + $\sum\limits_{s \in {{L}_{i}}} {{{K}_{{is}}}} {{D}_{s}}$ независимо от управлений лидеров, т.е. удовлетворяет и свойству ${{({\text{ЛВУ}})}_{1}}$.

Сформулируем ряд следствий из этой теоремы.

Так как все необходимые условия получены из свойства ${{({\text{ЛВУ}})}_{0}}$, то верно следующее утверждение.

Следствие 1. Свойства ${{({\text{ЛВУ}})}_{1}}$, $({\text{ЛВУ}})$ и ${{({\text{ЛВУ}})}_{0}}$ эквивалентны.

Замечание 3. Отметим также, что если перебирать все $(m \times n)$-матрицы Si, стабилизирующие пары матриц $({{A}_{i}},{{B}_{i}})$ ($i > {{l}_{0}}$), т.е. такие, что матрица ${{A}_{i}} + {{B}_{i}}{{S}_{i}}$ – гурвицева, а также все Ni, ${{\tilde {k}}_{i}}$, являющиеся решениями уравнений (7), (8) и все наборы матриц ${{{\text{\{ }}{{K}_{{is}}}{\text{\} }}}_{{s \in {{L}_{i}}}}}$, ${{K}_{{is}}} \in {{R}^{{m \times n}}}$, удовлетворяющие соотношению $\sum\limits_{s \in {{L}_{i}}} {{{K}_{{is}}}} = {{N}_{i}} - {{S}_{i}}$, то по формуле

${{u}_{i}}({{{\text{x}}}_{1}}, \ldots ,\;{{{\text{x}}}_{i}}) = {{S}_{i}}({{{\text{x}}}_{i}} + {{D}_{i}}) + \sum\limits_{s \in {{L}_{i}}} {{{K}_{{is}}}} ({{{\text{x}}}_{s}} + {{D}_{s}}) + {{\tilde {k}}_{i}}$
получим все управления, при которых выполняется свойство ${{({\text{ЛВУ}})}_{1}}$ (естественно, при выполнении условий 3) и 4) теоремы).

Далее заметим, что в случае, когда у формации больше одного лидера, то в теореме 1 условие 1) является следствием условий 2) и 4), так как можно взять ${{S}_{i}} = {{N}_{i}}$ и, кроме того, при этом можно взять ${{K}_{{is}}} = 0 \in {{R}^{{m \times n}}}$ для всех $i > {{l}_{0}}$, $s \in {{L}_{i}}$. Тогда получаем следующее

Следствие 2. Для линейной формации, имеющей больше одного лидера, выполняется свойство линейной внутренней устойчивости тогда и только тогда, когда выполнены следующие условия:

1) ${{A}_{i}} = {{A}_{1}}$ для всех $i = 1,2,\; \ldots ,\;{{l}_{0}}$ и матрица A1– гурвицева;

2) выполняются условия условий 2) и 3) теоремы 1.

При этом если Ni и ${{\tilde {k}}_{i}}$ удовлетворяют усло-вию 2) теоремы 1, то управления ведомых, обеспечивающие внутреннюю устойчивость, можно выбрать в виде

${{u}_{i}}({{{\text{x}}}_{i}}) = {{N}_{i}}({{{\text{x}}}_{i}} + {{D}_{i}}) + {{\tilde {k}}_{i}},\quad i > {{l}_{0}},$
независимо от состояний ведущих и управлений лидеров.

Замечание 4. Рассмотрим самую простую среди формаций, имеющих больше одного лидера. Она состоит из трех объектов: 1-й и 2-й – лидеры и ведущие для 3-го. По следствию 2 необходимыми и достаточными условиями ее линейной внутренней устойчивости будут следующие (учитывая, что для этой формации ${{D}_{1}} = {{D}_{2}} = 0$, $p(3)$ = 2 и ${{D}_{3}} = {{d}_{{32}}}$):

1) ${{A}_{2}} = {{A}_{1}}$ и матрица A1 – гурвицева;

2) существуют решения уравнений ${{B}_{3}}{{N}_{3}}$ = A1A3 и ${{B}_{3}}{{\tilde {k}}_{3}} = {{A}_{3}}{{d}_{{32}}}$;

3) ${{d}_{{31}}} = {{d}_{{32}}}$.

Однако очевидно, что (необходимое) условие ${{d}_{{31}}} = {{d}_{{32}}}$ невыполнимо и бессмысленно для реальных объектов. Причем можно показать, что подобная коллизия будет возникать в любой формации, имеющей больше одного лидера. Еще одним странным (противоречащим интуиции) фактом является возможность обеспечить внутреннюю устойчивость (в случае ее наличия) управлением, не зависящим от состояний ведущих (как в следствии 2).

Таким образом, данная концепция внутренней устойчивости не подходит для многолидерных формаций. Причиной этого, видимо, является то, что мы пытаемся выдержать заданное расположение ведомого объекта одновременно относительно всех его ведущих (т.е. пытаемся угнаться сразу за несколькими зайцами). Возможно, более приемлемой в данном случае может оказаться концепция внутренней устойчивости из [5], в которой используется некоторое агрегированное отклонение i-го объекта от идеального состояния вида ${{e}_{i}} = \sum\limits_{j \in {{L}_{i}}} {{{P}_{{ij}}}} {{{\text{z}}}_{{i,j}}}$, где Pij – некоторые матрицы-проекторы, и требуется выполнение оценки вида (2) для вектора, составленного из ei, $i > {{l}_{0}}$, а не всей матрицы z.

Далее рассмотрим случай формаций, для которых условие 3) теоремы 1 отсутствует. Это формации, граф которых является входящим деревом [15, с. 235], т.е. слабо связным бесконтурным графом, не имеющим полуконтуров, со стоком. Согласно [15, с. 236, теорема 16.4′] у формации с таким графом будет один лидер и, кроме того, у всех остальных объектов будет единственный ведущий, т.е. ${{L}_{i}} = {\text{\{ }}p(i){\text{\} }}$ при $i > 1$. В этом случае условие 3) теоремы 1 всегда выполнено и набор матриц ${{{\text{\{ }}{{K}_{{is}}}{\text{\} }}}_{{s \in {{L}_{i}}}}}$ состоит из одной матрицы ${{K}_{{ip(i)}}}$. Поэтому получаем следующее

Следствие 3. Для линейной формации, граф которой является входящим деревом, выполняется свойство линейной внутренней устойчивости тогда и только тогда, когда выполняются условия 1), 2) теоремы 1.

При этом если Si, Ni и ${{\tilde {k}}_{i}}$ удовлетворяют условиям 1), 2) теоремы 1, то управления ведомых, обеспечивающие внутреннюю устойчивость, можно выбрать в виде

$\begin{gathered} {{u}_{i}}({{{\text{x}}}_{i}},{{{\text{x}}}_{{p(i)}}}) = {{S}_{i}}({{{\text{x}}}_{i}} + {{D}_{i}}) + \\ + \;({{N}_{i}} - {{S}_{i}})({{{\text{x}}}_{{p(i)}}} + {{D}_{{p(i)}}}) + {{{\tilde {k}}}_{i}},\quad i > 1, \\ \end{gathered} $
независимо от управления лидера.

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

  1. Balch T., Arkin R.C. Behavior-Based Formation Control for Multirobot // IEEE Transactions on Robotics and Automation. 1998. V. 14. № 6. P. 926–939.

  2. Tanner H.G., Pappas G.J. Formation input-to-state stability // Proc. 15th IFAC World Congr. Autom. Control. Barcelona, 2002. P. 1512–1517.

  3. Tanner H.G., Kumar V., Pappas G.J. Stability properties of interconnected vehicles // Proc. 15th International Symposium on Mathematical Theory of Networks and Systems. South Bend, Indiana, 2002. P. 1–12. Paper number 4615–2. [CD-ROM].

  4. Tanner H.G., Pappas G.J., Kumar V. Input-to-state Stability on Formation Graphs // Proc. 41st IEEE Conference on Decision and Control. Las Vegas, NV, 2002. P. 2439–2444.

  5. Tanner H.G., Pappas G.J., Kumar V. Leader-to-formation stability // IEEE Transactions on Robotics and Automation. 2004. V. 20. № 3. P. 443–455.

  6. Oh K.K., Park M.C., Ahn H.S. A survey of multi-agent formation control // Automatica. 2015. V. 53. P. 424–440.

  7. Lü J., Chen F., Chen G. Nonsmooth leader-following formation control of nonidentical multi-agent systems with directed communication topologies // Automatica. 2016. V. 64. P. 112–120.

  8. Sontag E.D. Smooth stabilization implies coprime factorization // IEEE Trans. Automat. Control. 1989. V. 34. № 4. P. 435–443.

  9. Sontag E.D., Wang Y. On characterizations of the input-to-state stability property // Systems & Control Letters. 1995. V. 24. № 5. P. 351–359.

  10. Дашковский С.Н., Ефимов Д.В., Cонтаг Э.Д. Устойчивость от входа к состоянию и смежные свойства систем // АиТ. 2011. № 8. С. 3–40.

  11. Васильев С.Н., Козлов Р.И., Ульянов С.А. Анализ координатных и других преобразований моделей динамических систем методом редукции // Труды института математики и механики УрО РАН. 2009. Т. 15. № 3. С. 38–55.

  12. Васильев С.Н., Козлов Р.И., Ульянов С.А. Устойчивость многорежимных формаций // ДАН. 2014. Т. 455. № 3. С. 269–274.

  13. Ul’yanov S., Maksimkin N. Formation path-following control of multi-AUV systems with adaptation of reference speed // Mathematics in Engineering, Science and Aerospace. 2019. V. 10. № 3. P. 487–500.

  14. Халил X.К. Нелинейные системы. М.: РХД, 2009.

  15. Харари Ф. Теория графов. М.: Мир, 1973.

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

Инструменты

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