Cylindrical Projections of Occupied Diffusions

1. Énoncé du problème

Le papier aborde l'ingérabilité computationnelle de la simulation des diffusions occupées, une classe de processus stochastiques dépendants du chemin. Ces processus sont modélisés par des équations différentielles stochastiques (EDS) où les coefficients dépendent de la mesure d'occupation de l'historique du processus.

Les dynamiques sont données par :
$$\begin{cases} d O_t = \delta_{X_t} \lambda(O_t, X_t) dt \\ d X_t = b(O_t, X_t) dt + \sigma(O_t, X_t) dW_t \end{cases}$$
où :

• $X_t \in \mathbb{R}^d$ est le processus d'état.
• $O_t$ est une mesure finie sur $\mathbb{R}^d$ représentant le temps passé par le processus dans diverses régions (le flot d'occupation).
• La dépendance des coefficients ($b, \sigma, \lambda$) par rapport à $O_t$ introduit une dépendance au chemin de dimension infinie.

Le défi central : Bien que le système soit markovien sur l'espace produit des mesures et de l'espace euclidien ($\mathcal{M}(\mathbb{R}^d) \times \mathbb{R}^d$), l'espace d'état est de dimension infinie. La simulation directe est impossible car on ne peut ni stocker ni mettre à jour la mesure complète $O_t$ à chaque pas de temps. Les méthodes existantes peinent à approximer cet objet de dimension infinie tout en contrôlant l'erreur.

2. Méthodologie : Projections cylindriques

Les auteurs proposent une méthode d'approximation numérique appelée projections cylindriques pour réduire la mesure d'occupation de dimension infinie à un système de dimension finie.

A. Structure cylindrique

La méthode repose sur l'identification d'une mesure finie $o$ avec une séquence de ses intégrales par rapport à une famille dénombrable de fonctions tests $(f_k)_{k \in \mathbb{N}}$.

• Norme cylindrique : Une norme est définie sur l'espace des mesures en utilisant une famille séparante de fonctions $(g_j)$ :
  $$\|o\| := \left( \sum_{j=0}^\infty |o(g_j)|^2 \right)^{1/2}$$
• Projection : Pour un niveau de troncature $K$, la mesure d'occupation $O_t$ est projetée sur un vecteur de dimension finie $Z_t^K \in \mathbb{R}^{N_K+1}$ :
  $$Z_t^K := (O_t(f_0), \dots, O_t(f_{N_K}))$$
  où $(f_k^K)$ est une partition de l'unité subordonnée à un recouvrement de l'espace d'état.

B. Le système approché

L'EDS originale est approchée par un système de dimension finie $(X_t^K, Z_t^K)$ :
$$\begin{cases} d Z_t^K = f^K(X_t^K) \lambda^K(Z_t^K, X_t^K) dt \\ d X_t^K = b^K(Z_t^K, X_t^K) dt + \sigma^K(Z_t^K, X_t^K) dW_t \end{cases}$$
Ici, les coefficients sont « relevés » vers l'espace des mesures via un opérateur de reconstruction $o^K(z) = \sum z_k \delta_{x_k}$. Cela transforme le problème en une EDS standard de dimension finie qui peut être simulée à l'aide de schémas classiques comme Euler-Maruyama.

3. Contributions clés

1. Cadre théorique : Introduction des projections cylindriques comme outil rigoureux pour approximer les dynamiques dépendantes du chemin pilotées par des mesures d'occupation.
2. Vitesses de convergence forte : Les auteurs prouvent que l'approximation converge fortement vers le processus réel. Plus précisément, pour un horizon temporel fixe $T$ :
   $$\mathbb{E}\left[ \sup_{0 \le t \le T} |X_t - X_t^K|^2 \right]^{1/2} \le \frac{C}{K}$$
   Cela établit une vitesse de convergence linéaire par rapport au niveau de troncature $K$.
3. Analyse de l'erreur faible : Le papier étend l'analyse à la convergence faible (espérances de fonctionnelles), montrant que pour des fonctionnelles lipschitziennes $\phi$ :
   $$|\mathbb{E}[\phi(O_T, X_T)] - \mathbb{E}[\phi(o^K(Z_T^K), X_T^K)]| \le \frac{C}{K}$$
   Ceci est crucial pour la tarification des produits dérivés financiers.
4. Gestion des coefficients non bornés : Les auteurs traitent les cas où les coefficients peuvent croître linéairement (condition de Lipschitz) en utilisant des techniques de localisation (temps d'arrêt) et en dérivant des vitesses de convergence qui dépendent de facteurs logarithmiques si la croissance est non bornée.

4. Résultats numériques et exemples

La méthode est validée par des simulations Euler-Maruyama sur trois exemples distincts :

• Diffusion auto-attractive de Cranston–Le Jan :
  • Un processus attiré par sa propre moyenne glissante.
  • Une solution forte explicite existe, permettant une comparaison directe des erreurs.
  • Résultat : L'erreur numérique diminue linéairement avec $K$, confirmant le taux théorique $O(1/K)$.
• Diffusion auto-interagissante de Raimond :
  • Une diffusion multidimensionnelle où la dérive dépend de la direction des positions passées.
  • Comme aucune solution explicite n'existe, la méthode compare la projection contre une référence haute résolution ($K=75$).
  • Résultat : Une convergence est observée, limitée uniquement par le bruit d'échantillonnage et l'erreur de discrétisation.
• Modèle de volatilité occupée locale (LOV) :
  • Un modèle financier où la volatilité dépend du flot d'occupation exponentiel dans le temps.
  • Résultat : La méthode simule avec succès une volatilité dépendante du chemin, démontrant son applicabilité à des produits dérivés financiers complexes.

5. Importance et applications

• Faisabilité computationnelle : Le papier fournit la première méthode pratique et computationnellement traitable pour simuler des diffusions occupées générales, qui étaient auparavant considérées comme ingérables en raison de leur dimension infinie.
• Tarification des produits dérivés financiers : L'analyse de l'erreur faible soutient directement l'utilisation de méthodes de Monte Carlo pour la tarification d'options dépendantes du chemin (par exemple, les options asiatiques) dans des modèles à volatilité dépendante du chemin (LOV).
• Généralisabilité : L'approche par projection cylindrique est polyvalente et peut être appliquée à d'autres processus à valeurs mesures, y compris les dynamiques de type McKean-Vlasov et les systèmes avec interactions retardées.
• Efficacité : Contrairement aux EDS dépendantes du chemin générales où la dimension de l'état croît avec les pas de temps, la projection cylindrique maintient une dimension finie fixe $K$, offrant des avantages significatifs en termes de stockage et de calcul.

Conclusion

Tissot-Daguette et Zhang comblent avec succès le fossé entre l'existence théorique des diffusions occupées et leur simulation pratique. En introduisant les projections cylindriques, ils réduisent un problème de dimension infinie à un problème de dimension finie avec des bornes d'erreur rigoureuses ($O(1/K)$), permettant la simulation de systèmes auto-interagissants complexes et la tarification de produits dérivés financiers sophistiqués.