История кафедры исследования операций |
Содержание:
Становление кафедры (И. В. Романовский)
Предыстория и моя работа на кафедре (А. М. Вершик)
Быстрые алгоритмы дискретной оптимизации (И. В. Романовский)
Дискретный гармонический анализ (В. Н. Малозёмов)
Гарантированный поиск на графах (Т. В. Абрамовская)
Дифференциальный игры (Н. Н. Петров)
Становление кафедры (И. В. Романовский)
|
Предыстория и моя работа на кафедре
|
I. ПредысторияВ начале 60-х гг. Леонид Витальевич Канторович, переехавший к тому времени из Ленинграда в Новосибирск, продолжал активные действия по внедрению преподавания математической экономики и линейного программирования в программы обучения студентов ЛГУ и других вузов. Был организованы 6-й курс и, позже, кафедра кибернетики на экономическом факультете ЛГУ, а также чтение первых курсов на мат-мехе. В то время ещё не было ни достаточного числа нужных специалистов, ни готовых программ.В 1962 г. на кафедру вычислительной математики мат-меха, которую организовал в начале 50-х гг. Л.В. и которую возглавлял его ученик М.К.Гавурин было зачислено несколько молодых математиков, знакомых с предметом, для преподавания этого цикла математических курсов, в том числе я. Так случилось, что хотя я закончил университет по кафедре математического анализа и писал диссертацию по функциональному анализу, мне пришлось подробно познакомиться и учиться этим предметам у одного из моих учителей Л.В и у его ближайшего сотрудника Г.Ш.Рубинштейна. Чтение курсов и ведение семинаров по линейному программированию, комбинаторике, программированию и т.д. ставилось в задачи этих новых членов кафедры вычислительной математики. За несколько лет (до 1970 г.) по специальности математическая экономика мат-мех закончило несколько десятков выпускников. Формирование специальной кафедры ускорило постороннее обстоятельство — а именно внезапное решение властей образовать в ЛГУ прикладной факультет (примат), привлекая для этого часть сотрудников мат-меха. Профессора и преподаватели мат-меха неодобрительно встретили это решение, но он было проведено в жизнь, и лишь немногие преподаватели согласились перейти на новый факультет. Одновременно стал вопрос о том, что мат-мех должен сохранить свои прикладные направления и, поэтому, в частности, были созданы новые кафедры — кибернетики (зав. В.А.Якубович) и Исследования Операций (зав. М.А.Гавурин). Название последней не вполне удачно; выбор его следовал западной традиции объединять вместе ряд очень разных направлений, нацеленных, как сказали бы сейчас, на менеджмент и приложения. На вновь образованные кафедры мат-меха были переведены математики разных кафедр и, конечно, те, кто работал на вычислительной кафедре по соответствующим специальностям, в том числе и я, стали членами кафедры ИО. II. Моя работа на кафедре1. Разработка и чтение в течение почти 20 лет (1974-1992) основного курса кафедры «ЭКСТРЕМАЛЬНЫЕ ЗАДАЧИ» для студентов-математиков третьего курса. Курс совершенно оригинальный. Он включал основы вариационного исчисления, линейного и оптимального программирования, некоторых разделов комбинаторики и заканчивался изложением симплекс-метода. Курс более всего опирался на функциональный анализ, на работы Л.В.Канторовича и некоторые идеи автора. Он совершенствовался по мере его чтения. Через определённое время я организовал запись курса, её осуществляли последовательно Т.Кулаковская, А.Грибов, А.Барвинок. Впоследствии, А. Барвинок выпустил уже в США учебник, основанный в значительной мере на моих лекциях и их записей.2. Цикл работ с моим аспирантом П.Спорышевым о статистическом подходе к симплекс-методу. Основываясь на геометрическом (грасмановом) описании общего выпуклого многогранника; используя естественную (гауссову) меру, удалось оценить среднее число шагов симплекс-метода в типичной задаче линейного программирования. Работа получила распространение и известность, так как была одной из первых работ в сделавшимся сейчас популярном и важном изучении асимптотики числа граней разных размерностей в случайном многограннике. (начало 80-х гг.). Публикации:
3. Цикл работ (моих, и с аспирантом В.Гершковичем, и аспиранткой О.Граничиной) по неголономной геометрии и по неклассическому вариационному исчислению, неголономным потокам, виду геодезических и др. (начало 80-х гг.).
4. Работы с А.Барвинком по алгебро-геометрическому подходу к классификации массовых задач оптимизационного характера и о NP-полноте. В частности, обобщении задачи о назначении на многогранники алгебраического происхождения. (середина 80-х гг.)
5. Работы по бесконечномерному программирования и аппроксимации. (конец 60-х гг.)
6. Исследования метрики Канторовича и ее применения (конец 80 - начало 90-х гг.)
7. Работы на близкие темы.
III. Аспиранты кафедры, защитившие диссертации под моим руководствомВ.Темельт (ГДР) 1964 г. П.В.Спорышев (1984), В.Г. Гершкович (1984), А.И.Барвинок (1986), О.А.Граничина (1988), У.Бальде (Гвинея) (1991). Дипломантов около сотни. |
Быстрые алгоритмы дискретной оптимизации
|
Начиная с середины прошлого века, в математике и ее приложениях неуклонно растет роль дискретных моделей.
В этих моделях важную роль играют задачи оптимизации.
Заметный вклад в исследование этих вопросов сделали ученые нашего университета.
Отметим цикл задач раскроя и укладки, ведущий свое начало, в частности, от ленинградских работ Л.В.Канторовича и В.А.Залгаллера, известных математиков, которые в 1949-1950 гг. развили и реализовали на Вагоностроительном заводе в Ленинграде методы вычисления оптимальных раскроев листового материала.
Они не использовали тогда вычислительных машин.
Одни из первых программ, реализующих такие методы были созданы в ЛГУ, в лаборатории исследования операций НИИММ.
Сейчас в России есть несколько научных центров, в которых эта работа продолжается (в частности, в Уфе и Петрозаводске).
Отметим, что программа управления гофрокартонным оборудованием, разработанная Петрозаводским госуниверситетом, оказалась чрезвычайно эффективной и сейчас установлена на двух десятках крупных отечественных производств.
В связи с празднованием 100-летия со дня рождения Л.В.Канторовича профессор И.В.Романовский и ведущий научный сотрудник СПбЭМИ РАН А.А.Корбут (1934–2013) подготовили к печати 3-е издание книги Л.В.Канторовича и В.А.Залгаллера «Рациональный раскрой промышленных материалов». В это издание включен большой обзор современного состояния задачи оптимального раскроя и родственной проблематики, написанный Э.А.Мухачевой (1930–2011) и ее сотрудниками из уфимской школы. Естественно, что в обзоре отражена и роль сотрудников и аспирантов нашего университета. Задачи дискретной оптимизации характеризуются большими объёмами данных, поэтому для их решения разрабатываются быстрые алгоритмы в самых различных областях науки и приложений. Систематически понятие быстрого алгоритма стало использоваться в связи с быстрым преобразованием Фурье (БПФ). Открытие БПФ привело к созданию новой математической дисциплины «дискретного гармонического анализа», ориентированного, прежде всего, на задачи цифровой обработки сигналов. В этой области на кафедре исследования операций имеются серьёзные достижения: разработаны параметрические варианты БПФ, построены на базе БПФ вейвлет-пакеты. Опубликована книга: Малоземов В. Н., Машарский С. М. Основы дискретного гармонического анализа. СПб.: Изд-во «Лань», 2012. 304 с. Дискретные периодические сплайны нашли применение в геометрическом моделировании при построении замкнутых кривых и поверхностей по ограниченной информации. В частности, при построении периодических кривых Безье явно используется быстрое преобразование Фурье. Важным прикладным направлением является теория расписаний. Задачи теории расписаний имеют много практических приложений в системах обслуживания и оптимизации транспортного движения, планировании производства и организации работы многопроцессорных компьютерных комплексов. Практические задачи составления расписаний имеют большую размерность и ни один из известных в теории подходов к решению таких дискретных задач не гарантируют получения точного решения за приемлемое время. Разработка быстрых приближенных алгоритмов, которые позволяют получать достаточно близкие к оптимальным решениям результаты за приемлемое время, представляет большой теоретический и практический интерес. На кафедре исследования операций накоплен значительный опыт в разработке приближенных подходов к решению задач составления расписаний. Это генетические алгоритмы, методы локального поиска, метод, основанный на алгоритме муравьиной колонии, и гибридные методы, которые позволяют сочетать лучшие возможности разных подходов и учитывать специфику конкретной задачи. На проблематику наших исследований в этой области повлияла прикладная работа, которую в 2006 г. кафедра выполняла по договору СПбГУ с фирмой BHP Billton, Австралия. Нами было составлено несколько быстро работающих программ нахождения формы разрабатываемого открытым способом рудного месторождения, основанных на алгоритме построения максимального потока в сети А. Голдберга и Б. Черкасского, но была ориентирована на вариантные расчеты. В качестве возможного продолжения сотрудники фирмы предложили задачу составления укрупненного расписания разработки месторождения с интересной и мало встречавшейся ранее целевой функцией — взвешенным доходом от добытого рудного сырья, которой они занимаются и сами. Над этой очень трудной задачей мы продолжаем думать, рассматривая различные подготовительные модели. В 1998 г. в Петрозаводском университете по материалам лекций И.В.Романовского была издана книга «Субоптимальные решения». В ней излагался новый подход к приближенному решению задач дискретной оптимизации. С тех пор этот метод несколько раз применялся, но распространения не получил. Предполагается заняться оценкой трудоемкости этих методов (для этого сейчас в литературе есть достаточная методическая основа), а также рассмотреть вопрос о реализации таких методов на многопроцессорной вычислительной технике. |
Дискретный гармонический анализ
|
1. С конца 1980-х годов на волне широкого интереса к практическим задачам цифровой обработки сигналов (с ними нас познакомил проф. А. А. Ланнэ) на кафедре начались исследования в области, которая теперь называется «дискретный гармонический анализ». Основным объектом анализа являются дискретные периодические сигналы – комплекснозначные N-периодические функции целочисленного аргумента. Основным действием – разложение сигнала по ортогональным базисам. Основные базисы – система сдвигов единичного N-периодического импульса и экспоненциальный базис. Вектор коэффициентов разложения сигнала по экспоненциальному базису называется дискретным спектром Фурье. Цифровая обработка сигнала сводится к преобразованию его спектра. Сам преобразованный сигнал восстанавливается с помощью формулы обращения. 2. Первые результаты в области дискретного гармонического анализа были получены М. Г. Бером и В. Н. Малоземовым [1, 2]. Они разработали методы оптимального в некотором смысле восстановления дискретных периодических сигналов по их значениям (заданным либо точно, либо с некоторой погрешностью) на крупной сетке. В 1997 году В. Н. Малоземов м А. А. Третьяков предложили новый подход к быстрому вычислению дискретного преобразования Фурье в случае, когда длина периода N равна степени двойки [3]. Подход основан на построении рекуррентной последовательности ортогональных базисов в пространстве дискретных периодических сигналов. В качестве начального базиса берется система сдвигов единичного N-периодического импульса, финальным является экспоненциальный базис. Позднее выяснилось, что все промежуточные базисы имеют блочную структуру. Сигналы в каждом блоке различаются лишь сдвигом аргумента. Это позволило сформировать новые ортогональные базисы вейвлетного типа, в которые входят по одному блоку из каждого промежуточного базиса рекуррентной последовательности. В частности, так можно получить дискретный базис Хаара Отмеченный факт приводит к неожиданному следствию: в процессе вычисления дискретного спектра Фурье автоматически вычисляется и дискретный спектр Хаара [4]. В рекуррентной последовательности ортогональных базисов, упомянутой выше, каждый сигнал k-го базиса линейно зависит от двух сигналов (k – 1)-го базиса. В работе [5] было замечено, что при замене всех коэффициентов этих линейных зависимостей на единицу получается новая последовательность ортогональных базисов, причем финальным будет дискретный базис Уолша. В более общем виде вопрос о построении рекуррентных последовательностей ортогональных базисов с помощью замены коэффициентов в рекуррентных соотношениях изучался в работе [6]. 3. Дискретный гармонический анализ существенно зависит от арифметических свойств длины периода N, точнее, от его разложений на множители. В работах [7, 8] результаты, полученные в [3, 4, 5] при N, равном степени двойки, обобщаются на случай, когда N равно степени произвольного натурального числа n > 2. Соответствующие обобщения оказались возможными благодаря тому, что удалось разобраться с источником рекуррентных соотношений для базисов. Если раньше такие соотношения, по существу, предъявлялись, то теперь они выводятся на основе различных факторизаций матрицы Фурье. Эффективные расчетные формулы получаются путем использования индексной техники, когда при умножении разреженной матрицы на вектор убираются все операции с нулевыми элементами матрицы. 4. В работах [9, 10] В. Н. Малоземов и А. Б. Певный ввели дискретные периодические сплайны и изучали их фундаментальные свойства. Базисный сплайн (B-сплайн) определен с помощью циклической свертки (он принимает только целые неотрицательные значения). Общий дискретный сплайн представляет собой линейную комбинацию с комплексными коэффициентами сдвигов B-сплайна. Общий дискретный сплайн можно определить также как линейную комбинацию сдвигов дискретной функции Бернулли соответствующего порядка. Последнее определение позволяет доказать свойство минимальной нормы для дискретного интерполяционного сплайна. С помощью дискретных периодических сплайнов получено также эффектное решение задачи об оптимальном сглаживании дискретных периодических данных. В работе [11] построен ортогональный базис в пространстве дискретных периодических сплайнов и получены первые результаты в области дискретного сплайн-гармонического анализа. 5. С 2004 по 2014 годы пол руководством В. Н. Малоземова работал семинар по дискретному гармоническому анализу и геометрическому моделированию. Сайт семинара: http://www.dha.spb.ru. На сайте, в частности, в качестве препринтов публиковались тексты докладов участников семинара. Эти доклады стали основой для двухтомника [12, 13], опубликованного в 2014 году и содержащего 1188 страниц. Исследования на семинаре велись широким фронтом. Выделим четыре основных цикла работ и характерные публикации в каждом цикле:
В 2012 г. в издательстве «Лань» вышла из печати книга «Основы дискретного гармонического анализа» [19]. Её расширенный английский вариант был опубликован в 2020 г. [20].
|
Гарантированный поиск на графах
|
Одним из важных направлений исследований на кафедре в течение всего времени её существования являлась так называемая теория гарантированного поиска на графах.
Типичная задача гарантированного поиска на графах, на неформальном уровне, может быть сформулирована следующим образом. На графе располагаются две группы игроков. Первая группа — назовём её участников преследователями — ставит своей задачей «переловить» всех участников второй группы, которых будем называть убегающими. Впрочем, из дальнейшего станет ясно, что поставленную задачу достаточно решить для одного убегающего, что и будет предполагаться в нижеследующих рассуждениях. Цель команды преследователей — построить «программу» действий, которая обеспечивала бы им поимку «невидимого» убегающего при любом его поведении. Иными словами, поимка должна быть гарантирована даже в том случае, когда выбранная преследователями программа становится известной убегающему «до начала поиска». В каждой формализации задачи поиска фиксируются динамические возможности участников, условие успешного завершения поиска (условие поимки), арена поиска (комбинаторный или топологический граф), и ставится следующая задача: найти наименьшее число преследователей, необходимое для успешного завершения поиска. Это число называют поисковым числом графа с соответствующим прилагательным, которое зависит от принятой формализации. Столь естественная задача вполне могла быть поставлена в XVIII веке, во времена Эйлера, но первые содержательные результаты появились лишь в 70 – 80 годы прошлого века в рамках теории конфликтного управления с неполной информацией. В монографии Р. Айзекса [1], содержащей основы теории дифференциальных игр, рассматривается одна (неформально поставленная) задача преследования–уклонения — проблема «Принцесса и Чудовище», в которой участники не получают никакой информации о поведении друг друга. Пионером теории гарантированного поиска на кафедре исследования операций (и, по всей видимости, в Советском Союзе) был Николай Николаевич Петров. Будучи специалистом в области дифференциальных игр, он ставит задачу [2], созвучную с проблемой «Принцессы и Чудовища», но не допускающую смешанных стратегий. В основу поведения игроков положен «минимаксный принцип», траектории игроков - непрерывны, конфликт разворачивается на одномерном CW–комплексе. В то же время и независимо аналогичная задача ставится американским математиком Т. Парсонсом [3] (также на топологическом графе). В дальнейшем ученик Н. Н. Петрова П. А. Головач доказал, что обе формализации эквивалентны, так что в обоих случаях речь идёт об одном и том же поисковом числе. Более того, П. А. Головач свёл проблему нахождения поискового числа к некоторой дискретной задаче на соответствующем комбинаторном графе. Под руководстовм Н. Н. Петрова задачами гарантированного поиска на графах занимались выпускники кафедры исследования операций Т. В. Абрамовская, П. А. Головач, А. Б. Зеленевская, В. О. Капилевич, С. В. Лунегов, С. А. Старостина, М. А. Тетерятникова, И. Туре, Ф. В. Фомин и А. В. Чуманова. Выпускница кафедры Е. Э. Ржевская проводила исследования в этой области под руководством Т. В. Абрамовской. В 2010 году, к сорокалетию кафедры исследования операций, был подготовлен обзор результатов [4], полученных выпускниками кафедры на эту тему. Список литературы в упомянутой обзорной статье содержит 100 наименований, сорок из которых - за авторством перечисленных сотрудников и выпускников кафедры исследования операций.
|
Дифференциальный игры
|
В теории дифференциальных игр Николаем Николаевичем Петровым были получены следующие результаты:
|