УДК 336.763
АНАЛИТИЧЕСКОЕ ВЫРАЖЕНИЕ СТОИМОСТИ СВОПЦИОННОГО КОНТРАКТА
В МОДЕЛИ ХАЛЛА-УАЙТА
О.А. Бельснер, Д.Н. Жабин
Томский политехнический университет E-mail: [email protected]
Исследована неарбитражная модель временной структуры процентной ставки с непрерывным временем, представленная Хал-лом (Hull) и Уайтом (White) в 1990 г. Рассматривается возможность применения этой модели для свопциона, предлагается вывод аналитического выражения для его оценки. Модель является наиболее распространенной и часто используемой в настоящее время при оценивании стоимости финансовых инструментов, производных от процентных ставок. При выводе выражения для свопциона применялись методы финансовой математики и теории опционов. Для проверки адекватности полученных результатов производится апробация предложенного выражения стоимости цены свопционного контракта на основе рыночной информации, основные расчеты и калибровка выражения реализованы в пакете Mathematica.
В настоящее время в банковской практике все большее внимание уделяется таким производным финансовым инструментам как своп, представляющий собой обмен на некоторый срок парой активов между контрагентами. Внебиржевая природа свопов дает практически неограниченные возможности для формирования контрактов, устраивающих обе стороны. На сегодняшний день на мировых внебиржевых и биржевых рынках преобладают процентные инструменты, в частности процентные свопы и основанные на них опционные контракты (свопционы).
Определение цен опционов является одним из важнейших результатов научных исследований в области современной финансовой математики, так как цены динамично меняются во времени и принятие на себя обязательства продать или купить что-то в будущем по оговоренной цене связано с определенным уровнем риска для обеих сторон. Модели ценообразования, которые были разработаны для опционов (в частности модели, Блэка-Шоулза [1]), успешно используются финансистами и трейдерами для выявления равновесны цен производных финансовых инструментов.
Определение стоимости синтетических инструментов, к которым принадлежит и свопцион, является более сложной задачей в силу того, что необходимо учитывать свойства лежащего в основе инструмента, в данном случае, свопа. Являясь инструментом внебиржевого рынка, а, следовательно, обладая определенной гибкостью и способностью приспосабливаться к запросам сторон, участвующих в сделке, свопцион характеризуется отсутствием каких-либо стандартных условий, что осложняет процесс определения их равновесной цены.
Исследование данного вида финансовых дери-вативов и выявление его особенностей представляется актуальным в связи с широким распространением и ростом рынков операций «своп» как в зарубежных странах, так и в Российской Федерации (так, 26 сентября 2002 г. Банк России ввел в действие механизм рефинансирования кредитных организаций с использованием сделок «своп»).
Предположение о том, что процентные ставки являются постоянными или известными функциями от времени [2] верны лишь в случае краткосрочных контрактов, тогда как для долгосрочных необходимо более точное моделирование поведения процентной ставки. Нами рассматривается неарбитражная однофакторная модель процентной ставки, в основе которой лежит принцип моделирования поведения процентной ставки. Предположим, что лишь один параметр является случайным, так как невозможно реалистично спрогнозировать поведение ставки в некоторый момент времени в будущем, естественно сделать предположение о том, что она является случайной величиной. Ставка процента, поведение которой описывается при помощи модели, отражает доход по облигации с бесконечным сроком погашения [3].
Традиционные модели процентной ставки, используемые в финансовой математике, называются равновесными. Они используются для выявления соотношений между различными показателями в экономике, однако существенным недостатком является тот факт, что структура процентной ставки -конечный результат данных моделей. При оценке де-ривативов необходимо использовать модели, соответствующие рыночной структуре процентной ставки [2]. Неарбитражные модели, позволяющие прогнозировать эволюцию структуры процентной ставки на основе информации об ее предыдущем поведении, были разработаны именно с этой целью. Одно-факторные модели процентной ставки позволяют оценить кривую дохода, основываясь на информации о процентной ставки и параметров выбранной модели. При использовании такого рода моделей необходимо определить критерии выбора параметров и степень доверия полученным результатам [2].
Выбор параметров, основанный на временных рядах, приводит к тому, что результат будеттеорети-ческим. В лучшем случае данная кривая совпадает с рыночной. Естественно, возникает вопрос: из каких предположений исходить при выборе параметров и что является более достоверным результатом: полученная теоретическая кривая или реальные
рыночные цены? При выборе, однако, не следует забывать о том, что облигации и свопы являются высоколиквидным инструментами [4].
В работе основное внимание было уделено выражениям для оценки стоимости облигации и выражениям для оценивания стоимости опционного контракта.
В течение последних пятнадцати лет предпринималось большое количество попыток описать движение кривой дохода с использованием однофак-торной модели. При этом используется традиционный подход, предполагающий существование модели краткосрочной процентной ставки, вывод кривой текущего дохода и вариантов ее изменения, при этом параметры модели максимально приближаются к реальным рыночным значениям [3]. Примерами такого подхода могут служить работы Дотана (Donathan) (1978 г.) [4], Кокса (Cox), Ингерсолла (Ingersoll) и Росса (Ross) (1985 г.) [4]. Исследователями были разработаны альтернативные подходы. Один из них был предложен Халлом и Уайтом (Hull/White) в 1990 г. [3, 4]. Ими рассматривались две однофакторные модели краткосрочной процентной ставки, которые могут быть применены как к волатильности дохода облигации, по которой процент уплачивается по истечении срока вместе с погашением основного долга, так и к структуре текущей процентной ставки. Нами рассматривается модель временной структуры процентной ставки. Следует отметить, что в исследуемой модели учитывается эффект возвращения к среднему (mean reversion): если процентная ставки отклоняется от некоторого своего среднего значения, то в дальнейшем проявляется тенденция возвращения этой ставки процента к среднему значению [4].
Считается, что цены активов изменяются случайным образом под воздействием так называемой эффективной рыночной гипотезы. Существует несколько различных формулировок данного предположения, однако в каждой из них выделяются два основных момента [2]:
1. Предшествующие события полностью отражаются на «сегодняшней» цене актива, что не несет дополнительной информации.
2. Рынки немедленно реагируют на появление любой информации об активах.
Таким образом, моделирование цен активов равнозначно моделированию новой информации, оказывающей воздействие на цену данного актива. Принимая во внимание предположения, сделанные выше, можно сделать вывод о том, что непредвиденные изменения цены актива схожи с марковским процессом [4].
При каждом изменении цены актива ожидается определенная прибыль, определяемая как частное от деления разности цен на первоначальное значение стоимости [3].
Модель Халла-Уайта является наиболее распространенной и часто используемой в настоящее вре-
мя для оценки стоимости финансовых инструментов, производных от процентных ставок (будущее значение процентной ставки в коротком периоде имеет логарифмически нормальное распределение). Для оценки стоимости финансовых инструментов, производных от процентных ставок, используются модели временной структуры процентных ставок с непрерывным временем. Временная структура процентных ставок определяется доход-ностями облигаций с нулевыми купонами при различных сроках до погашения.
Цена облигации в момент времени t в модели Халла-Уайта оценивается из уравнения [4]:
P(t,T) = A(t,T)e-B(tT)r(t\ где
1 - e
- a (T-t )
B (t,T ) = -
a
lnA(t,T) = in-BitJ) 9lnP(0t)
P(0, t )
dt
1 -аТ -Ш\2/ 2аг
-ТТ(е - е ) (е - !)> 4а
где Р(1,Т) - цена облигации, дата погашения которой наступает в момент времени Т, Б(1, Т) - стоимость (в момент времени ¡) облигации с нулевым купоном, погашаемой через (Т—) лет, Д/,Т) - номинал облигации.
Стохастическое дифференциальное уравнение для определения величины изменения процентной ставки облигации имеет вид [4]:
ёг (/) = [в^) - а(/)г (/)]Ж + (/),
где ст(/,Т) - волатильность Р(1), ф) - безрисковая процентная ставка в коротком периоде в момент времени ¡, dz(t) - винеровский процесс.
Функция в(/) выбирается таким образом, чтобы модель соответствовала начальной временной структуре (зависимость процентной ставки от срока действия финансового инструмента) [4].
Функции я(/) и ст(/) являются параметрами волатильности и выбираются, исходя из предположения о том, что они должны соответствовать рыночным ценам ряда активно торгуемых процентных опционов.
Цена в начальный момент времени опциона покупателя с датой окончания в момент времени Т на бескупонную облигацию, срок погашения которой наступает в момент времени д, составляет [4]: сь = ЬР(0,8)Ы (к) - ХР(0, Т) N (к -ар),
где Ь - номинал облигации, X- цена исполнения, ар - стандартное отклонение логарифма цены облигации в момент времени Т,
к = ±Ы >Р(0- + 0,, ар Р(0, Т) X 2
=-(1 - e-a(s-T))11 e
2a
Цена опциона продавца на облигацию составляет [4]:
pb = XP(0, T)N(-h + op) - LP(0, s)N(-h),
h = -Lln LP(0J.S) + О, N(z) =1 f,
О p
P(0, T) X
! dx.
Sw = E L
P(0, tt)(Fo N (dl) - RxN (d 2)),
dl =
d 2 =
ln(F0 /RX) +ct2T/2 ln(F0 /RX)-ct2T/2
oJT
F = E(R (t)), t e [T, T + n]
= d l -^л/f,
l F l - P(T,s) —F = ——
да
e p itt+да
Данные формулы эквивалентны выражениям модели Блэка-Шоулза [1] для оценки стоимости опционного контракта.
Рассмотрим, как модель Халла-Уайта может быть использована при оценки опционов на облигации.
Напомним, что опционный контракт дает право, но не обязывает купить или продать базовый инструмент по некоторой цене - цене исполнения в определенную будущую дату (дату истечения срока) или до ее наступления. Опцион на облигацию -опцион на покупку или продажу определенной облигации по определенной цене в определенную дату в будущем. Кроме покупки и продажи на внебиржевом рынке, опционы на облигации часто могут быть встроены в облигации при их эмиссии, что делает выпускаемые облигации более привлекательными как для эмитента, так и для потенциального покупателя [4].
Перейдем к выводу аналитического выражения для свопциона в модели Халла-Уайта. Напомним, что под свопционом понимается опцион на заключение процентного своп-соглашения в определенный момент времени в будущем при определенной фиксированной ставке [5]. Процентный своп может быть представлен в виде обмена облигаций с фиксированной процентной ставкой на облигацию с плавающей ставкой процента (цена такой облигации будет равна ее номиналу в начальный момент времени при заключении своп соглашения) [2].
Таким образом, свопцион может быть представлен в виде опциона на облигацию с фиксированной процентной ставкой и ценой исполнения равной номиналу (то есть бескупонной облигации). Стоимость свопциона сроком п лет, в таком случае, составит [4]:
здесь з=1+п - время до окончания своп-соглашения.
Так как свопцион, наделяющий владельца правом осуществлять платежи по плавающей процентной ставке, является опционом продавца на облигацию с фиксированной ставкой, цена исполнения которого равна плавающей ставке, можем получить выражение для свопциона в модели Халла-Уайта.
Запишем выражение для выплат по опциону продавца на своп, выплаты по которому производятся в момент времени Т+1/т, на момент времени Т. Предположим, что 8=1/т - временной интервал между начислениями процентов, тогда выражение для выплат можем записать в виде: ¿8
Ps (l + RS)
max(R - Rx ,0),
После несложных преобразований перейдем к записи:
PS = max
L - L(l + SRx ) 0
(1)
- стоимость в момент времени T бес-
1+Я8
¿(1 + ^8 )
1 + Я8
купонной облигации с доходом Х(1+Л15) в момент времени Т+8. Таким образом, выражение (1) характеризует доход на момент времени Т по опциону «пут» на своп, выплаты по которому производятся в момент времени Т+1/т, а цена исполнения равна Ь.
Таким образом, представив свопцион как портфель опционов на бескупонные облигации, запишем выражение для свопциона в модели Халла-Уайта:
тп
Бы = ¿X Р(0, ^+1)М (-И + ар) -
1=1
-Р(0, ^ )(1 + 8ЯХ) М(-И), = Т + т, (2) к = 1 1п(1 + 8*х) Р(0, ^) + 8
P(0, t-)
sp=Sо -* ЦЧа
N (z) ="/= f<
! dx.
где Ь - номинал, т - частота начисления процентов (в год), - форвардная ставка свопа, Р(0,/;) - бескупонная облигация со сроком погашения в момент времени 1==Т+1/т, /=(1,...,тп), Т- срок погашения свопа, Ях - фиксированная ставка процента.
Значение может быть найдено из уравнения [4]:
Вследствие необходимой корректной оценки стоимости инструмента необходимо определение параметров а и о, т.е. проведение калибровки модели (определение параметров модели) [4].
В представленной работе производилась калибровка модели к кэплету, далее рассчитывалась стоимость свопционного контракта, производилось сравнение с рыночной, и был сделан вывод о состоятельности формулы. При проведении калибровки модели использовалась информация о кэ-плетах на ставку LIBOR (London Interbank Offered Rate, ставку процента, по которому большинство крупных международных банков финансируют
i=l
свою деятельность) за 2002 г. Ставка исполнения составляет 4,96 %.
Напомним, что контракт, определяющий верхнее граничное значение процентной ставки, называется кэпом. Кэп состоит из серии кэплетов, представляющих собой контракты со сроками погашения платежа в момент времени Т+тна выплату разницы между рыночным значением ставки ЯТ в момент времени Т и значением ставки Х в течение периода т, основная сумма при этом составляет Ь. Главным условием является то, что данная разность должна быть положительна, то есть стоимость (ра-уо$) кэплета в момент Т+т составляет [2]:
шахЬ[(Л—Х),0].
На практике калибровка модели Халла-Уайта осуществляется путем выбора значений параметров, характеризующих поведение процентной ставки и волатильности таким образом, чтобы они соответствовали рыночным ценам опционов. Эмпирические значения параметра а принадлежат интервалу [0, 0,1], значение же параметра а принадлежат интервалу [0,01, 0,03] [5]. При калибровке начальные значения параметров а и а принимались равными 0,1 и 0,01 соответственно.
Параметры определяются при помощи метода наименьших квадратов таким образом, чтобы сум-
ма квадратов отклонений цен, полученных по модели Халла-Уайта, и рыночных цен опционов, была минимальной. В работе при выводе выражения применяли методы финансовой математики и теории опционов. Для проверки адекватности полученных результатов использовалась рыночная информация, основные расчеты и калибровка выражения были произведены с использованием пакета Mathematica.
В результате калибровки модели были получены следующие значения параметров: (а—>0,108114, о"—>0,0112018}. Зная параметры модели, можно рассчитать цену свопциона, дающего держателю право осуществлять платежи при ставке процента 6,2 % согласно условиям 3-летнего свопа, действие которого начинается через 5 лет. Волатильность процентной ставки свопа составляет 20 %. Платежи осуществляются раз в полугодие, основная сумма долга составляет 100 USD.
Применив (2), получим значение стоимости свопционного контракта - 2,0035. Рыночная стоимость данного свопционного контракта составляет 2,0038 (данные предоставлены компанией «Эконо-физика-Томск»). Таким образом, полученное выражение позволяет оценить стоимость свопционного контракта, при этом отклонение модельной цены от соответствующей рыночной не превышает 1 %.
СПИСОК ЛИТЕРАТУРЫ
1. Black F., Scholes M. The pricing of options and Corporate Liabilities // Journal of political economy. - 1973. - № 6. - P. 637-659.
2. Wilmott P. Derivatives: the theory and practice of financial engineering. - N.Y.: John Willey & Sons, 2000. - 742 p.
3. Hull J., White A. One-factor Interest-Rate Models and the Valuation of Interest-Rate Derivative Securities // The Journal of Financial and Quantitative analysis. - 1993. - № 2. - P. 235-240.
4. Hull J. Options, futures and other derivatives. - Toronto: University of Toronto, 2002. - 680 p.
5. Marshall J.F., Kapner K.P. Understanding swap finance. - OH: South-Western, 1990. - 250 p.
6. Keith R. An introduction to derivatives. - N.Y.: John Willey & Sons, 2000. - 198 p.
УДК 533.6.011
ИТЕРАЦИОННЫЙ МЕТОД РЕШЕНИЯ ОДНОМЕРНЫХ ДИВЕРГЕНТНЫХ УРАВНЕНИЙ ГАЗОВОЙ ДИНАМИКИ
В.М. Галкин
Томский политехнический университет E-mail: [email protected]
Предлагается итерационный метод решения дифференциальных уравнений, имеющих дивергентный вид и описывающих одномерное стационарное течение с переходом через скорость звука. Метод основан на использовании априорной информации о монотонном возрастании числа Маха вдоль сопла. Сравнение с другими методами проводилось на точном решении, а также при расчете двухфазного течения.
Введение
В [1] был предложен и далее в [2] уточнен итерационный метод решения одномерных стационарных уравнений газовой динамики, который за счет
использования априорной информации о монотонном изменении числа Маха внутри рассматриваемой области по скорости сходимости на порядок превосходит метод установления на основе схемы Маккормака [3] и имеет простой алгоритм. В