автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.18, диссертация на тему:Численные методы расчета безарбитражных цен американских опционов в математических моделях финансовых рынков

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

Ромаданова Мария Михайловна

ЧИСЛЕННЫЕ МЕТОДЫ РАСЧЕТА БЕЗАРБИТРАЖНЫХ ЦЕН АМЕРИКАНСКИХ ОПЦИОНОВ В МАТЕМАТИЧЕСКИХ МОДЕЛЯХ ФИНАНСОВЫХ РЫНКОВ

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

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

005043396

1 7 МАЙ 2012

Великий Новгород - 2012

005043396

Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования "Санкт-Петербургский государственный архитектурно-строительный университет" на кафедре прикладной математики и информатики.

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

профессор Белопольская Яна Исаевна

Официальные оппоненты: доктор физико-математических наук

Смородина Наталья Васильевна (Санкт-Петербургский государственный университет, профессор)

кандидат физико-математических наук, старший научный сотрудник Солев Валентин Николаевич (Санкт-Петербургское отделение Математического института имени В.А. Стеклова РАН)

Ведущая организация: Федеральное государственное бюджетное

учреждение науки Институт проблем машиноведения РАН

Защита диссертации состоится 28 мая 2012 года в 13.30 на заседании диссертационного совета Д 212.168.04 при Новгородском государственном университете имени Ярослава Мудрого по адресу 173003, г. Великий Новгород, ул. Большая Санкт-Петербургская, д. 41.

С диссертацией можно ознакомиться в библиотеке Новгородского государственного университета имени Ярослава Мудрого.

Автореферат разослан " 2.& " апреля 2012г.

Ученый секретарь диссертационного совета Д 212.168.04,

кандидат физико-математических наук, доцент

орг

Токмачев М.С.

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

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

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

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

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

Для достижения этой цели необходимо решить следующие задачи:

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

2) вывести уравнения для расчета безарбитражных цен американских опционов в этих моделях;

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

4) на основе развитых в п. 3 методов создать алгоритмы расчета справедливых цен американских опционов;

5) разработать эффективные программы, реализующие эти алгоритмы, и проанализировать результаты численных экспериментов.

Методы исследования. Для решения поставленных задач использовались методы и результаты классической и современной теории уравнений в частных производных, теории стохастических уравнений, теории марковских цепей и марковских диффузионных и скачкообразных процессов. Программирование осуществлялось в среде МАТЬАВ и на языке С++.

На защиту выносятся следующие результаты:

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

2) численные методы расчета цен американских опционов ' в модифицированных моделях финансового рынка, основанные на вероятностных представлениях решений рассматриваемых задач;

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

Научная новизна диссертационного исследования состоит в следующем:

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

2) разработаны новые модификации численных методов построения решения краевых задач для ряда уравнений в частных производных, основанные на вероятностных представлениях соответствующих решений;

3) показано, что численная схема, построенная на основе вероятностного подхода, к решению задачи со свободной границей в модели Блэка-Шоулса дает результаты, близкие к результатам, полученным другими численными методами. Показана также эффективность разработанных в диссертации численных схем при нахождении цен американских опционов в более сложных моделях таких, как модели с переключениями, аффинные модели и модели типа Мертона;

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

параболического и эллиптического типа, а также для некоторого класса интегро-дифференциальных уравнений.

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

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

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

Апробация работы. Основные положения диссертационной работы докладывались и обсуждались на ряде научно-практических конференций:

- на 61-й, 62-й и 63-й Международной научно-технической конференции молодых ученых (аспирантов, докторантов) и студентов "Актуальные проблемы современного строительства" СПбГАСУ (2008, 2009, 2010гг.);

- на 66-й, 67-й и 68-й научной конференции профессоров, преподавателей, научных работников, инженеров и аспирантов университета СПбГАСУ (2009, 2010, 2011гг.);

- на Российской школе-конференции "Математика, информатика, их приложения и роль в образовании", Москва, РУДН, 14-18 декабря 2009г.;

- на международной конференции "Modern Stochastics: Theory and Applications II", Киев, Украина, 7-11 сентября 2010г.;

- на научной конференции преподавателей, аспирантов и студентов Новгородского государственного университета им. Ярослава Мудрого, Великий Новгород, март-апрель 2011г.;

- на международной конференции "Third Northern Triangular Seminar", Санкт-Петербург, 11-13 апреля 2011г.

Работа поддержана грантом для студентов, аспирантов, молодых ученых, молодых кандидатов наук вузов и академических институтов расположенных на территории Санкт-Петербурга. Диплом серия ПСП № 10577, номер гранта 2.1/02-06/016 (2010 год). Тема проекта: "Развитие численных методов, основанных на вероятностных подходах к решению задач оптимального управления и задач со свободной границей для параболических уравнений".

Также работа выполнена при финансовой поддержке СПбГАСУ, проект № 4-Ф-11 (2011-2012 гг.). Тема проекта: "Развитие вероятностных методов решения краевых задач для нелинейных задач математической физики и финансовой математики".

Публикации. Основные результаты диссертационной работы опубликованы в работах [1]-[Ю], в том числе 3 работы в научных журналах, рекомендованных перечнем ВАК.

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

КРАТКОЕ СОДЕРЖАНИЕ ДИССЕРТАЦИИ

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

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

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

А)^М= вир теТ[г,:г]

в)^м= ^ зиР ^Тг^-'Цу)

5*(т),теТ[,>г] Ь

где х^,» момент первого выхода процесса 5 на границу оптимального

исполнения 5 (б), бе [г,Т], которая определяется в процессе решения задачи, либо соотношениями С), О) вида

+ LF-rF = 0, t<T, s>S*(t),

D)

dF

at

F{t,s)=K-s, t<T, s < S*(t),

C)\F(t,s)>K-s, t<T,seR+, dF

- + LF-rF< 0, t <T, se R+, dt

F(T, s) = Ф(я) s max(A" - s, О) = [K - ,

(где L — генератор марковского процесса S(t), a S (t) интерпретируется как

неизвестная свободная граница), )р

- + LF-rF< 0, в {0,T)xR+, at

F(t,s)>0{s),

^ + Z.F-rFj(F-®)=0, в (0,Г)хй+,

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

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

Параграф 1.2 посвящен различным вариантам постановки задачи нахождения цены американского опциона в классической модели Блэка-Шоулса. Эти постановки используются в дальнейшем для создания различных вычислительных схем и алгоритмов, позволяющих находить численные решения.

Модель Блэка-Шоулса (Б-Ш) — это модель финансового рынка, состоящего из двух активов - безрискового и рискового, цены В(9) и S(9) которых удовлетворяют следующей системе уравнений

dB(Q) = rB(Q)dd, B(f) = 1, rfS'(8) = 5(0)He + o<Av(e)], S(t) = s, 0<t<Q<T, где w(i) — винеровский процесс, заданный на вероятностном пространстве

Расчет справедливой цены для американского колл-опциона в модели Б-Ш приводит к необходимости решать задачу Коши

~ + LBSF-rF = 0, F{T,s)=[s-K]+, at

г „ №12 2 д Р где = гх —- н—о 5 ——, а при рассмотрении пут-опциона или колл-

™ 2 э52

опциона на акции с выплатой дивидендов возникает необходимость решать задачу со свободной границей [1]. Здесь приведён также явный вид для цены бессрочного американского опциона в Б-Ш-модели.

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

В параграфе 1.3 рассматриваются постановки задачи о нахождении цен американских опционов в классической модели Мертона (М-модель), в которой цена базового рискового актива ."?(;) имеет как диффузионную, так и скачкообразную компоненту и удовлетворяет стохастическому уравнению

йК(е) = 5(9-) (г-5)Л9 + скЦ0)+ |(ег-1 , <>(/) =

** -I

где 0<?<9йТ, - пуассоновская случайная мера на с

параметром Ev(dQ,dz) = ¿&п(ск) и = V (¿/0,£/г)- £/0я(<&) -

компенсированная пуассоновская случайная мера.

В этом случае роль оператора Ь играет оператор Ьу [2] вида

г г 1 2 2^ Э^ ЬмР = -а я ——+ +

Эд

Эг

В параграфе 1.4 рассматривается рынок, на котором цены рисковых базовых активов (т.е. акций различных компаний) могут зависеть не только от функционирования самих компаний, но и от общего состояния рынка или от рейтинга компании. Другими словами, делается предположение, что динамика цен рисковых активов задается СДУ вида

. = к = \,...,с1, (1)

где _у(б)е {1,...,/,...,М} - марковская цепь с конечным числом состояний и генератором

и

П,ФМ= 2С/м(Я)[Ф(в,И1)-Ф(Я,/)]. (2)

т=1

При этом цена американского опциона с контрактной функцией Ф(л\/) удовлетворяет либо задаче Коши (для колл-опииона)

= 1 = \,...,М, (3)

т-1

Мт,з) = Ф/Ы, (4)

где

7=1 21=и>| М»;

либо соответствующей задаче со свободной границей. Вероятностное представление решения задачи (3), (4) задается теоремой 1.3.

Теорема 1.3. Пусть //(г,,у)г/(м,/) (0<?<7\ .геЛ^, /е {1.....-

классическое решение задачи (3), (4). Тогда оно допускает вероятностное представление вида

(5)

где процесс 5(0) удовлетворяет (1) и >>(9) - скачкообразный процесс с генератором П 5 вида (2).

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

Наконец, в параграфе 1.5 волатильность цены базового актива моделируется решением некоторого СДУ. С этой целью рассматривается аффинная модель финансового рынка (обобщенная модель Хестона), в которой динамика рискового актива относительно мартингальной меры Q задается следующей системой стохастических дифференциальных уравнений

&{Ъ) = Б{ъ)(гс}д + 4АЩ^(&1{д)), 5(/) = ^>0, 0<Г<13<7\ (6) Л(г>) = к, (О, - у(г>))^ + О] ¿Ле2(г)), 1>(/) = г>0. (7)

Здесь й^^е/г1, ^Ией1 - 2-винеровские процессы, = рсШ, г, Э], С[, А и С и коэффициент корреляции р -заданные константы. Будем считать, что г, к^, 0], — положительны, коэффициент корреляции | р | < 1, а константы Л и С удовлетворяют условию

(Ло,)2 <2к1(Л61+С), гарантирующему положительность С для 1<-&<Т.

Обозначим ¿j линейный оператор, действующий в пространстве Cq (л) по формуле

ds dv

2 Э52 2 Dv 3s3v

и пусть 5 (i,v) делит область G = {(i,s, v) с R+ xR+xR+} на две части

,00), где оптимально исполнение опциона, и

Q = [0,г]х[ S (t, v),o° ]х[0,°=>), где исполнение опциона не оптимально.

Задача нахождения безарбитражной цены американского пут-опциона в рассматриваемой модели, сформулированная как задача со

свободной границей, имеет вид -» р

— + LiF-rF = 0,{t,s,v)eQ, (8)

at

Л г->

—- + L\F — rF <0, (t,s,v)e N, (9)

dt

где Q, N задаются соотношениями

IF(t,s,v)= K-s, если S*(t,v)>s, тоесть seN,

F(t,s,v)> K-s, если S*(t,v)<s, тоесть seQ,

и S (i,v) - неизвестная функция, подлежащая определению в процессе решения задачи. Соответствующие начальные и граничные условия имеют вид lim = lim F(t,s,v) = K, lim F(i,s,v) = 0, (10)

lim F{t,s,v) = K-S*{t,v), lim fe[0,r]. (11)

Полезным при этом оказывается следующее утверждение. Лемма 1.7. Пусть двумерный случайный процесс S{d), v(-ö) (0 < t <Т) удовлетворяет системе стохастических дифференциальных уравнений (6)-(7). Тогда случайные процессы >'("0"), заданные соотношениями

S(ü) = exp(x(ü) + qy(i})), v(ü) = ^[(l + a2c2)y(^)-c], удовлетворяют

стохастическим уравнениям

dy{&) = к(0 - + о4уЩ dw2 (tf), y{t) = у,

где и'] (г), н'2 (?) - независимые виперовские процессы, а(у) = г-ак(0-у) —

2(1-р2)

к = кь в-(е1^ + С)(г-р2), а = Ао^1-р2 , а= Р 2 •

Лсг^-р )

При этом замена переменных

з = Кех™У, = V), (12)

ф(5<) = тах(К - 5,0) = К ш|| - ех+ау, о)= К[ф(х,.у)]+ сводит задачу (8)-(9) к задаче с постоянными коэффициентами по переменной

х и исключает смешанные производные.

Вторая глава посвящена развитию численных методов решения задачи Коши и задачи со свободной границей для параболических и интегро-дифференциальных уравнений и применению их к решению задачи нахождения справедливой цены американского опциона в Б-Ш-модели и М-модели.

В параграфах 2.1 и 2.2 ищется справедливая цена американского пут-опциона в Б-Ш-модели в постановке Б) как с помощью разностной схемы Бренана, так и с помощью метода подвижной границы, чтобы в дальнейшем сопоставить полученные численные решения с решениями, полученными с помощью численных методов, основанных на вероятностных подходах. Здесь же приведены явные выражения для цен бессрочных опционов американского типа.

В параграфе 2.3 рассматривается решение линейной двойственной задачи методом штрафов как в Б-Ш-модели, так и в М-модели. Поскольку в М-модели оператор Ьи уже не является эллиптическим оператором, то явное решение не удается получить даже для бессрочного опциона. Для получения численных решений в рамках метода штрафов линейная двойственная задача Б) сводится к нелинейному уравнению вида

где р{Р,Ф) - штрафное слагаемое и А = Ьвз или А = Ьм. Далее это

уравнение решается с помощью разностной схемы.

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

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

дифференциальных уравнений, что является важной и актуальной задачей. В этой главе используются методы и результаты работ Kappa [3], а также Боярченко и Левендорского [4, 5].

В параграфе 3.1 рассматривается задача нахождения цены F (s) бессрочного американского опциона на акции с выплатой дивидендов в Б-Ш-

модели. Цена F (s) задается соотношением F(s) = sup [ е_гс Ф (,S(x)) ],

теТ

которое в безразмерных координатах приобретает вид f(x)=f{x,h) = Ex[e~n>'^Y{4))]. Здесь фМ^-Цке*) и

Г(0)=х+ г-8-~а2 (в-/)+o[w(0)-w(i)] - диффузионный процесс с

генератором £r, h = 1п[^(т*)]— 1пК и xh = - момент оптимальной остановки процесса У (б). При вычислении f(x) важную роль играют

резольвентный оператор Щ процесса У(в), ~ Е

причем Rjt (с/ — Ly )ф = ф, и тесно связанные с ним операторы Zy, Z+, Z~ вида

ОО Q

О

2

где р , р~ - корни квадратного уравнения -—р2 -

¡е-1% {Y{t))dt Lo

2

■2 . (

P+g = 0.

S °2 r — o--

2

ч y

Используя соотношение 2Г/= 2*2 / = г~г+/, представляющее собой частный случай факторизации Винера-Хопфа [6] и справедливое для любой ограниченной измеримой функции /, доказывается следующее утверждение.

Лемма 3.1. Пусть функция g{x) = {r-Ly)^ непрерывна, монотонна (не убывает или не возрастает), изменяет знак и удовлетворяет некоторым дополнительным условиям. Тогда /{х,И) допускает представление

/М) = /іМ) = если < 0 <

= ЬЬ^г'^-^^+еХх), если ?(-~)>0>я(оо).

Используя теорему Дынкина и связи между введенными выше операторами можно найти оптимальное значение И как корень уравнения

т{х)=г я(х)= ек\х + и~)= 0 и явное выражение для функции /і(х,/г) л (х, ь)=г~1 (г+1{ь> =г-хе[і{11 „)(х+и+)т(х+и+),

I вида

где и± - экспоненциально распределенные случайные величины со значениями в [0,°°) и (-°°,0] соответственно с параметрами

V г.

Отметим еще одну важную интерпретацию операторов Zj-, Z+ и Z . Пусть Т - экспоненциально распределенная случайная величина со средним г~\ не зависящая от процесса V(t). Плотность распределения Т имеет вид ге~н и, следовательно, ZYg{x) = E[g{x + Y{T))]. Рассмотрим, наряду с процессом 7(i), также процессы F(f)= sup K(s) и Y(t) = inf Y(s). Известно,

что F(r) - это экспоненциально распределенная случайная величина на положительной полуоси со средним а У(Т) - это экспоненциально

1

распределенная случайная величина на отрицательной полуоси со средним —. Таким образом, действие операторов Z+ и Z~ можно представить в виде

Z+g(x)=E[g(x + r(r))] = rE Z-g{x) = E{g{x + Y{T))] = rE

\e-rtg{x + Y{t))dt

\e'rt g{x + Y{t))dt

Эта интерпретация и факторизация Винера-Хопфа для процесса Леви позволяют получить соответствующие выражения для цены бессрочного опциона и границы оптимального исполнения в модели Мертона.

Параграфы 3.2 и 3.3 посвящены расчету безарбитражных цен опционов американского типа с конечным сроком действия Т в Б-Ш-модели и М-модели. Как и в случае бессрочных опционов, рассматриваемая задача сводится к задаче со свободной границей для параболического уравнения с постоянными коэффициентами. Для решения этой задачи используется метод рандомизации Kappa [3], который представляет собой следующую трехшаговую процедуру.

Шаг 1. Рассмотрим американский опцион с конечным сроком действия Т и соответствующий ему американо-канадский опцион с датой исполнения х, где т - это экспоненциально распределенная случайная величина со средним

Ет = Х~1=Т. Пусть р\ (.y) = £[ е-гтФ(5(х)) ] - цена этого опциона.

Шаг 2. Найдем цену pn{s) американо-канадского опциона с датой

исполнения т", где т" - это сумма п независимых экспоненциально

распределенных случайных величин тк с параметрами Хк= — . При этом

п

распределение случайной величины т" = имеет вид

к=1

-h.

1 J («-1)! и = Т.

Шаг 3. Зададим p"{t,s) с помощью некоторой рекурентной процедуры, основанной на задании p"(s) и перейдем к пределу при о, чтобы

получить цену исходного американского опциона, имеющего силу до момента Г.

Для того, чтобы объяснить смысл третьего шага, заметим, что распределение случайной величины т" сходится к дельта функции Дирака, сосредоточенной в точке ЯГ1 = Т. При этом, если для непрерывной функции g(t) определить g„ (т) по формуле

(ntY-1

то получим lim gn{T)=g(T),

что и лежит в основе метода рандомизации

п—><*>

Kappa.

В Б-Ш-модели процедура нахождения цены американского опциона с конечным сроком действия имеет следующий вид. Пусть Т - фиксированная дата исполнения опциона. Разделим интервал [о,Г]

на п подинтервапов

0 = i0 <...<;„ =Г и обозначим A = tk+l-tk и q = А-1 + г. Аппроксимацию

границы hk оптимального исполнения на интервале {tk,tk+1) и

аппроксимацию /к функции f(tk,x) построим, используя обратную

индукцию. Для к = п положим fn{x) = §(x), а для к = п-\,п-2,... зададим /к как решение дискретизованного по времени варианта задачи в постановке С) (после перехода к безразмерным координатам) на интервале [tk,tk+J с оптимальной границей исполнения, определяемой так, чтобы максимизировать fk (х), т.е. рассмотрим задачу

(q-LY)fk{x) = A-'fk+\x)t x>hk, (]3)

fk(x) = ФМ, x<hk, (14)

а

г--

2

V

/х + — /хх - генератор марковского процесса У(/).

гдеІ7/= ^ 2

Отметим, что задачу (13), (14) можно получить в результате применения преобразования Лапласа-Карсона к задаче С) на интервале 1 ]что в

точности соответствует рандомизации момента исполнения опциона,

рассматриваемого на интервале

На интервале граница оптимального исполнения - это константа

Нк и процедура определения /к и Ик состоит в следующем: при к = п

зададим /" (х) = ф(х), а при к - п -1,п - 2,... пусть

/к{х) = Ея

О

где х^* - это оптимальный момент остановки, максимизирующий величину /к[х). В соответствии с теоремой Дынкина /к(х) удовлетворяет задаче (13),

(14), и х* — это первый момент попадания процесса У(?) в интервал ( ].

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

Теорема 3.5. Пусть функция g(x) = (г-1у )ф удовлетворяет некоторым

дополнительным условиям. Тогда для к-п — \,п — 2,...,0 и /к = — ф справедливы следующие утверждения:

1) функция / =2+(д-1/*+1-(г-£г)ф) является неубывающей функцией, имеющей единственный ноль при х = Ик, т.е. Ик - корень уравнения у^(дг) = 0;

2) хк — оптимальный момент остановки;

3) 7к ;

4)/к=/кя;

5) Функция не убывает и обращается в ноль при х < Ик.

Наконец, в 3.4 и 3.5 параграфах разработаны аналогичные алгоритмы для нахождения безарбитражных цен американских опционов, как бессрочных, так и с конечным сроком действия, в модели с переключениями и в обобщённой модели Хестона. В этом случае задача отыскания цены американского опциона сводится к решению задачи Коши или задачи со свободной границей для системы уравнений.

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

Разработанные в диссертации численные схемы и вычислительные алгоритмы, построенные на их основе, были реализованы в виде комплекса программ. Реализация этих схем, проведение численных экспериментов и полученные результаты представлены в диссертации в приложении В. При этом в приложении А приведен ряд основных фактов из теории процессов Леви и теории стохастических дифференциальных уравнений, которые были использованы при проведении исследования.

В качестве примера приведем здесь результаты расчета цены F(t,s,v) американского пут-опциона с контрактной функцией Ф(5) = в

обобщенной модели Хестона, которая задаётся уравнениями (6)-(7), со следующими параметрами: /- = 0,09, Щ=2, 9j=0,05, О\=0,2, А = 1,1, С = 0, р = -0Д, время действия опциона Г = 0,5 и цена исполнения X = 100. Программа выполнена в MATLAB.

БО

20

40

60

ЗХ 0.5

Рис.1. Цена американского пут-опциона

г,

0.2

0.5 0

Рис. 2. Цена американского пут-опциона

фиксированном значении волатильности).

^l(M) = F(i,.s,v)|v=0146 (при

F2(i,v) = F(f,i,v)|s=100 (при

фиксированной цене акции).

0.2 О

Рис. 3. Свободная граница S (t,v).

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ

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

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

3. Методы, разработанные в диссертации, были протестированы на модели Блэка-Шоулса, и было показано, что они дают результаты, совпадающие с результатами известных численных методов.

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

Список публикаций по теме диссертации

Публикации в журналах, входящих в перечень ВАК ведущих периодических изданий

1. Белопольская, Я. И. Вероятностный подход к задаче со свободной границей и расчет цен американских опционов / Я. И. Белопольская, М. М. Ромаданова // Записки научных семинаров ПОМИ. - СПб., 2010. -Т. 384. - С. 40-77.

2. Белопольская, Я. И. Вероятностные представления решений краевых задач для систем параболических уравнений / Я. И. Белопольская, М. М. Ромаданова // Вестник гражданских инженеров. - 2011. - № 4 (29). -С. 149-155.

3. Ромаданова, М. М. Расчет цен бессрочных американских опционов в обобщенной модели Хестона / М. М. Ромаданова // Вестник гражданских инженеров. -2012. - № 2 (31). - С. 350-354.

Публикации в других изданиях

4. Ромаданова, М. М. Расчет безарбитражных цен американских опционов / М. М. Ромаданова // Актуальные проблемы современного строительства : сб. материалов 61-й Международной научно-технической конференции молодых ученых / Санкт-Петербургский гос. архит.-строит. ун-т. — СПб., 2008. - Ч. IV. -С. 106-114.

5. Ромаданова, М. М. Расчет цен американских опционов в модели со скачками / М. М. Ромаданова // Актуальные проблемы современного строительства : сб. материалов 62-й Международной научно-технической

конференции молодых ученых / Санкт-Петербургский гос. архит.-строит. ун-т-СПб., 2009. - В 5 ч. Ч. I. - С. 117-122.

6. Белополъская, Я. И. Расчет цен американских опционов в моделях со скачками / Я. И. Белополъская, М. М. Ромаданова // Математика, информатика, их приложения и роль в образовании : тезисы докладов Российской школы конференции. - М.: РУДН, 2009. - С. 27-30.

7. Ромаданова, М. М. Расчет цен американских опционов в моделях финансовых рынков с диффузией и скачками / М. М. Ромаданова // Доклады 67-й научной конференции профессоров, преподавателей, научных работников, инженеров и аспирантов университета / Санкт-Петербургский гос. архит,-строит. ун-т. - СПб., 2010. - В 5 ч. Ч. IV. - С. 128-133.

8. Ромаданова, М. М. Расчет цен американских опционов в модели Хестона / М. М. Ромаданова // Актуальные проблемы современного строительства: 63-я Международная научно-техническая конференция молодых ученых / Санкт-Петербургский гос. архит.-строит. ун-т. - СПб., 2010. — В 3 ч. Ч. I.-C. 127-132.

9. Belopolskaya, Ya. American option pricing in various Financial market models / Ya. Belopolskaya, M. Romadanova // Abstracts of International Conference: Modern Stochastics: Theory and Applications II. - Kyiv: Taras Shevchenko National University of Kyiv, Ukraine, 2010. - P. 52.

10. Romadanova, M. American option pricing in stochastic volatility models / M. Romadanova // Programme and Abstracts of Third Northern Triangular Seminar. - St.Petersburg, 2011. - P. 17-19.

Список литературы

1. Detemple, J. American-Style Derivatives: Valuation and Computation / J. Detemple // New York: Chapman & Hall / CRC Financial Mathematics Series. -2005.-248 p.

2. Cont, R. Financial modelling with Jump Processes / R. Cont, P. Tankov // Chapman & Hall CRC Press. - 2009. - 543 p.

3. Carr, P. Randomization and the American put / P. Carr // Rev. Financ. Stud. -1998. - № 11 (3). - P. 597-626.

4. Boyarchenko, S. I. Non-Gaussian Merton-Black-Scholes Theory / S. I. Boyarchenko, S. Z. Levendorskii // Singapore: World Scientific. - 2002. -398 p.

5. Boyarchenko, S. I. Irreversible decisions under uncertainty: optimal stopping made easy / S. I. Boyarchenko, S. Z. Levendorskii // Berlin: Springer-Verlag. -2007.-283 p.

6. Sato, K. Levy processes and infinitely divisible distributions / K. Sato // Cambridge: Cambridge University Press. - 1999. - 486 p.

Ромаданова Мария Михайловна

ЧИСЛЕННЫЕ МЕТОДЫ РАСЧЕТА БЕЗАРБИТРАЖНЫХ ЦЕН АМЕРИКАНСКИХ ОПЦИОНОВ В МАТЕМАТИЧЕСКИХ МОДЕЛЯХ ФИНАНСОВЫХ РЫНКОВ

Подписано к печати 19.04.2012. Формат 60x84/16. Бумага офсетная. Печать офсетная. Гарнитура Times New Roman. Усл. печ. л.1. Тираж 100 экз. Заказ № 389.

Отпечатано в ЗАО «Новгородский технопарк». 173003, Великий Новгород, ул. Б. Санкт-Петербургская, 41. Тел. (816 2) 73-76-76.

61 12-1/860

САНКТ-ПЕТЕРБУРГСКИЙ ГОСУДАРСТВЕННЫЙ АРХИТЕКТУРНО-СТРОИТЕЛЬНЫЙ УНИВЕРСИТЕТ

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

Рол-

Ромаданова Мария Михайловна

ЧИСЛЕННЫЕ МЕТОДЫ РАСЧЕТА БЕЗАРБИТРАЖНЫХ ЦЕН АМЕРИКАНСКИХ ОПЦИОНОВ В МАТЕМАТИЧЕСКИХ МОДЕЛЯХ ФИНАНСОВЫХ РЫНКОВ

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

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

Научный руководитель:

доктор физико-математических наук,

профессор Белопольская Яна Исаевна

Санкт-Петербург - 2012

Оглавление

Введение 4

1 Модели финансовых рынков и задача отыскания безарбитражных цен американских опционов 25

1.1 Финансовые рынки и арбитраж, цены американских опционов . 25

1.2 Модель Блэка-Шоулса........................30

1.3 Модель с диффузией и скачками..................36

1.4 Модель с переключениями .....................39

1.5 Обобщенная модель Хестона....................47

2 Численные методы расчета цен американских опционов, основанные на решении уравнений в частных производных 52

2.1 Численный метод решения системы конечномерных неравенств 53

2.2 Метод подвижной границы.....................55

2.3 Метод штрафов............................62

2.3.1 Модель Блэка-Шоулса....................62

2.3.2 Модель финансового рынка с диффузией и скачками . . 70

3 Аналитические и численные методы расчета цен американских опционов, основанные на вероятностном подходе 75

3.1 Аналитические методы расчета цен бессрочных американских опционов................................75

3.1.1 Цена бессрочного американского онциона в модели Блэка-Шоулса ........................76

3.1.2 Цены бессрочных американских опционов в моделях типа Мертона..........................82

3.2 Расчет цен американских опционов с конечным сроком действия в Б-Ш-модели.........................86

3.3 Расчет цен американских опционов в модели с диффузией и скачками ...............................90

3.4 Расчет цены американского опциона в модели с переключениями 92

3.4.1 Безарбитражная цена бессрочного американского опциона в модели с переключениями..............94

3.4.2 Цена американского опциона с конечным сроком исполнения .............................97

3.5 Расчет цен американских опционов в модели Хестона......101

3.5.1 Расчет цены бессрочного американского иут-опциона . . 102

3.5.2 Американский пут-опцион с конечным сроком исполнения в модели Хестона . ...................109

Заключение 118

Литература 119

Приложение 126

А Необходимые сведения из теории стохастических дифференциальных уравнений 126

В Аналитические и численные результаты 133

Введение

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

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

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

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

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

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

Обзор литературы. Финансовая математика - одна из бурно развивающихся областей математики, находящаяся на стыке теории случайных процессов, теории уравнений в частных производных, численного анализа и ряда других. Теория арбитража или теория расчета справедливых цен финансовых инструментов составляет существенную часть этой области. К сожалению, на русском языке литература на эту тему не слишком обширна. Здесь следует упомянуть ряд монографий [1]- [6], представляющих собой введение в эту обширную область, а также некоторое количество статей например, [7]. С другой стороны, область эта активно развивается и количество монографий вышедших в мире и посвященных различным задачам теории арбитража исчисляется десятками, а число статей - сотнями. Упомянем здесь лишь наиболее близкие к теме диссертации монографии [8|- [18].

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

чи об оптимальной остановке для случайного процесса .5"(£), описывающего динамику базового актива, на который пишется опцион, по этой мере. Решение соответствующей задачи позволяет определить справедливую цену рассматриваемого опциона. Напомним, что справедливую цену называют также мартингальной или риск-нейтральной ценой.

В модели Блэка-Шоулса [19] (Б-Ш-модели) в такой постановке задача впервые была решена в работе Мертона [20], где был использован подход, предложенный Маккином [21] для решения соответствующей задачи для уравнения теплопроводности. Распространение этого подхода на более сложные современные модели рассматривалось в большом количестве работ, например, [22]- [26].

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

Описанные выше задачи можно решать как с использованием методов теории уравнений в частных производных, например, используя дискретизацию задачи или метод Фурье, так и с использованием чисто вероятностных вычислительных методов , в частности метода Монте Карло. Хотя в диссертационной работе в чистом виде последние методы не используются, однако, частично их идеи используются в главе 3.

1) Метод Фурье. Весьма эффективным при решении рассматриваемых задач является использование преобразования Фурье [31]- [34].

В частности, пусть Р{1, з) - цена европейского опциона с контрактной функцией , такой, что существует преобразование Фурье — ещ>(гхх)'ф{х)(1х, где Ф{х) = Ф(еа:) и г принадлежит некоторой полосе (зависящей от выбора Ф) в комплексной плоскости. Тогда

где У = 1п 5 + г(Т —£), ^ = и + т и (рх&, и) — Е схр(гиХ(1)) - характеристическая функция процесса Функция (рх&,и) при и £ К представляет собой преобразование Фурье плотности распределения случайного процесса Х{{) = 1п 3(1) и во многих случаях допускает аналитическое продолжение при и н-> г Е С регулярное в некоторой полосе ¿>х.

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

2) Метод Монте- Карло. Безарбитражная цена 5) американского опциона с контрактной функцией Ф(з) сроком жизни Т — I может быть определена соотношением

и, следовательно, её вычисление связано с вычислением среднего на траектории случайного процесса. Таким образом, для того чтобы воспользоваться методом Монте Карло, нужно, во-первых, уметь численно решать стохастические уравнения, моделирующие цены базовых активов и, кроме того, приближенно вычислять условные средние. Один из популярных алгоритмов метода Монте Карло может быть описан следующим образом. Пусть для простоты г = 0 и опцион может быть исполнен лишь в дискретные моменты времени £ = ¿о < ¿1 < ¿2 < • • • < ¿т = Т (иногда такое условие составляет часть контракта и такие контракты называют бермудскими). Пусть = £>(£,:) - значение случайного процесса в момент £,:, полученное, например, с помощью численного решения соответствующего стохастического дифференциального уравнения. Обратимся к формулировкам вычисления цены американского опциона, приведенным выше, и заметим, что они естественно связаны с методом динамического программирования. Пусть к(в, в) - зависящая от времени контрактная функция, в е [£, Т], к(Ь, в) = 0 и обозначает цену опциона в момент ^ при условии, что й'(^) = 5 и контракт не был исполнен ранее. Нас интересует цена контракта в момент ^ = ¿о - Определим ее с помощью

следующих рекурентных соотношений

Ъп(з) = М«), (1)

Е^(з) = тах{/г?_1(5), - §]}, г = 1,... ,т. (2)

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

= (3)

и наоборот, сопоставляя величину /¿(з) каждой паре {и, в) так. чтобы Р]п(,з) — Кт(в), мы определим момент остановки соотношением

т = гшп{г е {1,2,..., т] : Мй) > ^(«(й))}. (4)

Область исполнения £, определяемая значением цены ^(в) в г-й момент исполнения, - это множество вида {5 : /г,;(в) > ^-(б*)}, а область непрерывности С представляет собой ее дополнение. Момент остановки можно также задать как первый момент, в который цепь 5г- попадает в область исполнения. Цена С .¡{в) американского опциона в области С задается соотношением

<ЗД = = в], г = 1,...,т - 1 (5)

и может быть вычислена с помощью рекурентных соотношений

Ст = О, С;0) = ^[тах{^+1(5ш),(7г+1(^+1)}|5г = в], г = 0,...,т-1.

(6)

Цена опциона совпадает с величиной СДз) = СЬ(5о)-

Значения ^ определяют цену С^з) в области продолжения с помощью соотношений = = 5] и наоборот значения С^з) поз-

воляют задать Ц(в) = тах{/г,; (й),С\(з)},г = 1 ,...,т, так, что, зная

Сг(з) , мы получим

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

Для достижения этой цели необходимо решить следующие задачи:

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

2) вывести уравнения для расчета безарбитражных цен американских опционов в этих моделях:

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

4) на основе развитых в п. 3 методов создать алгоритмы расчета справедливых цен американских опционов;

5) разработать эффективные программы, реализующие эти алгоритмы, и проанализировать результаты численных экспериментов.

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

Методы исследования. Для решения поставленных задач использовались методы и результаты классической и современной теории уравнений в частных производных, теории стохастических уравнений, теории марковских цепей и марковских диффузионных и скачкообразных процессов. Программирование осуществлялось в среде МАТЬАВ и на языке С .

Научная новизна диссертационного исследования состоит в следующем:

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

2) разработаны новые модификации численных методов построения решения краевых задач для ряда уравнений в частных производных, основанные на вероятностных представлениях соответствующих решений:

3) показано, что численная схема, построенная на основе вероятностного подхода, к решению задачи со свободной границей в модели Блэка-Шоулса дает результаты, близкие к результатам, полученным другими численными методами. Показана также эффективность разработанных в диссертации численных схем при нахождении цен американских опционов в более сложных моделях таких, как модели с переключениями, аффинные модели и модели типа Мертона:

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

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

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

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

На защиту выносятся следующие результаты:

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

2) численные методы расчета цен американских опционов в модифицированных моделях финансового рынка, основанные на вероятностных представлениях решений рассматриваемых задач;

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

Апробация работы. Основные положения диссертационной работы докладывались и обсуждались �