close

Вход

Забыли?

вход по аккаунту

?

Управление беспилотным самолетом методом динамического программирования в теории графов.

код для вставкиСкачать
Известия ТулГУ. Технические науки. 2016. Вып. 11. Ч. 2
УДК 536.21:27.35:25
УПРАВЛЕНИЕ БЕСПИЛОТНЫМ САМОЛЕТОМ МЕТОДОМ
ДИНАМИЧЕСКОГО ПРОГРАММИРОВАНИЯ В ТЕОРИИ ГРАФОВ
Е.А. Пегачкова, Е.Л. Кузнецова, Ю.А. Горбунова
Рассматривается задача управления беспилотным самолетом, которому необходимо выполнить заданный маневр с минимальными затратами топлива. Используя математическую модель полета, описывающую движение самолета в канале тангажа, был составлен орграф, в котором в качестве длины дуги выступают затраты
топлива на выполнение одного из двух маневров: увеличение скорости при постоянной
высоте или набор высоты при фиксированной скорости. К рассматриваемой задаче
применяется обобщение классического принципа оптимальности и метода динамического программирования в терминах теории графов, то есть минимальный расход топлива является кратчайшим путем на ориентированном графе.
Ключевые слова: оптимальное управление, ориентированный граф, динамическое программирование.
Рассматриваемая задача полета летательного аппарата в вертикальной плоскости с минимальным расходом топлива носит прикладной характер. Предлагается упрощенный способ решения поставленной задачи,
дающий приемлемый результат с минимальной погрешностью.Составлена
упрощенная модель полета самолета в вертикальной плоскости, а также
предполагается, что для набора требуемых значений по высоте и скорости,
беспилотный летательный аппарат совершает ряд последовательных шагов, причем на каждом этапе происходит только увеличение высоты или
только увеличение скорости. Решение задачи рассматривается в классе непрерывно-дискретных систем. Используя математическую модель полета
строится орграф, в котором в качестве длины дуги выступают затраты топлива на выполнение одного из двух маневров: увеличение скорости при
постоянной высоте или набор высоты при фиксированной скорости. Вершины графа определяются парой координат, равных количеству шагов от
вершины, отвечающей начальным значениям, по осям скорости и высоты
до рассматриваемой вершины. Обобщением сформулированной задачи является одновременное увеличение скорости беспилотного самолета и набора им высоты, что приближает задачу к реальному полету, тогда перемещение по графу происходит по восходящей диагонали графа. К рассматриваемой задаче применяется обобщение классического принципа оптимальности и метода динамического программирования в терминах теории
графов, то есть минимальный расход топлива является кратчайшим путем
на ориентированном графе.Полученную на орграфе серию наборов высоты
и увеличения скорости далее необходимо подставить в соответствующие
системы дифференциальных уравнений, перейдя при этом к непрерывной
системе. Моделирование полета производится с помощью программного
352
Информатика, вычислительная техника, обработка и защита информации
комплекса, позволяющего строить траектории движения самолета и вычислять расход топлива на выполнение заданного маневра для различных
моделей и параметров беспилотных самолетов.
Постановка задачи. Рассматривается задача поиска оптимального
управления беспилотным самолетом (Aerosonde). Необходимо осуществить набор скорости и высоты полета до требуемых значений с минимальными затратами топлива. Задача рассматривается в классе непрерывнодискретных систем. Ограничимся движениями самолета с приращением
высоты на постоянной скорости и ускорениями на постоянной высоте.
Движение летательного аппарата описывается системой дифференциальных уравнений:
⋅ !
⋅⋅
− $ %&' (;
=
"
0⋅
!
⋅4
*
5
= 57,3 /12,3
−
67% (8 ;
"
@
(1)
9
=
:
%&'
(;
;
= : 67% (;
"
= −<сек ;
= = ( + ?,
где:– скорость; B– высота; C – время; ?– угол атаки; =– угол тангажа;
(– траекторный угол; сDE – коэффициент подъемной силы при = = 0; сD4 –
коэффициент производной первого порядка по углу атаки; с;E – коэффициент минимального лобового сопротивления; G– тяга двигателя; <сек – секундный расход топлива при максимальной тяге; Hℎ – плотность воздуха;
J– плотность крыла; K– размах крыла; L– коэффициент Освальда.
Приведем некоторые технические характеристикиAerosonde: размах
крыла – 2,8956 м; площадь крыла – 0,55 м2; коэффициент минимального
лобового сопротивления – 0,0434; коэффициент подъемной силы при нулевом угле тангажа – 0,23; коэффициент подъемной силы производной
первого порядка по углу атаки – 5,6106; масса самолета при полных баках
– 13,5 кг; масса самолета при пустых баках – 8,5 кг.
Граничные условия:
:CE = 20м/с; BCE = 50м; :O = 50м/с; BO = 1000м.
Требуется найти оптимальное управление, позволяющее выполнить
заданный маневр с минимальным расходом топлива:
V
RG = S
Gтек
<
UC → <&'.
G"T; сек "T;
353
Известия ТулГУ. Технические науки. 2016. Вып. 11. Ч. 2
Алгоритм решения задачи. Будем рассматривать задачу в классе
дискретных систем, в плоскости обобщенных координат высота – скоV ;и высоты
рость. Задавая интервалы изменения скорости ∆: =
Y
∆B =
, разбиваем полученные интервалы на [ и \ шагов, соответственно. На рис. 1 представлена схема такого разбиения.
9V 9
Z
Рис. 1. Схема разбиения
Далее строим взвешенный орграф], Г , каждая вершина которого
однозначно определяется парой координат, равных числу шагов от вершины, отвечающей начальным значениям, по осям скорости : и высоты B к
рассматриваемой вершине. Вес дуги – затраты топлива на выполнение самолетом элементарных операций по увеличению скорости на постоянной
высоте и по набору высоты при постоянной скорости. Расход топлива определяется из решения двух отельных систем, учитывая допущение о том,
что расход топлива считается от полного бака, т.к. изначально нет информации о количестве топлива при попадании в заданную точку на графе.Рассмотрим отдельно каждый маневр.
Набор высоты (вычисление весов вертикальных дуг графа): самолет находится одинаково ориентированным в пространстве и скорость и
угол тангажа постоянны, т.е. : = 67'%C,
= 0, = = 67'%C. Поэтому сис
тема (1) преобразуется к виду:
G = cx0 +
_
*
*
=
^ _
* ` a⋅b⋅c J ⋅ Hℎ
^ _
* `⋅
9
"
= : %&' ( ;
"
= −<сек ,
354
d
+ <$ %&' ( ;
!
5
− 67% ( ; @
(2)
где в качестве состояния выступает четверка =, (, B, < , а управление
производится тягой G.
Набор скорости (вычисление весов горизонтальных дуг графа): самолет летит на постоянной высоте с неизменным траекторным углом, т.е.
B = 67'%C, ( = 0, а значит система (1) принимает вид:
Информатика, вычислительная техника, обработка и защита информации
=
?=
!
⋅⋅ ⋅
fG − f6;E +
"
e
!
"⋅5
⋅⋅ ⋅
⋅4 ⋅
;
g ⋅ J ⋅ Hℎ ⋅
g ;@
(3)
;
= :;
"
= −<сек ,
где вектором состояния является :, h, < , а управление производится углом атаки ?.
Таким образом, получили взвешенный орграф специального вида,
правый верхний угол которого представлен на рис. 2. В соответствии с методом динамического программирования определим функциюi:, B , аналогичную функции Беллмана, как длину кратчайшего пути (минимальный
расход топлива) от предпоследней вершины до конечной вершины, при
этом ij :Y , BZ = 0.
a⋅b⋅c d
Рис. 2. Взвешенный орграф специального вида
Выпишем несколько начальных шагов решения. Если из данной
вершины путь единственный, то сразу вычисляем функцию i:, B как
сумму веса дуги при переходе от предыдущей вершины к следующей и ранее вычисленного расхода топлива:
ij
e :Y
e , BZ = U:Y
e , BZ → :Y , BZ + ij :Y , BZ ;
ij
e :Y , BZ
e = U:Y , BZ
e → :Y , BZ + ij :Y , BZ ;
ij
d :Y
d , BZ = U:Y
d , BZ → :Y , BZ
e + ij
e :Y , BZ
e ;
ij
d :Y , BZ
d = U:Y , BZ
d → :Y , BZ
e + ij
e :Y , BZ
e .
Если существует несколько путей, то выбирается тот, при котором
суммарный расход топлива будет минимальным:
355
Известия ТулГУ. Технические науки. 2016. Вып. 11. Ч. 2
ij
d :Y
e , BZ
e = <&'{U:Y
e , BZ
e → :Y , BZ
e + ij
e :Y , BZ
e ;
U:Y
e , BZ
e → :Y
e , BZ + ij
e :Y
e , BZ }.
Аналогичным образом последовательно вычисляем значения функции во всех вершинах орграфа, при этом расчеты ведутся, переходя от указанной выше «диагонали» к следующим нижестоящим. Последней вершиной, в которой будет вычислено значение функции, является начальная
вершина :E CE ; BE CE .Оптимальное управление строится, начиная с
вершины :E CE ; BE CE , с использованием условно-оптимальных траекторий, полученных на предыдущем шаге.
Полученную на орграфе серию наборов высоты и увеличения скорости подставляем в соответствующие системы дифференциальных уравнений, перейдя при этом к непрерывной системе. Построив траектории
движения самолета, вычисляем расход топлива на выполнение заданного
маневра.
Численное решение. Рассмотрим решение задачи с помощью программного комплекса, позволяющего задавать параметры летательного аппарата, граничные значения скорости и высоты, а также выбирать количество интервалов разбиения. Для наглядности приведем разбиение плоскости обобщенных координат 4х4.Числа между вершинами графа указывают
затраты топлива в килограммах на выполнение операций по увеличению
скорости или набору высоты. При этом считается, что при выполнении заданного маневра самолету не разрешается снижать скорость и высоту полета, т.е. изменение параметров полета должно проводиться на графе по
стрелкам, направленным вправо при ускорении или вверх при подъеме. На
рис. 3 красным цветом выделены условно-оптимальные управления, а синим – представлен оптимальный путь, вычисленный по приведенному выше алгоритму при следующих параметрах: угол тангажа = = 10°, начальная скорость :CE = 20м/с, конечная скорость :CE = 20м/с, текущая
высота BCE = 50м, конечная высотаBO = 100м, максимальное значение тяги двигателя G"T; = 70Н, секундный расход топлива при максимальной тяге <сек = 0,0002кг/с.
Расход топлива при таком режиме составляет 0,1301 кг. После построения полученной траектории в классе непрерывных систем расход в
силу введенного допущения незначительно сокращается до 0,1299 кг. Ниже на рис. 4 – 7 приведены графики полученных траекторий.
На рис. 6 – 7 “скачки” отражают, что переключение с режима набора высоты на набор скорости не происходят мгновенно и естественно такая
схема полета является субоптимальной и не отражает реальный полет, но
полученный таким образом расход топлива является верхней границей, т.е.
минимальные затраты топлива на выполнение требуемого маневра должны
получиться заведомо меньшими значения, полученного при решении в
классе дискретных систем.
356
Информатика, вычислительная техника, обработка и защита информации
Рис 3. Интерфейс программы
Рис.
Рис. 4. Зависимость скорости от времени
Рис. 5.. Зависимость высоты от времени
357
Известия ТулГУ.
ТулГУ Технические науки. 2016. Вып. 11. Ч. 2
Рис. 6. Зависимость траекторного угла от времени
Рис. 7. Зависимость угла атаки от времени
Заключение. Рассмотренная задача управления полетом беспилотного самолета в вертикальной плоскости до достижения заданных конечных условий с учетом минимизации расхода топлива решается в классе
непрерывно-дискретных
дискретных систем,
систем с учетом определенных допущений обоснован переход к дискретной системе. Полет самолета в работе разбивается
на отдельные маневры: увеличение скорости или набор высоты полета до
требуемых значений, расход топлива выступает в качестве длины дуги
орграфа и требуется найти кратчайший путь.Для решения задачи была разработана программа
программа, позволяющая рассчитывать оптимальный расход топлива как для дискретной
дискретной, так и для непрерывной системы.
Работа выполнена при финансовой поддержке гранта РФФИ16-3800079 мол_а, гранта Президента Российской Федерации МД
МД-4560.2015.8.
Список литературы
1. Горбунова Ю.А.,
Ю.А Пегачкова Е.А. Задача управления беспилотным
самолётом методом динамического программирования в теории графов //
XLII Международная молодёжная научная конференция «Гагаринские
«
358
Информатика, вычислительная техника, обработка и защита информации
чтения – 2016»: сборник тезисов докладов в 4 т. М.: Московский авиационный институт (национальный исследовательский университет), 2016.
Т. 1. С. 467-468.
2. Вентцель Е.С. Элементы динамического программирования. М.:
Наука, 1964. 175 с.
3. Беллман Р., Дрейфус С. Прикладные задачи динамического программирования. М.: Наука, 1965. 458 с.
4. Лебедев А.А., Чернобровкин Л.С. Динамика полета беспилотных
летательных аппаратов. М.: Оборонгиз, 1973. 615 с.
Пегачкова Елена Александровна, канд. физ.-мат. наук, доц., pegachkova@mail.ru, Россия, Москва, Московский авиационный институт (Национальный исследовательский университет),
Кузнецова Екатерина Львовна, д-р физ.-мат. наук, проф., lareyna@mail.ru,
Россия, Москва, Московский авиационный институт (Национальный исследовательский университет),
Горбунова Юлия Антоновна, студент, julia_gorbunova@list.ru, Россия, Москва, Московский авиационный институт (Национальный исследовательский университет)
CONTROL OF THE PILOTLESS AIRCRAFT BY METHOD OF DYNAMIC
PROGRAMMING IN THE THEORY OF COUNTS
E.A. Pegachkova, E.L. Kuznetsova, Yu.A. Gorbunov
The task of control of the pilotless aircraft which needs to execute the set maneuver
with the minimum costs of fuel is considered. Using the mathematical flight model describing
movement of the plane in the channel of pitch the orgraf in which fuel costs for accomplishment of one of two maneuvers act as length of an arch was constituted: increase in speed with
a fixed height or ascent at the fixed speed. Generalization of the classical principle of an optimality and method of dynamic programming in terms of the theory of counts is applied to the
considered task, that is the minimum consumption of fuel is in the shortest way on the
oriented count.
Key words: optimum control, the oriented count, dynamic programming.
Pegachkova Elena Aleksandrovna, candidate of physical and mathematical sciences,
docent, pegachkova@mail.ru, Russia, Moscow, Moscow Aviation Institute (National research
university),
Kuznetsova Ekaterina Lvovna, doctor of physical and mathematical sciences, professor, lareyna@mail.ru, Russia, Moscow, Moscow Aviation Institute (National research university),
Gorbunov Yulija Antonovna, student, julia_gorbunova@list.ru, Russia, Moscow,
Moscow Aviation Institute (National research university)
359
Документ
Категория
Без категории
Просмотров
16
Размер файла
371 Кб
Теги
методов, беспилотник, управления, графов, самолетов, программирование, теория, динамическое
1/--страниц
Пожаловаться на содержимое документа