Main menu

Метод Карра-Мадана: оценка опционов через быстрое преобразование Фурье (FFT)

Метод Карра-Мадана, предложенный в 1999 году, произвел революцию в вычислительных финансах, позволив оценивать целые спектры опционов за доли секунды с использованием аппарата быстрого преобразования Фурье (FFT). Этот подход особенно важен для моделей прыжковой диффузии и процессов Леви, где плотность вероятности цены неизвестна, но существует аналитическое выражение для ее характеристической функции.

Математическая суть метода заключается в переходе из пространства цен (spatial domain) в пространство частот (frequency domain). Пусть $x = ln(S_T)$ — логарифм будущей цены актива, а $q(x)$ — его риск-нейтральная плотность вероятности. Характеристическая функция определяется как преобразование Фурье плотности: $phi(u) = mathbb{E}[e^{iux}] = int_{-infty}^{infty} e^{iux} q(x) dx$. Для многих передовых моделей (Variance Gamma, CGMY, Heston) функция $phi(u)$ известна в точной алгебраической форме. Цена колл-опциона $C(K)$ выражается через интеграл от выплат по этой плотности, однако прямое преобразование Фурье функции выплат не интегрируемо (так как выплата $max(S_T - K, 0)$ растет до бесконечности при $S_T o infty$).

Гениальное решение Питера Карра и Дилипа Мадана состояло во введении экспоненциального демпфирующего фактора $alpha > 0$. Они рассмотрели модифицированную цену $c_T(k) = exp(alpha k) C(K)$, где $k = ln(K)$. Эта функция квадратично интегрируема, и для нее можно строго вычислить преобразование Фурье, которое элегантно выражается через известную характеристическую функцию $phi(u)$. Обратное преобразование Фурье позволяет получить исходную цену опциона: $C(K) = frac{exp(-alpha k)}{pi} int_0^infty e^{-ivk} psi(v) dv$, где $psi(v)$ — аналитически известная функция, зависящая от $phi(v - (alpha + 1)i)$.

Для численного расчета этого интеграла применяется алгоритм Быстрого преобразования Фурье (FFT). Если дискретизировать интеграл на сетке из $N$ точек (где $N$ — степень двойки, например, 4096), FFT позволяет вычислить интеграл одновременно для $N$ различных страйков за время порядка $O(N log N)$ операций вместо классического $O(N^2)$. Это означает, что аналитик получает цену сразу для всей "улыбки" опционов практически мгновенно. Выбор параметра демпфирования $alpha$ и шага сетки требует тонкой математической настройки для балансировки ошибок усечения (Truncation Error) и ошибок дискретизации (Aliasing Error). Внедрение FFT сделало возможным ежедневную промышленную калибровку многопараметрических стохастических моделей в любом инвестиционном банке.

Список литературы:
1. Carr P., Madan D. Option valuation using the fast Fourier transform. — Journal of Computational Finance, 1999.
2. Lewis A.L. Option Valuation under Stochastic Volatility: with Mathematica Code. — Finance Press, 2000.
3. Конт Р., Танков П. Финансовое моделирование с использованием процессов со скачками. — М.: ФИЗМАТЛИТ, 2004.
Эксперты: Питер Карр, Дилип Мадан, Алан Льюис.

Оценить
(0 votes)
Вверх

Соц. сети