2025-11-10T02:49:03.701673

Efficient simulation of prices for European call options under Heston stochastic-local volatility model: a comparison of methods

cai, Li
The Heston stochastic-local volatility model, consisting of a asset price process and a Cox--Ingersoll--Ross-type variance process, offers a wide range of applications in the financial industry. The pursuit for efficient model evaluation has been assiduously ongoing and central to which is the numerical simulation of CIR process. Different from the weakly convergent noncentral chi-squared approximation used in 25, this paper considers two strongly convergent and positivity-preserving methods for CIR process under Lamperti transformation, namely, the truncated Euler method and the backward Euler method. It should be noted that these two methods are completely different. The explicit truncated Euler method is computationally effective and remains robust under high volatility, while the implicit backward Euler method provides high computational accuracy and stable performance. Numerical experiments on European call options are presented to show the superiority of different methods.
academic

Эффективное моделирование цен европейских опционов колл в модели Хестона со стохастической-локальной волатильностью: сравнение методов

Основная информация

  • ID статьи: 2509.24449
  • Название: Efficient simulation of prices for European call options under Heston stochastic-local volatility model: a comparison of methods
  • Авторы: Мэн Цай, Тяньцзе Ли
  • Аффилиация авторов: Школа статистики и математики, Школа экономики, Центральный университет финансов и экономики
  • Классификация: q-fin.CP (вычислительные финансы)
  • Дата публикации: 16 октября 2025 г.
  • Ссылка на статью: https://arxiv.org/abs/2509.24449

Аннотация

В данной работе исследуется модель Хестона со стохастической-локальной волатильностью (HSLV), состоящая из процесса цены актива и процесса дисперсии типа Кокса-Ингерсолла-Росса, которая широко применяется в финансовой индустрии. В отличие от слабо сходящегося приближения нецентральным хи-квадратом, используемого в литературе 25, в данной работе рассматриваются два метода, которые обеспечивают сильную сходимость и сохранение положительности процесса CIR при преобразовании Лампертиса: метод усеченного Эйлера и метод обратного Эйлера. Явный метод усеченного Эйлера вычислительно эффективен и остается устойчивым при высокой волатильности, тогда как неявный метод обратного Эйлера обеспечивает высокую вычислительную точность и стабильность. Численные эксперименты на примере европейских опционов колл демонстрируют преимущества различных методов.

Исследовательский контекст и мотивация

Проблемный фон

  1. Требования к модели: Финансовая инженерия требует достижения баланса между совершенной калибровкой и корректной динамикой; модели со стохастической-локальной волатильностью (SLV) предоставляют решение этой проблемы
  2. Ограничения существующих моделей:
    • Чистые модели стохастической волатильности имеют ограниченные параметры и не могут идеально подогнать подразумеваемую волатильность для всех сроков и цен исполнения
    • Модели локальной волатильности, хотя и могут калибровать все текущие рыночные опционы, имеют проблемы с прогнозированием будущего поведения волатильности, что приводит к уплощению улыбки форвардной волатильности

Основные вызовы

Численное моделирование процесса CIR является ключевой сложностью при оценке модели, поскольку:

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

Исследовательская мотивация

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

Основные вклады

  1. Предложение двух новых численных методов: метода усеченного Эйлера и метода обратного Эйлера, оба обладают сильной сходимостью и сохранением положительности
  2. Теоретический анализ: доказательство сильной сходимости метода усеченного Эйлера (порядок сходимости 1/2)
  3. Комплексное численное сравнение: систематическое сравнение производительности различных методов через эксперименты по оценке цены европейских опционов колл
  4. Практические рекомендации: предоставление рекомендаций по выбору метода для различных сценариев применения

Подробное описание методов

Определение задачи

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

Архитектура модели

Определение модели HSLV

dS_t = rS_t dt + \sqrt{V_t} S_t \sigma(t, S_t)\left(\rho dW_t + \sqrt{1-\rho^2} d\tilde{W}_t\right) \\ dV_t = \kappa(\theta - V_t)dt + \gamma\sqrt{V_t} dW_t \\ dW_t d\tilde{W}_t = \rho dt \end{cases}$$ где: - $S_t$: процесс цены актива - $V_t$: процесс дисперсии (динамика CIR) - $\sigma(t,s)$: функция локальной волатильности - $\kappa, \theta, \gamma$: параметры процесса CIR #### Преобразование Лампертиса Для преодоления численных трудностей процесса CIR вводится преобразование $L_t = \sqrt{V_t}$: $$dL_t = \frac{1}{2}\kappa\left(\frac{\theta}{L_t} - L_t\right)dt + \frac{1}{2}\gamma dW_t$$ ### Основные численные методы #### 1. Метод усеченного Эйлера Определение функции усечения: $\pi_\tau(x) = (b\tau^{1/4}) \vee x$ Численная схема: $$L^N_{t_{n+1}} = L^N_{t_n} + \frac{\kappa}{2}\tau\left(\frac{\theta}{\pi_\tau(L^N_{t_n})} - \pi_\tau(L^N_{t_n})\right) + \frac{1}{2}\gamma\Delta W_n$$ **Характеристики**: - Явная схема, вычислительно эффективна - Гарантирует положительность через усечение - Остается устойчивой при высокой волатильности #### 2. Метод обратного Эйлера $$L^N_{t_{n+1}} = L^N_{t_n} + \frac{\kappa}{2}\tau\left(\frac{\theta}{L^N_{t_{n+1}}} - L^N_{t_{n+1}}\right) + \frac{1}{2}\gamma\Delta W_n$$ **Характеристики**: - Неявная схема, требует итеративного решения - Обеспечивает высокую вычислительную точность - Стабильная производительность, но высокие вычислительные затраты ### Технические инновации 1. **Сильная сходимость**: в отличие от слабо сходящегося метода нецентрального хи-квадрата, гарантирует сходимость на уровне траекторий 2. **Проектирование сохранения положительности**: обеспечение положительности дисперсии через математические преобразования и методы усечения 3. **Взаимодополняемость методов**: два метода имеют различные преимущества в эффективности и точности, применимы к различным сценариям ## Экспериментальная установка ### Параметры модели - **Базовые параметры Хестона**: $\gamma = 0.95$, $\kappa = 1.05$, $\rho = -0.315$, $\theta = 0.0855$, $S_0 = 1$, $r = 0$, $v_0 = 0.0945$, $T = 5$ - **Параметры коррекции HSLV**: введение параметра коррекции $p = 0.25$ для корректировки различных параметров ### Показатели оценки - Процент относительной ошибки - Время вычисления - Устойчивость при различных рыночных условиях ### Методы сравнения - Стандартный метод Эйлера - Метод AES (Andersen Exact Simulation) - Метод усеченного Эйлера - Метод обратного Эйлера ### Дизайн экспериментов 1. **Расчет условного математического ожидания**: использование 2D-COS метода в качестве эталона, 20 бинов, 10 000 траекторий, шаг 0.001 2. **Анализ чувствительности параметров**: тестирование производительности при различных значениях $\theta$ и $p$ 3. **Сравнение цен исполнения**: три случая $K = 70\%, 100\%, 150\%$ ## Результаты экспериментов ### Основные результаты #### Сравнение ошибок при различном числе шагов | Число шагов N | K=70% | K=100% | K=150% | |-------|--------|---------|---------| | **Метод усеченного Эйлера** | | 5 | 13.49%(0.24%) | 14.08%(0.24%) | 13.96%(0.85%) | | 40 | 4.66%(0.53%) | 3.86%(0.74%) | 4.80%(1.35%) | | **Метод обратного Эйлера** | | 5 | 15.49%(0.28%) | 15.97%(0.10%) | 17.32%(1.58%) | | 40 | 0.06%(0.66%) | 0.46%(0.68%) | 1.44%(1.15%) | #### Сравнение времени вычисления (40 шагов) - Метод AES: 0.1114 сек - Метод Эйлера: 0.0939 сек - Метод усеченного Эйлера: 0.0733 сек - Метод обратного Эйлера: 137.9476 сек ### Ключевые выводы 1. **При малом числе шагов**: метод усеченного Эйлера обычно имеет наименьшую ошибку 2. **При большом числе шагов**: метод обратного Эйлера обеспечивает наивысшую точность 3. **Устойчивость параметров**: метод усеченного Эйлера показывает медленный рост ошибки при увеличении $\theta$, демонстрируя превосходную устойчивость 4. **Вычислительная эффективность**: метод усеченного Эйлера вычислительно самый быстрый, метод обратного Эйлера требует значительно больше времени из-за необходимости итеративного решения ### Абляционные эксперименты #### Влияние параметра волатильности $\theta$ Эксперименты показывают, что с увеличением $\theta$ (моделирование экстремальных рыночных условий): - Ошибки методов AES и стандартного Эйлера значительно возрастают - Метод усеченного Эйлера остается устойчивым с низким темпом роста ошибки - Это указывает на большую применимость метода усечения при стресс-тестировании и моделировании хвостовых рисков #### Влияние параметра коррекции $p$ С увеличением $p$ метод обратного Эйлера постоянно показывает лучшие результаты, но имеет более высокую чувствительность к размеру шага. ## Связанные работы ### Основные направления исследований 1. **Калибровка SLV моделей**: пионерская работа Гюйона и Анри-Лабордера (2011) 2. **Численные методы для процесса CIR**: включая методы конечных разностей, методы частиц Монте-Карло и др. 3. **Применение глубокого обучения**: использование генеративных состязательных сетей для калибровки (Кучиеро и др., 2020) ### Преимущества данной работы По сравнению с существующими работами, данная статья предоставляет методы с сильной сходимостью и сохранением положительности, заполняя пробел слабо сходящихся методов и предоставляя более надежные численные инструменты для практического применения. ## Заключение и обсуждение ### Основные выводы 1. **Взаимодополняемость методов**: метод усеченного Эйлера подходит для крупномасштабного моделирования и управления рисками в реальном времени; метод обратного Эйлера подходит для стресс-тестирования и сценариев с высокими требованиями к точности 2. **Преимущества устойчивости**: метод усеченного Эйлера показывает отличную производительность в экстремальных рыночных условиях 3. **Компромисс точность-эффективность**: выбор подходящего размера шага и метода требует учета конкретных требований приложения ### Ограничения 1. Высокие вычислительные затраты метода обратного Эйлера 2. Выбор размера шага требует компромисса между точностью и эффективностью 3. Для некоторых экстремальных комбинаций параметров стабильность методов требует дальнейшей проверки ### Будущие направления 1. Разработка стратегий адаптивного размера шага 2. Исследование параллельной реализации 3. Расширение на более сложные варианты SLV моделей ## Глубокая оценка ### Преимущества 1. **Теоретическая строгость**: предоставление полных доказательств сходимости 2. **Высокая практическая ценность**: оба метода имеют преимущества, предоставляя гибкий выбор для практического применения 3. **Комплексные эксперименты**: систематическое сравнение производительности различных методов при различных условиях 4. **Ясное изложение**: точное математическое выражение, четкая логическая структура ### Недостатки 1. **Глубина теоретического анализа**: относительно меньше теоретического анализа метода обратного Эйлера 2. **Проверка на реальных рыночных данных**: отсутствие проверки на реальных рыночных данных 3. **Анализ вычислительной сложности**: отсутствие детального теоретического анализа сложности ### Влияние 1. **Академический вклад**: предоставление новых методов сильной сходимости для численного моделирования процесса CIR 2. **Практическая ценность**: прямое применение в оценке цены опционов и управлении рисками финансовых учреждений 3. **Воспроизводимость**: подробное описание методов, легко реализуемо ### Применимые сценарии 1. **Количественная торговля**: высокочастотная оценка цены опционов и стратегии хеджирования 2. **Управление рисками**: стресс-тестирование и оценка хвостовых рисков 3. **Калибровка модели**: оценка параметров SLV моделей и проверка ## Библиография Статья цитирует 25 связанных работ, охватывающих теорию SLV моделей, численные методы для процесса CIR, оценку цены опционов и другие ключевые области, предоставляя прочную теоретическую основу для исследования.