WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!


На правах рукописи

Гасникова Евгения Владимировна

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

Специальность 05.13.18 - математическое моделирование, численные методы и комплексы программ

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

Москва – 2012

Работа выполнена на кафедре анализа систем и решений Московского физико-технического института (государственного университета)

Научный консультант: кандидат физико-математических наук, доцент Меньшиков Иван Станиславович

Официальные оппоненты: Доктор физико-математических наук, профессор Бекларян Лёва Андреевич, Центральный экономико-математический институт РАН, главный научный сотрудник кандидат физико-математических наук Швецов Владимир Иванович, Институт системного анализа РАН, ведущий научный сотрудник

Ведущая организация: Институт прикладной математики им. М.В. Келдыша РАН

Защита состоится декабря 2012 года в часов на заседании диссертационного совета Д 212.156.05 при Московском физико-техническом институте (государственном университете) по адресу: 141700, г. Долгопрудный Московской обл., Институтский пер., д. 9, ауд. 903 КПМ.

С диссертацией можно ознакомиться в библиотеке Московского физикотехнического института (государственного университета).

Автореферат разослан ноября 2012 г.

Ученый секретарь диссертационного совета Федько О.С.

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

Актуальность темы. Работы Е.Т. Джейнса, А. Дж. Вильсона, И. Пригожина, Г. Хакена положили начало традиции использования термодинамического формализма для анализа различных макросистем, встречающихся в экономике, биологии, социальной области и т.д. В России систематические исследования в этом направлении велись Л.Н. Розоноэром, а в последнее время В.П. Масловым, Ю.С. Попковым и другими исследователями. В большинстве работ рассматривались случаи стационарных макросистем, которые уже находятся в равновесии.

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

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

Большой класс макросистем основывается на эволюционной динамике, приводящей к равновесию, которое может быть проинтерпретировано как равновесие Нэша. Возможность предложить и исследовать разумную динамику для той или иной стационарной макросистемы в ряде случаев позволяет также эффективно определить равновесие этой макросистемы.

Целью диссертационной работы является исследование равновесия макросистем, обладающих марковской динамикой, в основе которой лежат принципы стохастической химической кинетики.

Для этого решаются следующие задачи:

отыскание достаточных условий существования и единственности равновесия;

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

изучение наследования свойств динамической макросистемы при различных скейлингах (изменениях масштаба), в частности, исследование перестановочности предельных переходов по времени и размерности макросистемы;

изучение теоретико-игровых аспектов концепции равновесия динамических макросистем;

оценки скорости сходимости к равновесию и плотности концентрации инвариантной меры в равновесии;

анализ конкретных макросистем, рассматриваемых при математическом моделировании транспортных потоков; при ранжировании web-страниц; в экономических и социальных моделях.

Методы исследования. В работе используются: эргодические теоремы для марковских процессов, формула Стирлинга, метод Дарвина–Фаулера, современные варианты неравенства Чигера, понятие дискретной кривизны Риччи, второй метод Ляпунова, термодинамический формализм, аппарат производящих функций, метод внутренних штрафных функций, метод стохастического квазиградиентного спуска.

Научная новизна выносимых на защиту результатов состоит в том, что впервые:

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

выявлен характер связи функции Ляпунова макросистемы с инвариантной мерой этой макросистемы;

показано, что при достаточно естественных условиях время сходимости к равновесию оказывается малым;

приведена эволюционная интерпретация конструкции равновесия Нэша– Вардропа, отвечающая равновесному распределению потоков на графе транспортной сети, до этого использовавшаяся только стационарно. Это позволило, в том числе, усилить парадокс Браесса (1968), а также объяснить эвристический способ отбора Бар-Гира–Швецова (2010) единственного равновесия Нэша–Вардропа в случае, когда имеет место неединственность равновесий.

Практическая ценность работы. Диссертационная работа, в основном, носит теоретический характер, хотя во многом мотивированна конкретными приложениями.

В диссертации в рамках макросистемного подхода предложена модель, обобщающая известную на практике энтропийную модель расчета матриц корреспонденций. Эта известная модель была использована для анализа данных по транспортным потокам г. Москвы. Предложенная в диссертации модель позволила более адекватно описать имеющиеся данные. Сейчас предложенная модель тестируется в Научно-исследовательском и проектном институте Генерального плана города Москвы, как структурный блок общей четырехстадийной транспортной модели.

Диссертация выполнялась при поддержке Лаборатории структурных методов анализа данных в предсказательном моделировании, МФТИ, грант правительства РФ договор 11G34.31.0073, а также при поддержке грантов РФФИ 11-01-00494-а, мол_а_вед 12-01-33007.

Апробация и публикации. По материалам диссертации опубликовано печатных работ, в том числе, три – в изданиях, рекомендованных ВАК РФ [6, 12, 13].

Личный вклад соискателя в работы с соавторами состоит в доказательстве теорем о сходимости динамических макросистем к равновесию; исследовании ряда свойств равновесий макросистем (скорость сходимости, плотность концентрации инвариантной меры); исследовании макросистем, имеющих практическую ценность; систематизации ранее известных алгоритмов; разработке численных алгоритмов; разработке соответствующего комплекса программ.

Результаты работы докладывались, обсуждались и получили одобрение специалистов на 16 конференциях (в том числе, пяти международных), научных семинарах по экспериментальной экономике Вычислительного центра РАН (20092012), семинаре-конференции РАН под рук. акад. В.В. Козлова (2010), на семинарах кафедр математических основ управления и анализа систем и решений ФУПМ МФТИ (2009-2011), на семинарах по транспортным потокам в Независимом Московском Университете (2011-2012).

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

КРАТКОЕ СОДЕРЖАНИЕ ДИССЕРТАЦИИ Во введении сформулированы цели работы, обосновывается её актуальность, дан обзор работ, связанных с темой диссертации, приводятся основные результаты и их предпосылки.

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

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

Задача 1. Обобщение энтропийной модели А.Дж. Вильсона расчета матрицы корреспонденций.

Пусть в некотором городе имеется n районов, Li 0 – число жителей i -го nn района, Wj 0 – число работающих в j -м районе. При этом M j L W, – i i1 jобщее число жителей города. Обозначим через xij t 0 – число жителей, живущих в i-м районе и работающих в j-м в момент времени t. Со временем жители могут только меняться квартирами, поэтому во все моменты времени t n n xij t 0, x t Li, x t Wj, i, j 1,..., n, 1 n2 M. (A) ij ij j1 iОпишем имеющиеся стимулы к обмену. С одной стороны, работать далеко от дома плохо из-за транспортных издержек, но с другой, люди все-таки ездят на работу далеко. Это можно объяснить тем, что близко от дома не удается найти хорошую работу. Будем считать, что компромисс между этими двумя неприятностями описывается эффективной функцией затрат R l l lnl 2, где l 0 – расстояние от дома до работы, а 0, 0, 0 – настраиваемые параметры модели.

Теперь опишем саму динамику. Пусть в момент времени t 0 r-й житель живет в k-м районе и работает в m-м, а s-й житель живет в p-м районе и работает в q-м.

Тогда k,m; p,q t t o t – есть вероятность того, что жители с номерами r и s ( 1 r s M ) поменяются квартирами в промежутке времени t,t t.

Вероятность обмена местами жительства зависит только от мест проживания и работы обменивающихся:

k,m; p,q t k,m; p,q M exp R lkm R lpq R lpm R lkq 0, суммарные затраты суммарные затраты до обмена после обмена где коэффициент 0 характеризует интенсивность обменов.

Согласно эргодической теореме для марковских цепей:

n, n xij i1, j1 A lim P xij t xij, i, j 1,...,n t n, n def 1 n, n Z exp lnlij lij xij xij ! p xij i1, j1, i1, jгде статсумма Z находится из условия нормировки получившейся “пуассоновской” вероятностной меры. Отметим, что стационарное распределение n, n p xij i1, j1 удовлетворяет условию детального равновесия:

xkm 1 xpq 1 p x11,..., xkm 1,..., xpq 1,..., xpm 1,..., xkq 1,..., xnn k,m; p,q n, n xpmxkq p xij i1, j1 p,m; k,q.

n, n При M 1 распределение p xij i1, j1 сконцентрировано на множестве (A) в n, n * окрестности наиболее вероятного значения xij i1, j1, которое находится, как решение задачи энтропийно-линейного программирования (см. теорему 1 ниже):

n, n n, n n, n ln p xij i1, j1 xij ln xij e lij lnlij xij min.

n, n xij i1, j1 A i1, j1 i1, jРешение этой задачи можно представить как xij exp iL W lnlij lij, j n n где множители Лагранжа (двойственные переменные) iL i1 и W j j определяются из равенств системы (A). На практике мы имеем информацию о n, n n Li,Wi i1 и lij i1, j1. Из системы ограничений (A), мы найдем n, n n, n n xkm Li,Wi i1; lij i1, j1.

k1,mТакой способ расчета матрицы корреспонденций с 0 в литературе часто называют (энтропийной) гравитационной моделью:

xij AiBjLWj lij exp lij, i n n где Ai i1 и Bj j1 определяются из соотношений:

n n Ai B Wj lij exp lij , Bj A Li lij exp lij .

j i j1 i1 В классической статической энтропийной модели (второй половины 60-ых годов XX века) А.Дж. Вильсона 0. Однако по Москве наблюдается заметно лучшее соответствие реальных данных рассматриваемой модели, если допускать 0, и оптимально его подбирать.

Приведем теперь оценку концентрации стационарного распределения n, n n, n * p xij i1, j1 в окрестности наиболее вероятного значения xij i1, j1. Для этого, * * * предварительно, введем такую перестановку ij , что xij1 xij2 ... xij, а n затем введем обозначение nk t xijk t, k 1,...,n2.

Утверждение. q 0,..., n2 q 0, Tq,n M ln M : t Tq,n qn nk t P 1 , k 1,...,q 0.99.

* xijk M Рассмотрим теперь общий подход, частным случаем которого является изложенная выше модель.

Предположим, что некоторая макросистема может находиться в различных состояниях, характеризуемых вектором n с неотрицательными целочисленными компонентами. Будем считать, что в системе происходят случайные превращения (химические реакции). Пусть n n , , J – все возможные типы реакций (конечный набор). Введем, следуя М.А. Леонтовичу (1935), интенсивность реакции:

1 i i , n , n n M K ni ... ni i 1, i:i где K 0 – константа реакции (в общем случае K n – функция, см. например, теорему 2). При этом часто считают ni t M. Таким образом , n – i вероятность осуществления в единицу времени перехода n n . На макроуровне все это соответствует принципам химической кинетики (закон действующих масс Гульдберга–Вааге, 1864).

Определение. Условием унитарности (условием Штюкельберга–Батищевой– Пирогова) называется следующее соотношение:

jj 0: K K . (ШБП) j j jj : , J : , J Следующая теорема уточняет результаты а) В.В. Веденяпина (1999), б) В.А.

Малышева, С.А. Пирогова, А.Н. Рыбко (2004).

Теорема 1. а) ,n t ,n 0 (inv) тогда и только тогда, когда вектор ортогонален каждому вектору семейства .

, J б) Пусть выполняется условие (ШБП). Тогда * i i 1. мера n in e ni!, где i i*M, а – произвольное решение i (ШБП), будет инвариантной относительно предложенной стохастической марковской динамики.

2. на множестве (inv) эта мера экспоненциально быстро концентрируется, с ростом M, в окрестности наиболее вероятного состояния, которое и принимается за положение равновесия макросистемы.

3. задача поиска наиболее вероятного макросостояния асимптотически эквивалентна задаче максимизации энтропийного функционала n (воспользовались формулой Стирлинга n! 2n n e 1 1 ):

ln (n) ln ni i 1 на множестве, задаваемом условием (inv).

n i i Теорема 2 обобщает результаты В. Вайдлиха (2000) и др. на случай, когда рассматривается более общая схема, чем модель миграции населения.

Теорема 2. Пусть , J ; i,i 0,1, , i i ii K n exp ui ni 1 ui ni i ni 0.

, u i: i i:i Тогда справедливо утверждения п. б) теоремы 1 с ni n exp Ui ni ni!, Ui ni 2 , n ni ln ni Ui ni.

u i i i 1 i Путь, по которому мы до сих пор приходили к равновесию макросистемы, состоял из следующей последовательности предельных переходов: lim и lim.

t M Однако для качественного понимания того как эволюционирует макросистема бывает полезно осуществить последовательность переходов в обратном порядке.

Пусть теперь dim n и число реакций J не зависят от числа агентов M, и в начальный момент времени t 0 для любого i существует предел ci 0 lim ni 0 M 0. Тогда в произвольный момент времени t 0 для M п.н.

любого i существует не случайный предел ci t lim ni t M. Описанный выше M предельный переход называется каноническим скейлингом.

В результате такого скейлинга приходим к «динамике квазисредних» (терминология В. Вайдлиха):

dci j, c i i K c . (ДК) cj dt j , J Это технически нетривиальное утверждение следует из результатов ТроттераКуртца (1986).

Систему (ДК) можно получить и по-другому. А именно, как приближенную динамику средних ci t E ni t M . Приближенную в том смысле, что при выводе (ДК) используется приближение: F ci t E ni t M для F «достаточно хороших» функций F (например, полиномов). Это верно в случае пикообразного распределения ni t.

Можно показать (следуя Батищевой–Веденяпину, 2001), что если выполняются условия (ШБП), то траектория (ДК) сходится к неподвижной точке. Какой именно, зависит, вообще говоря, от «точки старта»; но можно сказать и точнее: к той единственной неподвижной точке из семейства неподвижных точек, которая принадлежит аффинному многообразию (inv), инвариантному относительно (ДК).

Для этого, вводится (минус) энтропия:

H c ln ci i 1, c i i и показывается, что она является функцией Ляпунова для системы (ДК). Обратим внимание, что инвариантная мера n (при каноническом скейлинге) “породила” функцию Ляпунова H c. Подобные закономерности наблюдаются для рассматриваемых моделей и без условия (ШБП).

Теорема 3. Пусть инвариантная вероятностная мера представляется в виде:

n M exp M H n M 1, M 1, где H c строго вогнутая функция.

Тогда H c – функция Ляпунова системы (ДК).

Естественно теперь задаться вопросом: А что будет, если условия (ШБП) не выполняются, однако система (ДК) имеет на внутренности пересечения неотрицательного ортанта и инвариантного аффинного многообразия (inv) единственную неподвижную точку (пример такой макросистемы имеется в докторской диссертации С.А. Пирогова)? Оказывается, имеет место Теорема 4. Если эта точка экспоненциально глобально устойчива, то 1) все инварианты (законы сохранения) системы (ДК) определяются (inv); 2) около положения равновесия инвариантная мера будет экспоненциально быстро концентрироваться (с ростом M ); 3) скорость сходимости к равновесию (mixing time) оценивается как M ln M ; 4) элементы корреляционной матрицы случайного вектора n t равномерно ограничены по времени; 5) предельные переходы lim и lim перестановочны: lim lim limM lim.

M t M t t Результаты главы 1 и ряда других работ наталкивают на гипотезу: аттрактор динамической системы (ДК), который может быть сколь угодно сложным множеством (например, в приложениях типичны случаи предельных циклов, нескольких положений равновесия, и даже хаотических аттракторов), является таким множеством, в малой окрестности которого на больших временах с большой вероятностью будет пребывать рассматриваемая макросистема. Эта гипотеза для многих наиболее интересных и важных на практике случаев обоснована в диссертации с помощью техники доказательства теоремы 3.

Здесь возможны разные варианты поведения. Например, в случае двух притягивающих равновесий динамика системы “притянется” к каждому из них с вероятностями, зависящими от точки старта, и далее система может пребывать в малой окрестности равновесия длительное время, до большой флуктуации, которая сможет перебросить её в другое положение равновесия (Вентцель– Фрейдлин, 1979). Другой тип поведения: постоянно циркулировать (перемешиваться) по аттрактору. Возможно, конечно, и сочетание отмеченных типов поведения.

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

Рассматривается общая задача энтропийно–линейного программирования:

m ln xk e max ; A : x q Tx 0l, F x d Gx 0w. (ЭЛП) x k m xA kДля отыскания единственного решения (ЭЛП) было предложено семейство двойственных барьерно-мультипликативных итерационных алгоритмов, которое содержит алгоритмы Брэгмана–Шелейховского, MART, GISM, Ю.С. Попкова и др.:

lw n n n x exp p , k t g q k 1,...,m;

pk qk p1 q n1 p hp qp m xk p 1,...,l;

n n (ИП) p t pk , k1 m n1 n n n q q q hq dq g xk q 1,...,w, qk , k1 l,m wm, где шаги T tpk p,k1, G gqk q,k1, 0, 0 – достаточно маленькие, а дважды l w гладкие функции hp , hq – монотонно возрастают и равны нулю в p1 qнуле. Например, часто выбирают hp y ln 1 y qp, p 1,...,l ; hq y y dq, q 1,...,w.

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

Глобальная сходимость итерационного процесса была установлена при следующих предположениях:

1) z 0m : q Tz 0l, d Gz 0w;

2) Строки матриц T и G линейно независимы в совокупности.

При этом было замечено, что предположение 2 не всегда выполняется в прикладных задачах (см., например, модель ранжирования web-страниц Дж.А.

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

Итерационные алгоритмы (ИП) особенно эффективны:

для сильно разреженных матриц T и G ;

при l w m.

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

Текст соответствующих компьютерных программ, написанных на МаtLab, имеется в приложении.

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

Исследование возникающих макросистем позволяет глубже понять концепцию равновесия Нэша–Вардропа.

Предварительно опишем модель Бэкмана равновесного распределения потоков на графе транспортной сети.

Пусть транспортная сеть города представлена ориентированным графом V, E, где V – узлы сети (вершины), E V V – дуги сети (рёбра графа).

Пусть W w i, j :i, j V – множество корреспонденций, т.е. возможных пар «исходный пункт» – «цель поездки»; p v1,v2,...,vm – путь из v1 в vm, если vk,vk1 E, k 1,...,m 1, m 1; Pw – множество путей, отвечающих корреспонденции wW, то есть если w i, j, то Pw – множество путей, начинающихся в вершине i и заканчивающихся в j ; P Pw – совокупность wW всех путей в сети ; xp [автомобилей/час] – величина потока по пути p, x xp : p P ; ye [автомобилей/час] – величина потока по дуге e :

1, e p ye x xp, где ep e p ;

ep pP 0, e ye – удельные затраты на проезд по дуге e (как правило, возрастающие, выпуклые, гладкие функции).

Рассмотрим теперь Gp x – затраты временные или финансовые на проезд по пути p. Естественно считать, что Gp x ye x ep, где e ye – удельные e eE затраты на проезд по дуге e. Как правило, предполагают, что это – возрастающие, гладкие функции от ye. В приложениях часто требуется учитывать также затраты на прохождения вершин графа, которые могут зависеть от величин всех потоков через рассматриваемую вершину.

Пусть также известно, сколько перемещений в единицу времени dw осуществляется согласно корреспонденции wW. Тогда вектор x, характеризующий распределение потоков, должен лежать в допустимом множестве: X x 0 : xp dw, wW.

pPw Это множество может иметь и более сложный вид, если дополнительно учитывать, например, конечность пропускных способностей рёбер (ограничения сверху на ye ).

Рассмотрим игру, в которой каждому элементу wW соответствует свой, достаточно большой ( dw 1), набор однотипных “игроков”, осуществляющих передвижение согласно корреспонденции w. Чистыми стратегиями игрока служат пути, а выигрышем – величина Gp x. Игрок “выбирает” путь следования p Pw, при этом, делая выбор, он пренебрегает тем, что от его выбора также “немного” зависят Pw компонент вектора x и, следовательно, сам выигрыш Gp x. Можно показать, что отыскание равновесия Нэша–Вардропа x* X (макро описание равновесия) равносильно решению задачи нелинейной комплементарности (принцип Вардропа): для любых wW, pPw выполняется x* Gp x* min Gq x* 0.

p qPw Заметим, что рассматриваемая нами игра принадлежит к классу, так называемых, потенциальных игр. В нашем случае это означает, что существует epxp pP такая функция x e z dz, что x xp Gp x для любого eE p P. Оказывается, что x* X – равновесие Нэша–Вардропа тогда и только тогда, когда оно доставляет минимум x на множестве X.

Введем теперь в эту статическую модель динамику. Пусть на шаге n игрок, осуществляющий перемещение согласно корреспонденции w, выбрал путь p.

Обозначим количество таких игроков через xp n. Тогда на шаге n 1 он с вероятностью 1 n выбирает тот же путь, а с вероятностью использует n смешанную стратегию: выбирает путь p Pw с вероятностью w Probw n 1 n max xp n,1 n exp Gp x n T Zn, wW, p w где Zn определяется из условия нормировки. Множитель max xp n,1 n характеризует желание имитировать, т.е. использовать стратегию большинства, а также надежность использования этой стратегии (при одинаковых затратах на разных путях «надежнее» выбрать тот путь «вчерашнего дня», который использовало большее число игроков – относительная флуктуация, характеризующая риски, будет меньше). Параметр характеризует «консерватизм» («ленивость»), чем меньше , тем более консервативный игрок;

«температура» T характеризует отношение к риску («горячность»), чем больше температура, тем более «горячий игрок», склонный к более рискованным действиям.

Как показали разнообразные численные эксперименты, часто вполне разумно выбирать 1 n. При таком выборе наблюдается сходимость к равновесию n n Нэша–Вардропа при наиболее общих условиях относительно T (вне зависимости от точки старта). Стоит также обратить внимание на высокую эффективность предложенной процедуры «нащупывания равновесия». Иначе говоря, на предложенный итерационный процесс можно смотреть просто как на эффективный способ численного нахождения равновесия Нэша–Вардропа.

Предложенную схему можно трактовать как стохастическую динамику наилучших ответов в эволюционной (популяционной) игре, при этом имеется много общего с концепциями «quantal response equilibria» (используется похожая рандомизация) и «minority games» (наблюдаются похожие колебания около положения равновесия). Однако наиболее близким к предложенному итерационному процессу является эффективный приближенный вероятностный алгоритм Григориадиса–Хачияна (см. также Фройнда–Шапире), являющийся предтечей современных вариантов метода стохастического зеркального спуска.

Имеет место следующее утверждение, при доказательстве которого активно используется недавняя работа Немировского–Юдитского–Шапиро и др. по стохастическому зеркальному спуску.

Утверждение. Пусть T 0 – достаточно мало, , . Тогда n n n1 nп.н.

x n x 0 – одно из равновесий.

x* n Если равновесие x* единственно, то x* x 0 x*.

Случай Случай Рис. Случай 1: x124 x134 3. Случай 2: x124 x1324 x134 2.

Полное время в пути T 83 мин Полное время в пути T 92 мин Пример (парадокс Браеса, 1968). Пусть корреспонденция d14 6 (тысяч [автомобилей/час]). Вес ребра e (удельные затраты на проезд по этому ребру) есть время движения по ребру (в минутах) e ye, если поток через ребро есть ye (тысяч [автомобилей/час]). Значения функций e ye приведены на рис. 1.

Напомним, что потоки по ребрам можно рассчитать исходя из потоков по путям, например, в случае 2 на рис. 1: y24 x124 x1324. Оба равновесия Нэша–Вардропа (в случаях 1 и 2) являются «притягивающими» положениями равновесия описанной выше динамики.

Кажется удивительным, но в парадоксе Браеса в случае строительства “паразитной” дороги участники движения, действуя согласно описанной выше динамике (каждый преследует свои интересы), действительно “приходят” к единственно возможному “плохому” равновесию Нэша–Вардропа, которое в данном примере не оптимально по Парето. Такие ситуации уже давно наблюдались в реальной жизни (например, в Германии и Голландии). Пару лет назад был поставлен эксперимент на базе парадокса Браеса с участием студентов в Лаборатории экспериментальной экономики МФТИ, который подтвердил, что после введения “паразитной” дороги водители–студенты приходили, причем довольно быстро, к “плохому” равновесию и пребывали там большую часть времени эксперимента.

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

В заключении приведены основные результаты диссертации.

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

ОСНОВНЫЕ РЕЗУЛЬТАТЫ ДИССЕРТАЦИИ 1. Получено новое достаточное условие, обеспечивающее существование и единственность равновесия динамической макросистемы. Показано, что при довольно естественных условиях время сходимости к равновесию достаточно мало.

2. Рассмотрена макросистема игровой природы (модель равновесного распределения потоков), когда равновесие этой макросистемы интерпретируется как равновесие Нэша. Введенная динамика позволяет в случае, если таких равновесий много, выбрать единственное. При этом процедура отбора основана на концепции равновесия макросистемы.

3. Выявлен характер связи функции Ляпунова макросистемы с инвариантной мерой этой макросистемы.

4. Для известной статической энтропийной модели расчета матриц корреспонденций предложено обобщение, учитывающее транспортную доступность мест работы, а также динамическая интерпретация, основанная на стохастической динамике и описывающая стремление людей жить поближе к месту работы.

5. Получено обобщение модели миграции населения В. Вайдлиха на случай более общих правил, описывающих возможности мигрирования.

6. Разработаны эффективные вычислительные алгоритмы поиска равновесий динамической макросистемы, сводящиеся к задачам энтропийно-линейного программирования. Полученные алгоритмы являются обобщением ранее известных алгоритмов Брэгмана–Шелейховского, MART, GISM, Ю.С.

Попкова и др. Разработан соответствующий комплекс программ.

СПИСОК ПУБЛИКАЦИЙ АВТОРА ПО ТЕМЕ ДИССЕРТАЦИИ 1. Гасникова Е.В. Двойственные мультипликативные алгоритмы для задачи энтропийно-линейного программирования // Труды 50-й научной конференции МФТИ. Современные проблемы фундаментальных и прикладных наук. – М., Долгопрудный: Изд-во МФТИ, 2007 г. Ч. 7. Т. 1. С.

106-109.

2. Гасникова Е.В. Барьерно-мультипликативные алгоритмы для двойственных задач // Тезисы докладов конференции молодых ученых и специалистов «Информационные технологии и системы ИТиС'08». 29 сентября – 3 октября 2008 года, изд-во ИППИ РАН, г. Геленджик. С. 423-429.

3. Гасникова Е.В. Барьерно-мультипликативные алгоритмы для задач, двойственнх к задачам оптимизации энтропийно подобных функционалов на выпуклых множествах простой (аффинной) структуры // Труды 51-й научной конференции МФТИ. Современные проблемы фундаментальных и прикладных наук. – М., Долгопрудный: Изд-во МФТИ, 2008 г. Ч. 7. Т. 1. С. 2125.

4. Гасникова Е.В. Поиск равновесий макросистем с помощью мультипликативнобарьерных модификаций итерационных алгоритмов Коши и Ньютона // Тезисы докладов международной конференции «Математическая теория систем» МТС-2009» 6 – 30 января 2009 года, М.: Изд-во МФТИ. С. 150-154.

5. Гасникова Е.В. Двойственные мультипликативно-барьерные алгоритмы отыскания вырожденных равновесий макросистем // Тезисы докладов международной конференции «Современные проблемы математики, механики и их приложений» 30 марта – 02 апреля 2009 года МГУ, М.: МАКС Пресс, ВМиК МГУ. С. 317.

6. Гасникова Е.В. Двойственные мультипликативные алгоритмы для задач энтропийно-линейного программирования // ЖВМ и МФ Т.49 № 3 Изд-во МАИК «Наука/Интерпериодика», 2009 г. С. 453-464.

7. Гасников А.В., Гасникова Е.В. О принципе максимума энтропии в сетевых моделях // Тезисы докладов IV Всероссийской научной конференции «Математическое моделирование развивающейся экономики и экологии» (ЭКОМОД-2009), Киров, 6 – 12 июля 2009. Киров: Изд-во ВятГУ, 2009. С. 38.

8. Гасников А.В., Гасникова Е.В. О принципе максимума энтропии в модели (А.Дж. Вильсона) расчета матрицы корреспонденций // Сборник трудов IV Всероссийской научной конференции «Математическое моделирование развивающейся экономики» (ЭКОМОД-2009), Киров, 6-12 июля 2009. Киров:

Изд-во ВятГУ, 2009. С. 110-121.

9. Гасников А.В., Гасникова Е.В., Дорн Ю.В. О некоторых примерах конечных однородных эргодических марковских процессов с огромным числом состояний // Труды 52-ой научной конференции МФТИ, – М., Долгопрудный:

Изд-во МФТИ, 2009. Ч. 7. Т. 1. С. 19-22.

10. Буслаев А.П., Гасников А.В., Гасникова Е.В., Холодов Я.А., Яшина М.В.

Возможная динамика в модели Дж. Вардропа. Труды VI международной конференции по исследованию операций. М.: МАКС Пресс. 19-23 октября 2010. С. 22-23.

11. Гасникова Е.В. О равновесиях макросистем. Труды VI международной конференции по исследованию операций. М.: МАКС Пресс. 19-23 октября 2010. С. 139-140.

12. Гасников А.В., Гасникова Е.В. О возможной динамике в модели расчета матрицы корреспонденций (А.Дж. Вильсона) // Труды МФТИ, Т. 2. № 4(8).

Изд-во МФТИ, 2010. С. 45-54.

13. Гасникова Е.В, Дорн Ю.В. О стохастической марковской динамике, приводящей к равновесию Нэша – Вардропа в модели распределения потоков // Труды МФТИ, Т. 2. № 4(8). Изд-во МФТИ, 2010. С. 55-57.

14. Гасникова Е.В. О социально-экономических приложениях уравнений стохастической химической кинетики, динамике средних и динамике, полученной в результате скейлинга // Труды 53-й научной конференции МФТИ, – М., Долгопрудный: Изд-во МФТИ, 2010. Ч. 7. Т. 1. С. 88-90.

15. Гасникова Е.В. О стохастической марковской динамике наилучших ответов, приводящей к равновесию Нэша–Вардропа // Труды 53-й научной конференции МФТИ. – М., Долгопрудный: Изд-во МФТИ, 2010. Ч. 7. Т. 1. С.

106-107.

16. Гасникова Е.В., Ишманов М.С., Колесников А.В., Нагапетян Т.А. О скорости сходимости к равновесию макросистемы // Тезисы докладов VI Всероссийской научной конференции «Математическое моделирование развивающейся экономики и экологии» (ЭКОМОД-2011), Киров, 27 июня – 3 июля 2011.

Киров: Изд-во ВятГУ, 2011, С. 35.

17. Гасникова Е.В., Нагапетян Т.А. О достаточных условиях существования равновесия макроситсемы // Тезисы докладов VI Всероссийской научной конференции «Математическое моделирование развивающейся экономики и экологии» (ЭКОМОД-2011), Киров, 27 июня – 3 июля 2011. Киров: Изд-во ВятГУ, 2011, С. 36.

18. Аввакумов С.Я., Гасников А.В., Гасникова Е.В. О сходимости к равновесию Нэша-Вардропа // Тезисы докладов Международной научной конференции «Динамические системы, нелинейный анализ и их приложения», Армения (г.

Ереван): Изд-во МАКС Пресс, 10 июля – 17 июля 2011. С. 47-48.

19. Гасников А.В., Гасникова Е.В., Колесников А.В., Нагапетян Т.А. О концепции равновесия макросистемы и её приложении к модели равновесного распределения потоков // Труды Международной научной конференции “лет ИППИ РАН”, – М.: Изд-во ИППИ РАН, 25 июля – 29 июля 2011. Москва, 2011. С. 53-61.

20. Gasnikov A.V., Gasnikova E.V., Nagapetyan T.A. On the convergence to one of the Nash-Wardrop's equilibriums // Traffic and granular flow’11. International conference. Moscow, 28 September – 1 October 2011. Изд-во МТУСИ. P. 296298.

21. Гасникова Е.В., Аввакумов С.Я. О равновесиях в моделях распределения потоков // Труды VI Всероссийской научной конференции «Математическое моделирование развивающейся экономики и экологии» (ЭКОМОД-2011), Киров, 27 июня – 3 июля 2011. Киров: Изд-во ВятГУ, 2011. С. 111-118.

22. Гасникова Е.В. О новых условиях существования равновесия макросистемы // Труды VI Всероссийской научной конференции «Математическое моделирование развивающейся экономики и экологии» (ЭКОМОД-2011), Киров, 27 июня – 3 июля 2011. Киров: Изд-во ВятГУ, 2011. С. 100-110.

23. Гасникова Е.В. Достаточные условия существования равновесия макросистем // Труды 54-й научной конференции МФТИ, – М., Долгопрудный: Изд-во МФТИ, 2011. Ч. 7. Т. 1. С. 64-65.

24. Гасникова Е.В. Концепция равновесия макросистемы в модели распределения потоков // Труды 54-й научной конференции МФТИ, – М., Долгопрудный:

Изд-во МФТИ, 2011. Ч. 7. Т. 1. С. 88-89.

25. Gasnikova E.V., Nagapetyan T.A. About new dynamical interpretations of entropic model of correspondence matrix calculation and Nash–Wardrop's equilibrium in Beckmann's traffic flow distribution model. Ninth International Conference on Traffic and Granular Flow, 28 September – 1 October 2011. Moscow: Springer, 2012. arXiv:1112.1628v26. Gasnikov A.V., Gasnikova E.V. A conception of equilibrium of a macrosystem in application to the traffic flow distribution in large transport networks // International Conference "Instabilities and Control of Excitable Networks: From macro – to nano-systems" Dolgoprudny, Moscow Area, RUSSIA, May 25 – 30, 2012. P. 16.

27. Гасников А.В., Гасникова Е.В., Петрашко Д.И. Макросистемный подход к ранжированию web-страниц // «Информационные технологии и системы – 2012» 35-я конференция молодых ученых и специалистов ИППИ РАН, ПреМоЛаб, СТРАДО. г. Петрозаводск: Изд-во ИППИ РАН. С. 423-429.

28. Gasnikov A.V., Gasnikova E.V. Stochastic optimization in the model of correspondences matrix calculation and traffic flow distribution // 21st International Symposium on Mathematical Programming (ISMP 2012). Springer: Berlin, August 19 – 24, 2012. P. 217.

29. Gasnikov A.V., Gasnikova E.V. Stochastic subgradient barrier-multiplicative descent for entropy optimization // International conference OPTIMA-2012. September 23 – 30, 2012. Costa da Caparica, Portugal: Из-во ВЦ РАН. P. 91-92.

30. Гасникова Е.В. О возможной динамике в модели расчета матриц корреспонденций // Приложение к монографии «Введение в математическое моделирование транспортных потоков» Под ред. А.В. Гасникова. М.:

МЦНМО, 2012. C. 250-273.

Гасникова Евгения Владимировна Моделирование динамики макросистем на основе концепции равновесия Автореферат Подписано в печать 02.11.2012.

Печать трафаретная Усл.п.л. – 1,Заказ № 4893 Тираж: 100 экз.

Типография «11-й ФОРМАТ» ИНН 77263309115230, Москва, Варшавское ш., (499) 788-78-www.autoreferat.ru







© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.