Модель волатильности SABR

редактировать

В математических финансах модель SABR - это стохастическая волатильность модель, которая пытается уловить улыбку волатильности на рынках производных финансовых инструментов. Название расшифровывается как «стохастический альфа, бета, ро », ссылаясь на параметры модели. Модель SABR широко используется практикующими специалистами в финансовой отрасли, особенно на рынках процентных производных. Его разработали Патрик С. Хэган, Дип Кумар, Эндрю Лесневски и Дайана Вудворд.

Содержание

  • 1 Динамика
  • 2 Асимптотическое решение
  • 3 SABR для отрицательных ставок
  • 4 Арбитражная задача в формула подразумеваемой волатильности
  • 5 Расширения
  • 6 Моделирование
  • 7 См. также
  • 8 Ссылки
  • 9 Внешние ссылки

Dynamics

Модель SABR описывает один форвардный F {\ displaystyle F}F , такой как LIBOR форвардный курс, курс форвардного свопа или форвардная цена акции. Это один из рыночных стандартов, используемых участниками рынка для определения волатильности. Неустойчивость форварда F {\ displaystyle F}F описывается параметром σ {\ displaystyle \ sigma}\ sigma . SABR - это динамическая модель, в которой оба F {\ displaystyle F}F и σ {\ displaystyle \ sigma}\ sigma представлены стохастическими переменными состояния, эволюция которых во времени задается следующей системой стохастические дифференциальные уравнения :

d F t = σ t (F t) β d W t, {\ displaystyle dF_ {t} = \ sigma _ {t} \ left (F_ {t} \ right) ^ {\ beta} \, dW_ {t},}{\ displaystyle dF_ {t} = \ sigma _ {t} \ left (F_ {t } \ righ t) ^ {\ beta} \, dW_ {t},}
d σ T = α σ td Z t, {\ displaystyle d \ sigma _ {t} = \ alpha \ sigma _ {t} ^ {} \, dZ_ {t}, }d \ sigma _ {t} = \ alpha \ sigma _ {t} ^ {{}} \, dZ_ {t},

с заданными значениями нулевого времени (наблюдаемыми в настоящее время) F 0 {\ displaystyle F_ {0}}F_{0}и σ 0 {\ displaystyle \ sigma _ {0}}\ sigma _ {0} . Здесь W t {\ displaystyle W_ {t}}W_{t}и Z t {\ displaystyle Z_ {t}}Z_{t}- это два коррелированных винеровских процесса с коэффициентом корреляции - 1 < ρ < 1 {\displaystyle -1<\rho <1}-1 <\ rho <1:

d W td Z t = ρ dt {\ displaystyle dW_ {t} \, dZ_ {t} = \ rho \, dt}{\ displaystyle dW_ {t} \, dZ_ {t} = \ rho \, dt}

Постоянные параметры β, α {\ displaystyle \ beta, \; \ alpha}\ beta, \; \ alpha удовлетворяют условиям 0 ≤ β ≤ 1, α ≥ 0 {\ displaystyle 0 \ leq \ beta \ leq 1, \; \ alpha \ geq 0}0 \ leq \ beta \ leq 1, \; \ alpha \ geq 0 . α {\ displaystyle \ alpha}\ alpha - параметр волатильности, подобный волатильности. ρ {\ displaystyle \ rho}\ rho - мгновенная корреляция между базовым активом и его волатильностью. α {\ displaystyle \ alpha}\ alpha , таким образом, контролирует высоту предполагаемого уровня волатильности банкомата. Корреляция ρ {\ displaystyle \ rho}\ rho управляет наклоном предполагаемого перекоса, а β {\ displaystyle \ beta}\ beta контролирует его кривизну.

Вышеупомянутая динамика является стохастической версией модели CEV с параметром асимметрии β {\ displaystyle \ beta}\ beta : фактически, он сводится к модели CEV, если α = 0 {\ displaystyle \ alpha = 0}\ alpha = 0 Параметр α {\ displaystyle \ alpha}\ alpha часто называют вольволом, и его значение заключается в логнормальной волатильности параметра волатильности σ {\ displaystyle \ sigma}\ sigma .

Асимптотическое решение

Мы рассматриваем Европейский опцион (скажем, колл) на форвард F {\ displaystyle F}F страйк на K {\ displaystyle K}K , который истекает T {\ displaystyle T}Tлет спустя. Значение этой опции равно подходящей дисконтированной ожидаемой величине выплаты max (FT - K, 0) {\ displaystyle \ max (F_ {T} -K, \; 0)}{\ displaystyle \ max (F_ {T} -K, \; 0)} в соответствии с распределением вероятностей процесса F t {\ displaystyle F_ {t}}F_t.

За исключением особых случаев β = 0 {\ displaystyle \ beta = 0}\ beta = 0 и β = 1 {\ displaystyle \ beta = 1}\ beta = 1 , выражение в закрытой форме для этого распределения вероятностей неизвестно. Общий случай может быть решен приближенно с помощью асимптотического разложения по параметру ε = T α 2 {\ displaystyle \ varepsilon = T \ alpha ^ {2}}\ varepsilon = T \ alpha ^ {2} . В типичных рыночных условиях этот параметр невелик, и приблизительное решение на самом деле довольно точное. Также важно то, что это решение имеет довольно простую функциональную форму, очень легко реализуется в компьютерном коде и хорошо подходит для управления рисками больших портфелей опционов в режиме реального времени.

Решение удобно выразить в терминах подразумеваемой волатильности σ impl {\ displaystyle \ sigma _ {\ textrm {impl}}}{\ displaystyle \ sigma _ {\ textrm {impl}}} варианта. А именно, мы переводим цену опциона по модели SABR в форму формулы оценки модели Блэка. Тогда подразумеваемая волатильность, которая представляет собой значение параметра логнормальной волатильности в модели Блэка, которое заставляет ее соответствовать цене SABR, приблизительно определяется следующим образом:

σ impl = α log ⁡ (F 0 / K) D (ζ) {1 + [2 γ 2 - γ 1 2 + 1 / (средний F) 2 24 (σ 0 C (средний F) α) 2 + ρ γ 1 4 σ 0 C (средний F) α + 2 - 3 ρ 2 24] ε}, {\ Displaystyle \ sigma _ {\ text {impl}} = \ alpha \; {\ frac {\ log (F_ {0} / K)} {D (\ zeta)}} \; \ left \ {1+ \ left [{\ frac {2 \ gamma _ {2} - \ gamma _ {1} ^ {2} + 1 / \ left (F _ {\ text {mid}} \ right) ^ {2} } {24}} \; \ left ({\ frac {\ sigma _ {0} C (F _ {\ text {mid}})} {\ alpha}} \ right) ^ {2} + {\ frac {\ rho \ gamma _ {1}} {4}} \; {\ frac {\ sigma _ {0} C (F _ {\ text {mid}})} {\ alpha}} + {\ frac {2-3 \ rho ^ {2}} {24}} \ right] \ varepsilon \ right \},}{\ displaystyle \ sigma _ {\ text {impl}} = \ alpha \; {\ frac {\ log (F_ {0} / K)} {D (\ zeta)}} \; \ left \ {1+ \ left [{\ frac {2 \ gamma _ {2} - \ gamma _ {1} ^ {2} + 1 / \ left (F _ {\ text {mid}} \ right) ^ {2 }} {24}} \; \ left ({\ frac {\ sigma _ {0} C (F _ {\ text {mid}})} {\ alpha}} \ right) ^ {2} + {\ frac { \ rho \ gamma _ {1}} {4}} \; {\ frac {\ sigma _ {0} C (F _ {\ text {mid}})} {\ alpha}} + {\ frac {2-3 \ rho ^ {2}} {24}} \ right] \ varepsilon \ right \},}

где для ясности мы установили C (F) = F β {\ displaystyle C \ left (F \ right) = F ^ {\ beta}}C \ left (F \ right) = F ^ { \ beta} . Значение F mid {\ displaystyle F _ {\ text {mid}}}F _ {{{\ text {mid }}}} обозначает удобно выбранную среднюю точку между F 0 {\ displaystyle F_ {0}}F_{0}и K {\ displaystyle K}K (например, среднее геометрическое F 0 K {\ displaystyle {\ sqrt {F_ {0} K}}}{\ sqrt {F_ {0} K}} или среднее арифметическое (F 0 + K) / 2 {\ displaystyle \ left (F_ {0} + K \ right) / 2}\ left (F_ {0} + K \ right) / 2 ). Мы также установили

ζ = α σ 0 ∫ KF 0 dx C (x) = α σ 0 (1 - β) (F 0 1 - β - K 1 - β), {\ displaystyle \ zeta = {\ frac {\ alpha} {\ sigma _ {0}}} \; \ int _ {K} ^ {F_ {0}} {\ frac {dx} {C (x)}} = {\ frac {\ alpha} {\ sigma _ {0} (1- \ beta)}} \; \ left (F_ {0} {} ^ {1- \ beta} -K ^ {1- \ beta} \ right),}{\ displaystyle \ zeta = {\ frac {\ alpha} {\ sigma _ {0}}} \; \ int _ {K} ^ {F_ {0}} {\ frac {dx} {C (x)}} = {\ frac {\ alpha} {\ sigma _ {0} (1- \ beta)}} \; \ left (F_ {0} {} ^ {1- \ beta } -K ​​^ {1- \ beta} \ right),}

и

γ 1 = C '(F mid) C (F mid) = β F mid, {\ displaystyle \ gamma _ {1} = {\ frac {C' (F _ {\ text {mid}})} {C (F _ {\ text {mid}})}} = {\ frac {\ beta} {F _ {\ text {mid}}}} \ ;,}{\displaystyle \gamma _{1}={\frac {C'(F_{\text{mid}})}{C(F_{\text{mid}})}}={\frac {\beta }{F_{\text{mid}}}}\;,}
γ 2 = C ″ (F mid) C (F mid) = - β (1 - β) (F mid) 2, {\ displaystyle \ gamma _ {2} = {\ frac {C '' (F _ {\ text {mid}})} {C (F_ {\ text {mid}})}} = - {\ frac {\ beta (1- \ beta)} {\ left (F _ {\ text {mid}} \ right) ^ {2}}} \ ;,}{\displaystyle \gamma _{2}={\frac {C''(F_{\text{mid}})}{C(F_{\text{mid}})}}=-{\frac {\beta (1-\beta)}{\left(F_{\text{mid}}\right)^{2}}}\;,}

Функция D (ζ) {\ displaystyle D \ left (\ zeta \ right)}D \ left (\ zeta \ right) , входящая в приведенную выше формулу, определяется как

D (ζ) = log ⁡ (1 - 2 ρ ζ + ζ 2 + ζ - ρ 1 - ρ). {\ Displaystyle D (\ zeta) = \ log \ left ({\ frac {{\ sqrt {1-2 \ rho \ zeta + \ zeta ^ {2}}} + \ zeta - \ rho} {1- \ rho }} \ right).}{\ displaystyle D (\ zeta) = \ log \ left ({\ frac {{\ sqrt {1-2 \ rho \ zeta + \ zeta ^ {2}}} + \ zeta - \ rho} {1- \ rho}} \ right).}

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

σ impl n = α F 0 - KD (ζ) {1 + [2 γ 2 - γ 1 2 24 (σ 0 C (F mid) α) 2 + ρ γ 1 4 σ 0 C (F mid) α + 2 - 3 ρ 2 24] ε}. {\ displaystyle \ sigma _ {\ text {impl}} ^ {\ text {n}} = \ alpha \; {\ frac {F_ {0} -K} {D (\ zeta)}} \; \ left \ {1+ \ left [{\ frac {2 \ gamma _ {2} - \ gamma _ {1} ^ {2}} {24}} \; \ left ({\ frac {\ sigma _ {0} C ( F _ {\ text {mid}})} {\ alpha}} \ right) ^ {2} + {\ frac {\ rho \ gamma _ {1}} {4}} \; {\ frac {\ sigma _ { 0} C (F _ {\ text {mid}})} {\ alpha}} + {\ frac {2-3 \ rho ^ {2}} {24}} \ right] \ varepsilon \ right \}.}{\ displaystyle \ sigma _ {\ text {impl}} ^ {\ text {n}} = \ alpha \; {\ frac {F_ {0} -K} {D (\ zeta)}} \; \ left \ {1+ \ l eft [{\ frac {2 \ gamma _ {2} - \ gamma _ {1} ^ {2}} {24}} \; \ left ({\ frac {\ sigma _ {0}) C (F _ {\ text {mid}})} {\ alpha}} \ right) ^ {2} + {\ frac {\ rho \ gamma _ {1}} {4}} \; {\ frac {\ sigma _ {0} C ( F _ {\ text {mid}})} {\ alpha}} + {\ frac {2-3 \ rho ^ {2}} {24}} \ right] \ varepsilon \ right \}.}

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

SABR для отрицательных ставок

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

d F t = σ t (F t + s) β d W t, {\ displaystyle dF_ {t} = \ sigma _ {t} (F_ {t} + s) ^ {\ бета} \, dW_ {t},}{\ displaystyle dF_ {t} = \ sigma _ {t} (F_ {t} + s) ^ {\ beta} \, dW_ {t},}
d σ t = α σ td Z t, {\ displaystyle d \ sigma _ {t} = \ alpha \ sigma _ {t} \, dZ_ {t},}{\ displaystyle d \ sigma _ {t} = \ alpha \ sigma _ {t} \, dZ_ {t },}

для некоторого положительного сдвига s {\ displaystyle s}s . Поскольку сдвиги включаются в рыночные котировки и существует интуитивно понятная мягкая граница того, как могут стать отрицательные ставки, смещенный SABR стал лучшей рыночной практикой для адаптации отрицательных ставок.

Модель SABR также может быть изменена для покрытия отрицательных процентных ставок следующим образом:

d F t = σ t | F t | β d W T, {\ Displaystyle dF_ {t} = \ sigma _ {t} | F_ {t} | ^ {\ beta} \, dW_ {t},}{\ displaystyle dF_ {t} = \ sigma _ {t} | F_ {t} | ^ {\ beta} \, dW_ { t},}
d σ t = α σ td Z t, {\ displaystyle d \ sigma _ {t} = \ alpha \ sigma _ {t} \, dZ_ {t},}{\ displaystyle d \ sigma _ {t} = \ alpha \ sigma _ {t} \, dZ_ {t },}

для 0 ≤ β ≤ 1/2 {\ displaystyle 0 \ leq \ beta \ leq 1/2}{\ displaystyle 0 \ leq \ beta \ leq 1/2} и свободное граничное условие для F = 0 {\ displaystyle F = 0}F = 0 . Доступны его точное решение для нулевой корреляции, а также эффективное приближение для общего случая.

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

Арбитражная проблема в формуле подразумеваемой волатильности

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

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

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

Расширения

Модель SABR может быть расширена, если предположить, что ее параметры зависят от времени. Однако это усложняет процедуру калибровки. Усовершенствованный метод калибровки зависящей от времени модели SABR основан на так называемых «эффективных параметрах».

Моделирование

Поскольку процесс стохастической волатильности следует за геометрическим броуновским движением, его точное моделирование несложно. Однако моделирование процесса форвардного актива - нетривиальная задача. Обычно рассматриваются схемы моделирования на основе Тейлора, такие как Эйлер – Маруяма или Мильштейн. Недавно были предложены новые методы для почти точного моделирования модели SABR методом Монте-Карло. Обширные исследования модели SABR недавно были рассмотрены в Cui et al. Для нормальной модели SABR (β = 0 {\ displaystyle \ beta = 0}\ beta = 0 без граничных условий в F = 0 {\ displaystyle F = 0}F = 0 ) известен метод моделирования в закрытой форме.

См. также

Ссылки

Внешние ссылки

Последняя правка сделана 2021-06-06 02:37:49
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте