УДК 007.52:629.3.05 DOI 10.18522/2311-3103-2022-1-121-139
В.Х. Пшихопов, М.Ю. Медведев
СРАВНИТЕЛЬНЫЙ АНАЛИЗ ЦЕНТРАЛИЗОВАННОГО И ДЕЦЕНТРАЛИЗОВАННОГО АЛГОРИТМОВ ДВИЖЕНИЯ СТРОЕМ БЛА МУЛЬТИКОПТЕРНОГО ТИПА*
Развитие робототехнических комплексов делает актуальным их групповое применение для решения различных задач. Эффективность выполнения задач обнаружения и определения координат объектов группой роботов существенно зависит от точности поддержания заданного строя. В этой связи практический интерес представляет задача определения алгоритмов планирования движения, обеспечивающих наибольшую точности поддержания заданного строя. Данная статья посвящена исследованию точности поддержания строя группой БЛА мультикоптерного типа с использованием централизованного алгоритма планирования движения и децентрализованного алгоритма. В централизованном алгоритме используется ведущий БЛА, который передает свои координаты ведомым БЛА. На основании полученных координат и заданной структуре строя ведомые БЛА планируют свое движение. В децентрализованной систем соседние БЛА группы передают свои координаты друг другу, на основании чего планируется движение отдельного БЛА. Точность исследуется в зависимости от погрешностей навигационной системы и частоты обновления данных о положении ведущего или соседних БЛА. Полагается, что БЛА группы в дискретные моменты времени определяют свои координаты, используя внешнюю навигационную систему. Централизованный и децентрализованный алгоритмы отрабатываются одинаковой системой управления движением. Алгоритмы исследуются в данной статье методами численного моделирования. В процессе моделирования учитываются модели кинематики, динамики и исполнительных механизмов, а также модели формирования погрешностей навигационной системы. Показано, что децентрализованный алгоритм группового планирования движения обеспечивает более высокую точность по сравнению с централизованным алгоритмом. Однако техническая реализация децентрализованного алгоритма более сложна с точки зрения организации системы групповой связи. В централизованной систем должна быть реализована передача данных от ведущего БЛА ведомым. В децентрализованной системе требуется реализовать сетевую связь.
Группа БЛА; формирование строя; централизованный алгоритм; децентрализованный алгоритм; групповое управление.
V.Kh. Pshikhopov, M.Yu. Medvedev
COMPARATIVE ANALYSIS OF CENTRALIZED AND DECENTRALIZED ALGORITHMS FOR THE MOVEMENT OF MULTICOPTER-TYPE UAVS
The development of robotics makes their group application relevant for solving various tasks. The effectiveness of performing the tasks of detecting and determining the coordinates of objects by a group of robots significantly depends on the accuracy of maintaining a given formation. In this regard, the task of determining motion planning algorithms that ensure the greatest accuracy of maintaining a given formation is of practical interest. This article is devoted to the study of the accuracy of maintaining the formation of a multicopter-type UAV group using a centralized motion planning algorithm and a decentralized algorithm. The centralized algorithm uses a master UAV, which transmits its coordinates to the slave UAVs. Based on the coordinates obtained and the given framework of the formation, the guided UAVs plan their movement. In a decentralized system, neighboring UAV groups transmit their coordinates to each other, on the basis of which the movement of a separate UAV is planned. The accuracy of the control system is investigated depending on the errors of the navigation system and the frequency of updating data on the position of the leading or neighboring UAVs. It is assumed that the group's UAVs determine their
* Работа при поддержке РНФ, проект № 22-29-00533.
coordinates in discrete moments of time using an external navigation system. Centralized and decentralized algorithms are worked out by the same motion control system. The algorithms are investigated in this article by numerical modeling methods. In the process of simulation, models of kinematics, dynamics and actuators are taken into account, as well as models for the formation of errors in the navigation system. It is shown that the de-centralized algorithm of group motion planning provides higher accuracy compared to the centralized algorithm. However, the technical implementation of a decentralized algorithm is more complicated from the point of view of organizing a group communication system. In a centralized system, data transmission from the master UAV to the slave should be implemented. In a decentralized system, it is required to implement network communication.
UAVs group; formation control; centralized algorithm; decentralized algorithm; group control.
Введение. Поддержание заданного строя актуально при решении задач движения колонной [1], слежения за целью [2], ретрансляции связи [3], поиска объектов [4], проведения георазведки [5]. Перспективной задачей является создание на базе группы БЛА многопозиционных радиолокационных или оптических систем [6-9]. Такие системы, обладая преимуществами многопозиционных систем, обеспечивают быстрое развертывание в заданном направлении и перемещение в требуемую область. Вместе с тем, в таких группах высокие требования предъявляются к системам управления и планирования движением. Если система навигации и управления обеспечивают высокую точность поддержания взаимного положения элементов группы в пространстве и синхронизированы между собой, что такая группа БЛА способная выполнять задачи радиолокации с высокой точностью. В этой связи в данной статье рассматриваются два подхода к пространственной координации - с выделением ведущего БЛА и с децентрализованной структурой.
В задаче создания многопозиционной радиолокационной системы на базе БЛА строй представляет собой линию, либо несколько линий, следующих друг за другом. В этой связи имеет смысл рассматривать метод виртуальных структур для формирования строя. При этом виртуальная структура может формироваться относительно ведущего БЛА [10-12] или децентрализованными методами [13-15]. При выделении ведущего БЛА его положение является точкой, относительно которой строится виртуальная структура. Ведущий БЛА передает свои координаты другим членам группы, которые рассчитывают собственные желаемые положения. При использовании децентрализованных методов каждый БЛА определяет координаты соседних БЛА и рассчитывает текущее желаемое положение.
В система с ведущим БЛА достаточно передавать только его координаты остальным членам группы. Это снижает требования к системе групповой связи. В режиме передачи данных работает только система связи ведущего БЛА. При большом расстоянии между отдельными БЛА группы может потребоваться ретрансляция координат ведущего БЛА. Также отметим, что в централизованной системе ведущий БЛА движется по заданной программе, не реагируя на текущее состояние остальных членов группы.
В централизованной систем задачей ведомых БЛА является слежение за точкой, положение которой определяется ведущим БЛА. В работе [1] рассмотрена задача слежения за роботом-лидером на основе оценок взаимного положения и ориентации, получаемых с помощью бортовой видеокамеры. Слежение за лидером осуществляется на основе уравнений кинематики материальной точки на плоскости [16]. Предложены алгоритмы слежения за лидером в системе координат робота-лидера и оценивания линейной и угловой скоростей робота-лидера. Алгоритм оценивания базируется на расширенном наблюдателе [17], обеспечивающем оценку скоростей и внешних возмущений. Алгоритм слежения базируется на методе линеаризации с обратной связью [18]. Предложен метод снижения ошибки слежения за лидером, учитывающий ограничения визуальной системы определения вза-
имного положения и ориентации. Однако использование уравнений кинематики материальной точки и упрощенных уравнений динамики ограничивает результаты статьи [1], особенно, при высоких требованиях к точности слежения за лидером, в трехмерных средах, а также в режимах маневрирования реальных роботов.
В статье [2] рассмотрена задача слежения за лидером группой роботов, движущихся заданным строем. При этом построение осуществляется за заданное конечное время, т.е. задача решается в терминальной постановке [19, 20]. За это время лидер должен оказаться в геометрическом центре группы. Все роботы группы описаны уравнениями кинематики материальной точки. При решении задачи делается предположение о том, что координаты и ориентации лидера группы известны только части ведомых роботов. Вместе с тем следует отметить, что задача решена без учета динамики движения.
В статье [21] также рассматривается задача движения за лидером на плоскости. Модель роботов представлена уравнениями кинематики и динамики материальной точки. Изучается проблема слежения за лидером в условиях, когда координаты и ориентация робота-лидера передаются через систему связи с ошибками. Предложен адаптивный динамический регулятор, использующий аппроксимирующую сеть радиальных базисных функций. Как в ранее рассмотренных статьях, основным недостатком является рассмотрение уравнений материальной точки, что не позволяет учесть ограничений, связанных с особенностями физической реализации роботов.
Также известны другие работы по групповому слежению за лидером, которые учитывают возможность не измерять ориентацию робота [22], обеспечивать устойчивость и робастность гладких траекторий движения [13, 23], ограничивать передачу информации о расстоянии между лидером и ведомыми [24]. Во всех указанных работах использованы уравнения кинематики материальной точки в двумерной среде, что не позволяет полноценно оценить возникающие в системе погрешности.
Аналогичная ситуация имеет место с методами децентрализованного управления. Так в работе [12] рассматриваются агенты, описываемые упрощенными уравнениями динамики материальной точки, для которых предлагается децентрализованный алгоритм, базирующийся на понятии виртуальных лидеров, которые рассчитываются каждым агентом в отдельности и являются его целевым положением.
В работе [14] представлены алгоритмы децентрализованной системы управления, в которых учитываются уравнения кинематики и динамики движения твердого тела. Указанные алгоритмы не предполагают выделение лидера или предварительно заданного строя. Роботы автоматически распределяются в заданной области и перестраиваются при обнаружении препятствий. Однако, используются сложные алгоритмы оптимизации, включая триангуляцию Делоне, минимаксную оптимизацию и нейронные сети.
В работе [15] рассмотрена группа роботов в 3-Б пространстве, которые описаны уравнениями твердого тела. Группа автоматически распределяется в прямоугольной области, а затем движется, поддерживая взаимное положение и перестраиваясь при встрече с препятствиями. Для создания формации используется триангуляция Делоне. Обход препятствий осуществляется с использованием неустойчивых по расстоянию до препятствий режимов движения, основанных на использовании третьей теоремы Ляпунова [17, 25]. Это позволяет не использовать сложных в вычислительном плане нейронных сетей на борту робота.
В работе [26] рассмотрен метод децентрализованного управления движением группы однородных роботов при перемещении группы грузов. Роботы описаны уравнениями кинематики и упрощенной динамики материальной точки. Упрощения состоят в том, что существует преобразование, приводящее исходную систему
к линейному виду. Предлагаются алгоритмы группового управления, получаемые в результате решения уравнения Риккати [27]. В итоге роботы подразделяются на подгруппы, которые формируют строй, осуществляющий перемещение грузов.
Проведенный обзор позволяет сделать выводы о том, что при разработке систем группового управления используются достаточно простые модели движения, не учитывающие динамические особенности мобильных роботов. В данной статье проводится сравнительный анализ методов централизованного и децентрализованного управления группой БЛА при решении задачи движения строем. При этом учитываются погрешности навигационных систем и запаздывания в каналах связи.
II. Математическая модель БЛА и формулировка задачи. Рассматривается БЛА мультикоптерного, которые описываются уравнениями кинематики и динамики твердого тела и исполнительных механизмов [28-30]
У = Я (у) х, (1)
X = М - 1 (К5 + Га), (2)
Т58 = Кии - 5, (3)
где у - вектор линейных и угловых координат БЛА в неподвижной системе координат размерностью 6; x - вектор линейных и угловых скоростей в связанной системе координат размерностью 6; R(y) - матрица кинематики; 5 - вектор координат исполнительных механизмов; Fd - вектор прочих сил и моментов, действующих на БЛА; M - матрица инерционных параметров; и - вектор управлений; ^ - матрица связи управляющих воздействий и переменных состояния исполнительных механизмов БЛА; ^ - матрица постоянных времени исполнительных механизмов БЛА.
В случае, если угловое положение описывается углами Эйлера [29], то матрица кинематики R(y), которая определяет связь между неподвижной и связанной системами координат, имеет вид:
А 03х3]
.Озхз Аа
R =
(4)
=
A=
cos у cos д sind
— sin у cos д —
— COS Y|/ Sin $ COS у +
+ sin у sin у cosd cosy cos у sin у + + sin у sin д cos у
cos у sin $ sin у + + sin у cosy — cos ß sin у
COS Y|/ COS У —
— sin у sin$ siny.
cosy cos$
siny — tan дcos у
siny cos$
cosy tan i9 sin y.
(5)
(6)
(7)
Матрица инерционных параметров сбалансированного БЛА имеет вид
М =
т 0 0 0 0 - -myT
0 m 0 0 0 0
0 0 m myT 0 0
0 0 myT Jx -Jxy —Jxz
0 0 0 -Jxy Jy —Jyz
—туТ 0 0 -Jxz Jyz Jz
(8)
где m - масса БЛА; Jx, Jy, Jz, Jxy, Jxг, Jyz - моменты инерции относительно осей связанной системы координат БЛА; ут - отклонение масс БЛА относительно начала координат связанной системы по оси OY.
Вектор сил и моментов Fd, действующих на БЛА, включает составляющие
— + F; + Fh,
(9)
где Fg - вектор силы тяжести; FJ - вектор инерционных сил вращающегося БЛА; Fh - вектор аэродинамических сил и моментов.
Вектор силы тяжести определяется в неподвижной системе координат имеет компоненты [0, -mg, 0]T. С помощью матрицы A (6) можно найти его компоненты в связанной системе координат K
— sin$ - cosд cosy cos д sin у ут cos д sin у 0
ут sin д
Вектор инерционных сил вращающегося БЛА равен
-m(azVy - (ayVz)-
Fa =
mg.
(10)
F] =
m((oxVz - (ÜzVx) m(axVy - G>yVx)
(Jy ~ Jz)®>z
(Jz -Jx)®z®x . (Jx ~ Jy)®x®y J
(11)
Вектор аэродинамических сил и моментов, действующих на БЛА, определяется выражением [31]
Fh = 0,5р72
Sex Scy Scz Umx Umy Urriv
(12)
где - характерная площадь БЛА; U - характерный объем БЛА; V - модуль вектора скорости движения БЛА относительно среды; ох, су, ох, шх, ту, тх - аэродинамические коэффициенты.
Если углы атаки а и скольжения р малы, то аэродинамические коэффициенты приближенно рассчитываются по следующим формулам [32]
= сх0 + + с£а2 + cßß2, а|а| + c®zg
cßß + czß2ß|ß|+czfflyräy,
су — Суо + Суа + Су
mx — mxß + m'
шу = ШуР + m^ßlßl + ШуУЮу,
ш2 = ш20 + ш^а + ш^2а|а| + шу2ю2, где ®ж = ax/V; ®у = ®у/7; ®z =
Матрица K5 зависит от типа БЛА. Для квадрокоптера, обладающего четырьмя винтами, вектор 5 координат исполнительных механизмов и указанная матрица определяются выражениями [33, 34]:
„ß2f
(13)
(14)
(15)
(17)
(18) (19)
К-: — Ке
kvlvsin (<PV)
-к
kvlvsin (ф )
kvlvsin (cpv)
К
-kvlvsin (<pv)
-kvlv sin (<pv) ~bv -kvlvsin (ф )
-kvlvsin (<pv)
К
kvlvsin (ф ) .
(20)
5 =
9 0 0 0 0
af 1 0 0 0
Юз , Kf = 0 0 0 1 0 0 0 0
со|. 0 0 1 0
0 0 0 1
(21)
где Д, - расстояние от начала системы координат, связанной с БЛА до пропеллеров; фv - угол, под которым размещены винты; кп Ь„ - аэродинамические коэффициенты винтов, - скорости вращения винтов.
Матрица постоянных времени двигателей является диагональной, т.к. винты БЛА работают независимо друг от друга, т.е.
г51 0 0 0
О Г51 0 0
ТS =
'81
о
1 81-1
(22)
Ки =
(23)
о о
0 0 о
Матрица коэффициентов передачи исполнительных механизмов также является диагональной в силу того, что винты БЛА работают независимо друг от друга:
ки1 0 оо
О ки2 0 0 0 0 К з о 0 0 0 киАЛ
Структура модели навигационной системы представлена на рис. 1. Она включает в себя блок вычисления координат и модуль параметров навигационной системы.
Входными переменными модели навигационной системы являются: 1) точные значения курса у, угла крена у, угла тангажа 9 и путевого угла ур БЛА в неподвижной системе координат; 2) точные значения угловых скоростей ах, ау, ах относительно осей связанной с БЛА системы координат; 3) точные значения угловых ускорений ах, ау , а2 относительно осей связанной с БЛА системы координат; 4) точные значения координат х, у, и г БЛА в неподвижной системе; точные значения линейных скоростей Ух, Уу, У БЛА в подвижной системе координат и скоростей X, у , г БЛА в неподвижной системе координат; 5) точные значения линейных ускорений БЛА в связанной системе координат и линейных ускорений
Х,у,г в неподвижной системе координат; 6) признак наличия сигналов внешней навигационной системы ПСНС; 7) точное время 1
Выходными переменными модели навигационной системы являются: 1) измеренные значения курса уи, угла крена уи, угла тангажа 9и и путевого угла ури БЛА в неподвижной системе координат; 2) измеренные значения угловых скоро-ю2и относительно осей связанной с БЛА системы координат; 3) изме-
стеи юх
ю„
относительно осей связанной с
ренные значения угловых ускорений ахи, а уи БЛА системы координат; 4) измеренные значения координат хи, уи, и ги БЛА в неподвижной системе; 5) измеренные значения линейных скоростей Ухи, Ууи, Уи БЛА в связанной с ним подвижной системе координат и скоростей Хи, уи, ¿и БЛА в неподвижной системе координат; 6) измеренные значения линейных ускорений Ухи, Ууи, Ц,и БЛА в связанной системе и линейных ускорений Хи, уи, ги в неподвижной системе координат; 7) время проведения текущего измерения координат /и.
Рис. 1. Структура модели навигационной системы
Измеренные углы Эйлера у, у, 9 и путевой угол ур на выходе модуля вычисления координат определяются по следующим формулам:
уи = у + (ту) + (С) , (24)
в и = в + Ев Д „ 1 (^) , (25)
Уи = у + Яуйу 1 (ту), (26)
Урл — Ч'р -^грр 1 С^гфф ) + Д ур 2 ( 0 , (27)
где Еу1, Е9, Е9у, Еур - СКО погрешностей измерения навигационной системы; RуЬ Яэ1, Яу1, Яур1 - функции генерации случайной помехи; Rу2, Яур2 - функции генерации систематической погрешности при отсутствии сигналов внешней навигации; ту, т9, ту, тур - параметры, определяющие время корреляции случайной составляющей измеряемых сигналов.
Значение константы Еу определяется погрешностью навигационной системы. Данная погрешность, в зависимости от класса используемой навигационной системы, может составлять от нескольких угловых минут до нескольких градусов.
Функция Rу1(ту) моделируется как коррелированный случайный шум [35], полученный путем пропускания белого шума через инерционное звено первого порядка [36]. Уравнение, описывающее такой шум, имеет вид
^ = -с).
аь у ъу) (28)
Ду1 =
где - переменная состояния инерционного звена первого порядка; £у - случайный сигнал типа белый шум.
Функция Rу2 генерирует нарастающую погрешность определения курса в соответствии с выражением:
Д у2 ( С) = £у2£:, ¿Д С < С < (I + 1 ) Д С, I = 0 , 1 , 2 . . . (29)
где Еу2 - параметр, определяющий скорость нарастания погрешности, Д1 - интервал коррекции навигационной системы от внешних сигналов.
Угол тангажа описывается уравнением (25), в котором отсутствует составляющая, нарастающая во времени. Функции Д ^ 1 (тв) генерируется с помощью уравнения вида (28), отличия состоят только в другом значении времени корреляции тв и СКО Е9. Аналогичным образом формируется модель измерения угла крена (26), для которой используется уравнение, аналогичное (28), с параметрами ту и Яу:
Путевой угол определяется аналогично углу рысканья, с использованием параметров Еи Еад2 в уравнениях (28), (29). Отметим, что путевой угол является суммой углов рысканья и скольжения [37].
Угловые скорости измеряются непосредственно датчиками угловых скоростей, основными погрешностями которых являются дрейф нуля и нестабильность коэффициента передачи [38]. Указанные погрешности моделируются с помощью следующих соотношений
ихи = их + Дих' (30)
ЮуИ = Юу+Д юу, (31)
ю2и = + , (32)
Rr,
D
^fflx 2' R «У 2
(s) =Г^шу 2' R «йМ = ~Л«
гюх2 Люх2' * ' ®у2\" ) „^и у2' ' -)
где СюхЬ Сюх2' СюуЬ Сюу2' Сгаг1' Сш22 - белые шумы.
Акселерометры измеряют так называемые кажущиеся ускорения. В этой связи показания акселерометров описываются следующими выражениями:
%И = % - 9 ъ 1 п $и + Д „х, (33)
^И = % - 9 соъ Уи с о ъ )9и + Д „у, (34)
Ити = К + 9 ъ 1пуи с о ъ ^ + Д„г. (35)
где ЯУу, ^Уг - случайные погрешности дрейфа нуля и неточностей коэффициентов передачи акселерометров и перекосов его осей.
Ускорения в неподвижной системе координат определяются расчетным путем на основе следующего выражения
Х0и "Кси"
Уои = Л "уи
-20и-
(36)
где матрица А определяется выражением (6).
Скорости БЛА в связанной системе координат вычисляются путем интегрирования ускорений в соответствии с выражениями
V =V (¿At) + f(i+1)At
"хи 'ihviui-/ ~ Jj
■>iAt
dx , i = 0 , 1 , 2 . . .
i-
\t
Vm = Vzll(iAt) + ,zll.
Vyll = Vyll(iAt) + J.(A+1)AtK,„dx
J ¡At Уи
r(i+l)&t
(37)
(38)
(39)
Ити ^ .
Скорости БЛА в неподвижной системе координат определяются в соответст вии с выражением
(40)
где матрица А , как и ранее, вычисляется в соответствии с выражением (6).
Координаты БЛА определяются путем интегрирования скоростей, согласно выражениям:
Х0и хи
Уои = А ^уи
.¿„И. ^и
х0и = х0и (¿А 0 + Х0иХ0'
Уои=Уои ( iA 0 + L, Уо и*0'
20и = Z0H(iAt) + Г.'
Ji&t (i + l)At . ¡At 20их0-
(41)
(42)
(43)
Структура модели системы связи представлена на рис. 2.
S, Модель блока
передачи данньк
s. Модель блока
приема данньк
Рис. 2. Структура модели системы связи
Модель системы связи состоит из модели блока передачи данных и модели блока приема данных. На вход модели блока передачи данных поступает сообщение Sb которое содержит следующую информацию: 1) полезное сообщение S1 и его размер D1 в килобайтах; 2) время появления сообщения tx; 3) текущие координаты БЛА, передающего сообщение.
Выходом модели блока передачи данных является сообщение S2, которое содержит информацию: 1) сообщение S2 и его размер D2 в килобайтах; 3) время начала передачи сообщения t2ii; 3) время окончания передачи сообщения t2s; 4) координаты передающего БЛА в момент начала передачи сообщения.
Сообщение S2 дублирует содержание сообщения Si при моделировании. Размер D2 сообщения S2 вычисляется в соответствии с выражением
D2 = KSDU (44)
где Ks - коэффициент, задающий степень увеличения размера передаваемого сообщения за счет служебной информации и кодирования. Принято Ks = 2.
Время начала передачи сообщения вычисляется в соответствии с выражением
t2H = t1+ Atb (45)
где At1 - время необходимое для подготовки к передаче данных.
Время окончания передачи сообщения t2s определяется выражением
t2K = t2H + D2/7nep. (46)
Входом блока приема данных является сообщение S3, которое равно сообщению сообщения S2, но время его приема вычисляется в соответствии с выражением
t3n = t2K + r/Kiep-
Ставится задача исследования точности поддержания строя группой из Nr БЛА при движении с постоянной скоростью в заданном направлении, в зависимости от точности навигационной системы, частоты получения данных о собственном положении г-го БЛА и частоты обновления данных о положении других БЛА группы, для следующих подходов поддержания строя: 1) централизованный подход, заключающийся на передаче координат лидера ведомым БЛА; 2) децентрализованный подход, использующий обмен координатами между соседними БЛА.
III. Алгоритмы формирования строя. Будем считать, что в начальный момент времени БЛА группы построены в линию на расстоянии гш. Рассмотрим вначале вариант движения с использованием лидера. Схема информационного обмена данными в такой группе представлена на рис. 3. Согласно этой схеме из пункта наземного управления поступает задание для БЛА-лидера. Задание включает в себя направление и скорость полета группы, а также требуемое расстояние, на котором БЛА должны находится друг от друга.
Рис. 3. Информационный обмен в системе с лидером
Определение собственных координат каждого БЛА определяется с помощью сигналов внешней навигационной системы. Кроме того, БЛА-лидер передает каждому ведомому БЛА свои координаты (х1, у1, х1), на основании которых каждый БЛА определяет свое желаемое положение в строю, используя выражения:
(47)
(48)
(49)
(50)
х" = V0t, х" = хь i = 2 , Nr
У1
о _
у0, i = l,Nr
_
z0 + (i~ 1 К, i = l,Nr
Vi
= у0, i = 2 , Nr
-xt
Рп = У? -у
к
где г0 - заданное положение первого БЛА; р - заданное желаемое направление движения группы БЛА. При этом принято, что система координат привязана к положению лидера в начальный момент времени [39], а ось ОХ направлена вдоль желаемого движения группы.
Проекции желаемой скорости БЛА группы в связанной системе координат вычисляются в соответствии с выражением
У°=АТУоРп1/| | рп | | , 1=х,уА (51)
где У0 - заданная скорость движения группы; матрица кинематики А определяется в соответствии с выражением (6), а направляющий вектор рп и его элементы рп вычисляются по формуле
(52)
Желаемый угол тангажа БЛА-лидера определяется выражением ( - V* - еа2), V 7Х° - Ух1 > еаг,
= Ь (7Х° - Ух1 + ей2), V 7Х° - Ух1 < —ей2, I 0, -ей2< 7Х° - Ух1< еаг,
где - зона нечувствительности, вводимая для фильтрации высокочастотных измерительных шумов [40]; Ух - проекция скорости БЛА на ось ОХ связанной системы координат; - параметр регулятора.
При этом угол тангажа дополнительно ограничен значениями ±дт ах.
Аналогичным образом, формируются требования к углу крена
" V* - ей2),\/ У2° - У21 > ей2, уо = \ ку СУ» - Ух1 + ей2), ЧУ° - У Л < - ей2, (54)
0, —ей2< У21 — У21< ей2,
(53)
где У2 I - проекция скорости БЛА на ось 02 связанной системы координат; ^0 - целевой угол крена; ку - параметр регулятора.
Угол крена, как и угол тангажа, дополнительно ограничен значениями ±/п ах. Ограничения на углы крена и тангажа связаны с особенностями управления. При превышении указанными углами некоторых максимальных значений, БЛА мультикоптерного типа падает [41-43].
Желаемая скорость полета в установившемся режиме определяется проекцией скорости , которая равна
УУ1 = V". (55)
Рассмотрим теперь вариант движения с роевыми принципами формирования группы [44, 45]. Схема информационного обмена данными в такой группе представлена на рис. 4. Согласно этой схеме, из пункта наземного управления поступает задание всех БЛА. Задание включает в себя направление и скорость полета группы, а также требуемое расстояние, на котором БЛА должны находится друг от друга.
Как и ранее, определение собственных координат каждого БЛА определяется с помощью сигналов внешней навигационной системы. Кроме того, все БЛА обмениваются своими координатами с соседними БЛА группы.
навигация £-\-^ Блд.3
Координаты _БЛА
БЛА-4
Рис. 4. Информационный обмен в роевой системе
Высота движения каждого БЛА для роевого принципа построения группы, как и для движения с лидером, задается выражением (48). Желаемое направление движения для всех БЛА группы вычисляется в соответствии с выражением (50).
Желаемая координата по оси 02 вычисляется в соответствии выражением
г° = 1 + гш, г" = 0, I = 2ТЖГ. (56)
Желаемая координата по оси 0Х в случае роевого принципа формирования строя, определяется выражениями
х° = У0С + ку (х2—х^, (57)
Х °1г = + кУ (ХУг _ 1 — XNr^' (58)
х° = 70 г + — , (59)
где .
Желаемые углы дифферента, крена и рысканья вычисляются в соответствии с выражениями (53) и (54), а вектор желаемых скоростей БЛА определяется выражением (51). Таким образом, при роевом управлении строем, каждый БЛА вычисляет свое положение по оси 02 и оси 0Х относительно соседних БЛА, а не БЛА-лидера.
Система управления движением строится на базе метода позиционно-траектор-ного управления. Описание процедуры синтеза алгоритма управления движением БЛА, а также алгоритмы адаптации системы управления, представлены в работах [15, 28, 43]. Общий вид алгоритма управления имеет следующую форму
A2R
M~1(KS5 + Fd) =
A2R + T2Vtr + №
(60)
0 0 0 1 0 0
¿2 = 0 0 0 0 1 0 , ¿3 =
.0 0 0 0 0 1. -t.
где 7\, Т2, Т3 - матрицы регулятора; Д - матрица, элементы которых являются производными по времени элементов матрицы Я;
4V = А2у +
Л4 = [0 1 0 0 0 0], А5 = [-7»], УР = А4у + А5.
IV. Результаты исследования алгоритмов формирования строя. Сравнение алгоритма управления с лидером и роевого алгоритма проведено с использованием численного моделирования.
Моделирование проведено при следующих параметрах БЛА и системы управления: начальные положения БЛА (0; 100; 0) - лидер, (0; 100; 50), (0; 100; 100) и (0; 100; 150); желаемое направление полета ц» = 0; желаемая высота полета у» = 100 м; расстояние между БЛА гш = 50 м; желаемая скорость полета 7» = 3 м/с; параметры регуляторов углов крена и тангажа ку = 0.1, = 0.1; значения углов тангажа и крена ограничены значениями +0,3 рад; максимальная скорость полета БЛА 4 м/с; инерционные параметры БЛА т=5 кг, Jx=5, Jy=0.1, Jz=5; максимальная тяга винтов 100 Н. Постоянная времени замкнутой системы управления 2.0 с.
На рис. 5 представлены результаты моделирования движения группы из 4-х БЛА с лидером при интервале обновления навигационных сигналов 1 с и максимальной погрешностью навигации 6,0 м.
Coordinate у
а) траектории БЛА
10 20 30 -И 50 60 70
б) высоты полета БЛА
Рис. 5. Результаты моделирования группы БЛА с лидером
В табл. 1 представлены погрешности поддержания заданного строя группой БЛА с лидеров при различных частотах обновления навигационных сигналов и различных погрешностях навигации.
Из табл. 1 следует, что с ростом погрешности навигационной системы растет и погрешность поддержания заданного положения в строю группой БЛА с лидером. При этом наибольшая точность обеспечивается по каналу управления высотой, т.к. она поддерживается независимо друг от друга. Наибольшая погрешность наблюдается в каналах управления по координатам ОХ и 02, что связано с необходимостью отслеживать группой положения лидера. При этом погрешность поддержания положения по оси ОХ максимальна, т.к. необходимо отслеживать изменяющуюся координату. В качестве погрешности навигационной системы выбрано средне-квадратичное отклонение [46], а в качестве погрешности системы управления - усредненное по реализациям значение от максимальных ошибок в процессе полета.
Таблица 1
Зависимость погрешностей системы управления с лидером от точности и частоты обновления навигационных данных
СКО НС, м 0,5 1,0 1,5 2,0
Ошибка по оси OY Дйав=0,1 с 0,73 1,07 1,53 2,6
Дйав=0,2 с 0,7 0,8 2,13 2,77
Дйав=0,5 с 0,7 0,9 2,17 2,67
Дйав=1,0 с 0,87 1,43 1,93 3,6
Ошибка по оси OX Дйав=0,1 с 1,7 3,2 4,3 6,8
Дйав=0,2 с 1,93 3,1 3,9 7,17
Дйав=0,5 с 1,67 2,87 5,27 6,77
Дйав=1,0 с 2,07 3,0 4,97 6,57
Ошибка по оси OZ Дйав=0,1 с 0,27 1,87 2,9 4,5
Дйав=0,2 с 0,33 1,97 2,33 3,77
Дйав=0,5 с 0,47 1,8 2,83 4,57
Дйав=1,0 с 0,47 1,97 3,2 4,17
На основе данных табл. 1 можно отметить, что погрешность системы группового управления с лидером слабо зависит от интервала обновления навигационных данных до того момента, пока он не приблизится к постоянной времени замкнутой системы. После этого система управления резко теряет качество и устойчивость. Так, в моделируемом примере потеря устойчивости наблюдается при частоте обновления навигационных данных около 1,5 с, что примерно соответствует постоянной времени замкнутой системы, равной 2,0 с.
В табл. 2 представлены результаты моделирования роевой системы управления, в которой желаемые положения БЛА по осям 0Х и 02 задаются выражениями (56)-(59).
Сравнивая данные, представленные в табл. 1 и 2 можно заметить, что наиболее существенное отличие роевой системы от системы группового управления с лидером - более высокая точность позиционирования относительно оси 0Х. В частности, для СКО навигационной системы 2,0 м, роевая система управления повышает точность позиционирования в среднем в 3,5 раза. Вместе с тем, погрешность позиционирования по высоте и по оси 02 существенно не изменяется. Это связано с тем, что роевой принцип позиционирования применен именно для позиционирования по фронту (по оси 0Х).
Таблица 2
Зависимость погрешностей роевой системы управления от точности и частоты обновления навигационных данных
СКО НС, м 0,5 1,0 1,5 2,0
Ошибка по оси OY Дйав=0,1 с 0,83 1,0 2,0 2,28
Дйав=0,2 с 0,83 0,97 1,75 2,23
Дйав=0,5 с 0.7 10,3 1,8 2,8
Дйав=1,0 с 0,9 1,34 3,33 3,0
Ошибка по оси OX Дйав=0,1 с 0,45 0,6 1,33 1,15
Дйав=0,2 с 0,5 1,07 1,46 1,54
Дйав=0,5 с 0,29 1,0 1,2 2,17
Дйав=1,0 с 0,32 0,7 1,7 2,4
Ошибка по оси OZ Дйав=0,1 с 0,47 1,58 2,85 3,6
Дйав=0,2 с 0,5 1,5 2,24 3,2
Дйав=0,5 с 0,35 1,57 2,3 4,3
Дйав=1,0 с 0,4 1,57 3,2 3,92
Можно отметить, что повышение точности поддержания строя требует более развитой системы связи или наличия на борту БЛА средств, позволяющий с достаточной точностью определять положения соседних БЛА группы. В случае группы с лидером, свои координаты остальным членам группы передает только лидер. В роевом методе БЛА передают свои координаты соседям. Совместная работа систем связи БЛА требует организации специальной сети связи [47], что в некоторых средах может создавать технические проблемы, например, под водой [48].
Отметим, что в данной работе не рассматриваются режимы значительных отклонений от желаемых положений. Это связано с тем, что в таких режимах требуется решать задачи планирования движения [49-52] и устранения коллизий.
Заключение. В данной статье проведено сравнение двух подходов к управлению движением группы БЛА при движении строем - управление с лидером и децентрализованное управление. Сравнение проведено по точности поддержания строя, которая имеет важное значение при выполнении ряда практических задач. Результаты моделирования показали, что система децентрализованного управления позволяет примерно в 3,5 раза уменьшить погрешность поддержания строя по сравнению с системой управления с лидером. Вместе с тем, реализация системы децентрализованного управления требует, чтобы соседние БЛА передавали друг другу информацию о своем местонахождении. В системе с лидеров передачу своих координат группе осуществляет только лидер, а ведомые БЛА находятся в режиме молчания. Таким образом, в ситуации, когда имеется возможность организовать сеть связи, применение децентрализованного метода управления представляется более перспективным.
Работа выполнена при поддержке гранта Российского научного фонда № 22-29-00533, выполняемого в АО "Научно-конструкторское бюро робототехники и систем управления".
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Arteaga-Escamilla C.M., Castro-Linares R., Alvarez-Gallegos J. Leader-follower formation with reduction of the off-tracking and velocity estimation under visibility constraints // International Journal of Advanced Robotic Systems. - 2021. - Vol. 18 (610).
2. Sun F., Li H., Zhu W., Kurths J. Fixed-time formation tracking for multiple nonholonomic wheeled mobile robots based on distributed observer // Nonlinear Dynamics. - 2021. - Vol. 106.
- P. 3331-3349.
3. DongX., Yu B., Shi Z., Zhong Y. Time-Varying Formation Control for Unmanned Aerial Vehicles: Theories and Applications // IEEE Transactions on Control Systems Technology. - Vol. 23 (1).
- P. 340-348.
4. Pack D.J., DeLima P., Toussaint G.J., York G. Cooperative control of UAVs for localization of intermittently emitting mobile targets // IEEE Transactions on Systems, Man, and Cybernetics, Part B (Cybernetics). - 2009. - Vol. 39 (4). - P. 959-970.
5. Безрук Г.Г., Мартынова Л.А., Саенко И.Б. Динамический метод поиска антропогенных объектов в морском дне с использованием автономных необитаемых подводных аппаратов // Тр. СПИИРАН. - 2018. - Т. 3 (58). - С. 203-226.
6. Шепета А.П., Ненашев В.А. Точностные характеристики определения координат объектов в двухпозиционной системе малогабаритных бортовых РЛС // Информационно-управляющие системы. - 2020. - Т. 2. - P. 31-36.
7. Nenashev V.A., Shepeta A.P., Kryachko A.F. Fusion Radar and Optical Information in MultiPosition on-Board Location Systems // 2020 Wave Electronics and its Application in Information and Telecommunication Systems (WECONF). - 2020. - P. 1-5. - DOI: 10.1109/ WEC0NF48837.2020.9131451.
8. Nenashev V.A., Khanykov I.G. Formation of Fused Images of the Land Surface from Radar and Optical Images in Spatially Distributed On-Board Operational Monitoring Systems // Journal of. Imaging. - 2021. - Vol .7 (251). - https://doi.org/10.3390/jimaging7120251.
9. Ненашев В.А., Ханыков И.Г. Формирование комплексного изображения земной поверхности на основе кластеризации пикселей локационных снимков в многопозиционной бортовой системе // Информатика и автоматизация. - 2021. - Т. 20, № 2. - С. 302-340.
10. Lewis M.A., Tan K.-H. High Precision Formation Control of Mobile Robots Using Virtual Structures // Autonomous Robots. - 1997. - No. 4. - P. 387-403.
11. Tan K-H., Lewis M. Virtual Structures for High-Precision Cooperative Mobile Robotic Control // Proc. ofthe IEEE/RSJ Intern. Conf. on Intelligent Robots and Systems. - 1996. - P. 132-139.
12. Морозова Н.С. Виртуальные формации и виртуальные лидеры в задаче о движении строем группы роботов // Вестник Санкт-Петербургского университета. Серия 10. Прикладная математика // Информатика. Процессы управления. - 2015. - № 1. - C. 135-149.
13. Эндо Т., Маэда Р., Мацуно Ф. Анализ устойчивости роя гетерогенных роботов с ограниченным полем зрения // Информатика и автоматизация. - 2020. - Т. 19 (5). - С. 942-966.
14. Гайдук А.Р., Мартьянов О.В., Медведев М.Ю., Пшихопов ВХ., Хамдан Н., Фархуд А. Нейросетевая система управления группой роботов в неопределенной двумерной среде // Мехатроника, автоматизация, управление. - 2020. - Т. 21 (8). - С. 470-479.
15. ПшихоповВ.Х., МедведевМ.Ю. Групповое управление движением мобильных роботов в неопределенной среде с использованием неустойчивых режимов // Тр. СПИИРАН.
- 2018. - Вып. 60. - C. 39-63.
16. Medvedev M., Pshikhopov V., Gurenko B., Hamdan N. Path planning method for mobile robot with maneuver restrictions // Proc. of the International Conference on Electrical, Computer, Communications and Mechatronics Engineering (ICECCME) 7-8 October 2021, Mauritius. 10.1109/ICECCME52200.2021.9591090.
17. Carlos A., Hebertt S., Jos'e A. Stability of Active Disturbance Rejection Control for Uncertain Systems: a Lyapunov Perspective // International Journal of Robust Nonlinear Control. - 2017.
- Vol. 27. - P. 4541-4553. - DOI: 10.1002/rnc.3812.
18. Воротников В.И., Вохмянина А.В. Метод линеаризующей обратной связи в задаче управления по части переменных при неконтролируемых помехах // Тр. СПИИРАН.
- 2018. - Вып. 6 (61). - С. 61-93.
19. Фельдбаум А.А. О распределении корней характеристического уравнения системы регулирования // Автоматика и телемеханика. - 1948. - № 4. - С. 253-279.
20. Finaev V.I., Medvedev M.Yu., Pshikhopov V.K., Pereverzev V.A., Soloviev V.V. Unmanned Powerboat Motion Terminal Control in an Environment with Moving Obstacles // Mekhatronika, Avtomatizatsiya, Upravlenie. - 2021. - Vol. 22 (3). - P. 145-154.
- https://doi.org/10.17587/mau.22.145-154.
21. Park B.-S., Yoo, S.-J. Adaptive Secure Control for Leader-Follower Formation of Nonholonomic Mobile Robots in the Presence of Uncertainty and Deception Attacks // Mathematics. - 2021. - Vol. 9. - https://doi.org/10.3390/math9182190.
22. Hirata-Acosta J., Pliego-Jiménez J., Cruz-Hernádez C., Martínez-Clark R. Leader-Follower Formation Control of Wheeled Mobile Robots without Attitude Measurements // Applied Sciences. - 2021. - Vol. 11 (12), 5639. - https://doi.org/10.3390/app11125639.
23. Maghenem M., Loria A., Panteley E. Cascades-based leader-follower formation tracking and stabilization of multiple nonholonomic vehicles // IEEE Transactions on Automatic Control.
- Vol. 65 (8). - P. 3639-3646.
24. Wang Z., Wang L., Zhang H., Chen Q., Liu J. Distributed regular polygon formation control and obstacle avoidance for non-holonomic wheeled mobile robots with directed communication topology // IET Control Theory and Applications. - 2020. - Vol. 14 (9). - P. 1113-1122.
25. Sun J., Chen J. A survey on Lyapunov-based methods for stability of linear time-delay systems // Frontiers of Computer Science. - 2017. - Vol. 11. - P. 555-567. - https://doi.org/ 10.1007/s11704-016-6120-3.
26. Hu J., BhowmickP., Lanzon A. Group Coordinated Control of Networked Mobile Robots With Applications to Object Transportation // IEEE Transactions on Vehicular Technology. - 2021.
- Vol. 70 (8). - P. 8269-8274. - DOI: 10.1109/TVT.2021.3093157.
27. Arnold W.F., Laub A.J. Generalized Eigenproblem Algorithms and Software for Algebraic Riccati Equations // Proceedings of the IEEE. - 1984. - Vol. 72 (12). - P. 1746-1754.
28. Pshikhopov V., Medvedev M. Multi-Loop Adaptive Control of Mobile Objects in Solving Trajectory Tracking Tasks // Automation and Remote Control. - 2020. - Vol. 81 (11). - P. 2078-2093.
- https://doi.org/10.1134/S0005117920110090.
29. Stephen T. Thornton, Jerry B. Marion. Classical Dynamics of Particles and Systems. Brooks Cole. - 5th ed. - 2003. - 672 p. - ISBN-10: 0534408966.
30. Пшихопов ВХ., Медведев М.Ю., Гайдук А.Р., Нейдорф Р.А., Беляев В.Е., Федоренко Р.В., Костюков В.А., Крухмалев В.А. Система позиционно-траекторного управления роботизированной воздухоплавательной платформой: математическая модель // Мехатроника, автоматизация и управление. - 2013. - № 6. - С. 14-21.
31. Gotten F., Finger D.F., Havermann M. et al. Full configuration drag estimation of short-to-medium range fixed-wing UAVs and its impact on initial sizing optimization // CEAS Aeronautical Journal. - 2021. - Vol. 12. - P. 589-603. - https://doi.org/10.1007/s13272-021-00522-w.
32. Milne-Thomson L.M. Theoretical aerodynamics. - Courier Corporation, 2012.
33. Li X., Qi G., Zhang L. Time-varying formation dynamics modeling and constrained trajectory optimization of multi-quadrotor UAVs // Nonlinear Dynamics. - 2021. - Vol. 106. - P. 3265-3284.
- https://doi.org/10.1007/s11071-021-06788-3.
34. Medvedev M., Pshikhopov V. Path Planning of Mobile Robot Group Based on Neural Networks // Lecture Notes in Artificial Intelligence. - 2020. - P. 51-62. - https://doi.org/ 10.1007/978-3-030-55789-8_5.
35. Ren X.X., Yang G.H. Noise covariance estimation for networked linear systems under random access protocol scheduling // Neurocomputing. - 2021. - Vol. 455 (30). - P. 68-77.
36. Golnaraghi F., Kuo B.C. Automatic Control Systems. - 9-th ed. - JOHN WILEY & SONS, 2010. - 944 p.
37. DieboldF. Elements of Forecasting (Fourth ed.). - Thomson/South-Western, 2007. - 366 p.
38. El-Sheimy N., Hou H., Niu X. Analysis and modeling of inertial sensors using Allan variance // IEEE Transactions on Instrumentation and Measurement. - 2008. - Vol. 57 (1). - P. 140-149.
39. Consolini L., Morbidi F., Prattichizzo D., et al. Leader-follower formation control of nonholonomic mobile robots with input constraints // Automatica. - 2008. - Vol. 44 (5).
- P. 1343-1349.
40. Cocetti M., Tarbouriech S., Zaccarian L. High-Gain Dead-Zone Observers for Linear and Nonlinear Plants // IEEE Control Systems Letters. - 2019. - Vol. 3 (2). - P. 356-361. - DOI: 10.1109/LCSYS.2018.2880931.
41. Liu Y., Chen C., Wu H. et al. Structural stability analysis and optimization of the quadrotor unmanned aerial vehicles via the concept of Lyapunov exponents // International Journal of Advanced Manufacturing Technology. - 2018. - Vol. 94. - P. 3217-3227. - https://doi.org/ 10.1007/s00170-016-9311-z.
42. Omürlü V.E., Büyük^ahin U., Artar R et al. An experimental stationary quadrotor with variable DOF // Sadhana. - 2013. - Vol. 38. - P. 247-264. https://doi.org/10.1007/s12046-013-0132-6.
43. Пшихопов В.Х., Медведев М.Ю., Гуренко Б.В. Алгоритмы терминального управления подвижными объектами мультикоптерного типа // Мехатроника, автоматизация и управление. - 2019. - Т. 20, № 1. - С. 44-51. - 10.17587/mau.20.44-51.
44. Bayindir L. A review of swarm robotics tasks // Neurocomputing. - 2016. - https://doi.org/ 10.1016/ j.neucom.2015.05.116.
45. Shi Y. Particle swarm optimization: Developments, applications and re- 23 sources // Evolutionary Computing, 2001 IEEE International Conference on. - 2001. - Vol. 1. - P. 81-86.
46. Bruce P.C., Bruce A.G. Practical Statistics for Data Scientists. - 1-st ed. O'Reilly Media. USA, 2016.
47. Zhou P., FangX., Fang Y., He R., Long Y. and Huang G. Beam Management and Self-Healing for mmWave UAV Mesh Networks // IEEE Transactions on Vehicular Technology. - 2019. -Vol. 68 (2). - P. 1718-1732. - DOI: 10.1109/TVT.2018.2890152.
48. Li N., Cürüklü B., Bastos J., Sucasas V., Fernandez J.A.S., Rodriguez J. A Probabilistic and Highly Efficient Topology Control Algorithm for Underwater Cooperating AUV Networks // Sensors. - 2017. - Vol. 17, 1022. - https://doi.org/10.3390/s17051022.
49. Kostjukov V., Medvedev M., Pshikhopov V. Method for Optimizing of Mobile Robot Trajectory in Repeller Sources Field // Informatics and Automation. - 2021. - Vol. 20 (3). - P. 690-726.
50. Медведев М.Ю., Костюков В.А., Пшихопов ВХ. Оптимизация движения мобильного робота на плоскости в поле конечного числа источников-репеллеров // Тр. СПИИРАН.
- 2020. - Вып. 19 (1). - С. 43-78. - https://doi.org/10.15622/10.15622/sp.2020.19.L2.
51. Sánchez-Ibáñez J.R., Pérez-del-Pulgar C.J., García-Cerezo A. Path Planning for Autonomous Mobile Robots: A Review // Sensors. - 2021. - Vol. 21. - P. 7898. - https://doi.org/ 10.3390/s21237898.
52. Mac T.T., Copot C., Tran D.T., De Keyser R. Heuristic approaches in robot path planning: A survey // Robotics and Autonomous Systems. - 2016. - Vol. 86. - P. 13-28.
REFERENCES
1. Arteaga-Escamilla C.M., Castro-Linares R., Alvarez-Gallegos J. Leader-follower formation with reduction of the off-tracking and velocity estimation under visibility constraints, International Journal of Advanced Robotic Systems, 2021, Vol. 18 (610).
2. Sun F., Li H., Zhu W., Kurths J. Fixed-time formation tracking for multiple nonholonomic wheeled mobile robots based on distributed observer, Nonlinear Dynamics, 2021, Vol. 106, pp. 3331-3349.
3. DongX., Yu B., Shi Z., Zhong Y. Time-Varying Formation Control for Unmanned Aerial Vehicles: Theories and Applications, IEEE Transactions on Control Systems Technology, Vol. 23 (1), pp. 340-348.
4. Pack D.J., DeLima P., Toussaint G.J., York G. Cooperative control of UAVs for localization of intermittently emitting mobile targets, IEEE Transactions on Systems, Man, and Cybernetics, PartB (Cybernetics), 2009, Vol. 39 (4), pp. 959-970.
5. Bezruk G.G., Martynova L.A., Saenko I.B. Dinamicheskiy metod poiska antropogennykh ob"ektov v morskom dne s ispol'zovaniem avtonomnykh neobitaemykh podvodnykh apparatov [Dynamic Method of Searching Anthropogenic Objects in Use of Seabed with Autonomous Underwater Vehicles], Tr. SPIIRAN [SPIIRAS Proceedings], 2018, Vol. 3 (58), pp. 203-226.
6. Shepeta A.P., Nenashev V.A. Tochnostnye kharakteristiki opredeleniya koordinat ob"ektov v dvukhpozitsionnoy sisteme malogabaritnykh bortovykh RLS [Accuracy characteristics of object location in a two-position system of small onboard radars], Informatsionno-upravlyayushchie sistemy [Informatsionno-Upravliaiushchie Sistemy], 2020, Vol. 2, pp. 31-36.
7. Nenashev V.A., Shepeta A.P., Kryachko A.F. Fusion Radar and Optical Information in MultiPosition on-Board Location Systems, 2020 Wave Electronics and its Application in Information and Telecommunication Systems (WECONF), 2020, pp. 1-5. - DOI: 10.1109/ WEC0NF48837.2020.9131451.
8. Nenashev V.A., Khanykov I.G. Formation of Fused Images of the Land Surface from Radar and Optical Images in Spatially Distributed On-Board Operational Monitoring Systems, Journal of Imaging, 2021, Vol .7 (251). Available at: https://doi.org/10.3390/jimaging7120251.
9. Nenashev V.A., Khanykov I.G. Formirovanie kompleksnogo izobrazheniya zemnoy poverkhnosti na osnove klasterizatsii pikseley lokatsionnykh snimkov v mnogopozitsionnoy bortovoy sisteme [Formation of Fused Images of the Land Surface from Radar and Optical Images in Spatially Distributed On-Board Operational Monitoring Systems], Informatika i avtomatizatsiya [Informatics and Automation], 2021, Vol. 20, No. 2, pp. 302-340.
10. Lewis M.A., Tan K.-H. High Precision Formation Control of Mobile Robots Using Virtual Structures, Autonomous Robots, 1997, No. 4, pp. 387-403.
11. Tan K-H., Lewis M. Virtual Structures for High-Precision Cooperative Mobile Robotic Control, Proc. of the IEEE/RSJ Intern. Conf. on Intelligent Robots and Systems, 1996, pp. 132-139.
12. Morozova N.S. Virtual'nye formatsii i virtual'nye lidery v zadache o dvizhenii stroem gruppy robotov [Virtual formations and virtual leaders in the task of moving a group of robots in formation], Vestnik Sankt-Peterburgskogo universiteta. Seriya 10. Prikladnaya matematika // Informatika. Protsessy upravleniya [Vestnik of Saint Petersburg University Applied Mathematics. Computer Science. Control Processes], 2015, No. 1, pp. 135-149.
13. Endo T., Maeda R., Matsuno F. Analiz ustoychivosti roya geterogennykh robotov s ogranichennym polem zreniya [Stability Analysis of Swarm Heterogeneous Robots with Limited Field of View], Informatika i avtomatizatsiya [Informatics and Automation], 2020, Vol. 19 (5), pp. 942-966.
14. Gayduk A.R., Mart'yanov O.V., Medvedev M.Yu., Pshikhopov V.Kh., Khamdan N., Farkhud A. Neyrosetevaya sistema upravleniya gruppoy robotov v neopredelennoy dvumernoy srede [Neural network based control system for robots group operating in 2-d uncertain environment], Mekhatronika, avtomatizatsiya, upravlenie [Mekhatronika, Avtomatizatsiya, Upravlenie], 2020, Vol. 21 (8), pp. 470-479.
15. Pshikhopov V.Kh., Medvedev M.Yu. Gruppovoe upravlenie dvizheniem mobil'nykh robotov v neopredelennoy srede s ispol'zovaniem neustoychivykh rezhimov [Group control of autonomous robots motion in uncertain environment via unstable modes], Trudy SPIIRAN [SPIIRAS Proceedings], 2018, Issue 60, pp. 39-63.
16. Medvedev M., Pshikhopov V., Gurenko B., Hamdan N. Path planning method for mobile robot with maneuver restrictions, Proc. of the International Conference on Electrical, Computer, Communications and Mechatronics Engineering (ICECCME) 7-8 October 2021, Mauritius. 10.1109/ICECCME52200.2021.9591090.
17. Carlos A., Hebertt S., Jos'e A. Stability of Active Disturbance Rejection Control for Uncertain Systems: a Lyapunov Perspective, International Journal of Robust Nonlinear Control, 2017, Vol. 27, pp. 4541-4553. DOI: 10.1002/rnc.3812.
18. Vorotnikov V.I., Vokhmyanina A.V. Metod linearizuyushchey obratnoy svyazi v zadache upravleniya po chasti peremennykh pri nekontroliruemykh pomekhakh [Feedback Liniarization Method for Problem of Control of a Part of Variables in Uncontrolled Disturbances], Tr. SPIIRAN [SPIIRAS Proceedings], 2018, Issue 6 (61), pp. 61-93.
19. Fel'dbaum A.A. O raspredelenii korney kharakteristicheskogo uravneniya sistemy regulirovaniya [On the distribution of the roots of the characteristic equation of control systems], Avtomatika i telemekhanika [Automation and Remote Control], 1948, No. 4, pp. 253-279.
20. Finaev V.I., Medvedev M.Yu., Pshikhopov V.K., Pereverzev V.A., Soloviev V.V. Unmanned Powerboat Motion Terminal Control in an Environment with Moving Obstacles, Mekhatronika, Avtomatizatsiya, Upravlenie, 2021, Vol. 22 (3), pp. 145-154. Available at: https://doi.org/10.17587/mau.22.145-154.
21. Park B.-S., Yoo, S.-J. Adaptive Secure Control for Leader-Follower Formation of Nonholonomic Mobile Robots in the Presence of Uncertainty and Deception Attacks, Mathematics, 2021, Vol. 9. Available at: https://doi.org/10.3390/math9182190.
22. Hirata-Acosta J., Pliego-Jiménez J., Cruz-Hernâdez C., Martinez-Clark R. Leader-Follower Formation Control of Wheeled Mobile Robots without Attitude Measurements, Applied Sciences, 2021, Vol. 11 (12), 5639. - https://doi.org/10.3390/app11125639.
23. Maghenem M., Loria A., Panteley E. Cascades-based leader-follower formation tracking and stabilization of multiple nonholonomic vehicles, IEEE Transactions on Automatic Control, Vol. 65 (8), pp. 3639-3646.
24. Wang Z., Wang L., Zhang H., Chen Q., Liu J. Distributed regular polygon formation control and obstacle avoidance for non-holonomic wheeled mobile robots with directed communication topology, IET Control Theory and Applications, 2020, Vol. 14 (9), pp. 1113-1122.
25. Sun J., Chen J. A survey on Lyapunov-based methods for stability of linear time-delay systems, Frontiers of Computer Science, 2017, Vol. 11, pp. 555-567. Available at: https://doi.org/10.1007/s11704-016-6120-3.
26. Hu J., BhowmickP., Lanzon A. Group Coordinated Control of Networked Mobile Robots With Applications to Object Transportation, IEEE Transactions on Vehicular Technology, 2021, Vol. 70 (8), pp. 8269-8274. DOI: 10.1109/TVT.2021.3093157.
27. Arnold W.F., Laub A.J. Generalized Eigenproblem Algorithms and Software for Algebraic Riccati Equations, Proceedings of the IEEE, 1984, Vol. 72 (12), pp. 1746-1754.
28. Pshikhopov, V., Medvedev, M. Multi-Loop Adaptive Control of Mobile Objects in Solving Trajectory Tracking Tasks, Automation and Remote Control, 2020, Vol. 81 (11), pp. 2078-2093. - https://doi.org/10.1134/S0005117920110090.
29. Stephen T. Thornton, Jerry B. Marion. Classical Dynamics of Particles and Systems. Brooks Cole. 5th ed., 2003, 672 p. ISBN-10: 0534408966.
30. Pshikhopov V.Kh., Medvedev M.Yu., Gayduk A.R., Neydorf R.A., Belyaev V.E., Fedorenko R.V., Kostyukov V.A., Krukhmalev V.A. Sistema pozitsionno-traektornogo upravleniya robotizirovannoy vozdukhoplavatel'noy platformoy: matematicheskaya model' [Positional trajectory control system of a robotic aeronautical platform: mathematical model], Mekhatronika, avtomatizatsiya i upravlenie [Mekhatronika, Avtomatizatsiya, Upravlenie], 2013, No. 6, pp. 14-21.
31. Gotten F., Finger D.F., Havermann M. et al. Full configuration drag estimation of short-to-medium range fixed-wing UAVs and its impact on initial sizing optimization, CEAS Aeronautical Journal, 2021, Vol. 12, pp. 589-603. Available at: https://doi.org/10.1007/s13272-021-00522-w.
32. Milne-Thomson L.M. Theoretical aerodynamics. Courier Corporation, 2012.
33. LiX., Qi G., Zhang L. Time-varying formation dynamics modeling and constrained trajectory optimization of multi-quadrotor UAVs, Nonlinear Dynamics, 2021, Vol. 106, pp. 3265-3284. Available at: https://doi.org/10.1007/s11071-021-06788-3.
34. Medvedev M., Pshikhopov V. Path Planning of Mobile Robot Group Based on Neural Networks, Lecture Notes in Artificial Intelligence, 2020, pp. 51-62. Available at: https://doi.org/ 10.1007/978-3-030-55789-8 5.
35. RenX.X., Yang G.H. Noise covariance estimation for networked linear systems under random access protocol scheduling, Neurocomputing, 2021, Vol. 455 (30), pp. 68-77.
36. Golnaraghi F., Kuo B.C. Automatic Control Systems. 9-th ed. JOHN WILEY & SONS, 2010, 944 p.
37. DieboldF. Elements of Forecasting (Fourth ed.). Thomson/South-Western, 2007, 366 p.
38. El-Sheimy N., Hou H., Niu X. Analysis and modeling of inertial sensors using Allan variance, IEEE Transactions on Instrumentation and Measurement, 2008, Vol. 57 (1), pp. 140-149.
39. Consolini L., Morbidi F., Prattichizzo D., et al. Leader-follower formation control of nonholonomic mobile robots with input constraints, Automatica, 2008, Vol. 44 (5), pp. 1343-1349.
40. Cocetti M., Tarbouriech S., Zaccarian L. High-Gain Dead-Zone Observers for Linear and Nonlinear Plants, IEEE Control Systems Letters, 2019, Vol. 3 (2), pp. 356-361. DOI: 10.1109/LCSYS.2018.2880931.
41. Liu Y., Chen C., Wu H. et al. Structural stability analysis and optimization of the quadrotor unmanned aerial vehicles via the concept of Lyapunov exponents, International Journal of Advanced Manufacturing Technology, 2018, Vol. 94, pp. 3217-3227. Available at: https://doi.org/ 10.1007/s00170-016-9311-z.
42. Omürlü V.E., Büyük§ahin U., Artar R. et al. An experimental stationary quadrotor with variable DOF, Sadhana, 2013, Vol. 38, pp. 247-264. Available at: https://doi.org/10.1007/s12046-013-0132-6.
43. Pshikhopov V.Kh., Medvedev M.Yu., Gurenko B.V. Algoritmy terminal'nogo upravleniya podvizhnymi ob"ektami mul'tikopternogo tipa [Algorithms for terminal control of moving objects of multicopter type], Mekhatronika, avtomatizatsiya i upravlenie [Mekhatronika, Avtomatizatsiya, Upravlenie], 2019, Vol. 20, No. 1, pp. 44-51. 10.17587/mau.20.44-51.
44. Bayindir L. A review of swarm robotics tasks, Neurocomputing, 2016. Available at: https://doi.org/ 10.1016/j.neucom.2015.05.116.
45. Shi Y. Particle swarm optimization: Developments, applications and re- 23 sources, Evolutionary Computing, 2001 IEEE International Conference on, 2001, Vol. 1, pp. 81-86.
46. Bruce P.C., Bruce A.G. Practical Statistics for Data Scientists. 1-st ed. O'Reilly Media. USA, 2016.
47. Zhou P., FangX., Fang Y., He R., Long Y. and Huang G. Beam Management and Self-Healing for mmWave UAV Mesh Networks, IEEE Transactions on Vehicular Technology, 2019, Vol. 68 (2), pp. 1718-1732. DOI: 10.1109/TVT.2018.2890152.
48. Li N., Cürüklü B., Bastos J., Sucasas V., Fernandez J.A.S., Rodriguez J. A Probabilistic and Highly Efficient Topology Control Algorithm for Underwater Cooperating AUV Networks, Sensors, 2017, Vol. 17, 1022. Available at: https://doi.org/10.3390/s17051022.
49. Kostjukov V., MedvedevM., Pshikhopov V. Method for Optimizing of Mobile Robot Trajectory in Repeller Sources Field, Informatics and Automation, 2021, Vol. 20 (3), pp. 690-726.
50. Medvedev M.Yu., Kostyukov V.A., Pshikhopov V.Kh. Optimizatsiya dvizheniya mobil'nogo robota na ploskosti v pole konechnogo chisla istochnikov-repellerov [Optimization of mobile robot movement on a plane with finite number of repeller sources], Tr. SPIIRAN [SPIIRAS Proceedings], 2020, Issue 19 (1), pp. 43-78. Available at: https://doi.org/10.15622/ 10.15622/sp.2020.19.1.2.
51. Sánchez-Ibáñez J.R., Pérez-del-Pulgar C.J., García-Cerezo A. Path Planning for Autonomous Mobile Robots: A Review, Sensors, 2021, Vol. 21, pp. 7898. Available at: https://doi.org/ 10.3390/s21237898.
52. Mac T.T., Copot C., Tran D.T., De Keyser R. Heuristic approaches in robot path planning: A survey, Robotics and Autonomous Systems, 2016, Vol. 86, pp. 13-28.
Статью рекомендовал к опубликованию д.т.н., профессор С.Г. Капустян.
Медведев Михаил Юрьевич - НИИ робототехники и процессов управления Южного федерального университета; e-mail: medvmihal@sfedu.ru; г. Таганрог, Россия; тел.: 88634371694;
д.т.н.; в.н.с.
Пшихопов Вячеслав Хасанович - e-mail: pshichop@rambler.ru; д.т.н.; профессор; директор.
Medvedev Mikhail Yurjevich - R&D Institute of Robotics and Control Systems; e-mail:
medvmihal@sfedu.ru; Taganrog, Russia; phone: 88634371694; dr. of eng. sc.; leading researcher.
Pshikhopov Viacheslav Khasanovich - e-mail: pshichop@rambler.ru; dr. of eng. sc.; professor;
director.