Московский Государственный Университет
имени М. В. Ломоносова
Механико-математический факультет
На правах рукописи
УДК 519.21
Ярыкин Павел Николаевич
НЕЛИНЕЙНЫЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ
И АНАЛИЗ СИСТЕМ
ВЗАИМОДЕЙСТВУЮЩИХ ЧАСТИЦ
теория вероятностей и математическая статистика) (01.01.05АВТОРЕФЕРАТ
диссертации на соискание ученой степени кандидата физико-математических наук
Москва 2006
Работа выполнена на кафедре теории вероятностей Механико-математического факультета Московского государственного университета им. М. В. Ломоносова.
кандидат физико-математических наук,
Научный руководитель:
доцент А. Д. Манита доктор физико-математических наук,
Официальные оппоненты:
профессор Б. М. Гуревич кандидат физико-математических наук, с. н. с. С. В. Анулова Институт проблем передачи
Ведущая организация:
информации РАН
Защита диссертации состоится " " 2006 г. в 16 часов 15 минут на заседании диссертационного совета Д.501.001.85 в Московском государственном университете им. М. В. Ломоносова по адресу: ГСП-2, Москва, Ленинские горы, МГУ, Механико-математический факультет, аудитория 16–24.
С диссертацией можно ознакомиться в библиотеке Механико-математического факультета (главное здание, 14 этаж).
Автореферат разослан " " 2006 г.
Ученый секретарь диссертационного совета Д.501.001.85 в МГУ, Т. П. Лукашенко доктор физико-математических наук, профессор
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность темы.
1. Исторический контекст. Математическая теория (классических) стохастических дифференциальных уравнений (СДУ) возникла в 40-x годах XX века с работы Ито1, в которой было введено понятие стохастического дифференциала и его свойства. С тех пор теория классических СДУ была значительно развита, и в настоящее время их решения изучены для широкого круга условий2.
Широкий интерес представляют математические модели больших систем взаимодействующих броуновских частиц. С одной стороны, эти модели являются многомерными классическими СДУ. Но большая размерность сильно затрудняет получение конкретных результатов о поведении системы. В частности, интересны системы с парным взаимодействием типа среднего поля, когда влияние, оказываемое на частицу со стороны других частиц усреднено.
То есть речь идет о системе частиц, которая может быть описана системой СДУ N i,N i,N i,N j,N dXt = dWti + a Xt (1) + b X t, Xt dt N j= с начальными условиями X0. Здесь Wti iN i независимые стандартные броуновские движения, a(x) задает снос частицы (процесса) под воздействием внешнего поля и b(x, y) функция попарного взаимодействия частиц.
i Предполагается, что начальные условия X0 являются независимыми одинаково распределенными случайными величинами.
В 50-х годах при описании эволюции частиц разреженного газа (уравнеi,N ние Больцмана) Кац3 предположил, что частицы Xt являются асимптотически независимыми при больших N. В силу симметрии, маргинальные распределения всех частиц совпадают. Соответственно, распределение каждой частицы в момент времени t близко к эмпирической мере µN,X (dx) = t 1 Ito K. Stochastic integral. Proc. Imp. Acad. Tokyo. 1944, v. 20, p. 519–524.
2 Анулова С. В., Веретенников А. Ю., Крылов Н. В., Липцер Р. Ш., Ширяев А. Н. Стохастическое исчисление. Итоги науки и техн. Соврем. пробл. матем. Фундам. направления. ВИНИТИ, 1989. т. 49.
3 Kac M. Foundation of kinetic theory. Proc. Third Berkley Symp. on Math. Stat. and Prob. 1956, v. 3, p. 171–197.
а поведение выделенной частицы можно аппроксимироi=1 Xt (dx), вать решением нелинейного СДУ где мера µX является распределением Xt. Таким образом, нелинейные СДУ позволяют исследовать асимптотическое поведение физических систем.
Процессы со взаимодействием типа среднего поля и связанные с ними нелинейные случайные процессы образуют широкий класс случайных процессов и представляют научный интерес. В настоящее время, они исследуются многими авторами при различных предположениях4,5,6.
Нелинейные СДУ оказались существенно сложнее классических СДУ для изучения, и их решения были исследованы в меньшей степени. В значительной мере это связано с тем, что решение нелинейного СДУ не является марковским процессом, хотя уравнение (2) и является обобщением однородного по времени классического СДУ.
2. Актуальные результаты и задачи. До недавнего времени основные результаты по теме исследования были получены в работах Шнитмана7 и Тамуры8,9 . Шнитман получил результаты о пределе среднего поля для системы взаимодействующих частиц (1) в случае ограниченного липшицевого взаимодействия. Результаты Тамуры о существовании, единственности и устойчивости стационарного решения нелинейного СДУ (2) были получены в случае, когда взаимодействие быстро убывает с ростом расстояния и имеется сильное полиномиальное центростремительное внешнее поле. Под устойчивостью стационарного распределения понимается слабая сходимость реМалышев В. А., Манита А. Д. Фазовые переходы в модели синхронизации времени.
Теория вероятн. и примен., 2005, т. 50, в. 1, с. 150–158.
5 Karpelevich F. I., Rybko A. N. Thermodynamic Limit for the Mean Field Model of Simple Symmetrical Closed Queueing Network. Markov Processes and Related Fields. 2000, v. 6, p. 89–105.
6 Manita A., Shcherbakov V. Asymptotic analysis of particle system with mean-eld interaction. Markov Processes and Related Fields. 2005, v. 11, p. 489–518.
7 Sznitman A. S. Topics in propagation of chaos. Lect. Notes in Math., Springer, Berlin.
1989, v. 1464, p. 165–250.
8 Tamura Y. On asymptotic behaviors of the solution of a non-linear diusion equation.
J. Fac. Sci. Univ. Tokio Sect. IA, Math. 1984, v. 31, p. 195–221.
9 Tamura Y. Free energy and the convergence of distributions of diusion processes of McKean type. J. Fac. Sci. Univ. Tokio Sect. IA, Math. 1987, v. 34, p. 443–484.
шения (его распределения в момент времени t) с любыми начальными условиями к стационарному распределению при t +.
В последнее время появился ряд работ, развивающих вышеуказанные результаты. Как правило10,11, авторы предполагают наличие достаточно сильного центростремительного внешнего поля, подавляющего нелинейную составляющую сноса процесса. На этом фоне выделяются работы12,13, в которых результат типа результата Тамуры получен для выпуклого ядра попарного взаимодействия в предположении отсутствия внешнего поля. Получение подобных результатов для других видов взаимодействия в предположении отсутствия внешнего поля является актуальной задачей.
Также заметим, что во всех изученных случаях накладывались ограничения, которые приводили к единственности инвариантного распределения.
Интересной задачей является явное описание и изучение случая, когда имеется несколько стационарных распределений.
Цель работы. Целью данной работы является разностороннее изучение нелинейного случайного процесса Xt (решения нелинейного СДУ (2)) на вещественной прямой R и его связи с системой взаимодействующих частиц (решением системы СДУ (1)). Исследования ограничиваются случаем, когда внешнее воздействие на систему отсутствует, то есть a(x) 0, а ядро взаимодействия (x) состоит из двух компонент: линейно возрастающей силы притяжения и ограниченного липшицевого возмущения, то есть В части работы на ядро взаимодействия накладывается более сильное ограничение:
10 Veretennikov A. Yu. On ergodic measures for McKean-Vlasov stochastic equations. Isaac Newton Inst. for Math. Sci. 2003, preprint NI03066.
11 Carrillo J. A., McCann R. J., Villani C. Kinetic equilibration rates for granular media and relates equations: entropy dissipation and mass transportation estimates. Revista Matematica Iberoamericana. 2003, v. 19, p. 1–48.
12 Benachour S., Roynette B., Talay D., Vallois P. Nonlinear self-stabilizing processes I:
Existence, invariant probability, propagation of chaos. Stochastic Processes and Appl. 1998, v. 75, p. 173–201.
13 Benachour S., Roynette B., Vallois P. Nonlinear self-stabilizing processes II: Convergence to invariant probability. Stochastic Processes and Appl. 1998, v. 75, p. 203–224.
В рамках этого предположения решение нелинейного СДУ исследуется для всех значений параметра.
Еще одной целью работы является явное указание такого взаимодействия, когда СДУ (2) имеет несколько стационарных решений, а также изучение свойств соответствующего нелинейного процесса и стационарных распределений. Эта цель выполняется в вышеуказанных рамках, поскольку искомым взаимодействием оказывается взаимодействие с ядром (x) = x + sin x при достаточно больших.
Методы исследования. В работе применены методы теории вероятностей и случайных процессов, в частности, теория марковских случайных процессов, характеристические функции, преобразования стохастических дифференциалов, сходимость вероятностных мер. Также были использованы методы функционального анализа. Еще одним важным инструментом является функционал свободной энергии.
Научная новизна. Основные результаты работы являются новыми и состоят в следующем:
1) доказаны существование и единственность сильного решения нелинейного СДУ (2);
2) доказан предел среднего поля для системы частиц (решения системы СДУ (1)) в случае возмущения 1 (x) с константой Липшица < 1/4;
3) доказано существование стационарных решений нелинейного СДУ (2) и получены результаты об их количестве (в частности, при достаточно больших получена неединственность стационарного распределения);
4) исследована устойчивость найденных стационарных распределений.
Теоретическая и практическая ценность. Полученные в диссертации результаты носят теоретический характер. Работа относится к области случайных процессов. Ее результаты могут быть использованы для изучения больших систем частиц с парным взаимодействием типа среднего поля.
Апробация диссертации. Результаты диссертации докладывались автором на следующих конференциях и семинарах.
• Научно-исследовательский семинар Теория вероятностей и статистическая физика механико-математического факультета под руководством Оселедеца В. И. и Гуревича Б. М. МГУ, 2002 г.
• XXV Конференции молодых ученых механико-математического факультета МГУ им. М.В. Ломоносова. Название доклада: Среднеполевая аппроксимация для одной системы взаимодействующих частиц. (МГУ, апрель 2003 г.) • Научно-исследовательский семинар Теория вероятностей и статистическая физика механико-математического факультета под руководством Оселедеца В. И. и Гуревича Б. М. (МГУ, 2004 г.) • Научно-исследовательский семинар Вероятностные методы в биологии механико-математического факультета МГУ под руководством Малышева В. А. (МГУ, 2004 г.) • Ломоносовские чтения. Название доклада: Устойчивость нелинейного стохастического процесса, аппроксимирующего систему взаимодействующих частиц. (МГУ, апрель 2005 г.) • Научно-исследовательский семинар Добрушинской математической лаборатории ИППИ РАН под руководством Минлоса Р. А.( ИППИ РАН, • Большой Семинар кафедры теории вероятностей механико-математического факультета МГУ под руководством Ширяева А. Н. (МГУ, октябрь 2005 г.) Публикации. Результаты диссертации опубликованы в 5 работах автора [1–5], список которых приводится в конце автореферата.
Структура и объем работы. Диссертация состоит из введения, трех глав, дополнения и списка литературы. Общий объем работы составляет 85 страниц. Список литературы включает 45 наименований.
Поддержка. Исследования по теме диссертации частично были поддержаны Российским фондом фундаментальных исследований (гранты 99-01и 02-01-00945).
КРАТКОЕ СОДЕРЖАНИЕ ДИССЕРТАЦИИ
Во введении приводится исторический обзор развития теории нелинейных СДУ и систем взаимодействующих частиц. Там же сформулированы основные результаты данной работы, которые являются новыми.Первая глава посвящена предварительному изучению нелинейного СДУ для симметричного трансляционно инвариантного взаимодействия b(x, y) = = (x y) = b(x, y) вида (x) = x + 1 (x), где 1 (x) ограниченная функция с константой Липшица.
В параграфе §1.1 приводится определение решения нелинейного СДУ (3) и доказывается теорема о его существовании. Попутно, как необходимая часть доказательства, получается свойство конечности первого момента.
В параграфе §1.2 выводится ряд свойств найденного решения, преимущественно касающихся его гладкости и различных оценок его плотности.
Полученные оценки плотности будут активно использоваться в главе 3.
Теорема 1 является основным результатом первой главы.
Теорема 1. Уравнение (3) при начальном условии X0 с конечным первым моментом имеет единственное сильное решение Xt на промежутке t [0, ). Более того, математическое ожидание EXt существует и ограничено.
Решение нелинейного СДУ понимается в следующем смысле.
Определение. Пусть имеется (расширенная P-нулевыми событиями) фильтрация F = (Ft )t 0, винеровский процесс W = {Wt, t 0} относительно F и начальное условие X0 = X0, измеримое относительно F0. Сильным решением нелинейного СДУ (3) на промежутке [0, ) называется случайный процесс X = {Xt, t 0}, имеющий п. н. непрерывные траектории, согласованный с фильтрацией F и такой, что при подстановке его и семейства его распределений {µX, t 0} в левую и правую части формулы (3) при каждом t 0 получается равенство с вероятностью единица.
Найденное решение СДУ (3) обладает следующими свойствами:
• неизменность первого момента;
• абсолютная непрерывность распределения и гладкость его плотности;
• непрерывность по времени в слабой топологии;
• полиномиальная скорость убывания плотности (и ее производных по x) Вторая глава посвящена изучению системы случайных частиц Xt со взаимодействием типа среднего поля и ее связи с нелинейным случайным процессом Xt. При этом многочастичный процесс удовлетворяет следующей системе СДУ, ассоциированной со СДУ (3):
где X0 независимые одинаково распределенные случайные величины, независимые реализации стандартного броуновского движения.
Результаты данной главы применимы при < 1/4. В частности, доказывается, что изучаемый процесс Xt действительно является пределом среднеN Полученный результат также позволяет заключить, что стационарное распределение существует и единственно. Более того, любой процесс, удовлетворяющий нелинейному уравнению и имеющий второй момент, будет слабо сходиться при t + к стационарному распределению.
Основным результатом данной главы является теорема 3, которая заключается в следующем:
Теорема 3. Пусть EX0 и EX0 конечны. Положим X0 = X0 и Wt1 = Wt.
Тогда при < выполнено где c(t) является некоторым коэффициентом, зависящим от времени.
Ключевую роль в этой главе играет лемма, которая устанавливает равномерную по времени оценку расстояния между случайными процессами, близкими к нелинейному и многочастичному соответственно.
Точнее, в лемме рассматриваются процессы Yti,N и Yti,N, задаваемые формулами:
образование позволяет избежать размывания многочастичного случайного процесса в пространстве путем перехода к системе отсчета, привязанной к его центру масс. Аналогичное преобразование N независимых реализаций нелинейного случайного процесса применяется для их согласования с многочастичным процессом.
Оказывается, что выделенная частица Yti,N многочастичного процесса в подвижной системе координат и нелинейный процесс Yti,N в своей подвижной системе координат близки равномерно по времени, то есть Кроме того, такое же неравенство выполнено для разности нелинейного процесса Xt Xt и его модификации Yt1,N.
Похожее, но не равномерное по времени неравенство связывает Xt и Yt1,N :
Более того, процесс Yti,N является эргодическим, что позволяет поi= лучить теорему 4.
Теорема 4. Пусть EX0 и EX0 конечны и < 1/4. Тогда имеется слабая сходимость распределений где (одномерное) распределение стационарного решения СДУ (3).
В третьей главе изучается предельное поведение решения нелинейного уравнения при = x + sin x. Основные результаты этой главы относятся к случаю отрицательных и случаю больших положительных.
В случае 1/2 доказывается, что стационарное решение единственно и что решение с любым начальным условием слабо сходится к стационарному решению (при t ). При малых || получена также скорость этой сходимости.
Для больших положительных доказывается, что существует два устойчивых стационарных распределения и что остальные стационарные распределения неустойчивы. Здесь устойчивость стационарного распределения понимается в смысле сходимости решения стохастического дифференциального уравнения с близким начальным условием к соответствующему стационарному распределению (локальная устойчивость). Эта часть результатов получена в инвариантном (замкнутом относительно динамики процесса) подмножестве четных распределений Psym (R) множества вероятностных распределений P(R).
Глава разбита на 5 параграфов.
Параграф §3.1 содержит результаты о стационарных распределениях процесса. В нем находится в явном виде двухпараметрическое семейство распределений, которое содержит в себе все стационарные распределения.
Теорема 5. Любое стационарное решение уравнения (3) в классе вероятностных плотностей с нулевым средним имеет вид причем решение существует для любого. Для 1/2 решение единственно и соответствует паре вида (a, b)() = (a(), 0).
Более того, для всех R находятся (уже в неявном виде) значения параметров, задающие в классе четных распределений Psym (R) стационарное (или стационарные, если их несколько) распределение. Также доказывается единственность стационарного решения при < 0 (где 0 > 1 ) и существование нескольких стационарных решений при 0.
Заметим, что в классе Psym (R) стационарное распределение задается одним параметром a, поскольку параметр b равен 0.
Теорема 6. Рассмотрим стационарные распределения СДУ (3) в классе Psym (R). Тогда существуют 0 > 1 и a0 < 0 такие, что:
а) при < 0 существует ровно одно стационарное распределение, соответствующее значению a1 () > 0;
б) при = 0 существует два стационарных распределения, соответствующих значениям a1 () > 0 и a2 () = a0 ;
в) при > 0 существует три стационарных распределения, соответствующих значениям a1 () > 0, a2 () < a0 и a0 < a3 () < 0.
Основные результаты данной главы находятся в параграфах §3.4, §3. и содержатся в теоремах 12–14. При этом, теорема 12 аналогична классическим результатам о предельном поведении нелинейных процессов, в то время как теоремы 13 и 14 относятся к принципиально иному случаю, когда существует несколько стационарных решений СДУ (3).
Теорема 12. При 1/2 существует единственное стационарное решение СДУ (3) с заданным математическим ожиданием. При этом оно глобально устойчиво в классе распределений с тем же математическим ожиданием, то есть любое решение СДУ (3) с таким X0, что существуют EX0 и EX0, слабо сходится при t + к стационарному распределению с тем же математическим ожиданием.
Обозначим стационарные плотности из класса Psym (R), соответствующие параметру ai, через pai.
Теорема 13. Стационарная плотность pa3 неустойчива, то есть для любой окрестности в слабой топологии Upa3 Psym (R) существует p Upa3 такое, что решение СДУ (3) с начальным распределением p не сходится к стационарному решению с плотностью pa3.
Теорема 14. Стационарные плотности pa1 и pa2 локально устойчивы.
То есть существуют окрестности точек pa1 и pa2 в Psym (R) со слабой топологией такие, что решение СДУ (3) с начальным распределением из этих окрестностей сходится слабо при t + к соответствующему стационарному решению.
Заметим также, что ключевым понятием этого раздела является функционал свободной энергии В частности, результаты теорем 12–14 получены исходя из пространственно-временных свойств данного функционала, сформулированных в нескольких теоремах и леммах.
По своей сути функционал свободной энергии является аналогом функции Ляпунова в теории дифференциальных уравнений. Его основные свойства получены в параграфе §3.3. К ним относятся:
• значение функционала от распределения Xt убывает как функция от t (кроме случаев, когда одномерное распределение Xt является стационарным);
• функционал свободной энергии ограничен снизу;
• стационарные распределения СДУ (3) в точности соответствуют критическим точкам функционала свободной энергии;
• стационарное распределение СДУ (3) является (локально) устойчивым тогда и только тогда, когда оно соответствует (локальному) минимуму функционала свободной энергии.
В параграфе §3.2 при малых по модулю с помощью преобразования Фурье показана экспоненциальная скорость слабой сходимости решения СДУ (3) с произвольным начальным условием к его стационарному решению.
В дополнение вынесены вспомогательные замечания и доказательства.
Благодарность. Автор выражает искреннюю благодарность кандидату физико-математических наук, доценту Анатолию Дмитриевичу Маните за постановку задач, постоянное внимание, многочисленные ценные советы и помощь в работе.
СПИСОК РАБОТ АВТОРА ПО ТЕМЕ ДИССЕРТАЦИИ
[1] П. Н. Ярыкин. Устойчивость нелинейного стохастического процесса, аппроксимирующего систему взаимодействующих частиц. Теория вероятностей и ее применения, 2006, т. 51, вып. 2, с. 400–409.[2] П. Н. Ярыкин. Предельные свойства нелинейного стохастического процесса, аппроксимирующего систему взаимодействующих броуновских частиц. Деп. в ВИНИТИ 07.12.05 №1606–В2005, 2005, 41 с.
[3] П. Н. Ярыкин. Поведение нелинейного случайного процесса в окрестности его стационарных распределений. Успехи математических наук, 2006, т. 61, вып. 4, с. 199–200.
[4] П. Н. Ярыкин. Поведение стохастического процесса, описывающего систему частиц со взаимодействием среднего поля. Вестник Московского Университета, Сер. 1, 2004, № 2, с. 55–58.
[5] П. Н. Ярыкин. Среднеполевая аппроксимация для одной системы взаимодействующих частиц. Труды XXV Конференции молодых ученых (31 марта – 5 апреля 2003 г.), 2004, т. II, с. 251–255.