Санкт-Петербургский государственный университет
Фундаментальная математика и механика
Механика жидкости, газа и плазмы
Баринова Ольга Вячеславовна
«Частные решения уравнения Колмогорова-Чепмена и их связь с
уравнениями математической физики»
Выпускная квалификационная работа
Научный руководитель:
доктор физико-математических наук, профессор Мирошин.Р.Н.
Рецензент:
доктор физико-математических наук, профессор Халидов.И.A.
Санкт-Петербург
2016г.
SAINT-PETERSBURG STATE UNIVERSITY
Fundamental Mathematics and Mechanics
Mechanics of a liquid , gas and plasma
Barinova Olga
“Particular solutions of Chapman- Kolmogorov equation and their connection with
equations of mathematical physics”
Graduation Thesis
Scientific supervisor:
Professor, Doctor of Physics and Mathematics Miroshin Roman
Reviewer:
Professor, Doctor of Physics and Mathematics Khalidov Iscander
Saint-Petersburg
2016
Оглавление
Введение…………………………………………………………………………..4
1. Глава I. Случайные марковские процессы……………………………….........6
1.1. Процессы с непрерывной траекторией…...................................................8
1.2. Чисто разрывные процессы……………...................................................10
2. Глава II. Уравнения Колмогорова………………………………………........12
2.1. Частное решение уравнения Колмогорова-Чепмена для непрерывного
процесса. (Пример 1)………………………………………………………...12
2.1.1. Об оценке непрерывной функции распределения. (Лемма 1)…...12
2.1.2. О связи математической физики и плотности вероятности
перехода. (Утверждения 1, 2). ……………………………………………15
2.2. Частное решение уравнения Колмогорова-Чепмена для
скачкообразного процесса. (Пример 2)………………………..……………19
2.2.1. Об оценке интеграла. (Лемма 2)…………………………………....21
2.2.2. Разрывный процесс (Пример 3)……………………………………24
2.2.3. Об оценке разрывной функции распределения (Лемма 3)….........26
2.2.4. Дополнительная оценка интеграла (Лемма 4)…………………….28
2.2.5. Применение дробных производных в первом и втором уравнениях
Колмогорова. (Утверждение 3, Лемма 5)……………................................31
3. Заключение……………………………………………………………………34
4. Литература…………………………………………………………………….35
3
Введение
При изучении различного рода явлений мы нередко сталкиваемся с теми
процессами, течение которых заранее предсказать в точности невозможно.
Иными словами, это те процессы, физическая система которых с течением
времени переходит из одного состояние в другое случайным образом. Теория
случайных процессов в последние десятилетия бурно развивается, потому
как является востребованной в таких науках как механика, физика, биология,
экономика, социология и в других областях. Особое место среди случайных
процессов занимает так называемые марковские случайные процессы,
впервые введенные русским ученым А.А.Марковым в 1907г.
Отличительную черту марковского процесса можно сформулировать
следующим образом: Если именовать момент t настоящим, а моменты
прошлым,
(
то
зависимость
настоящем только от значений процесса в последний
процесса
в
перед настоящим
момент времени из прошлого. Итак, марковские процессы отличаются тем,
что нам нет необходимости знать "предысторию" состояния системы, а лишь
достаточно знать ее "настоящее", и тогда мы сможем предсказать "будущее"
состояние системы.
Настоящая работа включает в себя две главы, первая из которых
отведена теории марковских процессов с непрерывной и разрывной
траекториями,
вторая
—
исследованию
примеров,
описывающих
вышеупомянутые процессы, взятых из статьи [2]. В частности, в качестве
примеров, были рассмотрены три частных
уравнения
Колмогорова-Чепмена
и
по
ним
решения интегрального
получены
уравнения
Колмогорова. Среди примеров были взяты решения в виде интеграла без
особенности,
со
слабой
особенностью
и
сингулярный
интеграл.
Первостепенная цель данной работы — по имеющимся решениям уравнения
Колмогорова-Чепмена получить первое и второе уравнения Колмогорова,
используя асимптотику и приемы из математической физики. Среди них
4
уравнения
диффузии,
использованием
дробных
интегро-дифференциальные
производных, ранее
в
уравнения,
теории
с
марковских
процессов нам неизвестны, хотя дробные производные в прикладных науках
в последнее время нередко встречаются.
5
Глава I. Случайные марковские процессы.
Эта глава посвящена основным определениям, к которым мы будем
возвращаться в дальнейшем. Основополагающим источником литературы
является Б.В.Гнеденко [1]. Автор этого пособия - крупный специалист по
теории вероятности и математической статистике и ее приложениях в науке и
технике.
В этом и в последующих параграфах, мы будем рассматривать
исключительно марковские случайные процессы с непрерывным параметром
t, именуемым временем. Также мы будем предполагать, что множество
всевозможных состояний системы принадлежит множеству действительных
чисел.
Итак, дадим определение случайного марковского процесса с
непрерывным временем.
О п р е д е л е н и е 1.2. [1, 3]
Случайный процесс
называется марковским если выполняется
равенство для любой
где
плотность
,
распределения процесса.
В частности, при
Запись
можно преобразовать в запись без интегралов, учитывая
6
произвольность
:
или
где
плотность вероятности перехода процесса из
состояния
Таким
в момент времени
образом,
в состояние
марковский
начального распределения
процесс
в момент времени
характеризуется
функцией
и плотностью вероятности перехода
,
уравнение для которой имеет вид
где
Уравнение
называют обобщенным уравнением Маркова (в
терминологии Гнеденко [1]) или уравнением Колмогорова-Чепмена. В то же
время, уравнение, аналогичное
можно записать и для вероятности
перехода, учитывая соотношение
где
пространство состояний.
Тогда преобразованное уравнение
7
примет вид
К последнему уравнению
мы будем обращаться в дальнейшем.
1.1. Процессы с непрерывной траекторией.
Настоящий параграф посвящен процессам, траектория которых является
непрерывной. Будем говорить, что случайный процесс
непрерывен, если
за сколь угодно малые промежутки времени лишь с малой вероятностью
может получить заметные по величине приращения. Иными словами,
сформулируем определение, сославшись на [1].
О п р е д е л е н и е 1.2. [1]
Случайный процесс
непрерывен, если для
существует
Здесь
величина
вероятность того, что в момент
случайная
принимает значение, меньшее , если известно, что в момент
времени
процесс принимал значение, равное
Заметим, как и всякая функция распределения,
при любых
удовлетворяет следующим условиям:
В то же время функция
определена для
и для нее
выполняются дополнительные условия:
Также предположим существование частных производных первого и
второго порядков
8
Пусть они непрерывны при любых
Тогда для
существуют пределы:
На основе этих предположений сформулируем первую теорему [1].
Т е о р е м а 1.1. [1]
Если выполнены условия (1.1), (1.2) и (1.3), тогда функция распределения
удовлетворяет дифференциальному
уравнению в частных
производных второго порядка
Доказательство выводится из уравнения Колмогорова-Чепмена
свойств функции распределения
и
[1]. В то же время уравнение
(1.8) можно записать с использованием плотности вероятности перехода, при
ряде допущений, а именно:
1) Пусть существует
2)
непрерывные производные
9
Т е о р е м а 1.2. [1]
Если выполнены условия (1.5), (1.6) и (1.7), а также вышеизложенные
допущения 1), 2), то функция распределения
дифференциальному
уравнению в частных производных второго порядка:
Определив непрерывный случайный процесс, мы ознакомились с
уравнениями Колмогорова для непрерывного процесса. Наша ближайшая
цель определить разрывный процесс и выявить для них уравнения.
1.2. Скачкообразные марковские процессы.
Значимую роль играют процессы, изменение состояния системы
которых с течением времени меняется не непрерывно, а скачкообразно.
Случайный процесс
чисто разрывен, если за любой промежуток времени
траектория процесса, начавшаяся в точке x, остается равной с
вероятностью
, а с вероятностью
может претерпеть скачком некоторое изменение.
О п р е д е л е н и е 1.3. [1]
Пусть
условная вероятность распределения случайной величины
, при условии, что в момент t произошел скачок, а до скачка величина
принимала
значение
при
равное
х.
Тогда
функция
распределения
выражается через функции
следующим образом:
Отметим, что, как и для всякой функции распределения, так и для
10
выполняются равенства
Также будем иметь ввиду неотрицательность и непрерывность
относительно t, x.
Т е о р е м а 1.3. [1]
Функция
распределения
вида
(1.9)
удовлетворяет
двум
интегро-
дифференциальным уравнениям:
В завершении первой главы, хочется отметить, что нам удалось
познакомиться
с
различными процессами, а также с уравнениями,
описывающими те самые процессы. Ближайшими задачами являются
получить уравнения Колмогорова на примерах, основываясь на теоремы 1.1,
1.2, 1.3.
11
Глава II. Уравнения Колмогорова.
Определяющим уравнением в теории марковских процессов является
уравнение
Колмогорова-Чепмена.
условиях
из
этого
С.Н.Бернштейн
интегрального
уравнения
показал
можно
при
каких
получить
дифференциальное уравнение второго порядка параболического типа [4]. В
дальнейшем, О.В.Сарманов получил частное решение в виде ряда [5].
Известен общий вид решения уравнения Колмогорова, полученный в работе
[2]. Нам следует выяснить, к какому виду дифференциального уравнения
сводятся
некоторые
частные
решения
интегрального
уравнения
Колмогорова-Чепмена.
2.1.
Частное
решение
уравнения
Колмогорова-Чепмена
для
непрерывного процесса. (Пример 1).
В статье [2] представлены некоторые решения интегрального уравнения
Колмогорова-Чепмена. В этом параграфе будем рассматривать плотность
вероятности перехода, которая имеет следующий вид
Как показано в [2], плотность (2.1) является решением уравнения
Колмогорова-Чепмена
Рассмотрим, какому линейному уравнению
отвечает это решение. Следуем процедуре, описанной в предыдущей главе.
Глядя на решение (2.1), интуитивно понятно, что траектория является
непрерывной, но чтобы убедиться в этом, докажем следующую лемму.
2.1.2. Об оценке непрерывной функции распределения.
Л е м м а 1.
Марковский случайный процесс с плотностью вероятности перехода
12
имеет непрерывную траекторию, т.е. выполняются условия
Доказа ль
о
В
л (a)
ла м зам н
. По л ч о
ло
Для удобства разобьём
первый случай, когда
оо
а
чным
о о о н
ала
м
на два предельных интеграла и рассмотрим
анало
. Случай
нно
ф нкц я
ф нкц я ош бок ко о ая о
Анало
м нных:
оц
ля
ам
ош бок
яч
ол ча м
13
ч н
о олн
льная
з
льно
знач н
ля
ак м об азом
об а н
зам н
ло
ло
а
акж
а
ы олн но П о з о я
ы олня
я Ч о
бо ало ь
оказа ь.
П
к
й
м
оказа ль
ь к оказа ль
ло
я н
Ра мо
ак
ы но
м л чай
ж зам н
азобьём
на
омян
а н
ала
о о о
(b)
нахож
н я, а м нно
н
ф нкц й
ый ко а
ыш
о а ы аж н
14
ж
а акж
ля
ла м
ля
об
м
а
л
П
й
о а льно
м к нахож
хож м ы аж н
я ля
н
я коэфф ц
н ом
ой ф нкц
а ным н л
), ко о ая за а
я
м а м нно
Таким образом, мы получили сумму двух интегралов равную двум.
л
о а льно,
Л мма оказана.
В итоге нам удалось проверить условия на непрерывность с плотностью
вероятности
перехода
(2.1),
после
чего
сформулируем
некоторые
утверждения, основанные на теоремах из первой главы.
2.1.2. О связи математической физики и плотности вероятности
перехода.
У т в е р ж д е н и е 1.
Если условия (a) и (b) в лемме 1 выполнены, тогда функция распределения
15
и плотность вероятности перехода
удовлетворяют первому (прямому) уравнению Колмогорова
Действительно, случайный процесс с функциями (2.2) и (2.3) является
непрерывным, что показано выше в лемме. В общем виде из теоремы 1
первое уравнение Колмогорова для функции распределения имеет вид
или для плотности
где
Итак, подставив данные коэффициенты в вышеуказанные уравнения,
получаем (2.4) и
Полученные уравнения в частных производных второго порядка
параболического типа, являются уравнениями теплопроводности или же
уравнениями, описывающими явление диффузии.
Второе уравнение Колмогорова было получено физиками Фоккером и
Планком, в связи с развитием теории диффузии. Выясним, какой вид будет
иметь второе (или обратное) уравнение Колмогорова, где частные
16
производные первого и второго порядка берутся соответственно, не по t, а по
s и не по х, а по y.
У т в е р ж д е н и е 2.
Если выполнены условия (a) и (b) в лемме 1, а также выполняется (2.4),
тогда для непрерывного случайного процесса плотность распределения
вероятностей (2.3) удовлетворяет второму (или обратному) уравнению
Колмогорова
Действительно,
при
подстановке
коэффициентов
и
во второе уравнение Колмогорова в общем виде, упомянутом в
первой главе (теореме 2), получаем (2.5).
Также в данном параграфе покажем сведение к уравнению
, но уже
другой процедурой, а именно, методом моментов.
В [3] показано уравнение (2.6), эквивалентное уравнению Колмогорова
, и полученное методом моментов. Обозначив
, уравнение
можно записать следующим образом
где
коэффициенты интенсивности, определяющиеся равенством
Также
имеет
место
тождество:
А
математическое ожидание случайной величины А. Или в другой записи:
Ниже приведены вычисления моментов для случаев
17
Вычисленные моменты
подставим в выражение (2.7) для
отыскания коэффициентов интенсивности. Заметим, моменты при
равны нулю.
Поставим (2.8) в уравнение (2.6) и получим уравнение Колмогорова
18
2.2. Частное решение уравнения Колмогорова-Чепмена для
скачкообразного процесса. (Пример 2).
Перейдем к примерам, описывающим процесс с разрывной траекторией.
Так же, как и в предыдущем параграфе, мы будем рассматривать частные
решения Колмогорова Чемпена [2].
Имеем следующую плотность вероятности перехода
Пользуясь классическим определением функции распределения
и взяв интеграл от плотности, получаем
Используем первые члены ряда Тейлора для этих функций [6]:
19
Тогда имеем два случая
или в упрощенном виде:
С помощью функции Хэвисайда, т.е.
функцию распределения можно записать в виде
или
20
По
а ляя
ольз я ь
а н н
ой
ом
Ч
м на
олмо о о а
ф нкц
ол ча м
где
Убедимся в том, что последним интегралом в
можно пренебречь.
2.2.2. Об оценке интеграла.
Л е м м а 2.
Для интеграла имеет место следующая оценка
Доказательство.
Для удобства, обозначим
и проинтегрируем
по частям.
21
Преобразовывая последнее равенство, получаем
Или в другом виде:
Т.к
величина более высокого порядка, то первым слагаемым в
правой части предыдущего равенства можно пренебречь, а во втором
слагаемом мы ограничимся величиной , следовательно
22
Осталось доказать при
интеграл сходится к нулю.
Здесь, в (2.13) интеграл понимается в смысле главного значения по Коши.
Сделав замену
, положив
причем
имеем
следующую запись:
Ограничиваясь первым членом в разложении выражения в круглых
скобках по степеням p, получаем
Перейдем к вычислению предела c учетом вышеуказанного выражения (2.15)
При выполнении соотношений
, выражение (2.16) стремится к
нулю:
23
Лемма доказана.
В итоге мы получили, слагаемым
можно пренебречь в
.
где
2.2.2. Разрывный процесс. (Пример 3).
Рассмотрим еще одно из решений уравнения Колмогорова-Чепмена,
полученного в [2].
где
синус-преобразования Фурье от функций
24
н
ы ная ф нкц я
Как известно, функция распределения связана с плотностью вероятности
перехода следующим соотношением:
Подставляя (2.17) в (2.19) получим выражение для
Взяв интеграл по переменной у в (2.20) получаем функцию в виде:
Отметим основные свойства функции распределения
1) Неотрицательность;
2) при
, функция принимает значение, равное нулю;
3) при
значение функции стремится к единице.
Свойства 1)
3) нетрудно проверить. Неотрицательность
вытекает из неотрицательности
Заметим,
Осталось показать, что
что
Следовательно, упомянутая функция является характеристической (теорема
Пойа [3]), а это гарантирует неотрицательность синус-преобразования Фурье
(теорема Бохнера [3]).
Во втором свойстве легко убедиться, если в (2.21) вычислить интеграл
при
В самом деле,
25
И наконец, перейдем к третьему свойству. Положим
тогда
При помощи формул Эйлера, определим функции следующим образом:
Тогда можно применить лемму Эрдейи [7] к (2.22) и получаем (2.23).
мя я к б кон чно
где
о
(2.23) будет выполняться при
л нном
на
м
В итоге, нам удалось убедиться в вероятностном смысле решения (2.17),
проверив на свойствах функцию (2.21).
Теперь выведем новое соотношение для функции распределения
, но уже без интегралов.
2.2.3. Об оценке разрывной функции распределения.
Л е м м а 3.
Функция распределения
может быть записана в виде
26
где функция представима
т.е.
л
л
Доказательство.
Действительно, имеем выражение для функции
в котором h устремляем к нулю, причем
Возьмем интеграл по переменной z, и в то же время, оставим от
два
первых
о а м
м
члена
о л
асимптотики
при
условии,
что
ы аж н
Также (2.24) можно получить, используя лемму Эрдейи [7].
н
Воспользуемся [8.68] при
Подставим (2.25) в (2.24), получаем преобразованное выражение:
Как и в предыдущих примерах, мы будем следовать процедуре , описанной в
[1]. Для этого подставим последнее выражение для
27
в уравнение
Колмогорова Чепмена (4 ) и получаем, положив
и используя свойство
ф нкц
где
Следующая лемма направлена на оценку последнего интеграла в (2.26),
для возможного использования производной
особенность функции
, невзирая на
в точке
2.2.4. Дополнительная оценка интеграла.
Л е м м а 4.
При предположении
справедлива следующая оценка
где
Доказательство.
Пусть
тогда проинтегрировав (2.27) по частям,
получаем:
Преобразуем (2.28), сославшись на формулу (2.21)
28
Первый член последнего выражения порядка
осталось разобраться со
вторым. Преобразуем последнее выражение к виду:
Таким образом, учитывая оценки для интегралов
мы получили оценку
Лемма доказана.
Перейдем к выводу уравнения Колмогорова, но перед этим, определим
некоторые понятия.
О п р е д е л е н и е 2.1. [8]
В равенстве
функция
называется дробной производной от
Г.Харди и Д.Е.Литтльвуда [6] и обозначается
Согласно [6], при
т.е.
29
степени
в смысле
При
дробная производная
совпадает с обычной
Определяя дробную производную в смысле Гельдера [4] равенством
Видим, что она совпадает с правой частью (2.3018), но, в отличие от (2.32),
нелокальная (зависит от
восстановление
Дробная производная (2.29) допускает
, как по обычной производной.
О п р е д е л е н и е 2.2. [8]
Правосторонняя и левосторонняя дробные производные Римана-Лиувилля
порядка
от вещественной функции
определяется следующими
соотношениями
л
о о онняя
а о о онняя
Интегрируя в этих равенствах по частям, находим что
Таким образом, получаем
30
2.2.5. Применение дробных производных в первом и втором
уравнениях Колмогорова.
У т в е р ж д е н и е 3.
Вероятности
удовлетворяют
перехода
первому
и
второму
уравнениям
Колмогорова
соответственно:
где
Доказательство.
Обратимся к уравнению (2.21) при
Решение (2.33) подставим в уравнение Колмогорова-Чепмена и получим
31
Л е м м а 5.
При
имеет место следующая оценка:
Доказательство.
Воспользуемся
соотношением
запишем
(2.34) через плотность вероятности перехода, предварительно сделав замену
переменной
Используя последнее, возьмем интеграл (2.36)
л
где справедлива оценка
о а льно
справедлива и оценка
Аналогичными рассуждениями, мы получаем для вероятности перехода
второе уравнение Колмогорова:
32
Заключение
В настоящей работе установлена связь некоторых частных решений
билинейного
уравнения
Колмогорова-Чепмена
с
линейными
дифференциальным и интегро-дифференциальными уравнениями. Причем,
рассмотрены примеры случайных марковских процессов с непрерывной и
разрывной
траекторией.
Оказалось,
среди
интегро-дифференциальных
уравнений есть уравнения с сингулярными интегралами и уравнения с
дробными производными. При выводе использовались нестандартные
оценки интегралов посредством леммы Эрдейи и процедуры взятия
интеграла с помощью главного значения по Коши.
33
Литература
1. Гнеденко Б. В. Курс теории вероятностей. М.: Книжный дом
«ЛИБРИКОМ», 2011. 488 с.
2. Мирошин Р. Н. О некоторых решениях интегрального уравнения
Колмогорова-Чепмена // Вестн. С.-Петерб. Ун-та. Сер. 1: Математика,
механика, астрономия. 2007. Вып. 4. С. 22-29.
3. Мирошин Р. Н. О некоторых решениях интегрального уравнения
Колмогорова-Чепмена // Вестн. С.-Петерб. Ун-та. Сер. 1: Математика,
механика, астрономия. 2007. Вып. 4. С. 22-29.
4. Бернтштейн С. Н. О зависимостях между случайными величинами //
Собр. Соч. Т. 4.: Наука, 1964. С. 235-254.
5. Сарманов О. В. Исследование стационарных марковских процессов
методом разложения по собственным функциям // Труды Мат. Ин-та АН
СССР. Т. 60. М.: Наука, 1961. С. 238-259.
6. Двайт Г. Б. Таблицы интегралов и другие математические формулы//
Пер. с англ. М.: Наука, 1966. 228 с.
7. Федорюк М. В. Метод аеревала. М.:Наука, 1977. 368 с.
8. Нахушев А.М. Элементы дробного исчисления и их применение. М.:
Физмалит, 2003. 272 с.
9. Бейтмен Г., Эрдейи А. Таблицы интегральных преобразований // Пер. с
англ. Т. 1. М.: Наука, 1969. 344 с.
34
Отзывы:
Авторизуйтесь, чтобы оставить отзыв