Научный журнал
Научное обозрение. Технические науки
ISSN 2500-0799
ПИ №ФС77-57440

МОДЕЛИРОВАНИЕ ТРАЕКТОРИИ С ОГРАНИЧЕНИЯМИ НА КРИВИЗНУ

Нефедова В.А. 1 Дубанов А.А. 1
1 ФГБОУ ВО «Бурятский государственный университет имени Доржи Банзарова»
Предметом исследования в нашей статье является моделирование и расчет оптимальной траектории между преследующим и преследуемым объектом. Предлагается метод имитационного моделирования, позволяющий решать практические исследовательские задачи. В статье описаны этапы моделирования траектории движущегося тела, имеющей ограничение по кривизне, следует учесть то, что тело, движущееся по заданной траектории, начинает движение под начальным заданным углом и завершает его также под определенным углом. Наша задача имеет прикладное значение, так как рассматривается подход к организации маршрутной траектории полета летательного аппарата в вертикальной плоскости при наличии ограничений на переменные состояния. Основные результаты работы: введены математические модели движущихся тел, приведены расчеты для моделирования их траектории в каждый момент времени движения. Разработан программный код для математического моделирования траектории согласно поставленной задачи на основе аппарата аналитической геометрии, построение моделируемой траектории осуществлено в некотором классе функций. На основе этого описывается новый алгоритм планирования траектории. Данный алгоритм решения задачи является актуальным, так как результаты работы могут найти свое применение при решении прикладных задач робототехники, авиамоделирования, ракетного моделирования и др. Решение задачи в виде программного кода было произведено в среде программирования MathCAD.
моделируемая траектория
имитационное моделирование
траектория летательных аппаратов
ограничение по кривизне
преследующий объект
касательная к траектории
1. Айзекс Р. Дифференциальные игры. М.: Мир, 1967. 480 с.
2. Бурдаков С.В., Сизов П.А. Алгоритмы управлением движения мобильным роботом в задаче преследования // Научно-технические ведомости Санкт-Петербургского государственного политехнического университета. Информатика. Телекоммуникации. Управление. 2014. № 6 (210). С. 49–58.
3. Дубанов А.А. Конструирование базовых элементов формы и параметризация составных поверхностей // Информационные технологии. 2003. № 10. С. 42–51.
4. Ахметжанов А.Р. Динамические игры преследования на поверхностях: автореф. дис. ... канд. физ.-мат.наук. Москва, 2019. С. 28.
5. Соловьева Н.А. Групповое преследование в рекуррентных диффенциальных играх: 01.01.02-дифференциальные уравнения: дис. ... канд. физ.-мат.наук. Ижевск, 2016. С. 97.
6. Изместьев И.В., Ухоботов В.И. Задача преследования маломаневренных объектов с терминальным множеством в форме кольца // Итоги науки и техн. сер. Соврем. Мат. и ее прил. Темат. Обз. 2018. Т. 148. С. 25–31.
7. Кузьмина Л.И., Осипов Ю.В. Расчет длины траектории для задачи преследования // Вестник МГСУ. 2013. С. 20–26.
8. Саматов Б.Т. Задача преследования убегания при интегрально-геометрических ограничениях на управления преследователя // Автомат. и телемех. 2013. № 7. С. 17–28.
9. Романников Д.О. Пример решения минимаксной задачи преследования с использованием нейронных сетей // Сборник научных трудов НГТУ. 2018. № 2 (92). С. 108–116.
10. Пашко С.В. Эффективные стратегии преследования. Основанные на использовании функции Ляпунова // Доклады Национальной академии наук Украины. 2016. № 1. С. 26–33.
11. Пашко С.В. NP-трудность задач оптимизации коллективного преследования // Материалы 8-й Международной конференции по программированию UkrPROG’2014, Украина, Киев, 2014. С. 44–52.
12. Программный код в системе MathCAD представлен. [Электронный ресурс]. URL: http://dubanov.exponenta. ru/books.htm (дата обращение: 05.12.2019).
13. Анимированное изображение по результату работы программного кода. [Электронный ресурс]. URL: https://www.youtube.com/watch?v=SNKeL3pKLwo&feature=youtu.be (дата обращение: 05.12.2019).

В настоящей статье мы рассмотрим вопросы моделирования траекторий между преследующим и преследуемым объектами. На плоскости (или на поверхности) передвигаются два объекта, так как мы интерпретируем нашу задачу в задачу преследования, это участники задачи преследования. Объекты, преследующий и преследуемый, будут двигаться по определенной траектории преследования. Моделируемая траектория должна удовлетворять некоторым условиям. Одно из условий – это то, что минимальный радиус кривизны траектории не должен быть меньше некоторого порогового значения. Еще одно условие в рассматриваемой задаче – это то, что моделируемая траектория выходит под определенным углом из точки нахождения преследующего объекта и приходит в точку нахождения преследуемого объекта также под заданным углом. Будем рассматривать моделирование траектории в каждый момент времени, что предполагает точность вычислений. Теоретические положения для данной работы изложены в источниках [1, 2]. Некоторые особенности моделирования составных кривых взяты из источника [3].

Цель исследования: моделирование оптимальной траектории движения в каждый момент времени, при наличии ограничений на кривизну. Также, при заданных углах выхода преследующего объекта из точки нахождения и прихода в конечную точку преследуемого объекта под заданным углом. И описание нового алгоритма планирования заданной траектории в виде програм- много кода для применения нашего решения на практических задачах в областях робототехники, авиамоделирования, моделирования робототехнических комплексов и других моделях, где присутствуют элементы задачи преследования.

Материалы и методы исследования

Рассмотрим случай на горизонтальной плоскости проекций. Пусть преследуемый объект движется по траектории nefed01.wmf, а преследующий объект движется по траектории nefed02.wmf (рис. 1).

NEF1.tif

Рис. 1. Моделируемая траектория

В нашей задаче, которая решается в настоящей статье, моделирование траектории в каждый момент времени t. Моделируемая траектория может выходить из точки nefed03.wmf под любым углом. Но в программном коде, написанном по материалам статьи, мы приняли угол между касательной к моделируемой траектории и касательной к траектории nefed04.wmfравным 0, что не имеет принципиального значения в плане моделирования алгоритма. Также угол между касательной к моделируемой траектории и касательной к траектории nefed05.wmf в листинге программы принят равным nefed06.wmf. При этом моделируемая траектория имеет ограничения по своей кривизне. То есть радиус кривизны моделируемой кривой не может быть меньше некоторой пороговой величины. Реализация тестовых программ будет сделана в системе компьютерной математики MathCAD. Поскольку мы имеем дело с динамическими объектами, результат будет снабжен ссылкой на анимированные изображения.

При моделировании траекторий в задачах преследования, приняты во внимание работы, посвященные динамическим играм преследования [4], групповому преследованию [5], простому преследованию [6], расчету длин траекторий в задачах преследования [7].

Также использовались результаты, достигнутые в работе по играм преследования с интегрально-геометрическими ограничениями на функции управления преследователя [8].

Как возможность применения результатов, полученных в данной статье, рассматривалась работа моделирования задачи преследования в нейронных сетях [9].

Также приняты во внимание при написании статьи методы коллективного преследования, изложенных в работах [10, 11].

Результаты исследования и их обсуждение

Моделирование окружности у преследующего объекта

Мы уже говорили, что преследующий объект движется по траектории nefed07.wmf (при моделировании общего алгоритма угол между траекторией и ее касательной не имеет принципиального значения). Выполнение условия того, что угол между касательными к моделируемой траектории и к траектории nefed08.wmf равен 0, и условие ограничения по кривизне, приводят к построению окружности, изображенной на рис. 2.

NEF2a.tif

Рис. 2. Окружность к траектории преследующего объекта

На рис. 2 показано, что вдоль траектории преследующего объекта nefed09.wmf формируется локальный динамический базис (vf, nf), где vf – единичный вектор, сонаправленный с вектором скорости nefed10.wmf преследующего объекта. Вектор nf – есть вектор нормали к касательной к траектории nefed11.wmf.

Отсюда следует, центр моделируемой окружности Cf(t) формируется следующим образом: Cf(t) = Rf(t) + r?nf(t), где r – это минимальный пороговый радиус траектории преследующего объекта.

Центр окружности может формироваться иначе (рис. 3): nefed12.wmf. В таком случае (как на рис. 3) моделируемая траектория выходит под углом ∞ к траектории преследующего объекта nefed13.wmf.

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

NEF3a.tif

Рис. 3. Формирование окружности на траектории преследующего объекта

Моделирование окружности у преследуемого объекта

Рассмотрим формирование окружности при входе моделируемой траектории в траекторию преследуемого объекта (рис. 4). Здесь центр окружности формируется так: nefed14.wmf, где nefed15.wmf – это точка траектории преследуемого объекта, а r – минимальное значение радиуса моделируемой траектории. Здесь направление вектора nr(t) полностью определяется направлением вектора Vm(t) скорости моделируемой траектории в момент пересечения с траекторией преследуемого объекта nefed16.wmf. Векторы nr(t) и Vm(t) должны быть перпендикулярны. В листинге программы по материалам статьи используется перпендикулярность векторов Vm(t) и vr.

NEF4a.tif

Рис. 4. Моделирование окружности у преследуемого объекта

NEF5a.tif

Рис. 5. Линия сопряжения в моделируемой траектории

NEF6a.tif

Рис. 6. Сопряжение двух окружностей

Линии сопряжения между двумя окружностями

Рассмотрим вопрос о сопряжении между двумя окружностями одного радиуса передвигающимися по траекториям nefed17.wmf и nefed18.wmf (рис. 5).

Исходя из данных, представленных на рис. 5, к построению общей касательной между окружностями (Cf, r) и (Cr, r). Если мы перейдем в систему координат, образованную точкой Cf(t) и векторами (рис. 6):

nefed19.wmf, nefed20.wmf,

то задача построения линий сопряжения между двумя окружностями преобразуется к виду, изображенному на рис. 6, где nefed21.wmf. Как видно на рис. 6, существует четыре варианта линии сопряжения. Выбор варианта сопряжения оставим разработчику. Кроме того, существуют варианты расположения окружностей, которые изображены на рис. 7. В листинге программы, написанной по материалам статьи, мы использовали один вариант взаимного расположения окружностей и построили для демонстрации четыре линии сопряжения.

Для точек nefed22.wmf, nefed23.wmf, nefed24.wmf на рис. 6 существуют решения:

nefed25.wmf, nefed26.wmf,

nefed27.wmf,

nefed28.wmf,

nefed29.wmf, nefed30.wmf,

nefed31.wmf,

nefed32.wmf.

Но эти решения найдены для локальной системы координат (Cf, E1, E1) (рис. 5), чтобы перевести точку nefed33.wmf, принадлежащую отрезку nefed34.wmf, в мировую систему координат (О, Н1, Н2), необходимо выполнить преобразование по формуле

nefed35.wmf,

где h1, h2 – это векторы nefed36.wmf и nefed37.wmf, разложенные по векторам базиса (Cf, E1, E1), то есть путем преобразования уже в мировую систему координат nefed38.wmf и nefed39.wmf.

Заключение

Результаты моделирования траекторий

В настоящей статье предлагаются модели и алгоритмы моделирования траектории путем решения задачи преследования. По результатам, изложенным в данной статье, был написан программный код в системе компьютерной математики MathCAD. С текстом программы можно ознакомиться на следующем ресурсе [12]. Кроме того, по результатам работы программ изготовлено анимированное изображение, которое можно просмотреть в следующем источнике [13].

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

Предполагается, что данные модели и алгоритм решения поставленной задачи являются актуальными так как результаты работы могут быть использованы при решении прикладных задач робототехники, авиамоделирования, моделирования робототехнических комплексов и других моделей, где присутствуют элементы задачи преследования.


Библиографическая ссылка

Нефедова В.А., Дубанов А.А. МОДЕЛИРОВАНИЕ ТРАЕКТОРИИ С ОГРАНИЧЕНИЯМИ НА КРИВИЗНУ // Научное обозрение. Технические науки. – 2019. – № 6. – С. 38-43;
URL: https://science-engineering.ru/ru/article/view?id=1268 (дата обращения: 24.04.2024).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1,674