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 articolo: 2509.24449
- Titolo: Efficient simulation of prices for European call options under Heston stochastic-local volatility model: a comparison of methods
- Autori: Meng Cai, Tianze Li
- Affiliazione autori: Scuola di Statistica e Matematica, Scuola di Economia, Università Centrale di Finanza ed Economia
- Classificazione: q-fin.CP (Finanza computazionale)
- Data di pubblicazione: 16 ottobre 2025
- Link articolo: https://arxiv.org/abs/2509.24449
Questo articolo esamina il modello Heston di volatilità stocastica-locale composto da un processo di prezzo dell'asset e da un processo di varianza di tipo Cox-Ingersoll-Ross, che ha ampia applicazione nel settore finanziario. A differenza dell'approssimazione chi-quadrato non centrale a convergenza debole utilizzata nella letteratura 25, questo articolo considera due metodi che presentano convergenza forte e preservano la positività sotto la trasformazione di Lamperti del processo CIR: il metodo di Euler troncato e il metodo di Euler all'indietro. Il metodo di Euler troncato esplicito è computazionalmente efficiente e mantiene robustezza ad alta volatilità, mentre il metodo di Euler all'indietro implicito fornisce elevata precisione computazionale e stabilità delle prestazioni. Gli esperimenti numerici su opzioni call europee dimostrano i vantaggi dei diversi metodi.
- Esigenze del modello: L'ingegneria finanziaria richiede di raggiungere un equilibrio tra calibrazione perfetta e dinamiche corrette; il modello di volatilità stocastica-locale (SLV) fornisce una soluzione a questo proposito
- Limitazioni dei modelli esistenti:
- I modelli di volatilità stocastica pura hanno parametri limitati e difficilmente si adattano perfettamente alla volatilità implicita di tutti i termini e prezzi di esercizio
- Sebbene i modelli di volatilità locale possano calibrare tutte le opzioni del mercato attuale, presentano problemi nella previsione del comportamento futuro della volatilità, causando un appiattimento del sorriso della volatilità forward
La simulazione numerica del processo CIR è un punto critico nella valutazione del modello, poiché:
- Il termine di diffusione è illimitato e i coefficienti non sono globalmente Lipschitz continui
- Il metodo tradizionale di Euler-Maruyama può produrre varianze negative, prive di significato sia matematico che finanziario
- Le approssimazioni chi-quadrato non centrale esistenti possono solo raggiungere convergenza debole
Sviluppare metodi con convergenza forte e che preservino la positività, garantendo che le traiettorie della soluzione esatta e approssimata siano vicine tra loro, migliorando così la determinazione complessiva dei prezzi delle opzioni.
- Proposizione di due nuovi metodi numerici: il metodo di Euler troncato e il metodo di Euler all'indietro, entrambi con convergenza forte e preservazione della positività
- Analisi teorica: dimostrazione della convergenza forte del metodo di Euler troncato (ordine di convergenza 1/2)
- Confronto numerico completo: attraverso esperimenti di determinazione dei prezzi delle opzioni call europee, confronto sistematico delle prestazioni di diversi metodi
- Guida pratica: raccomandazioni sulla scelta dei metodi per diversi scenari applicativi
Risolvere il problema della determinazione dei prezzi delle opzioni call europee secondo il modello Heston di volatilità stocastica-locale, dove il nucleo è la simulazione efficiente e accurata del processo di varianza 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}$$
Dove:
- $S_t$: processo di prezzo dell'asset
- $V_t$: processo di varianza (dinamica CIR)
- $\sigma(t,s)$: funzione di volatilità locale
- $\kappa, \theta, \gamma$: parametri del processo CIR
#### Trasformazione di Lamperti
Per superare le difficoltà numeriche del processo CIR, si introduce la trasformazione $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$$
### Metodi numerici fondamentali
#### 1. Metodo di Euler troncato
Definire la mappatura di troncamento: $\pi_\tau(x) = (b\tau^{1/4}) \vee x$
Schema numerico:
$$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$$
**Caratteristiche**:
- Schema esplicito, computazionalmente efficiente
- Garantisce positività attraverso il troncamento
- Mantiene robustezza ad alta volatilità
#### 2. Metodo di Euler all'indietro
$$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$$
**Caratteristiche**:
- Schema implicito, richiede risoluzione iterativa
- Fornisce elevata precisione computazionale
- Prestazioni stabili ma costi computazionali elevati
### Punti di innovazione tecnica
1. **Convergenza forte**: A differenza del metodo chi-quadrato non centrale a convergenza debole, garantisce convergenza a livello di traiettoria
2. **Progettazione della preservazione della positività**: Garantisce che la varianza rimanga sempre positiva attraverso trasformazioni matematiche e tecniche di troncamento
3. **Complementarità dei metodi**: I due metodi hanno vantaggi reciproci in efficienza e precisione, adatti a diversi scenari
## Configurazione sperimentale
### Parametri del modello
- **Parametri Heston di riferimento**: $\gamma = 0.95$, $\kappa = 1.05$, $\rho = -0.315$, $\theta = 0.0855$, $S_0 = 1$, $r = 0$, $v_0 = 0.0945$, $T = 5$
- **Parametri di correzione HSLV**: Introduzione del parametro di correzione $p = 0.25$, regolazione di vari parametri
### Indicatori di valutazione
- Percentuale di errore relativo
- Tempo di calcolo
- Robustezza in diverse condizioni di mercato
### Metodi di confronto
- Metodo di Euler standard
- Metodo AES (Andersen Exact Simulation)
- Metodo di Euler troncato
- Metodo di Euler all'indietro
### Progettazione degli esperimenti
1. **Calcolo dell'aspettativa condizionata**: Utilizzo del metodo 2D-COS come riferimento, 20 bin, 10.000 traiettorie, ampiezza del passo 0.001
2. **Analisi di sensibilità dei parametri**: Test delle prestazioni con diversi valori di $\theta$ e $p$
3. **Confronto dei prezzi di esercizio**: Tre casi $K = 70\%, 100\%, 150\%$
## Risultati sperimentali
### Risultati principali
#### Confronto degli errori con diversi numeri di passi
| Numero passi N | K=70% | K=100% | K=150% |
|-------|--------|---------|---------|
| **Metodo di Euler troncato** |
| 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%) |
| **Metodo di Euler all'indietro** |
| 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%) |
#### Confronto dei tempi di calcolo (40 passi)
- Metodo AES: 0.1114 secondi
- Metodo di Euler: 0.0939 secondi
- Metodo di Euler troncato: 0.0733 secondi
- Metodo di Euler all'indietro: 137.9476 secondi
### Scoperte chiave
1. **Con pochi passi**: Il metodo di Euler troncato di solito presenta l'errore più basso
2. **Con molti passi**: Il metodo di Euler all'indietro ha la massima precisione
3. **Robustezza dei parametri**: Il metodo di Euler troncato mostra una crescita lenta dell'errore all'aumentare di $\theta$, dimostrando eccellente robustezza
4. **Efficienza computazionale**: Il metodo di Euler troncato è il più veloce; il metodo di Euler all'indietro richiede tempi significativamente più lunghi a causa della risoluzione iterativa
### Esperimenti di ablazione
#### Effetto del parametro di volatilità $\theta$
Gli esperimenti mostrano che all'aumentare di $\theta$ (simulando condizioni di mercato estreme):
- Gli errori dei metodi AES e Euler standard aumentano significativamente
- Il metodo di Euler troncato mantiene robustezza con un tasso di crescita dell'errore inferiore
- Ciò indica che il metodo di troncamento è più adatto per test di stress e modellazione del rischio di coda
#### Effetto del parametro di correzione $p$
All'aumentare di $p$, il metodo di Euler all'indietro mostra sempre le migliori prestazioni, ma presenta maggiore sensibilità all'ampiezza del passo.
## Lavori correlati
### Principali direzioni di ricerca
1. **Calibrazione del modello SLV**: Lavoro pioneristico di Guyon e Henry-Labordère (2011)
2. **Metodi numerici per il processo CIR**: Inclusi metodi alle differenze finite, metodi di particelle Monte Carlo, ecc.
3. **Applicazioni di apprendimento profondo**: Cuchiero et al. (2020) utilizzano reti generative avversarie per la calibrazione
### Vantaggi di questo articolo
Rispetto ai lavori esistenti, questo articolo fornisce metodi con convergenza forte che preservano la positività, colmando le lacune dei metodi a convergenza debole e fornendo strumenti numerici più affidabili per applicazioni pratiche.
## Conclusioni e discussione
### Conclusioni principali
1. **Complementarità dei metodi**: Il metodo di Euler troncato è adatto per simulazioni su larga scala e gestione del rischio in tempo reale; il metodo di Euler all'indietro è adatto per test di stress e scenari che richiedono elevata precisione
2. **Vantaggi di robustezza**: Il metodo di Euler troncato mostra eccellenti prestazioni in condizioni di mercato estreme
3. **Compromesso precisione-efficienza**: La scelta dell'ampiezza del passo e del metodo appropriato deve essere basata sulle esigenze specifiche dell'applicazione
### Limitazioni
1. Il costo computazionale del metodo di Euler all'indietro è relativamente elevato
2. La scelta dell'ampiezza del passo richiede un compromesso tra precisione ed efficienza
3. Per alcune combinazioni di parametri estreme, la stabilità dei metodi richiede ulteriore verifica
### Direzioni future
1. Sviluppo di strategie di ampiezza del passo adattive
2. Esplorazione di implementazioni parallele
3. Estensione a varianti più complesse del modello SLV
## Valutazione approfondita
### Punti di forza
1. **Rigore teorico**: Fornisce prove complete di convergenza
2. **Alto valore pratico**: I due metodi hanno vantaggi reciproci, offrendo scelte flessibili per applicazioni pratiche
3. **Esperimenti completi**: Confronto sistematico delle prestazioni di diversi metodi in varie condizioni
4. **Scrittura chiara**: Espressione matematica accurata, struttura logica chiara
### Carenze
1. **Profondità dell'analisi teorica**: L'analisi teorica del metodo di Euler all'indietro è relativamente limitata
2. **Validazione del mercato reale**: Mancanza di verifica con dati di mercato reali
3. **Analisi della complessità computazionale**: Mancanza di analisi teorica dettagliata della complessità
### Impatto
1. **Contributo accademico**: Fornisce nuovi metodi a convergenza forte per la simulazione numerica del processo CIR
2. **Valore pratico**: Direttamente applicabile alla determinazione dei prezzi delle opzioni e alla gestione del rischio nelle istituzioni finanziarie
3. **Riproducibilità**: La descrizione dei metodi è dettagliata e facile da implementare
### Scenari applicabili
1. **Trading quantitativo**: Determinazione dei prezzi delle opzioni ad alta frequenza e strategie di copertura
2. **Gestione del rischio**: Test di stress e valutazione del rischio di coda
3. **Calibrazione del modello**: Stima dei parametri del modello SLV e verifica
## Bibliografia
L'articolo cita 25 riferimenti correlati, che coprono i campi chiave della teoria del modello SLV, metodi numerici per il processo CIR, determinazione dei prezzi delle opzioni e altri lavori importanti, fornendo una base teorica solida per la ricerca.