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

Модель стохастической волатильности, используемая на рынках деривативов

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

Динамика

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

г Ф т = σ т ( Ф т ) β г Вт т , {\displaystyle dF_{t}=\sigma _{t}\left(F_{t}\right)^{\beta }\,dW_{t},}
г σ т = α σ т г З т , {\displaystyle d\сигма _{t}=\альфа \сигма _{t}^{}\,dZ_{t},}

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

г Вт т г З т = ρ г т {\displaystyle dW_{t}\,dZ_{t}=\rho \,dt}

Постоянные параметры удовлетворяют условиям . — параметр, подобный волатильности, для волатильности. — мгновенная корреляция между базовым активом и его волатильностью. Начальная волатильность контролирует высоту уровня подразумеваемой волатильности ATM . Как корреляция , так и контролирует наклон подразумеваемого перекоса. Волатильность волатильности контролирует его кривизну. β , α {\displaystyle \бета,\;\альфа} 0 β 1 , α 0 {\displaystyle 0\leq \beta \leq 1,\;\alpha \geq 0} α {\displaystyle \альфа} ρ {\displaystyle \ро} σ 0 {\displaystyle \сигма _{0}} ρ {\displaystyle \ро} β {\displaystyle \бета} α {\displaystyle \альфа}

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

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

Мы рассматриваем европейский опцион (скажем, колл) на форвард , выполненный в , который истекает через несколько лет. Стоимость этого опциона равна соответствующим образом дисконтированному ожидаемому значению выплаты при распределении вероятностей процесса . Ф {\displaystyle F} К {\displaystyle К} Т {\displaystyle Т} макс ( Ф Т К , 0 ) {\displaystyle \max(F_{T}-K,\;0)} Ф т {\displaystyle F_{t}}

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

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

σ импл = α бревно ( Ф 0 / К ) Д ( ζ ) { 1 + [ 2 γ 2 γ 1 2 + 1 / ( Ф середина ) 2 24 ( σ 0 С ( Ф середина ) α ) 2 + ρ γ 1 4 σ 0 С ( Ф середина ) α + 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\},}

где для ясности мы установили . Формула не определена, когда , поэтому мы заменяем ее пределом как , который получается путем замены множителя на 1. Значение обозначает удобно выбранную среднюю точку между и (например, среднее геометрическое или среднее арифметическое ). Мы также установили С ( Ф ) = Ф β {\displaystyle C\left(F\right)=F^{\beta }} К = Ф 0 {\displaystyle К=Ф_{0}} К Ф 0 {\displaystyle K\to F_{0}} бревно ( Ф 0 / К ) Д ( ζ ) {\displaystyle {\frac {\log(F_{0}/K)}{D(\zeta)}}} Ф середина {\displaystyle F_{\text{середина}}} Ф 0 {\displaystyle F_{0}} К {\displaystyle К} Ф 0 К {\displaystyle {\sqrt {F_{0}K}}} ( Ф 0 + К ) / 2 {\displaystyle \left(F_{0}+K\right)/2}

ζ = α σ 0 К Ф 0 г х С ( х ) = α σ 0 ( 1 β ) ( Ф 0 1 β К 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),}

и

γ 1 = С ( Ф середина ) С ( Ф середина ) = β Ф середина , {\displaystyle \gamma _{1}={\frac {C'(F_{\text{mid}})}{C(F_{\text{mid}})}}={\frac {\beta }{F_{\text{mid}}}}\;,}
γ 2 = С ( Ф середина ) С ( Ф середина ) = β ( 1 β ) ( Ф середина ) 2 , {\displaystyle \gamma _{2}={\frac {C''(F_{\text{mid}})}{C(F_{\text{mid}})}}=-{\frac {\beta (1-\beta )}{\left(F_{\text{mid}}\right)^{2}}}\;,}

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

Д ( ζ ) = бревно ( 1 2 ρ ζ + ζ 2 + ζ ρ 1 ρ ) . {\displaystyle D(\zeta )=\log \left({\frac {{\sqrt {1-2\rho \zeta +\zeta ^{2}}}+\zeta -\rho }{1-\rho }}\right).}

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

σ импл н = α Ф 0 К Д ( ζ ) { 1 + [ 2 γ 2 γ 1 2 24 ( σ 0 С ( Ф середина ) α ) 2 + ρ γ 1 4 σ 0 С ( Ф середина ) α + 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\}.}

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

Точность аппроксимации и степень арбитража можно дополнительно улучшить, если использовать эквивалентную волатильность в рамках модели CEV с той же самой для ценообразования опционов. [2] β {\displaystyle \бета}

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

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

г Ф т = σ т ( Ф т + с ) β г Вт т , {\displaystyle dF_{t}=\sigma _{t}(F_{t}+s)^{\beta }\,dW_{t},}
г σ т = α σ т г З т , {\displaystyle d\sigma _{t}=\alpha \sigma _{t}\,dZ_{t},}

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

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

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

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

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

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

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

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

Расширения

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

В качестве альтернативы Герреро и Орландо [7] показывают, что зависящая от времени локальная стохастическая волатильность (SLV) может быть сведена к системе автономных PDE, которые могут быть решены с использованием теплового ядра с помощью метода факторизации Вей-Нормана и алгебраических методов Ли. Явные решения, полученные с помощью указанных методов, сопоставимы с традиционным моделированием Монте-Карло, что позволяет сократить время численных вычислений.

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

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

Смотрите также

Ссылки

  1. ^ Хаган, Патрик С.; Кумар, Дип; Кесневски, Эндрю С.; Вудворд, Диана Э. (январь 2002 г.). «Управление риском улыбки» (PDF) . Уилмотт . Том 1. стр.  84–108 . Архивировано (PDF) из оригинала 30.04.2022 . Получено 30.04.2022 .
  2. ^ Чой, Джэхёк; Ву, Лисинь (июль 2021 г.). «Эквивалентная константа-эластичность-дисперсии (CEV) волатильности стохастической-альфа-бета-ро (SABR) модели». Журнал экономической динамики и управления . 128 : 104143. arXiv : 1911.13123 . doi : 10.1016/j.jedc.2021.104143. S2CID  235239799. SSRN  3495464. Получено 30.04.2022 .
  3. ^ Антонов, Александр; Коников, Майкл; Спектор, Майкл (2015-01-28). «Свободная граница SABR: естественное расширение до отрицательных ставок». SSRN  2557046.
  4. ^ Grzelak, Lech A.; Oosterlee, Cornelis W. (февраль 2017 г.) [2016-07-04]. «От арбитража к безарбитражной подразумеваемой волатильности». Journal of Computational Finance . 20 (3): 31– 49. doi :10.21314/JCF.2016.316. ISSN  1755-2850. SSRN  2529684. Получено 2022-04-30 .
  5. ^ Ле Флок, Фабьен; Кеннеди, Гэри (2016-08-15). "Методы конечных разностей для SABR без арбитража". Журнал вычислительных финансов . ISSN  1755-2850 . Получено 2022-04-30 .
  6. ^ Ван дер Стоп, Антон В.; Гржелак, Лех Александр; Остерли, Корнелис В. (2015-09-28). «Зависящая от времени модель FX-SABR: эффективная калибровка на основе эффективных параметров». Международный журнал теоретических и прикладных финансов . 18 (6): 1550042. doi :10.1142/S0219024915500429. SSRN  2503891. Получено 30 апреля 2022 г.
  7. ^ Герреро, Хулио; Орландо, Джузеппе (сентябрь 2021 г.). «Модели стохастической локальной волатильности и метод факторизации Вей-Нормана». Дискретные и непрерывные динамические системы — S . 15 (12): 3699– 3722. arXiv : 2201.11241 . doi :10.3934/dcdss.2022026. ISSN  1937-1632. S2CID  246295004 . Получено 30.04.2022 .
  8. ^ Лейтао, Альваро; Гржелак, Лех А.; Остерли, Корнелис В. (2017-04-10) [2016-04-13]. «Об эффективном многошаговом моделировании Монте-Карло модели SABR». Количественные финансы . 17 (10): 1549– 1565. doi : 10.1080/14697688.2017.1301676 . SSRN  2764908.
  9. ^ Cui, Zhenyu; Kirkby, Justin L.; Nguyen, Duy (2018-04-24). «Общая структура оценки для SABR и моделей стохастической локальной волатильности». Журнал SIAM по финансовой математике . 9 (2): 520–563 . doi :10.1137/16M1106572. S2CID  207074154.
  10. ^ Чой, Джэхёк; Лю, Ченру; Со, Бён Ки (2018-10-31). «Модель гиперболической нормальной стохастической волатильности». Журнал фьючерсных рынков . 39 (2): 186–204 . arXiv : 1809.04035 . doi : 10.1002/fut.21967. S2CID  158662660. SSRN  3068836. Получено 2022-04-30 .

Дальнейшее чтение

  • Хаган, Патрик; Лесневски, Эндрю; Вудворд, Диана (2005-03-22). "Распределение вероятностей в модели SABR стохастической волатильности" (PDF) . Архивировано из оригинала (PDF) 2021-03-08 . Получено 2022-04-30 .
  • Бартлетт, Брюс (февраль 2006 г.). "Хеджирование в рамках модели SABR" (PDF) . Уилмотт . Архивировано из оригинала (PDF) 2020-12-30 . Получено 2022-04-30 .
  • Хаган, Патрик; Лесневски, Эндрю (2008-04-30). "Модель рынка LIBOR со стохастической волатильностью в стиле SABR" (PDF) . Архивировано из оригинала (PDF) 2022-03-03 . Получено 2022-04-30 .
  • Хаган, Патрик С.; Кумар, Дип; Лесневски, Эндрю С.; Вудворд, Диана Э. (29.01.2014). «SABR без арбитража». Wilmott . Vol. 2014, no. 69. pp.  60–75 . doi :10.1002/wilm.10290 . Получено 30.04.2022 .
  • Облой, Ян (18 марта 2008 г.). «Настройте свою улыбку – поправка к Хагану и др.». arXiv : 0708.0998 [q-fin.CP].
  • Уэст, Грэм. "Краткое изложение подходов к модели SABR для производных инструментов на акции smile". Riskworx . Архивировано из оригинала 2015-09-14 . Получено 2022-04-30 .
  • Анри-Лабордер, Пьер (15.02.2005). «Объединение моделей BGM и SABR: краткий экскурс в гиперболическую геометрию». arXiv : physics/0602102 .
  • Джордан, Ричард; Тир, Чарльз (17.05.2011). «Асимптотические приближения моделей CEV и SABR». SSRN  1850709.
  • "Калибровка SABR". 2012-12-26. Архивировано из оригинала 2016-05-27 . Получено 2022-04-30 .
  • Антонов, Александр; Спектор, Майкл (2012-03-23). ​​«Расширенная аналитика для модели SABR». SSRN  2026350.
  • Антонов, Александр; Коников, Майкл; Спектор, Майкл (2019-05-02). Современная аналитика SABR: формулы и идеи для квантов, бывших физиков и математиков (Springer Briefs in Quantitative Finance) 1-е изд . doi : 10.1007/978-3-030-10656-0. ISBN 978-3-030-10655-3. ISSN  2192-7014. S2CID  182484805./ Пресс-релиз, Нью-Йорк, штат Нью-Йорк – 24 июня 2019 г.
  • Гулисашвили, Арчил; Хорват, Бланка; Жакье, Антуан (Джек) (2016-11-22). "Масса в нуле в некоррелированной модели SABR и асимптотика подразумеваемой волатильности". SSRN  2563510.
Retrieved from "https://en.wikipedia.org/w/index.php?title=SABR_volatility_model&oldid=1245080983"