автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.18, диссертация на тему:Аналитические и численные процедуры построения решений некоторых задач управления

кандидата физико-математических наук
Лебедев, Павел Дмитриевич
город
Екатеринбург
год
2009
специальность ВАК РФ
05.13.18
цена
450 рублей
Диссертация по информатике, вычислительной технике и управлению на тему «Аналитические и численные процедуры построения решений некоторых задач управления»

Автореферат диссертации по теме "Аналитические и численные процедуры построения решений некоторых задач управления"

□ □347

На правах рукописи УДК 517.977.58

ЛЕБЕДЕВ Павел Дмитриевич

Аналитические и численные процедуры построения решений некоторых задач управления

05.13.18 — математическое моделирование, численные методы и комплексы программ

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата физико-математических наук

1 5 ОПТ 20ПР

Екатеринбург - 2009

003479942

Работа выполнена в Институте математики и механики Уральского отделения Российской академии наук в отделе динамических систем.

Научный руководитель — кандидат физико-математических наук

Успенский Александр Александрович.

Официальные оппоненты

Ведущая организация

— доктор физико-математических наук, профессор

Тонков Евгений Леонидович

— кандидат физико-математических наук Дарьин Александр Николаевич

— ГОУ ВПО «Челябинский государственный университет», г. Челябинск.

Защита диссертации состоится "ZB"eOc,}nJ<SpjL 2009 г. в 15° часов на заседании диссертационного совета Д 212.286.10 по защите докторских и кандидатских диссертаций при ГОУ ВПО «Уральский государственный университет им. A.M. Горького» по адресу: 620000, г.Екатеринбург, пр.Ленина, 51, комн. 248.

С диссертацией можно ознакомиться в научной библиотеке ГОУ ВПО «Уральский государственный университет им. A.M. Горького»

Автореферат разослан "-2.S" ОЦЩЯ&^Х 2009 года

Ученый секретарь диссертационного совета доктор физико - математических наук,

профессор — В.Г. Пименов

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ1

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

Современный облик теории управления движением динамической системы сформировался в значительной степени под влиянием работ отечественных математиков JI.C. Понтрягина и H.H. Красовского, ставших основателями известных научных школ но теории управления. Весомые, основополагающие результаты были получены их коллегами и учениками — представителями московской научной школы Е.Ф. Мищенко, Р.В. Гамкрелидзе, В.Г. Болтянским, представителями уральской научной школы A.B. Куржанским, Ю.С. Осиповым, А.И. Субботиным, а также зарубежными учеными Р. Калманом, Р. Беллманом, Р. Айзексом, У. Флемингом, Ж-П. Обеном. Существенный прогресс в становлении и развитии теории управления связан также с именами Э.Г. Альбрехта, В.Д. Ба-тухтина, Р.Ф. Габасова, А.Я. Дубовицкого, С.Г. Завалищина, Ф.М. Кирилловой, A.A. Меликяна, A.A. Милютина, М.С. Никольского, A.A. Пет-росяна, Б.Н. Пшеничного, H.H. Субботиной, В.М. Тихомирова, E.JI. Тон-кова, В.Е. Третьякова, В.Н. Ушакова, А.Г. Ченцова, Ф.Л. Черноусько, A.A. Чикрия и многих других. Актуальность изучения управляемых систем обусловлена наличием многочисленных приложений в различных отраслях знания — в механике, робототехнике, оптике, экономике, биологии. Немаловажным побудительным мотивом в исследовании являются внутренние потребности, возникшие в математической теории управления динамическими системами в условиях конфликта и неопределенности, а также стремление исследователей привлечь для изучения динами-

1 Работа выполнена при финансовой поддержке гранта Президента РФ по поддержке ведущих научных школ № НШ-2640.2008.1, гранта РФФИ № 08-01-00587_а, федеральной программы Президиума РАН № 29

ческих задач конструкции из других разделов математики — негладкого анализа [8], дифференциальной геометрии [18], теории особенностей дифференцируемых отображений [2].

Осуществляемые в работе исследования проводятся в рамках концепции теории позиционных дифференциальных игр, развиваемой в научной школе H.H. Краеовского [10, 11]. Теория позиционных дифференциальных игр, обогащенная результатами его соратников, учеников и последователей, объединяет конструктивные методы решения широкого круга проблем от теорем существования и единственности решения до разработки и реализации вычислительных алгоритмов [24, 7,16]. Также используются конструкции отпочковавшейся от этой теории и получившей глубокое развитие в работах А.И. Субботина теории минимаксного решения уравнения в частных производных первого порядка [21].

Результатам последнего времени, связанным с проблемой построения (аналитического или приближенного) обобщенных решений уравнений типа Гамильтона-Якоби [21], предшествовали работы С.К. Годунова, Е. Hopf, P.D. Lax, других авторов. В 70-х годах прошлого века С.Н. Кружков [12], реализуя метод исчезающей вязкости, ввел определение фундаментального (обобщенного) решения уравнения в частных производных первого порядка типа эйконала посредством предельного перехода по параметру малости при старшей производной от решений соответствующих уравнений в частных производных второго порядка. Позднее, в 80-е годы, действуя аналогичным образом, М. Крэндалл и П.Л. Лионе ввели определение обобщенного решения уравнения в частных производных первого порядка, названного ими вязкостным решением. Несмотря на разницу минимаксного и вязкостного подходов, определяемый объект является одним и тем же. Отличительной чертой минимаксного подхода при изучении и численном построении обобщенных решений является активное вовлечение методов, средств и конструкций выпуклого и негладкого анализа. При этом используются наработки отечественных и зарубежных математиков [8]. Полезными с точки зрения разработки конструктивных подходов к построению обобщенных решений уравнений в частных производных первого порядка оказались мсто-

ды теории особенностей дифференцируемых отображений, разрабатываемые В.И. Арнольдом [2, 3], его коллегами, зарубежными авторами Т. Поетоном, И. Стюартом, Дж. Брусом, П. Джиблиным [5]. Средствами этой теории, в частности, формируются списки типичных особенностей каустик и волновых фронтов, предлагаются подходы к построению дискриминантных множеств.

Минимаксный подход вкупе с вязкостным подходом применяется для исследования обобщенных решений функциональных уравнений типа Гамильтона-Якоби, возникающих в задачах конфликтного управления динамическими системами, описываемыми дифференциальными уравнениями с дискретными и\или распределенными параметрами [13]. Исследуются вопросы существования и единственности обобщенных решений, выявляются условия оптимальности гарантированного результата управления в таких задачах в случае исходных данных, удовлетворяющих условию Липшица [7].

Задачами быстродействия занимались многие исследователи. Привлечение идей,результатов и конструкций указанных выше теорий позволяет в рамках минимаксного подхода разрабатывать аналитические и аппроксимационные процедуры построения функции оптимального результата для задач динамического управления по быстродействию [9]. В диссертации эти подходы распространены на задачи быстродействия для частного случая вектограммы скоростей, а также задачи геометрической оптики [20]. Разработанные процедуры позволяют проводить численно-аналитическое конструирование эволюции волновых фронтов [2|, обобщенного эйконала [20], позволяют численно строить функцию оптимального результата [1].

Геометрия волновых фронтов изучалась еще X. Гюйгенсом, в частности им был сформулирован принцип прямолинейного распространения света [20]. Негладкие особенности фронтов были классифицированы и изучены В.И. Арнольдом и его учениками еще в семидесятые годы XX века. Ими выделены так называемые «множества симметрии», в которых системы лучей, каустик и волновых фронтов имеют особенности [2, 3]. Ирландский математик П. Джиблин предложил способы использования

этих множеств в геометрической оптике и компьютерной графике [5]. Топология множеств симметрии изучена В.Д. Седых: выведены их эйлеровы характеристики в пространствах размерности до 6 включительно, количество и связь составляющих их гладких многообразий [19].

Понятие меры певыпуклости множества в произвольном евклидовом пространстве впервые предложено В.Н. Ушаковым в работе [22]. Эта мера имеет смысл угла и опирается на свойство проекций точки х на заг мкнутое множество М (т.е. ближайших к х в евклидовой метрике точек из М). Вводится понятие а-множества, обобщающее понятие выпуклого множества.

В семидесятые года XX века H.H. Красовский и А.И. Субботин [10] ввели понятие и-стабильных и и-стабильных функций, которые мажорируют и минорируют функцию цены. Функция цены является единственной функцией, которая одновременно u-стабильна и и-стабильна, а в точках дифференцируемости удовлетворяет дифференциальному уравнению в частных производных первого порядка (уравнению Айзекса-Беллмана). Эти свойства определяют одно и только одно обобщенное (минимаксное) решение уравнения Айзекса-Беллмана. Оно может быть определено несколькими эквивалентными способами, в том числе опираясь на субдифференциалы и супердифференциалы функции.

Множества негладкости функции цены дифференциальной игры или задачи управления имеют особый смысл с точки зрения построения оптимальных траекторий. Различные виды сингулярных поверхностей классифицированы Р. Айзексом в работе [1], и исследованы затем H.H. Красовский и А.И. Субботиньм [10]. В настоящее время эти работы продолжаются в ИММ УрО РАН [9].

Уравнение типа эйконала изучались в геометрической оптике с целью построения линий распространения света и волновых фронтов от источника на плоскости или в трехмерном пространстве. Эйконал называется еще оптическим путем волны и в общем случае является негладкой функцией. Один из подходов к его построению, как было отмечено выше, предложен С.Н. Кружковым [12]. В диссертации однако делается упор на построение минимаксного решения уравнения Айзекса-Беллмана для за-

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

Стабильные мосты в дифференциальных играх, которые рассматриваются в третьей главе, являются одной из главных тем научной школы H.H. Красовского на Урале. Выделение максимального стабильного моста в пространстве позиций — одна из основных и наиболее сложных задач, возникающих на пути построения решения дифференциальной игры. Задачи о сближении с целевым множеством в момент $ и на отрезке времени [to,i9] являются одними из наиболее важных в теории дифференциальных игр. Они связаны с многими крупными задачами оптимального гарантированного управления динамическими системами [24], в частности, — с задачей об оптимальном быстродействии для динамических систем, подверженных влиянию помех [10]. Кроме того, в рамки общей постановки таких задач укладываются многие конкретные дифференциальные игры [1].

Ведя исследования дифференциальных игр в рамках позиционного подхода, центральными элементами которого являются множества позиционного поглощения — максимальные u-стабильные мосты [10], В.Н. Ушаков, Х.Г. Гусейнов и A.M. Тарасьев [6] предложили инфини-тизимальные конструкции при построении стабильных мостов, которые используют методы негладкого анализа. Они позволяют свести установление совпадения максимальных м-стабильных мостов для стационарных систем к проверке относительно простых соотношений для векто-граммы скоростей в точках границы целевого множества [23].

Цель работы.

К основным целям диссертации относятся:

1. Изучение и характеризация свойств плоских невыпуклых множеств с негладкой границей посредством множеств симметрии;

2. Разработка и реализация алгоритмов аналитического и численного построения функции оптимального результата для одного класса задач быстродействия с невыпуклым целевым множеством;

3. Разработка и реализация алгоритмов численного построения эволюции волновых фронтов и обобщенного эйконала для среды с постоянным показателем преломления;

4. Выявление необходимых и достаточных условий совпадения стабильных мостов в задаче сближения в двух, вообще говоря, различных игровых постановках — в игре сближения «в момент» и в игре сближения «к моменту».

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

Исследования проводятся в рамках подхода, разрабатываемого в научной школе H.H. Красовского [10, 11] по управлению и дифференциальным играм. Один из способов опирается на построение нормалей к кривой, ограничивающей целевое множество. В другом случае узлы выбираются фиксированными, и для каждого из них находится аппроксимация проекций на целевое множество.

Нахождение меры невыпуклости множества сводится к поиску максимума функции нескольких переменных на объединении нескольких нульмерных и одномерных многообразий. За счет параметризации точек эта задача в свою очередь сводится к нахождению экстремума функции одной переменной на нескольких ограниченных либо неограниченных интервалах. Для некоторых множеств найдено аналитическое значение меры невыпуклости.

Проверка совпадения максимальных u-стабильных мостов в двух различных задачах о сближении осуществляется различными способами. Одни из них базируется на проверке условий стабильности для векто-грамм скоростей и конусов, аппроксимирующих целевое множество в точках его границы. Другие — на свойствах гамильтониана системы в точках границы дМ множества М и на описании множества в виде системы неравенств.

Вычисления производились в программном пакете MATLAB 6.1 [15], который позволяет использовать математические библиотеки для ускорения составления алгоритмов. В нем предусмотрена визуализация результатов, включая анимацию и трехмерную графику.

Научная новизна.

Данная работа является продолжением исследований задач управления, которые проводились и проводятся в настоящий момент в Институте математики и механики УрО РАН и, в том числе, в отделе динамических систем Института. Эта диссертация продолжает исследования по вычислению меры невыпуклости плоских множеств, построению обобщенных решений уравнений Гамильтона-Якоби, а также по нахождению условий совпадения максимальных ц-стабильных мостов в дифференциальных играх в момент и к моменту.

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

При решении второй задачи тоже используются множества симметрии (биссектрисы). С их помощью строятся волновые фронты — линии уровня функции оптимального результата в задаче быстродействия с круговой индикатрисой скоростей. На границе множества выделяются характеристические точки (псевдовершины), отвечающие за зарождение биссектрисы множества.

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

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

Практическая ценность работы обусловлена возможностью применения полученных в ней результатов в различных отраслях знания. Построения волновых фронтов и вычисление эйконала является важной задачей оптики и электродинамики [20]. В диссертации исследована геометрия волновых фронтов с круговой вектограммой для случая невыпуклого источника, имеющего кусочно-гладкую границу. Разработаны аналитические и численные алгоритмы эволюции волновых фронтов, предложены процедуры построения многообразий, «сотканных» из изломов волновых фронтов [2]. Обоснована формула минимаксного (обобщенного) решения задачи Дирихле для уравнения типа эйконала в случае изотропной среды при предположении, что краевое множество замкнуто, причем имеет кусочно-гладкую границу. Предложен конструктивный подход к построению минимаксного решения. Результаты, полученные в части совпадения максимальных и-стабильных в различных дифференциальных играх позволяют сводить решение более сложной задачи к решению более простой [14, 17].

Апробация работы. Главные положения диссертационной работы докладывались и обсуждались на российских и международных конференциях. Сделаны доклады

• Конференция «Теория управления и математическое моделирование». Ижевск, УдГУ, 3-8 июля 2006.

• Научный семинар Математическая теория оптимального управления и теория дифференциальных включений. М., МГУ. 2006.

• IX Международная Четаевская конференция «Аналитическая механика, устойчивость и управление движением», Иркутск, июнь 2007.

• XXII Международная конференция «Дифференциальные уравнения и смежные вопросы», посвященной памяти И.Г.Петровского, Москва, МГУ, 21-25 мая 2007.

• «Колмогоровские чтения. Общие проблемы управления и их приложения. Проблемы преподавания математики», Тамбов, октябрь 2007.

• Международная конференция Дифференциальные уравнения и топология. М.: МГУ. 2008.

• Международная конференция Дифференциальные уравнения и топология. М.: МГУ. 2008.

• Международная конференция «Дифференциальные уравнения. Функциональные пространства. Теория приближений», посвященная 100-летию со дня рождения C.JI. Соболева, Ин-т математики СО РАН. Новосибирск, октябрь 2008.

• Научный конференция-семинар «Теория управления и математическое моделирование», посвященная памяти профессора Н.В. Азбелева, Ижевск, май 2008.

• Международная конференция «Дифференциальные уравнения и топология». М.: МГУ. 2008.

• «Алгоритмический анализ неустойчивых задач». Международная конференция, посвященная 100-летию В.К. Иванова. Екатеринбург. ИММ УрО РАН, 2008.

Публикации. Основные результаты диссертации опубликованы в шестнадцати работах, приведенных в конце автореферата, шесть из них в изданиях, включенных в перечень ВАК. В работах, выполненных в соавторстве с научным руководителем, A.A. Успенским, и чл.-корр. РАН В.Н. Ушаковым, исследованы задачи оптимального управления и дифференциальные игры методами, предложенными в диссертации. Приведены примеры численного моделирования ряда задач и визуализации результатов.

Структура и объем диссертации. Диссертационная работа состоит из введения, пяти глав, списка литературы, включающего 117 названий, и приложения. Общий объем работы составляет 150 страниц.

СОДЕРЖАНИЕ ДИССЕРТАЦИИ

Диссертационная работа состоит из настоящего введения, списка сокращений и обозначений, трех глав, объединяющих 17 параграфов и списка литературы. Общий объем диссертации составляет 150 страниц, библиографический список включает 131 наименование, иллюстративный материал насчитывает 35 рисунков. Нумерация параграфов осуществляется в пределах каждой главы. Нумерация формул двойная: в первой позиции указывается номер параграфа, в котором приведена формула, во второй — порядковый номер формулы в этом параграфе. Такая же нумерация принята для определений, теорем, замечаний, примеров и рисунков.

Глава I диссертации посвящена изучению меры невыпуклости множества. В ней исследуются множества симметрии, волновые фронты, угловая характеристика точек относительно множества и их взаимосвязь. В параграфе 1.1 дается определение угловой характеристики точки ам (г) и меры невыпуклости а(М) замкунутого множества в пространстве К2. Для их вычисления вводятся новые понятия — биссектрисы Ь(М) и псевдовершины \У множества М С К2. Биссектриса Ь(М) является частным случаем множества симметрии с точки зрения геометрии множества М, для всех ее точек угловая характеристика отлична от нуля. Псевдовершины — характеристические точки границы дМ множества М, отвечающие за зарождение биссектрисы.

В параграфе 1.2 изучаются свойства биссектрисы и псевдовершин. Приводятся условия нахождения псевдовершин и формулы для построения гладких участков биссектрисы. Доказана теорема о том, что псевдовершины, лежащие на участках гладкости кривой Г, совпадающей с границей дМ множества М, есть точки локального максимума кривизны Г. Выписано уравнение, связывающее координаты двух точек, являющихся проекциями на множество М одной точки у биссектрисы Ь{М) (то есть являющимися ближайшими к ней в евклидовой метрике точками множества М).

Параграф 1.3 посвящен изучению свойств проекций точек биссектри-

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

Параграф 1.4 описывает строение биссектрис в окрестности точек прекращения. Выписаны координаты этих точек для разных типов, порождающих их псевдовершин. Доказана теорема о том, что если в псевдовершине W определен радиус кривизны R кривой Г, то точка прекращения V биссектрисы L(M) совпадает с центром кривизны Г в псевдовершине W. В случае негладкой псевдовершины W имеет место совпадение V и W. Если в псевдовершине определены предельные значения радиуса кривизны Г с двух сторон, то точка V есть линейная комбинация пределов центра кривизны в W слева и справа.

Параграф 1.5 посвящен нахождению достаточных условия гладкости биссектрисы и построению к ней касательной. Доказано, что если точка у € L(M) имеет ровно две проекции xi и хг, лежащие на гладких участках кривой Г, то в ней определена касательная П к биссектрисе L(M). При этом П совпадает с биссектрисой угла Zxi,y,x2-

В параграфе 1.6 предложены аналитические методы нахождения меры невыпуклости а(М) для случаев, когда М — подграфик кусочно-гладкой функции / ; ß н 1. Приведена и доказана теорема о мере невыпуклости подграфика выпуклой всюду дифференцируемой функции f(x):

а(М) = lim arctg/'(a:)— lim arctg/'(a;).

Z-t-rOO x—*—oo

Аналогично приведена теорема для подграфика функции f(x), дифференцируемой всюду, кроме точки для которой выполняется условия

1) lim f'(x) < lim f'(x)

, x<xm x—*x,, x>x*

2)f(x) — выпуклая вверх функция на интервалах (—оо, х*) и (х„, оо). Тогда мера невыпуклости множества М = hypo / равна

а(М) — lim arctg/'(:r) — lim arctg/'(a;).

X—*>x^x>x» z~*x,,x<x.

Параграф 1.7 содержит примеры построения биссектрис и вычисления меры невыпуклости для множеств с кусочно-гладкой границей. Рассмотрены множества, являющиеся подграфиками алгебраических, тригонометрических, показательных функций, а также их склеек. Для наиболее простых случаев мера невыпуклости находится аналитически по формулам параграфа 1.6. В остальных мера невыпуклости находится численными методами.

Глава II диссертации посвящена изучению свойств решения и = и(х) задачи быстродействия с круговой вектограммой скоростей и невыпуклым целевым множеством.

В параграфе 2.1 приводятся различные постановки задачи, в том числе в виде краевой задачи для уравнения эйконала

/ди\2 /ди\2 , „

«|г = о,

и уравнения в частных производных типа Гамильтона-Якоби

. ( ди ди\ , „ mm уц- + vi— ) +1 — 0, «/:|И<1\ OX dyj

u|r = 0.

Здесь ||i/|| = а/vf + i/| — норма вектора v = (fi,z/2). Краевое условие u|r = 0 определено на границе Г = дМ замкнутого множества М С R2.

Собственно задача быстродействия для динамической системы имеет вид:

х = Vi V = v2

где управление v — (t/j, г/2) стеснено ограничением ||i/|| ^ 1.

В параграфе 2.2 сформулирована и доказана теорема о том, что функция евклидового расстояния от точки до множества М является обобщенным (минимаксным) решением уравнения в частных производных Айзекса-Беллмана с краевым условием, определенным на границе множества М. Изучаются свойства биссектрисы множества, как множества негладкое™ функции оптимального результата и сингулярной (рассеивающей) линии для задачи управления.

Параграф 2.3 посвящен изучению асимптотики линий негладкости решения и(х). Доказано теорема о достаточных условиях существования асимптот биссектрисы. Если множество М — односвязное и невыпуклое, его дополнение £) = со М\М до выпуклой оболочки со М состоит из конечного числа компонент связности Д, г = 1,..., п и некоторая компонента связности Di — ограниченное множество, то множество дВ{ П Э(со М) является отрезком, и срединный перпендикуляр к нему т] — асимптота биссектрисы Ь(М). Одновременно показано, что данное условие не является необходимым. Возможно существование асимптот в частности в случае, когда дополнение со М\М — множество односвязное, но неогра-пичешюе.

В параграфе 2.4 приведен пример аналитического построения функции оптимального результата для задачи быстродействия с невыпуклым целевым множеством. А именно, когда множество М — подграфик функции /(ж) = х2. Биссектриса лежит на оси ординат

Ь(М) = {(х,у) : а; >0.5, у = 0}.

А функция оптимального результата и(х, у) вычисляется по формуле

и(х, у) = ^{х - хРУ + {у - Хр)2,

т2 (2у - l)3 ,

ТБ

V5?

arceos

z2 Tü

cos ■

2у-1

arceos

—З-у/бж

(2у -1)3/2 "3 Зл/бж

ж ^ 0, 2/ > Е(х)

■ cos-

(2t/-l)3/2 Т-

, х > 0, у > Е(х)

27 Ï63

Здесь Е(х) = - +

Параграф 2.5 содержит примеры численного конструирования функции оптимального результата и(х. у) в задаче быстродействия для множеств с кусочно гладкой границей на основе построения их биссектрис. В качестве целевых множеств взяты подграфики функций и одно множество, ограниченное эллиптической кривой. Выделены области их негладкости, в частности проиллюстрировано, что в точках биссектрисы L{M) функция и(х, у) суисрдифференцируема.

Глава III диссертации связана с изучением взаимосвязи дифференциальных игр сближения-уклонения «в момент» и «к моменту».

В параграфе 3.1 рассматривается конфликтно-управляемая система, поведение которой на конечном промежутке времени [i0, t>] описывается векторным дифференциальным уравнением.

х = f{t,x,u, v), ¡r[i0] = Xq, и 6 P, v e Q.

Здесь x — m-мерный фазовый вектор системы, и — управление первого игрока, v — управление второго игрока, Р и Q — компакты в евклидовых пространствах Rp и R9 соответственно.

Изучаются и сравниваются две игровые задачи о сближении системы с терминальным множеством M в фазовом пространстве [11]. В первой из них первому игроку требуется обеспечить с помощью позиционного управления попадание фазового вектора системы на M в конечный момент времени Во второй задаче требуется обеспечить с помощью позиционного управления попадание фазового вектора системы на M не

позже момента д.

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

В параграфе 3.3 приводятся теоремы о геометрических условиях совпадения максимальных и-стабильных мостов в двух задачах о сближении.

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

^ = В{х)и + С(х)ь,

где

х — (х\,х2) € К2 — фазовый вектор системы,

и£ Р, V 6 <Э — управления 1 и 2 игроков,

Р = <5 = {(21,2:2) : —1 ^ 11 ^ 1,12 = 0} ~ ограничения на управление,

В одном случае в качестве целевого множества М взят квадрат с вершинами Тх = (1,1); Т2 = (1, —1); Тз = (—1, -1); Т4 = (—1,1), и показаг но, что максимальные м-стабильные мосты в обеих задачах совпадают. В другом взят треугольник с вершинами Т1 = (0,1); Т2 = (—\/3/2, —0.5); Т3 = (\/3/2, —0.5). Тогда максимальные м-стабильные мосты не совпадают, а точнее мост в задаче о сближении «в момент» строго вложен в мост в задаче о сближении «к моменту».

СПИСОК ЦИТИРОВАННОЙ ЛИТЕРАТУРЫ

1. Айзеке. Р. Дифференциальные игры. М.: Мир, 1967. — 479 с.

2. Арнольд В.И. Особенности каустик и волновых фронтов. М.: ФАЗИС, 1996. - 334 с.

3. Арнольд В.И. Теория катастроф. М.: Изд.-во МГУ, 1983. — 80 с.

4. Бабич В.М., Булдырев B.C., Молотков И.А. Пространственно-временной лучевой метод: Линейные и нелинейные волны. М.: Изд-во Ленингр. ун-та, 1985. - 272 с.

5. Брус Дж., Джиблин П. Кривые и особенности. М.: Мир, 1988. — 262 с.

6. Гусейнов Х.Г., Ушаков В.Н. Дифференциальные свойства интегральных воронок стабильных мостов // Прикладная математика и механика. 1991. №55 (1), С. 72-78.

7. Дарьип А.Н. Синтез управлений при двойных и неоднородных ограничениях // Автореферат диссертации на соискание ученой степени кандидата физ.-мат. наук. 01.01.02. — дифференциальные уравнения. МГУ им. М.В. Ломносова. М. 2004. 15 с.

8. Демьянов В.Ф., Васильев Л.В. Недифференцируемая оптимизация. М.: Наука, 1981. - 384 с.

9. Исакова Е.А., Логунова Г.В., Пацко B.C. Построение стабильных мостов в линейной дифференциальной игре с фиксированным моментом окончания — В кн.: Алгоритмы и программы решения линейных дифференц. игр: Сб. ст. ИММ УНЦ АН СССР. Свердловск, 1984. С. 127-158.

10. Красовский H.H., Субботин А.И. Позиционные дифференциальные игры. М.: Наука, 1974. — 456 с.

11. Красовский H.H. Игровые задачи о встрече движений. М.: Наука, 1970. - 420 с.

12. Кружков С.Н. Обобщенные решения уравнений Гамильтона-Якоби типа эйконала //I. Математический сборник. 1975. Т. 98, вып. 3. С. 450-493.

13. Кряжимский A.B., Осипов /О. С. Дифференциально-разностная игра сближения с функциональным целевым множеством. // Прикл. математика и механика. 1973. Т. 37, выи. 1. С. 3-13.

14. Курэюанский A.B. Альтернированный интеграл Понтрягина в теории синтеза управлений // Тр. Мат. Ин-та им. Стеклова. 1999. Т. 224. С. 234-248.

15. Мэтьюз Д., Фипк К. Численнтые методы. Использование MATLAB (3-е издание). СПб.: Вильяме, 2001. — 720 с.

16. Панасенко Е.А., Тонкое Е.Л. Инвариантные и устойчиво инвариантные множества дифференциальных включений // Тр. МИАН. 2008. Т. 262, С. 202-221.

17. Понтрягип Л.С. О линейных дифференциальных играх. I // Докл. АН СССР. 1967. Т. 174, № 6. С. 1278-1281.

18. Рашевский П.К. Курс дифференциальной геометрии. М.: Едиториал, УРСС, 2003. - 432 с.

19. Седых В. Д. Соотношения между эйлеровыми характеристиками многообразий особенностей коранга 1 фронта общего положения / / Докл. РАН. - 2002. - т.383, № 6. - с. 735-39. - Библиогр.: 7 назв.

20. Слюсарев Г.Г. Геометрическая оптика. М.: Издательство Академии наук СССР, 1946. - 332 с.

21. Субботин А.И. Обобщенные решения уравнений в частных производных первого порядка. Перспективы динамической оптимизации. Москва-Ижевск: Институт компьютерных технологий, 2003. — 336 с.

22. Успенский A.A., Ушаков В.Н., Фомин А.Н. а-множества и их свойства// Ин-т математики и механики УрО РАН. — Екатеринбург, 2004. — 62 е.: 38 ил. — Библиогр.: 7 назв. — Рус — Деп. в ВИНИТИ 02.04.04, № 543-В2004

23. Ушаков В.Н. Процедуры построения стабильных мостов в дифференциальных играх.: Дис. д-ра физ.-мат. наук: 01.01.02. Свердловск, 1991. 308 с. / Ин-т математики и механики УрО АН СССР.

24. Чепцов А.Г. Об игровой задаче сближения в заданный момент времени // Мат. сб. 1976. Т. 99, № 3. С. 394-420.

ПУБЛИКАЦИИ АВТОРА ПО ТЕМЕ ДИССЕРТАЦИИ

Статьи, опубликованные в ведущих рецензируемых научных журналах:

1. Лебедев П.Д. Вычисление меры невыпуклости плоских множеств // Труды Института математики и механики. Том 13, № 3. 2007 г. С. 84-94.

2. Лебедев П.Д., Успенский A.A. Геометрия и асимптотика волновых фронтов // Известия вузов. Математика. №3 (550). 2008 г. С. 27-37.

3. Лебедев П.Д., Успенский A.A., Ушаков В.Н. Построение минимаксного решения уравнений типа эйконала // Труды Института математики и механики. Том 14, № 2. 2008 г. С. 182-191.

4. Лебедев П.Д., Успенский A.A. Алгоритмы построения функции оптимального результата в задаче с простой динамикой // Вестник Удмуртского университета. Математика. Механика. Компьютерные науки. Вып. 2. 2008 г. С. 152-154.

5. Успенский A.A., Лебедев П.Д. Условия трансверсальности ветвей решения нелинейного уравнения в задаче быстродействия с круговой индикатрисой // Труды Института математики и механики. Том 14, № 4. 2008. С. 82-99.

6. Успенский A.A., Лебедев П.Д. Процедуры вычисления меры невыпуклости плоского множества // Журнал вычислительной математики и математической физики. Т. 49, № 3. 2009. С. 431-440.

7. Успенский A.A., Лебедев П.Д. Построение функции оптимального результата в задаче быстродействия на основе множества симметрии // Автоматика и телемеханика. № 7. 2009. С. 50-57.

8. Ушаков В.Н., Гусейнов Х.Г., Латушкин Я.А., Лебедев П.Д. О совпадении максимальных стабильных мостов в двух игровых задачах о сближении для стационарных управляемых систем // Труды Института математики и механики. Том 15, № 3. 2009. С. 219—240.

Другие публикации:

9. Лебедев П.Д., Успенский A.A. К вопросу о геометрии волновых фронтов // Известия Института математики и информатики. Ижевск, УдГУ. 2006. Вып.З (37). С. 79-80.

10. Лебедев П.Д., Успенский A.A. К вопросу о геометрии волновых фронтов // Научный семинар Математическая теория оптимального управления и теория дифференциальных включений. М., МГУ. 200G. С. 40-41.

11. Лебедев П.Д., Успенский A.A. Исследование геометрии и асимптотики волновых фронтов в некоторых задачах управления // Тр. 9-й междунар. Четаевской конф. 2007. Т. 5. С. 224-236.

12. Лебедев П.Д., Успенский A.A. Аналитическое и численное конструирование функции оптимального результата для одного класса задач быстродействия // Прикладная математика и информатика. Труды факультета ВМК МГУ. №27. 2007 г. С. 65-79.

13. Лебедев П.Д., Успенский A.A. Численно-аналитические методы построения волновых фронтов в задачах управления и геометрической оптике // Вестник Тамбовского университета. 2007. Т. 12, вып. 4. С. 538539.

14. Лебедев П.Д., Успенский A.A. Построение обобщенного решения уравнения в частных производных первого порядка // Алгоритмический анализ неустойчивых задач. Тезисы докладов международной конференции, посвященной 100-летию В.К. Иванова. 2008 г. С. 240-241.

15. Успенский A.A., Лебедев П.Д. Аналитические и численные подходы к построению функции оптимального результата для задачи быстродействия ,// Международная конференция «Дифференциальные уравнения и топология». Тезисы Докладов. М.: МГУ. 2008. С. 409-410.

16. Ушаков В.Н., Латпушкин Я.А., Лебедев П.Д. Критерии совпадения максимальных стабильных мостов в двух игровых задачах о сближении // Международная конференция «Дифференциальные уравнения и топология». Тезисы докладов. М.: МГУ. 2008. С. 411-412.

17. Ушаков В.Н., Успенский A.A., Лебедев П.Д. Конструкции негладкого анализа при построении минимаксных решений уравнений в частных производных первого порядка // Международная конференция «Дифференциальные уравнения и топология». Тезисы Докладов. М.: МГУ. 2008. С. 412.

18. Успенский A.A., Лебедев П.Д. Численно-аналитический подход к построению функции оптимального результата для одного класса задач быстродействия // Международная конференция «Управление динамическими системами». ИПМех РАН. Тезисы докладов. 26-30 января 2009. С. 65.

Подписано в печать 1С.09.2009. Формат 60x90 1/16 Гарнитура Times. Печать ротационная трафаретная. Усл. печ. л. 5. Уч.-изд.л. 4,3 Тираж 150 экз. Заказ № 2152 Отпечатано в типографии «Артикул» 620026 г. Екатеринбург, ул. Декабристов, 20

Оглавление автор диссертации — кандидата физико-математических наук Лебедев, Павел Дмитриевич

Введение

Список публикаций автора по теме диссертации

Список докладов на конференциях по теме диссертации

1 Мера невыпуклости множества и алгоритмы ее вычисления

1.1 Основные понятия, определения и формулировки.

1.2 Построение псевдовершин и биссектрисы множества.

1.3 Расположение а-симметричных точек в окрестности псевдовершины

1.4 Крайние точки биссектрисы.

1.5 Гладкость биссектрисы

1.6 Процедуры отыскания меры певыпуклости а(М)

1.7 Примеры вычисления меры невыпуклости для плоских множеств.

2 Применение множеств симметрии при построении решений задач быстродействия и уравнения эйконала

2.1 Постановки краевых задач

2.2 Структура минимаксного решения уравнения в частных производных первого порядка типа эйконала.

2.3 Асимптотика биссектрисы.

2.4 Пример аналитического построения функции оптимального результата.

2.5 Примеры численного построения функции оптимального ре

3 Критерии совпадения максимальных стабильных мостов в двух игровых задачах о сближении для стационарных управляемых систем

3.1 Постановка задач о сближении.

3.2 Оператор стабильного поглощения и стабильные мосты в задачах о сближении

3.3 Критерии совпадения максимальных w-стабильных мостов для стационарных управляемых систем.

3.4 Аналитические критерии совпадения максимальных и-стабильных мостов для стационарных управляемых систем

Введение 2009 год, диссертация по информатике, вычислительной технике и управлению, Лебедев, Павел Дмитриевич

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

Предыстория и актуальность темы

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

Современный облик теории управления движением динамической системы сформировался в значительной степени под влиянием работ отечественных математиков JI.C. Поптрягина, Н.Н. Красовского, ставших основателями известных научных школ по теории управления. Весомые, основополагающие результаты были получены их коллегами и учениками — представителями московской научной школы Е.Ф. Мищенко, Р.В. Гамкрелидзе, В.Г. Болтянским, представителями уральской научной школы А.Б. Кур-жанским, Ю.С. Осиповым, А.И. Субботиным, а также зарубежными учеными Р. Калманом, Р. Беллманом, Р. Айзексом, У. Флемингом, Ж-П. (Эбеном. Существенный прогресс в становлении и развитии теории управления связан также с именами Э.Г. Альбрехта, В.Д. Батухтипа, Р.Ф. Габа-сова, А.Я. Дубовицкого, С.Г. Завалищина, Ф.М. Кирилловой, А.А. Мели-кяна, А.А. Милютина; М.С. Никольского, А.А. Петросяна, Б.Н. Пшеничного, Н.Н. Субботиной, В.М. Тихомирова, E.JI. Тонкова, В.Е. Третьякова, В.Н. Ушакова, А.Г. Ченцова, Ф.Л. Черноусько, А.А. Чикрия и многих других. Актуальность изучения управляемых систем'обусловлена наличием' многочисленных приложений в различных отраслях знания — в механике, робототехнике, оптике, экономике, биологии. Немаловажным побудительным мотивом в, исследовании являются внутренние потребности, возникшие в математической теории управления динамическими системами в условиях конфликта и неопределенности, а также стремление исследователей привлечь для изучения динамических задач конструкции из других разделов математики — негладкого анализа [26,27], дифференциальной геометрии [18], теории особенностей дифференцируемых отображений [7-9,106].

Осуществляемые в работе исследования проводятся в рамках концепции теории позиционных дифференциальных игр, развиваемой в научной школе Н.Н. Красовского [35-42]. Теория позиционных дифференциальных игр, обогащенная результатами его соратников, учеников и последователей, объединяет конструктивные методы решения широкого круга проблем от теорем существования и единственности решения-до разработки и реализации вычислительных алгоритмов [25,65,103,104,121-124]. Также используются конструкции отпочковавшейся от этой теории и получившей глубокое развитие в работах А.И. Субботина теории минимаксного решения уравнения в частных производных первого порядка [83-88,130,131].

Результатам последнего времени, связанным с проблемой построения (аналитического или приближенного) обобщенных решений уравнений типа Гамильтона-Якобй [90,91,118,119], предшествовали работы С.К. Годунова [22], Е. Hopf [115,116], P.D. Lax [120], других авторов. В 70-х годах прошлого века С.Н. Кружков [43-45], реализуя метод исчезающей вязкости, ввел определение фундаментального (обобщенного) решения уравнения в частных производных первого порядка типа эйконала посредством предельного перехода по параметру малости при старшей производной от решений соответствующих уравнений в частных производных второго порядка. Позднее, в 80-е годы, действуя аналогичным образом, М. Крэндалл и П.Л. Лионе [110,111] ввели определение обобщенного решения уравнения в частных производных первого порядка, названного ими вязкостным решением. Несмотря на разницу минимаксного и вязкостного подходов, определяемый объект является одним и тем же. Отличительной чертой минимаксного подхода при изучении и численном построении обобщенных решений является активное вовлечение методов, средств и конструкций выпуклого и негладкого анализа. При этом используются наработки отечественных и зарубежных математиков [26,27,74,105]. Значимость тематики обобщенных решений уравнений в частных производных первого порядка для современной математики иллюстрируется тем фактом, что в 1994 г. за работы по теории вязкостных решений П.-Л. Лионе был удостоен Филдсовской премии. Вопросами построения разностных операторов в рамках вязкостного подхода успешно занимались и занимаются П. Соуганидис [129], М. Бар-ди [107], М. Фалькопе [108]. Надо сказать, что полезными с точки зрения разработки конструктивных подходов к построению обобщенных решений уравнений в частных производных первого порядка оказались методы теории особенностей дифференцируемых отображений, разрабатываемые В.И. Арнольдом [3-8,106] , его коллегами, зарубежными авторами Т. По-стоном, И. Стюартом, Дж. Брусом, П. Джиблиным [17,70,109,112,113]. Средствами этой теории, в частности, формируются списки типичных особенностей каустик и волновых фронтов, предлагаются подходы к построению дискриминантных множеств.

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

V. ческими системами, описываемыми дифференциальными уравнениями с дискретными и\или распределенными параметрами [46]. Исследуются вопросы существования и единственности обобщенных решений, выявляются условия оптимальности гарантированного результата управления в таких задачах в случае исходных данных, удовлетворяющих условию Липшица [60,61].

Задачами быстродействия занимались многие исследователи, см. например [15,16,66]. Привлечение идей,результатов и конструкций указанных выше теорий позволяет в рамках минимаксного подхода разрабатывать I аналитические и аппроксимационные процедуры построения функции оптимального результата для задач динамического управления по быстродействию [30,31,117]. В диссертации эти подходы распространены на задачи быстродействия для частного случая вектограммы скоростей, а также задачи геометрической оптики [11]. Разработанные процедуры позволяют проводить численно-аналитическое конструирование эволюции волновых фронтов [14], обобщенного эйконала [19], позволяют численно строить функцию оптимального результата [15].

Геометрия волновых фронтов изучалась еще X. Гюйгенсом, в частности им был сформулирован принцип прямолинейного распространения света [28]. Негладкие особенности фронтов были классифицированы и изучены В.И. Арнольдом и его учениками еще в семидесятые годы XX века. Ими выделены так называемые «множества симметрии», в которых системы лучей, каустик и волновых фронтов имеют особенности [3-8,106]. Ирландский математик П. Джиблин предложил способы использования'этих множеств в геометрической оптике и компьютерной графике [17,112,113]. Топология множеств симметрии изучена В.Д. Седых: выведены их эйлеровы характеристики в пространствах размерности до 6 включительно, количество и связь составляющих их гладких многообразий [75,125,126]. Им же подробно изучены свойства невыпуклых множеств в плане особенностей их выпуклых оболочек [77-79].

Понятие меры невыпуклости множества в произвольном евклидовом пространстве впервые предложено В.Н. Ушаковым в работе [94]. Эта мера имеет смысл угла и опирается на свойство проекций точки х на замкнутое множество М (т.е. ближайших к х в евклидовой метрике точек из М). Вводится понятие «-множества, обобщающее понятие выпуклого множества.

В семидесятые года XX века Н.Н. Красовский и А.И. Субботин [35, 36] ввели понятие ii-стабильных и v-стабильных функций, которые мажорируют и минорируют функцию цены. Функция цепы является единственной функцией, которая одновременно гг-стабильна и ^-стабильна, а в точках дифференцируемое™ удовлетворяет дифференциальному уравнению в частных производных первого порядка (уравнению Айзекса-Белл-мана). Эти свойства определяют одно и только одно обобщенное (минимаксное) решение уравнения Айзекса-Беллмана. Оно может быть определено несколькими эквивалентными образами, в том числе опираясь на субдифференциалы и супердифференциалы функции.

Множества негладкости функции цены дифференциальной игры или задачи управления имеют особый смысл с точки зрения прохождения в них оптимальных траекторий. Различные виды сингулярных поверхностей классифицированы Р. Айзексом в работе [1], и исследованы затем Н.Н. Красовский и А.И. Субботиным [35]. В настоящее время эти работы продолжаются в ИММ УрО РАН.

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

С.Н. Кружковым [43]. В диссертации однако делается упор на построение минимаксного решения уравнения Айзекса-Веллмана для задачи быстродействия, которое в некотором смысле эквивалентно уравнению эйконала. Его решение для некоторых частных случаев изучались в частности в монографии [10]. Близкие задачи изучаются в отделе динамических систем ИММ УрО РАН [23,67,89]. Из зарубежных исследователей схожими задачами занимается, например, J.A. Sethian [127,128].

Стабильные мосты в дифференциальных играх, которые рассматриваются в третьей главе, являются одной из главных тем научной школы Н.Н. Красовского на Урале. Выделение максимального стабильного моста в пространстве позиций — одна из основных и наиболее сложных задач, возникающих на пути построения решения дифференциальной игры. Задачи о сближении с целевым множеством в момент д и па отрезке времени [£о,$] являются одними из наиболее важных в теории'дифференциальных игр. Они связаны с многими крупными задачами оптимального гарантированного управления динамическими системами, в частности, — с задачей об оптимальном быстродействии для динамических систем, подверженных влиянию помех [35,36,84]. Кроме того, в рамки общей постановки таких задач укладываются многие конкретные дифференциальные игры [1].

В отделе динамических систем ИММ УрО РАН найдено широкое практическое применение теории дифференциальных игр и задач оптимального управления. Одно из них — восстановления траектории самолета в процессе поступления замеров его геометрического положения [29,123,124].

Ведя исследования дифференциальных игр в рамках позиционного подхода, центральными элементами которого являются множества позиционного поглощения — максимальные ■и-стабильные мосты [35,36], В!Н. Ушаков, Х.Г. Гусейнов и A.M. Тарасьев [24,90-92,114] предложили инфините-зимальные конструкции при построении стабильных мостов, которые используют методы негладкого анализа. Они позволяют свести установление совпадения максимальных w-стабильных мостов для стационарных систем к проверке относительно простых соотношений для вектограммы скоростей в точках границы целевого множества. В диссертации на их основе разработаны и доказаны критерии совпадения максимальных стабильных мостов в задачах о достижении цели в фиксированный момент времени и на отрезке времени [96-101].

Цели работы

К основным целям диссертации относятся:

1. Изучение и характеризация свойств плоских невыпуклых множеств с негладкой границей посредством множеств симметрии;

2. Разработка и реализация алгоритмов аналитического и численного построения функции оптимального результата для одного класса задач быстродействия с невыпуклым целевым множеством;

3. Разработка и реализация алгоритмов численного построения эволюции волновых фронтов и обобщенного эйконала для среды с постоянным показателем преломления;

4. Выявление необходимых и достаточных условий совпадения стабильных мостов в задаче сближения в двух, вообще говоря, различных игровых постановах — в игре сближения «в момент» и в игре сближения «к моменту».

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

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

Нахождение меры невыпуклости множества сводится к поиску максимума функции нескольких переменных на объединении нескольких нульмерных и одномерных многообразий. За счет параметризации точек она в свою очередь сводится к нахождению экстремума функции одной переменной на нескольких ограниченных либо неограниченных интервалах. Для некоторых множеств найдено аналитическое значение меры невыпуклости.

Проверка совпадения максимальных w-стабильных мостов в двух различных задачах о сближении осуществляется различными способами. Одни из них базируется на проверке условий стабильности для вектограмм скоростей и конусов, аппроксимирующих целевое множество в точках его границы. Другие на свойствах гамильтониана системы в точках границы дМ множества М и на описании множества в виде системы неравенств.

Вычисления производились в программном пакете MATLAB 6.1 [2,32, 62,102], который позволяет использовать математические библиотеки для ускорения составления алгоритмов. В нем предусмотрена визуализация результатов, включая анимацию и трехмерную графику. Часть вычислений, относящаяся к главе III, выполнена А.Р. Матвийчуком в пакете С++.

Научная новизна

Данная работа является продолжением исследований задач управления, которые проводились и проводятся в настоящий момент в Институте математики и механики УрО РАН и, в том числе, в отделе динамических систем Института. Эта диссертация продолжает исследования по вычислению меры невыпуклости плоских множеств, построению обобщенных peшепий уравнений Гамильтона Якоби, а также по нахождению условий совпадения максимальных ад-стабильных мостов в дифференциальных играх в момент и к моменту.

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

При решении второй задачи тоже используются множества симметрии (биссектрисы). С их помощью строятся волновые фронты — линии уровня функции оптимального результата в задаче быстродействия с круговой индикатрисой скоростей. На границе множества выделяются характеристические точки (псевдовершины), отвечающие за зарождение биссектрисы множества.

Задача о совпадении максимальных 71-стабильных мостов в двух дифференциальных играх исследуется с применением инфинетизимальных конструкций, как то контингентные конуса, производные множества стабильного моста в прямом и обратном времени. Условие стабильности выписывается также для случая, когда целевое множество задано виде системы неравенств для непрерывно дифференцируемых функций и неравенства для функции дифференцируемой по направлениям. Используются свойства гамильтониана системы, в частности его выпуклость.

Теоретическая и практическая ценность

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

Практическая ценность работы обусловлена возможностью применения полученных в ней результатов в различных отраслях знания. Построения волновых фронтов и вычисление эйконала является важной задачей оптики и электродинамики [34,81]. В диссертации исследована геометрия волновых фронтов с круговой индикатрисой для случая невыпуклого источника, имеющего кусочно-гладкую границу. Разработаны аналитические и численные алгоритмы эволюции волновых фронтов, предложены процедуры построения многообразий, «сотканных» из изломов волновых фронтов [3]. Обоснована формула минимаксного (обобщенного) решения задачи Дирихле для уравнения типа эйконала в случае изотропной среды при предположении, что краевое множество замкнуто, причем имеет кусочно-гладкую границу. Предложен конструктивный подход к построению минимаксного решения. Результаты, полученные в части совпадения максимальных ii-стабильных в различных дифференциальных играх позволяют сводить решение более сложной задачи к решению более простой [47,48,63,68,69].

Структура, объем и краткое содержание

Диссертационная работа состоит из настоящего введения, списка сокращений и обозначений, трех глав, объединяющих 17 параграфов и списка литературы. Общий объем диссертации составляет 150 страниц, библиографический список включает 131 наименование, иллюстративный материал насчитывает 35 рисунков. Нумерация параграфов осуществляется в пределах каждой главы. Нумерация формул двойная: в первой позиции указывается номер параграфа, в котором приведена формула, во второй — порядковый номер формулы в этом параграфе. Такая же нумерация принята для определений, теорем, замечаний, примеров и рисунков.

Заключение диссертация на тему "Аналитические и численные процедуры построения решений некоторых задач управления"

Заключение

Диссертация посвящена теории оптимального управления и дифференциальных игр.

Основными результатами диссертации являются следующие.

1. Разработаны и реализованы алгоритмы построения множеств симметрии. Для двух классов множеств доказаны теореме об аналитическом выражении их меры невыпуклости. Предложены численные методы ее вычисления в остальных случаях.

2. Обоснована формула обобщенного (минимаксного) решения уравнения задачи Дирихле для в частных производных первого порядка типа Гамльтона-Якоби. Показана ее связь с решением уравнения эйконала. Предложены и реализованы алгоритмы численнрго построения решения и множеств его негладкости.

3. Исследована геометрия и эволюция волновых фронтов для одного класса задач быстродействия — линий уровня решения задачи Дирихле для соответствующего уравнения Айзекса-Беллмана. Для одного класса целевых множеств найдены достаточные условия существования асимптот у линий их негладкости.

4. Исследованы две игровые задачи о сближении системы с терминальным множеством М в фазовом пространстве: в момент и на отрезке времени. Получены необходимые и достаточные условия совпадения в них максимальных и-стабильных мостов, опирающиеся на унификацию вектограммы скоростей системы. Они записаны для случаев различной аппроксимации множества М инфинитезимальными конструкциям.

Библиография Лебедев, Павел Дмитриевич, диссертация по теме Математическое моделирование, численные методы и комплексы программ

1. Айзеке. Р. Дифференциальные игры. М.: Мир, 1967. — 479 с.

2. Ануфриев И.Е., Смирнов А.Б., Смирнова ЕЛ. MATLAB 7. СПб.: БХВ-Петербург, 2005. 1104 с.

3. Арнольд ВЛ. Особенности каустик и волновых фронтов. М.: ФАЗИС, 1996. 334 с.

4. Арнольд ВЛ. Теория катастроф. М.: Изд.-во МГУ, 1983. — 80 с.

5. Арнольд ВЛ. Особенности систем лучей. Успехи мат. наук. 1983. № 38 (2). С. 77-147.

6. Арнольд В Л. Инварианты и перестройки фронтов па плоскости // Труды мат. института им. Стсклова. 1995. № 209. С. 14-64.

7. Арнольд В Л. Лагранжевы многообразия с особенностями, асимптотические лучи и раскрытый ласточкин хвост / /Функциональный анализ и его приложения. 1981. № 15 (2). С. 1-14.

8. Арнольд В Л., Афраймович B.C., Ильяшенко Ю.С., Шильников Л.П. Теория бифуркаций. Москва, 1986: ВИНИТИ, том 5, 218 стр.

9. Арнольд ВЛ., Варченко АЛ., Гуссейн-Заде СЛ. Особенности дифференцируемых отображений. 1. Классификация критических точек, каустик и волновых фронтов. М., Наука. 1982. — 304 с.

10. Бабич В.М., Булдырев B.C., Молотков И.А. Пространственно-временной лучевой метод: Линейные и нелинейные волны. М.: Изд-во Лепингр. ун-та, 1985. — 272 с.

11. Бегунов Б.Н. Геометрическая оптика. Л.: Изд-во МГУ, 1966. — 211 с.

12. Благодатеких ВЛ. Введение в оптимальное управление. М: Высшая школа, 2001. 238 с.

13. Благодатских В.И., Филиппов А.Ф. Дифференциальные включения и оптимальное управление // Тр. Мат. ин-та им. Стеклова. 1985. Т. 169. С. 194-252.

14. Богаевский И.А. Перестройки фронтов в эволюционных семействах. // Тр. Мат. института им. Стеклова. 1995. № 209. С. 65-83.

15. Болтянский В.Г. Математические методы оптимального управления. М.: Наука, 1966. — 308 с.

16. Болтянский В.Г., Гамкрелидзе Р.В., Понтрягин Л. С. Теория оптимальных процессов. I. Принцип максимума В. Г. Болтянский, Изв. АН СССР-Сер. матем. 1960. № 24:1. С. 3-42

17. Брус Дж., Джиблин П. Кривые и особенности. М.: Мир, 1988. — 262 с.

18. Бюшгенс С.С. Дифференциальная геометрия. Государственное издательство технико-теоретической литературы, 1940. — 300 с.

19. Герцбергер М.М. Современная геометрическая оптика. М.: Изд-во иностр. литератууры, 1962. — 487 с.

20. Гилмор Р. Прикладная теория катастроф. Т. 1. М.: Мир, 1984. — 350 с.

21. Гилмор Р. Прикладная теория катастроф. Т. 2. М.: Мир, 1984. — 284 с.

22. Годунов С. К. О численном решении краевых задач для систем линейных обыкновенных дифференциальных уравнений // УМН. № 16:3 (99), 171-174. 1961.

23. Григорьева С.В., Пахотинских В.Ю., Успенский А.А., Ушаков В.Н. Конструирование решений в некоторых дифференциальных играх с фазовыми ограничениями // Математический сборник. 2005. Т. 196, № 4. С. 51-78.

24. Гусейнов Х.Г., Ушаков В.Н. Дифференциальные свойства интегральных воронок стабильных мостов // Прикладная математика и механика. 1991. №55 (1), С. 72-78.

25. Дарьин А. И. Синтез управлений при двойных и неоднородных ограничениях // Автореферат диссертации на соискание ученой степени кандидата физ.-мат. наук. 01.01.02. — дифференциальные уравнения. МГУ им. М.В. Ломносова. М. 2004. 15 с.

26. Демьянов В.Ф., Васильев JI.B. Недифференцируемая оптимизация. М.: Наука, 1981. — 384 с.

27. Демьянов В.Ф., Рубинов A.M. Основы негладкого анализа и квазидифференциальное исчисление. М.: Наука, 1990. — 432 с.

28. Залгаллер В.А. Теория огибающих. М.: Наука, 1975. — 100 с.

29. Камнева Л. В. Разрывная функция цены в игровых задачах быстродействия // диссертация кандидата физико-математических наук: 01.01.02 Екатеринбург. 2007. 120 с. РГБ ОД, 61:07-1/978.

30. Камнева Л. В. Достаточные условия стабильности для функции цены дифференциальной игры в терминах сингулярных точек // Приклад-пая математика и механика. 2003. Т. 67, вып. 3. С. 366-383.

31. Кетков Ю.Л., Кетков А.Ю., Шульц М.М. MATLAB 7: программирование, численные методы. СПб.: БХВ-Петербург, 2005. — 752 с.

32. Колокольцов В.Н., Маслов В.П. Задача Коши для однородного уравнения Беллмана // Доклады АН СССР. 1987. № 296 (4), С. 796-800. ,

33. Кравцов Ю.А., Орлов Ю.И. Геометрическая оптика неоднородных сред. М.: Наука, 1980. — 304 с.

34. Красовский Н.Н., Субботин А.И. Позиционные дифференциальные игры. М.: Наука, 1974. — 456 с.

35. Красовский Н.Н., Субботин А.И. О структуре игровых задач динамики // Прикладная математика и механика. 1971. № 35(1). С. 110-122.

36. Красовский Н.Н. Игровые задачи о встрече движений. М.: Наука, 1970. 420 с.

37. Красовский Н.Н. К задаче унификации дифференциальных игр // Докл. АН СССР. 1976. Т. 226, № 6. С. 1260-1263.

38. Красовский Н.Н. Унификация дифференциальных игр // Труды Института математики и механики УНЦ АН СССР. Свердловск, 1977. Вып. 24: Игровые задачи управления. С. 32-45.

39. Красовский Н.Н. О синтезе в дифференциальной игре // Прикл. математика и механика. 1986. Т. 50, № 6. С. 898-902.

40. Красовский Н.Н. О некоторых задачах управления. // Тр. Мат. ин-та РАН. 1999. Т. 224. С. 208-217.

41. Красовский Н.Н., Субботин А.И., Ушаков В.Н. Минимальная дифференциальная игра // ДАН СССР. 1972. Т. 206, № 2, С. 277-280.

42. Крулсков С.Н. Обобщенные решения уравнений Гамильтона-Якоби типа эйконала //I. Математический сборник. 1975. Т. 98, вып. 3. С. 450493.

43. Кружков С.Н. К методам построения обобщенных решений задачи Коши для квазилинейного уравнения первого порядка //УМН. 1965. Т. 20., № 6(126), С. 112-118.

44. Крусисков С.Н. Обобщенные решения нелинейных уравнений первого порядка со многими независимыми переменными //I. Математический сборник. Т. 70, вып. 112(3). С. 394-415. 1966.I

45. Кряжимский А.В., Осипов Ю.С. Дифференциально-разностная игра сближения с функциональным целевым множеством. // Прикл. математика и механика. 1973. Т. 37, вып. 1. С. 3-13.

46. Куржанский А. В. Альтернированный интеграл Поптрягина в теории синтеза управлений // Тр. Мат. Ин-та им. Стеклова. 1999. Т. 224. С. 234-248.

47. Куржанский А.Б. О синтезе управлений по результатам измерений // Прикл. математика и механика. 2004. Т. 68, вып. 4. С. 547-563.

48. Лебедев П.Д., Успенский А.А. К вопросу о геометрии волновых фронтов // Известия Института математики и информатики. Ижевск, УдГУ. 2006. Вып. 3 (37). С. 79-80.

49. Лебедев П.Д., Успенский А.А. К вопросу о геометрии волновых фронтов // Научный семинар Математическая теория оптимального управления и теория дифференциальных включений. М., МГУ. 2006. С. 4041.

50. Лебедев П.Д., Успенский А.А. Исследование геометрии и асимптотики волновых фронтов в некоторых задачах управления // Тр. 9-й между-нар. Четаевской конф. 2007. Т. 5. С. 224-236.

51. Лебедев П.Д., Успенский А.А. Численно-аналитические методы построения волновых фронтов в задачах управления и геометрической оптике // Вестник Тамбовского университета. 2007. Т. 12, вып. 4. С. 538-539.

52. Лебедев П.Д. Вычисление меры невыпуклости плоских множеств // Труды Института математики и механики. 2007. Т. 13, № 3. С. 84-94.

53. Лебедев П.Д., Успенский А.А. Аналитическое и численное конструирование функции оптимального результата для одного класса задач быстродействия // Прикладная математика и информатика. Труды факультета ВМК МГУ. 2007. № 27. С. 65-79.

54. Лебедев П.Д., Успенский А.А. Геометрия и асимптотика волновых фронтов // Известия высших учебных заведений. Математика. 2008. № 3 (550). С. 27-37.

55. Лебедев П.Д., Успенский А.А. Построение обобщенного решения уравнения в частных производных первого порядка // Алгоритмический анализ неустойчивых задач. Тезисы докладов международной конференции, посвященной 100-летию В.К. Иванова. 2008. С. 240-241.

56. Лебедев П.Д., Успенский А.А., Ушаков В.Н. Построение минимаксногоiрешения уравнений типа эйконала // Труды Института математики и механики. 2008. Т. 14, № 2. С. 182-191.

57. Лебедев П.Д., Успенский А.А. Алгоритмы построения функции оптимального результата в задаче с простой динамикой // Вестник Удмуртского университета. Математика. Механика. Компьютерные науки. 2008. Вып. 2. С. 152-54.

58. Лейхтвейс К. Выпуклые множества. М.: Наука, 1985. — 335 с.

59. Лукоянов Н.Ю. Уравнения Гамильтопа-Якоби для наследственных систем: минимаксное и вязкостное решения // Доклады РАН, 2008. Т. 418, № 3. С. 300-303.

60. Лукоянов Н.Ю. Стабильные функционалы в задачах управления с наследственной информацией // Вестник УдГУ. Математика". Механика. Компьютерные науки, Вып. 2, 2008. С.77-80.

61. Мэтъюз Д., Финк К. Численнтые методы. Использование MATLAB (3-е издание). СПб.: Вильяме, 2001. — 720 с.

62. Никольский М.С. Об альтернированном интеграле Л.С. Понтрягина // Мат. сб. 1981. Т. 116, № 1. С. 136-144.

63. Обен Ж.-П., Экланд И. Прикладной нелинейный анализ. М.: Мир, 1988. 512 с.

64. Панасенко Е.А., Тонкое Е.Л. Инвариантные и устойчиво инвариантные множества дифференциальных включений // Тр. МИАН. 2008. Т. 262, С. 202-221.

65. Паиасюк А.П. Уравнения динамики множеств достижимости в задачах оптимизации и управления в условиях неопределенности // Прикл. математика и механика. 1986. Т. 50, вып. 4. С. 596-604.

66. Пахотипских В.Ю., Успенский А.А., Ушаков В.Н. Конструирование стабильных мостов в дифференциальных играх с фазовыми ограничениями // Прикладная математика и механика. 2003. Т. 67, вып. 5. С. 771-783.

67. Понтрягин Л. С. О линейных дифференциальных играх. I // Докл. АН СССР. 1967. Т. 174, № 6. С. 1278-1281.

68. Понтрягин Л. С. О линейных дифференциальных играх. II // Докл. АН СССР, 1967. Т. 175, № 4. С. 764-766.

69. Постои Т., Стюарт И. Теория катастроф и ее приложения. М.: Мир, 1980. 607 с.

70. Пшеничный Б.Н. Структура дифференциальных игр // Докл. АН СССР. 1969. Т. 84, № 2. С. 285-287.

71. Пшеничный Б.Н. Выпуклый анализ и экстремальные задачи. М.: Наука, 1980. — 320 с.

72. Рашевский П.К. Курс дифференциальной геометрии. М.: Бдиториал, УРСС, 2003. 432 с.

73. Рокафеллар Р. Выпуклый анализ. М.: Мир, 1973. — 469 с.

74. Седых В.Д. Соотношения между эйлеровыми характеристиками многообразий особенностей коранга 1 фронта общего положения // Докл. РАН. 2002. - т.383, № 6. - С. 735-739. — Библиогр.: 7 назв.

75. Седых В.Д. Функциональные модули особенностей выпуклых оболочек многообразий коразмерности 1 и 2 // Мат.сб. 1982. Т. 119, № 2. С. 233247.

76. Седых В.Д. Особенности выпуклых оболочек // Сиб. мат. журнал. 1983. Т. 24, № 2. С. 158-175.

77. Седых В.Д. Выпуклые оболочки и преобразование Лежандра // Сиб. Мат.журнал. 1983. Т. 24, № 3. С. 122-134.

78. Седых В.Д. Стабилизация особенностей выпуклых оболочек, Мат.сб. Т. 135, № 4. 1988. С. 514-519.

79. Слюсарев Г.Г. Геометрическая оптика. М.: Издательство Академии наук СССР, 1946. 332 с.

80. Солимено О., Крозинъяни В., Ди Порто П. Дифракция и волноводнос распространение оптического излучения. М.: Мир, 1989. — 662 с.

81. Спесивцев JI.B., Ушаков В.Н. Приближённое построение ядра инвариантности в дифференциальных включениях // Журнал вычислительной математики и математической физики. 2005. Т. 45, № 4. С. 592-602.

82. Субботин А.И. Обобщенные решения уравнений в частных производных первого порядка. Перспективы динамической оптимизации. Москва-Ижевск: Институт компьютерных технологий, 2003. — 336 с.

83. Субботин А.И., Чепцов А.Г. Оптимизация гарантии в задачах управления // М.: Наука, 1981. 288 с.

84. Субботин А.И. Минимаксные решения уравнений с частными производными первого порядка // УМН. 1996. Т. 51, № 2(308). С. 105-138.

85. Суббот;ин А.И. Об одном свойстве субдифференциала // Матем. сб. 1991. № 182 (9). С. 1315 1330.

86. Субботин А.И., Ченцов А.Г. Итерационная процедура построения минимаксных и вязкостных решений уравнений Гамильтона-Якоби и ее обобщения // Труды МИАН. 1999. № 224. С. 311-334.

87. Субботин А.И. Минимаксные решения уравнений Гамильтона-Якоби. // Итоги науки и техн. Сер. Соврем, мат. и ее прил. Темат. обз. 1999. № 64. С. 222-231.

88. Субботина Н.Н. Метод характеристик Коши и обобщенные решения уравнений Гамильтона-Якоби-Беллмана // Доклады АН СССР. 1991. Т 320, № 3 С. 556-561.

89. Тарасьев A.M., Ушаков В.Н., Хрипунов А.П. Об одном вычислительном алгоритме решения игровых задач управления // Прикл. математика и механика. 1987. Т. 51, вып. 2. С. 216-222.

90. Тарасьев A.M., Ушаков В.Н. О построении стабильных мостов в минимаксной игре сближения-уклонения. Свердловск, 1983. 61 с. Деп. в ВИНИТИ. № 2454-83.

91. Успенский А.А. Аналитические методы вычисления меры невыпуклости плоских множеств // ИММ УрО РАН. — Екатеринбург, 2007. — 21 е.: 10 ил. Библиогр.: 9 назв. - Рус. - Деп. в ВИНИТИ 07.02.07, № 104-В2007

92. Успенский А.А., Ушаков В.Н., Фомин А.Н. си-множества и их свойства/ / Ин-т математики и механики УрО РАН. — Екатеринбург, 2004. 62 е.: 38 ил. - Библиогр.: 7 назв. - Рус - Деп. в ВИНИТИ 02.04.04, № 543-В2004

93. Успенский А.А., Лебедев П.Д. Аналитические и численные подходы к построению функции оптимального результата для задачи быстродействия // Международная конференция «Дифференциальные уравнения и топология». Тезисы Докладов. М.: МГУ. 2008. С. 409-410.

94. Ушаков В.Н. К вопросу стабильности в дифференциальных играх // Позиционное упр. с гарантированным результатом / УрО АН СССР.

95. Свердловск, 1988. С. 101-109.

96. Ушаков В.Н. Процедуры построения стабильных мостов в дифференциальных играх.: Дис. д-ра физ.-мат. наук: 01.01.02. Свердловск, 1991. 308 с. / Ин-т математики и механики УрО АН СССР.

97. Ушаков В.Н., Хрипунов А.П. О приближенном построении решений вигровых задачах управления // Прикл. математика и механика. 1997.1. Т. 61, № 3. С. 413-421.

98. Ушаков В.Н., Латушкин Я.А., Лебедев П.Д. Критерии совпадения максимальных стабильных мостов в двух игровых задачах о сближении // Международная конференция «Дифференциальные уравнения и топология». Тезисы Докладов. М.: МГУ. 2008. С. 411-412.

99. Ушаков В.Н., Латушкин Я.А. Критерии совпадения максимальных стабильных мостов в двух игровых задачах о сближении / / Тр. МИ АН им. Стеклова. 2008. Т. 262. С. 253-271.

100. Чен К., Джиблин П., Ирвинг A. MATLAB в математических исследованиях. М.: Мир, 2001. — 346 с.

101. Ченцов А.Г. Об игровой задаче сближения в заданный момент времени // Мат. сб. 1976. Т. 99, № 3. С. 394-420.

102. Ченцов А.Г. Об структуре игровой задачи сближения в заданный момент времени // Доклады АН СССР. 1975. Т. 16. С. 1404-1406.

103. Экланд И., Темам Р. Выпуклый анализ и вариационные проблемы. М.: Мир, 1979. 399 с.

104. Arnol'd V.I. Wave front evolution and equivariant. Morse lemma. Comm. Pure Appl. Math., no 29 (6), 1976. Pp. 557-582.

105. Bardi M. and Osher S. The nonconvex multi-dimensionsl Riemann problem for Hamilton-Jacobi equations // SIAM. J. Math. Anal. 1991. Vol. 22, No 2. Pp. 344-351.

106. Bardi M. and Falcone M. An Approximation Scheme for the minimax Time Function // SIAM J. Control Optimiz. Vol. 28. No 4, 1994. Pp. 950965.

107. Bruce J.W., Giblin P.J. and Gibson C.G. Symmetry sets. Proc. Roy. Soc. Edinburgh Sect. A, 101. 1983. Pp. 163-186.

108. Crandall M.G. and Lions P.-L. Viscosity solutions of Hamilton-Jacobi equations // Trans. Amer. Math. Soc. Vol. 277, No 1, 1983. Pp. 1-42.

109. Crandall M.G. and Lions P.-L. Two approximations of solutions of Hamilton-Jacobi equations // Math. Comput. Vol. 277, No 1,1984. P. 1-19.

110. Giblin P.J. Symmetry sets and medial axes in two and three dimensions. The Mathematics of Surfaces IX, (eds. Roberto Cipolla and Ralph Martin), Springer-Verlag. 2000. Pp. 306-321.

111. Giblin P. J., Kimia B.B. The local form and Transitions of Symmetry Sets, Medial Axes and Shocks. Int. J. Computer Vision, 54. 2003. Pp. 143-157.

112. Guseinov H.G., Subbotin A.I. and Ushakov V.N. Derivatives for Multivalued Mappings with Applications to Game-Theoretical Problems of Control // Problems Control Inform. Theory. 1985. Vol 14, no 6. Pp. 405419.

113. Hopf E. The partial differential equation ut + uux = fiuxx// Comm. Pure and Appl. Math., Vol. 3, 1950. P. 201—230.

114. Hopf E. Generalized solutions of non-linear equations of first order// J. Math, and Mech, Vol. 14, No 6, 1965. Pp. 951-973.

115. Kamneva L. V. Discontinous Value Function in Time-Optimal Differential Game// International Conference "Differental Equatijns ant Topology". MSU, Moscow. 2008. P. 259.

116. Krasovskii A. A. and Tarasyev A.M. Dynamic Optimization of Investments in the Economic Growth Models. Automation and Remote Control, Vol. 68, No. 10, P. 1765-1777, 2007.

117. Krasovskii A.A. and Tarasyev A.M. Conjugation of Hamiltonian Systems in Optimal Control Problems. Proceedings of the 17th World Congress, The Interna tional Federation of Automatic Control, P. 7784-7789, 2008.

118. Lax P.D. Weak solutions of nonlinear hyperbolic equations and their numerical computation // Comm. Pure Appl. Math. Vol. 7, No 1, 1954. P. 159-193.

119. Melik'jan A.A. Method characteristics for constructing of singular paths and manifolds in optimal control and differentials games // Differential Games — Developements in Modelling and Computation. 1991. Vol. 156. P. 81-90.

120. Melikyan A.A. Generalized Characteristics of the First Order PDEs: Applications in Optimal Control and Differential Games. Boston: Birkhauser. 1998. — 310 p.

121. Patsko VS., Botkin N.D., Kein V.M., Turova V.L. and Zarkh M.A. Control of Aircraft Landing in Windshear// Journal Opt. Theory Appl., No 83(2), 1994. P. 237-267.

122. Sedykh V.D. On the topology of symmetry sets of smooth submanifolds in Rk// Advanced Studies in Pure Mathematics 43, 2006. Singularity Theory and Its Applications, P. 401-419.

123. Sedykh V.D. Some invariants of convex manifolds, In Workshop, on Real and Complex Singularities (Sao Carlos, 1992), Mat. Contemp. 5. 1993. P. 187-198.

124. Sethian J.A. Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid Me. Cambridge University Press, 1999. 400 p.

125. Sethian J.A. Fronts Propagating with Curvature Dependent Speed: Algorithms Based on Hamilton-Jacobi Formulations / / Journal of Computational Physics, Vol 79, 1988. P. 12-49,

126. Souganidis P.E. Approximation schemes for viscosity solutions of Hamilton-Jacobi equations //J. of Different. Equat. Vol. 59, 1985. P. 143.

127. Subbotin A.I. Generalized Solutions of First-Order PDEs. The Dynamical Optimization Perspective. Boston: Birkhauser. 1995. 312 P. (System & Control: Foundation & Appl.)

128. Subbotin A.I. and Tarasyev A.M. Stability properties of the value function of a differential game and viscosity solutions of Hamilton-Jacobi equations // Problems Control Inform. Theory. 1986. 15 (6). P. 451-463.