Интеграл Ито. Стохастические дифференциальные уравнения
[27/0%]Приведите пример адаптированного, но не прогрессивно измеримого процесса.
Следующие две -алгебры (опииональных и предсказуемых множеств) играют в стохастическом анализе важную роль в самых разнообразных вопросах.
Определение 8. Опциональная -алгебра подмножеств есть -алгебра, порожденная всеми согласованными процессами , рассматриваемыми как отображения из в , траектории которых принадлежат пространству Скорохода . Процессы, которые измеримы относительно -алгебры , называются опциональными.
Докажите, что равносильным образом опциональная -алгебра может быть определена как -алгебра подмножеств в , порожденная стохастическими интервалами вида
где — марковские моменты (относительно потока ).
Определение 9. Предсказуемая -алгебра подмножеств есть -алгебра, порожденная всеми согласованными процессами , рассматриваемыми как отображения из в , траектории которых принадлежат пространству непрерывных на действительных функций. Процессы, которые измеримы относительно -алгебры , называются предсказуемыми.
Докажите что равносильным образом предсказуемая -алгебра может быть определена как -алгебра, порожденная всеми согласованными процессами, траектории которых лишь непрерывны слева на .
Имеются и другие альтернативные (и во многих рассмотрениях более удобные) определения, см., например, [26] и [46].
Проверьте, что -алгебра порождается каждым из следующих классов подмножеств:
с и множествами вида с ;
с и стохастическими интервалами
где — марковский момент;
с и множествами вида с .
Убедитесь, что между введенными -алгебрами имеют место следующие включения:
Напомним, что по теореме 4 всякий согласованный с фильтрацией непрерывный слева или непрерывный процесс является прогрессивно измеримым.
Пусть процесс является прогрессивно измеримым и — конечный марковский момент. Покажите, что случайная величина будет -измеримой (это свойство является едва ли не определяющим целесообразность введения понятия прогрессивной измеримости).
Проверьте, что простые (кусочно-непрерывные справа) функции вида
где величины — -измеримы, , являются опциональными.
Покажите, что простые (кусочно-непрерывные слева для ) функции вида
где величина — -измерима, — -измеримы, , являются предсказуемыми.
Проведите построение интеграла Ито, изложенное в § 2 и § 3, отправляясь от класса опциональных простых функций, фигурирующих в упражнении 7 (которым приписывается то же значение интеграла Ито, что и для предсказуемых простых функций из упражнения 8). Можно ли утверждать, что замыкание в данной совокупности простых функций, входящих в , приводит к пространству ?
Пусть — марковский момент такой, что при всех (T — положительная константа). Докажите, что , где .
Докажите, что если есть неслучайная функция из пространства Скорохода , то процесс , является гауссовским. Найдите его среднее и ковариационную функцию.
Обсудим теперь некоторые обобщения интеграла Ито. Интегрирование можно определить для следующего класса функций .
Определение 10. Будем говорить, что функция принадлежит классу , если прогрессивно измерима и
Идея построения интеграла для функции состоит в том, чтобы подобрать последовательность функций , для которых интеграл уже построен, и таких, что
Отсюда выводится, что последовательность является фундаментальной по вероятности. Следовательно, существует случайная величина, обозначаемая , такая, что . Предел записывают, по-прежнему, как или . Часто броуновское движение обозначают , и тогда пишут .
Докажите, что существует последовательность , описанная выше в (85). Докажите, что если , то существует п. н. непрерывная модификация процесса .
Отметим, что для процесс не обязан быть мартингалом, но является локальным мартингалом (т.е. процессом, для которого существует последовательность моментов остановки таких, что п.н. и для каждого "остановленные" процессы , где , являются мартингалами).
Рассмотрим для процесс
Докажите, что (по-прежнему, все дифференциальные соотношения означают сокращенную запись соответствующих интегральных).
Подробнее со стохастическим исчислением можно ознакомиться по книгам , . Кратные стохастические интегралы рассматриваются в [163].
Стохастический интеграл по винеровскому процессу (броуновскому движению) позволяет дать "явное" представление фрактального броуновского движения.
Определение 11. Случайньй процесс со значениями в называется автомодельным (самоподобным), если для каждого можно найти такое , что
Иначе говоря, для таких процессов изменение временно́й шкалы ( приводит к тому же самому результату, что и изменение фазовой шкаль ( . Если в определении (87) для любого параметр , то случайный процесс называется автомодельным процессом с показателем Харста Н. Величина называется статистической фрактальной размерностью случайного процесса .
Напомним, что фрактальное броуновское движение для определяется как центрированный гауссовский процесс с ковариационной функцией
При процесс является броуновским движением.
Докажите, что , где , есть автомодельный процесс (с показателем Харста ).
Процессы впервые рассматривались А. Н. Колмогоровым в 1940 г. (см. [33]) и были названы им спиралями Винера. Термин фрактальное (или дробное, "fractional") броуновское движение был введен Б. Мандельброттом и Й. Ван Нессом в 1968 г. в работе [164], и там же была получена следующая теорема.
Теорема 10 (Мандельбротт, Ван Несс). Для и фрактальное броуновское движение допускает представление
где нормирующая константа выбрана так, что (см., например, [86; т. 1, с. 281]), a — независимые стандартные винеровские процессы, причем для .
рактальное броуновское движение (и его дискретные аналоги) используются во многих важных моделях, в частности, для описания динамики финансовых индексов (см. [86]). Сложность изучения процессов с состоит в том, что за исключением случая (броуновского движения) и эти процессы не являются семимартингалами (т.е. не входят в важный класс процессов, для которых развито стохастическое исчисление, см. [46; гл. 4]). Если в представлении (88) вместо взять гёльдеровскую функцию (т. е. ) со значениями в , то получается случайный процесс, называемый мультифрактальным броуновским движением.
Докажите, что для процесса при
Пусть — пространство прогрессивно измеримых функций , удовлетворяющих (84) при .
Следующая теорема, основанная на понятии стохастического интеграла, описывает структуру броуновских функционалов.
Теорема 11 (Ито, Кларк). Пусть — естественная фильтрация броуновского движения и является измеримой случайной величиной. Справедливы следующие утверждения.
-
Если , то найдется такой стохастический процесс , что
-
Если , то представление (89) справедливо с некоторым процессом .
-
Если — положительная случайная величина (т.е. и , то найдется процесс такой, что допускает представление вида , где определяется согласно (86).
Объясните, почему величина , фигурирующая в теореме 11 , есть функционал от броуновского движения, т.е. , где и .
Теорема 11 (Ито, Кларк). Пусть — естественная фильтрация броуновского движения и является -измеримой случайной величиной. Справедливы следующие утверждения.
-
Если , то найдется такой стохастический процесс , что
-
Если , то представление выше справедливо с некоторым процессом .
-
Если — положительная случайная величина (т.е. ) и , то найдется процесс такой, что допускает представление вида , где .
Доказательства теоремы 11 в различных вариантах даны К. Ито, Дж. Кларком, Дж. Дубом и приведены во многих книгах, см., например, .
В качестве примера использования представления функционалов от броуновского движения в виде стохастических интегралов приведем следующую экстремальную задачу:
среди всех моментов остановки (относительно естественной фильтрации броуновского движения найти такой момент , для которого достигается
Можно представить себе модель, в которой колебания цен на акции в течение единичного периода времени (например, день) описываются броуновским движением (разумеется, это лишь интерпретация, так как цены не бывают отрицательными) и требуется выбрать момент для самой выгодной в среднем квадратическом смысле продажи имеющихся акций. Заметим, что дает смещенную оценку для , поскольку и (поясните в качестве простого упражнения). Поэтому можно обратиться к величине
Легко проверить (убедитесь в этом), что
Положим
Теорема 12 (Граверсен, Пешкир, Ширяев). Решение приведенной выше задачи дается формулой
где константа находится из уравнения
здесь и ф соответственно функция распределения и плотность стандартной нормальной величины. При этом .
Ключевую роль при доказательстве этого результата играет представление вида
где const (укажите эту константу) и
Широкое использование результатов и методов стохастического анализа связано с возможностью построения теории стохастического интегрирования по процессам более общей структуры, нежели броуновское движение.
Определение 12. Согласованный с фильтрацией процесс называется возрастающим, если для п.в. имеем , где — неубывающая функция по и для каждого .
Определение 13. Действительный случайный процесс принадлежит классу Дирихле ( ), если семейство величин равномерно интегрируемо, где обозначает совокупность всех конечных марковских моментов (относительно фильтрации ).
Справедлив следующий фундаментальный результат (см. [26, т. 1, с. 68]; сравните с теоремой 1 главы IV).
Теорема 13 (разложение Дуба-Мейера). Субмартингал класса (D) с траекториями типа càdlàg допускает и притом единственное (с точностью до стохастической неразличимости) разложение вида
где — возрастающий предсказуемый интегрируемый (, где ) процесс, а есть равномерно интегрируемый мартингал.
Рассмотрим для простоты непрерывный квадратично интегрируемый мартингал на фильтрованном вероятностном пространстве такой, что . Класс таких мартингалов обозначается . Будем предполагать, что фильтрация удовлетворяет обычным условиям.
Наметим схему построения интеграла по процессу . Пусть — возрастающий предсказуемый процесс, участвующий в разложении Дуба-Мейера субмартингала , обладающий тем свойством, что .
Как и в случае броуновского движения, интеграл
вообще говоря, не может быть определен потраекторно, как интеграл Лебега-Стилтьеса (можно определить в тривиальном случае при ).
Процесс назовем простым, если существуют строго возрастающая числовая последовательность с и , а также последовательность действительных случайных величин и некоторая константа , для которых при всех , причем и
Класс простых процессов обозначим и для положим
Далее доказывается, что определение (92) можно распространить (снова используя должные аппроксимации простыми процессами) с на более широкий класс процессов. Интересно отметить, что этот класс будет зависеть от того, является ли функция для п.в. абсолютно непрерывной относительно меры Лебега. Детальное построение и дальнейшие обобщения можно почерпнуть в книгах .
С помощью формулы Ито докажите, что для детерминированной функции из справедлива формула "интегрирования по частям":
Постройте пример случайной функции , для которой формула (93) не имеет места.
Пример 2. Рассмотрим стохастически возмущенное уравнение роста популяции
Точнее говоря, рассмотрим стох астическое дифференциальное уравнение вида
с начальным условием и , где константы .
Из формулы Ито непосредственно следует, что процесс
является сильным решением уравнения (95). Поскольку , то
Иначе говоря, при средний рост (убывание, если ) будет таким же, как в случае детерминированной модели (94).
Уравнение (95) играет важную роль в финансовой математике. Его решение (96) часто называют "геометрическим броуновским движением". Локальный снос характеризует скорость изменения среднего значения процесса , диффузию в финансовой литературе часто называют волатильностью. Видимо, П. Самуэльсон был первым, кто осознал важность геометрического броуновского движения для описания эволюции цен, используя для него термин "экономическое броуновское движение".
Определения и результаты, излагавшиеся в главе VIII, ценой некоторых усложнений естественным образом переносятся на многомерный случай. Так, уравнение (59) можно понимать как систему, в которой — векторная, — матричная функции, точнее говоря,
— -мерное броуновское движение (относительно фильтрации . Считаем также, что
Интеграл от векторной функции определяется как вектор из интегралов от компонент, а обозначает вектор-функцию с -й компонентой, равной
Определение прогрессивной измеримости очевидным образом переносится на многомерный случай. Теорема 13 остается в силе, если в условиях (61) и (62) модуль понимать как норму в , а модуль понимать как норму матрицы, например, считать, что .
При исследовании систем стохастических дифференциальных уравнений важнейшую роль играют многомерные варианты формулы Ито, т. е. формулы, описывающие стохастический дифференциал процесса , где , а — векторный процесс Ито. Другими словами, имеет стохастический дифференциал вида
где и — соответственно матричная и векторная прогрессивно измеримые функции, . Иначе говоря, пусть
где при любом , всех и
Теорема 14. Пусть функция такова, что , т. е. существуют непрерывные производные и . Тогда процесс , дде процесс Ито определяется согласно (100), (102), имеет стохастический дифференииал, который задается формулой
здесь и d , фигурирующие в (101), “перемножаются по правилу":
Доказательство формулы Ито см., например, в [26]. По поводу разнообразных обобщений этой формулы см., например, .
Приведем здесь лишь один результат в направлении упомянутых обобщений. Пусть — стандартное (одномерное) броуновское движение . Пусть функция является абсолютно непрерывной,
причем функция , т.е. для любого
Определим — квадратическую ковариацию процессов и следующим образом:
Здесь - обозначает предел по вероятности, а при каждом — римановские последовательности детерминированных моментов , т.е. для при каждом и при ; кроме того, для любого
Важно подчеркнуть, что поскольку процесс не является, вообще говоря, семимартингалом, нетривиален факт существования предела по вероятности, задающего . Один из результатов работы состоит именно в доказательстве существования этого предела.
Теорема 15 (Проттер, Фёльмер, Ширяев; [129]). При сделанных выше предположениях о функции имеет место формула
Пусть функция . Докажите, что
Тем самым, формула (103) переходит в формулу Ито.
Пусть . Докажите, что , где — локальное время (П. Леви) броуновского движения в точке 0 , определяемое следующим образом:
Таким образом, в данном случае формула (103) переходит в формулу Танака для броуновского движения.
Обобщению интеграла Ито на функции многих переменньх посвящена монография [163]. Приведем только один результат этого направления.
Теорема 16 (Ито). При всех и справедлива формула
где — полином Эрмита степени п, т.е.
Обратимся теперь к однородной диффузии, т.е. рассмотрим решение уравнения
где — -мерный винеровский процесс, а функции и удовлетворяют условиям теоремы 8 , которые в данном случае сводятся к одному требованию: существует такое, что
(для вектора обозначает евклидову норму и ), из которого вытекает, что для некоторого
Обозначим единственное сильное решение уравнения (106) при .
Докажите, что для каждого определенный выше процесс , является однородным марковским процессом.
Докажите, что при выполнении условий теоремы 8 единственность сильного решения (на ) будет в более широком классе процессов, у которых существует для .
Теорема 8. Рассмотрим стохастическое дифференциальное уравнение
где и определены на . Пусть выполнено условие Липшица: существует константа такая, что
и условие линейного роста: для некоторого
Пусть величина является -измеримой, причем . Тогда существует единственное сильное решение уравнения с -измеримым начальным условием такое, что для любого и функция ограничена на .
Покажите, что если условия теоремы 8 выполнены на полуинтервале , то ее утверждение также справедливо на этом промежутке, только в этом случае рассматривается класс процессов , для которых функция будет ограничена на каждом отрезке, вложенном в .
Теорема 8. Рассмотрим стохастическое дифференциальное уравнение
где и определены на . Пусть выполнено условие Липшица: существует константа такая, что
и условие линейного роста: для некоторого
Пусть величина является -измеримой, причем . Тогда существует единственное сильное решение уравнения с -измеримым начальным условием такое, что для любого и функция ограничена на .
С помощью формулы Ито докажите, что решением (одномерного) уравнения
является процесс
(продолжение упражнения 23). Докажите, что п.н. при . Таким образом, решение уравнения (109) представляет собой броуновский мост над отрезком с закрепленными концами и . Стандартный броуновский мост получается при и .
Из различного рода обобщений теоремы существования и единственности решения стохастического дифференциального уравнения приведем несколько неожиданный результат А.К. Звонкина (см., например, [86; с. 322]), утверждающий, что для существования сильного решения стохастического дифференциального уравнения
достаточно лишь измеримости по и равномерной ограниченности .
Докажите, что стохастическое дифференциальное уравнение
с "плохим" коэффициентом имеет и, притом, единственное сильное решение.
Определение 14. Стохастическое дифференциальное уравнение (59) с начальным условием , где — заданная мера на , имеет слабое решение на промежутке , если найдутся фильтрованное вероятностное пространство , броуновское движение и непрерывньй (п.н.) случайньй процесс такие, что и Р-п. н. для каждого выполнено равенство (60).
Разумеется, вместо отрезка можно рассматривать промежутки вида , где , или для .
Важно подчеркнуть, что в отличие от сильного решения, рассматриваемого на заданном фильтрованном вероятностном пространстве с заданным на нем броуновским движением, в определении слабого решения такие объекты (вероятностное пространство и броуновское движение) не фиксируются, а требуется лишь, чтобы они нашлись. Очевидно, сильное решение является и слабым.
Определение 15. Слабая единственность (слабого или сильного) решения уравнения (59) означает, что у любых двух решений (в рассматриваемом смысле) с одинаковыми начальными условиями совпадают распределения, т. е. совпадают конечномерные распределения.
Докажите, что уравнение
где , имеет по крайней мере два решения (но "слабых"). Кроме того, на некоторых вероятностных пространствах у этого уравнения может вовсе и не быть "сильного" решения.
Для выполнения последнего упражнения следует убедиться, что процесс
является броуновским движением на отрезке . Это утверждение обеспечивает (см., например, [83]) следующая
Теорема 17 (П. Леви). Пусть — -п.н. непрерывный квадратично интегрируемый мартингал, заданный на некотором фильтрованном вероятностном пространстве . Пусть также является мартингалом, т. е.
Тогда есть стандартное броуновское движение.
Заметим, что М. Барлоу [92] доказал, что уравнение (113) может не иметь сильного решения даже в случае непрерывных ограниченных функций .
Пусть выполнены условия теоремы 8. Докажите, что решение (слабое или сильное) уравнения (59) слабо единственно.
Упражнение 26 показывает, что могут существовать слабые решения, но не быть сильных. Поэтому естественно ожидать существования слабого решения при менее ограничительных условиях на коэффициенты уравнения (59).
Один из первых результатов в этом направлении (см. об этом [86; т. 1, с. 325]) формулируется следуюшим образом. Рассмотрим стохастическое дифференциальное уравнение
с начальным распределением таким, что
Если коэффициенты и являются непрерывными ограниченными функциями, то уравнение (115) имеет слабое решение.
Замечание 3. В предположении ограниченности и невырожденности коэффициента утверждение о существовании и единственности (по распределению) слабого решения уравнения (115) остается в силе, если потребовать лишь ограниченность и измеримость коэффициента . Приведенные результаты о слабых решениях также допускают обобщения на многомерный случай, на случай коэффициентов, зависящих от прошлого и т. д. (см., например, [86; т. 1, гл. III]).
Важную роль в теории стохастических дифференциальных уравнений (в частности, для построения слабых решений) играет теорема Гирсанова об абсолютно непрерывной замене меры, для формулировки которой введем необходимые обозначения.
Пусть — некоторое фильтрованное вероятностное пространство, — -мерное броуновское движение, . Пусть , где — прогрессивно измеримый -мерный случайный процесс такой, что
здесь и .
Образуем процесс , полагая
где .
Лемма 10 (см., например, [86; т. 2, с. 326]). Если выполнено условие Новикова:
то и введенный в (117) процесс будет равномерно интегрируемым мартингалом.
В силу положительности (-п.н.) и условия на можнозадать вероятностную меру , полагая
Теорема 18 (Гирсанов). Для введенных выше процессов и определим
Тогда является -мерным броуновским движением на фильтрованном пространстве .
Эта теорема, как уже было отмечено, дает возможность строить слабые решения стохастических дифференциальных уравнений вида
а также более общих уравнений (см., например, [86; гл. VII, § 36], [148; гл. 5]).
Наряду с интегралом Ито в ряде задач полезно обращение к "симметризованным" интегралам Стратоновича (или фиска-Стратоновича):
где — броуновское движение, а функция входит в определенный класс. Чтобы пояснить суть дела, скажем, что для некоторого класса функций интеграл (121) строится с помощью интегральных сумм вида
где .
Доводы в пользу трактовки стохастического дифференциального уравнения как интегрального уравнения на основе интеграла (121) можно прочитать в [169]. В частности, замена переменных в интеграле Стратоновича не приводит к появлению дополнительных членов второго порядка, возникающих в формуле Ито. Это используется при изучении стохастических дифференциальных уравнений на многообразиях (см. [128]). В то же время при изучении стохастических уравнений вида
требуются бо́льшие ограничения на коэффициенты, чем при изучении уравнений вида (60). Если функция дифференцируема по , то уравнение (122) равносильно следующему уравнению Ито:
В частности, если не зависит от , то обе интерпретации уравнения (60) совпадают.
Подчеркнем также, что интеграл Ито является мартингалом (например, при условиях теоремы 2), а интеграл Стратоновича этим свойством не обладает.
Для дальнейшего ознакомления с различными аспектами теории стохастических дифференциальных уравнений можно обратиться, например, к . Стохастические уравнения в частных производных изучаются в [196]. Квантовая диффузия рассматривается, например, в [88].
Заключительная часть этого раздела посвящена некоторым проблемам, относящимся к самостоятельным научным направлениям в рамках общей теории случайных процессов и ее приложений.
Исключительно важным направлением является теория фильтрации. Предположим, что на полезный сигнал налагается "шум" и требуется выделить (отфильтровать) этот сигнал из зашумленных наблюдений. Формализация так просто описанного явления нетривиальна. Мы обратимся к модельной ситуации, которая привела к построению очень красивой теории.
Пусть интересующий нас -мерный процесс , — полезный, но "скрытый" сигнал — описывается стохастическим дифференциальным уравнением (системой) (59) с начальным условием и пусть "наблюдаемый" -мерньй процесс , задается стохастическим уравнением
где . Считается, что все рассматриваемые случайные процессы определены на некотором вероятностном пространстве , причем -мерное броуновское движение не зависит от -мерного броуновского движения , фигурирующего в (59).
Одна из возможных здесь задач состоит в том, чтобы по наблюдениям дать оценку величины для . Обозначив , будем искать оптимальную оценку такую, что
Легко видеть, что . Мы будем предполагать, что все рассматриваемые случайные величины обладают конечным вторым моментом (это относится к компонентам векторов), а также, что выполнены условия, обеспечивающие существование решений используемых стохастических уравнений.
Следующий классический результат (см., например, [169]) для линейных систем явился одним из первых результатов, дающих удобные для практической реализации формулы, описывающие динамику оптимального фильтра.
Теорема 19 (фильтр Кальмана-Бьюси). Пусть
Пусть матрица обратима для всех и норма матрицы ограничена на каждом конечном промежутке изменения . Тогда решение задачи (125) удовлетворяет стохастическому дифференциальному уравнению
с начальным условием , а матричная функция
подчиняется матричному дифференциальному уравнению Риккати
Пример 3. Применив теорему 19 к определенному ниже частному случаю скалярных процессов, найдем . Пусть
Уравнение Риккати (127) для приобретает вид
откуда
Для из (126) находим, что
Отсюда, используя формулу Ито для диффузионных процессов, находим, что
следовательно,
Теории фильтрации посвящена, например, книга [160]. В качестве введения в данную область рекомендуем прочитать главу IV в [169].
Мы лишь укажем на еще одно важное обстоятельство, которому уделяется пристальное внимание последние десять лет. Если в упомянутой вьше задаче фильтрации "шум" рассматривался как помеха, от которой желательно избавиться, то имеются задачи (относящиеся к анализу нелинейных систем), в которых шум играет полезную роль. Речь идет о том, что некоторые системы под воздействием определенного случайного шума могут приобретать черты устойчивого поведения.
В последние годы большое внимание уделяется проблемам стохастической финансовой математики.
Здесь для моделей с непрерывным временем вместо разностных уравнений возникают стохастические дифференциальные уравнения. При этом ценой значительных технических усложнений удается получить аналоги результатов, известных для дискретных моделей. Мы не приводим переформулировок "дискретных" определений, данных в дополнении к главе IV, на случай непрерывного времени (см. [86]). Заметим, что широкую известность получили исследования в рамках модели Блэка-Мертона-Шоулса, используюшие для описания эволюции цен геометрическое броуновское движение. Работы М. Шоулса и Р. Мертона были отмечены в 1997 г. Нобелевской премией по экономике ( . Блэк скончался, а посмертно Нобелевские премии не присуждаются).
Итак, если теперь -рынок описывается стохастическими дифференциальными уравнениями
где — банковский счет, а — цены акций, то имеет место следующий замечательный результат.
Теорема 20 (формула Блэка-Шоулса). В модели (128) с постоянными справедливая стоимость стандартного опциона-колл европейского типа с платежной функцией определяется формулой
где
а — функция распределения стандартной нормальной величины. В частности, при и имеем
при этом , когда .
Имеются разные подходы к доказательству этой теоремы (см., например, [86]). Отметим, что наряду с мартингальным подходом, опирающимся на теорему Гирсанова, имеется подход (который и был использован Блэком и Шоулсом), состояший в получении при определенных условиях для цены совершенного хеджирования
так называемого фундаментального уравнения
с краевым условием . При этом величина (см. (129)).
Полное представление о проблемах стохастической финансовой математики можно получить, ознакомившись с .