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.
- 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
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.
- 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
- 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
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
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.
- 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é
- Analyse théorique : Preuve de la convergence forte de la méthode d'Euler tronquée (ordre de convergence 1/2)
- 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
- Orientation pratique : Recommandations de sélection de méthodes pour différents scénarios d'application
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.
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
| Nombre d'étapes N | K=70% | K=100% | K=150% |
|-------|--------|---------|---------|
| **Méthode d'Euler tronquée** |
| 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%) |
| **Méthode d'Euler rétrograde** |
| 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%) |
#### 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.