Simulation efficace des prix des options d'achat européennes sous le modèle de volatilité stochastique-locale de Heston : une comparaison de méthodes

Informations fondamentales

• ID de l'article : 2509.24449
• Titre : Efficient simulation of prices for European call options under Heston stochastic-local volatility model: a comparison of methods
• Auteurs : Meng Cai, Tianze Li
• Affiliations des auteurs : École de statistique et mathématiques, École d'économie, Université centrale des finances et de l'économie
• Classification : q-fin.CP (Finance computationnelle)
• Date de publication : 16 octobre 2025
• Lien de l'article : https://arxiv.org/abs/2509.24449

Résumé

Cet article étudie le modèle de volatilité stochastique-locale de Heston composé d'un processus de prix d'actif et d'un processus de variance de type Cox-Ingersoll-Ross, qui a des applications largement répandues dans le secteur financier. Contrairement à l'approximation du chi-carré non-central à convergence faible utilisée dans la littérature 25, cet article considère deux méthodes à convergence forte et préservant la positivité sous la transformation de Lamperti pour le processus CIR : la méthode d'Euler tronquée et la méthode d'Euler rétrograde. La méthode d'Euler tronquée explicite est efficace en calcul et robuste à forte volatilité, tandis que la méthode d'Euler rétrograde implicite offre une haute précision de calcul et une stabilité de performance. Les expériences numériques sur les options d'achat européennes démontrent les avantages des différentes méthodes.

Contexte de recherche et motivation

Contexte du problème

1. Besoins du modèle : L'ingénierie financière nécessite d'atteindre un équilibre entre un étalonnage parfait et une dynamique correcte ; le modèle de volatilité stochastique-locale (SLV) fournit une solution à cet égard
2. Limitations des modèles existants :
   • Les modèles de volatilité stochastique pure ont des paramètres limités et ne peuvent pas s'adapter parfaitement à la volatilité implicite pour tous les termes et prix d'exercice
   • Bien que le modèle de volatilité locale puisse étalonner toutes les options du marché actuel, il présente des problèmes dans la prédiction du comportement futur de la volatilité, ce qui entraîne un aplatissement du sourire de volatilité à terme

Défis fondamentaux

La simulation numérique du processus CIR est un point critique difficile de l'évaluation du modèle, car :

• Le terme de diffusion est non borné et les coefficients ne sont pas globalement Lipschitz continus
• La méthode d'Euler-Maruyama traditionnelle peut produire des variances négatives, ce qui n'a aucun sens mathématiquement et financièrement
• L'approximation du chi-carré non-central existante ne peut réaliser que la convergence faible

Motivation de la recherche

Développer des méthodes à convergence forte et préservant la positivité, garantissant que les trajectoires des solutions exactes et approximées sont proches les unes des autres, améliorant ainsi la tarification globale des options.

Contributions fondamentales

1. Proposition de deux nouvelles méthodes numériques : la méthode d'Euler tronquée et la méthode d'Euler rétrograde, toutes deux possédant la convergence forte et la préservation de la positivité
2. Analyse théorique : Preuve de la convergence forte de la méthode d'Euler tronquée (ordre de convergence 1/2)
3. Comparaison numérique complète : Comparaison systématique des performances de différentes méthodes par le biais d'expériences de tarification des options d'achat européennes
4. Orientation pratique : Recommandations de sélection de méthodes pour différents scénarios d'application

Détails des méthodes

Définition de la tâche

Résoudre le problème de tarification des options d'achat européennes sous le modèle de volatilité stochastique-locale de Heston, où l'élément clé est la simulation efficace et précise du processus de variance CIR.

Architecture du modèle

Définition du modèle HSLV

$\begin{cases} 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}$

Où :

• $S_t$ : Processus de prix d'actif
• $V_t$ : Processus de variance (dynamique CIR)
• $\sigma(t,s)$ : Fonction de volatilité locale
• $\kappa, \theta, \gamma$ : Paramètres du processus CIR

Transformation de Lamperti

Pour surmonter les difficultés numériques du processus CIR, on introduit la transformation $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$

Méthodes numériques fondamentales

1. Méthode d'Euler tronquée

Définition de l'application de troncature : $\pi_\tau(x) = (b\tau^{1/4}) \vee x$

Schéma numérique : $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$

Caractéristiques :

• Schéma explicite, calcul efficace
• Garantit la positivité par troncature
• Robustesse maintenue à forte volatilité

2. Méthode d'Euler rétrograde

$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$

Caractéristiques :

• Schéma implicite, nécessite une résolution itérative
• Offre une haute précision de calcul
• Performance stable mais coût de calcul élevé

Points d'innovation technique

1. Convergence forte : Contrairement à la méthode du chi-carré non-central à convergence faible, garantit la convergence au niveau des trajectoires
2. Conception préservant la positivité : Assure que la variance reste toujours positive par transformation mathématique et techniques de troncature
3. Complémentarité des méthodes : Les deux méthodes ont chacune des avantages en termes d'efficacité et de précision, adaptées à différents scénarios

Configuration expérimentale

Paramètres du modèle

• Paramètres de Heston de référence : $\gamma = 0.95$, $\kappa = 1.05$, $\rho = -0.315$, $\theta = 0.0855$, $S_0 = 1$, $r = 0$, $v_0 = 0.0945$, $T = 5$
• Paramètres de correction HSLV : Introduction d'un paramètre de correction $p = 0.25$, ajustant les différents paramètres

Indicateurs d'évaluation

• Pourcentage d'erreur relative
• Temps de calcul
• Robustesse dans différentes conditions de marché

Méthodes de comparaison

• Méthode d'Euler standard
• Méthode AES (Andersen Exact Simulation)
• Méthode d'Euler tronquée
• Méthode d'Euler rétrograde

Conception expérimentale

1. Calcul de l'espérance conditionnelle : Utilisation de la méthode 2D-COS comme référence, 20 bins, 10 000 trajectoires, pas de temps 0.001
2. Analyse de sensibilité des paramètres : Test des performances pour différentes valeurs de $\theta$ et $p$
3. Comparaison des prix d'exercice : Trois cas : $K = 70\%, 100\%, 150\%$

Résultats expérimentaux

Résultats principaux

Comparaison des erreurs pour différents nombres d'étapes

Comparaison des temps de calcul (40 étapes)

• Méthode AES : 0.1114 secondes
• Méthode d'Euler : 0.0939 secondes
• Méthode d'Euler tronquée : 0.0733 secondes
• Méthode d'Euler rétrograde : 137.9476 secondes

Découvertes clés

1. Avec peu d'étapes : La méthode d'Euler tronquée présente généralement l'erreur la plus faible
2. Avec beaucoup d'étapes : La méthode d'Euler rétrograde offre la meilleure précision
3. Robustesse des paramètres : La méthode d'Euler tronquée montre une croissance lente des erreurs lorsque $\theta$ augmente, démontrant une excellente robustesse
4. Efficacité de calcul : La méthode d'Euler tronquée est la plus rapide ; la méthode d'Euler rétrograde est significativement plus coûteuse en raison de la résolution itérative requise

Expériences d'ablation

Impact du paramètre de volatilité $\theta$

Les expériences montrent qu'avec l'augmentation de $\theta$ (simulation de conditions de marché extrêmes) :

• Les erreurs des méthodes AES et Euler standard augmentent significativement
• La méthode d'Euler tronquée reste robuste avec un taux de croissance d'erreur plus faible
• Cela indique que la méthode de troncature est plus appropriée pour les tests de stress et la modélisation du risque de queue

Impact du paramètre de correction $p$

Avec l'augmentation de $p$, la méthode d'Euler rétrograde offre toujours les meilleures performances, mais elle est plus sensible au choix du pas de temps.

Travaux connexes

Principaux axes de recherche

1. Étalonnage du modèle SLV : Travaux fondateurs de Guyon et Henry-Labordère (2011)
2. Méthodes numériques pour le processus CIR : Incluant les différences finies, les méthodes de particules Monte Carlo, etc.
3. Applications d'apprentissage profond : Cuchiero et al. (2020) utilisant les réseaux antagonistes génératifs pour l'étalonnage

Avantages de cet article

Par rapport aux travaux existants, cet article fournit des méthodes à convergence forte préservant la positivité, comblant les lacunes des méthodes à convergence faible et offrant des outils numériques plus fiables pour les applications pratiques.

Conclusions et discussion

Conclusions principales

1. Complémentarité des méthodes : La méthode d'Euler tronquée convient aux simulations à grande échelle et à la gestion des risques en temps réel ; la méthode d'Euler rétrograde convient aux tests de stress et aux scénarios nécessitant une haute précision
2. Avantages de robustesse : La méthode d'Euler tronquée offre des performances excellentes dans les conditions de marché extrêmes
3. Compromis précision-efficacité : Le choix du pas de temps et de la méthode appropriée dépend des besoins spécifiques de l'application

Limitations

1. Le coût de calcul de la méthode d'Euler rétrograde est relativement élevé
2. Le choix du pas de temps nécessite un compromis entre précision et efficacité
3. Pour certaines combinaisons de paramètres extrêmes, la stabilité des méthodes nécessite une vérification supplémentaire

Directions futures

1. Développer des stratégies de pas de temps adaptatif
2. Explorer les implémentations parallélisées
3. Étendre à des variantes plus complexes du modèle SLV

Évaluation approfondie

Points forts

1. Rigueur théorique : Fournit des preuves complètes de convergence
2. Valeur pratique élevée : Les deux méthodes ont chacune des avantages, offrant des choix flexibles pour les applications pratiques
3. Expériences complètes : Comparaison systématique des performances de différentes méthodes dans diverses conditions
4. Clarté de la rédaction : Expression mathématique précise et structure logique claire

Insuffisances

1. Profondeur de l'analyse théorique : L'analyse théorique de la méthode d'Euler rétrograde est relativement limitée
2. Validation sur marché réel : Manque de vérification avec des données de marché réelles
3. Analyse de complexité de calcul : Absence d'analyse théorique détaillée de la complexité

Impact

1. Contribution académique : Fournit de nouvelles méthodes à convergence forte pour la simulation numérique du processus CIR
2. Valeur pratique : Directement applicable à la tarification des options et à la gestion des risques dans les institutions financières
3. Reproductibilité : Description détaillée des méthodes, facile à mettre en œuvre

Scénarios d'application

1. Trading quantitatif : Tarification d'options haute fréquence et stratégies de couverture
2. Gestion des risques : Tests de stress et évaluation du risque de queue
3. Étalonnage du modèle : Estimation et validation des paramètres du modèle SLV

Références

L'article cite 27 références pertinentes, couvrant la théorie des modèles SLV, les méthodes numériques pour le processus CIR, la tarification des options et d'autres domaines clés, fournissant une base théorique solide pour la recherche.