Cylindrical Projections of Occupied Diffusions

1. Probleemstelling

Het artikel behandelt de computationele onberekenbaarheid van het simuleren van Occupied Diffusions, een klasse van pad-afhankelijke stochastische processen. Deze processen worden gemodelleerd door Stochastische Differentiaalvergelijkingen (OSDE's) waarbij de coëfficiënten afhankelijk zijn van het occupatiemaat van de geschiedenis van het proces.

De dynamiek wordt gegeven door:
$$\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}$$
waarbij:

• $X_t \in \mathbb{R}^d$ het toestandsproces is.
• $O_t$ een eindig maat is op $\mathbb{R}^d$ die de tijd vertegenwoordigt die het proces in verschillende gebieden heeft doorgebracht (de occupatiestroom).
• De afhankelijkheid van de coëfficiënten ($b, \sigma, \lambda$) van $O_t$ een pad-afhankelijkheid van oneindige dimensie introduceert.

De Kernuitdaging: Hoewel het systeem Markoviaans is op het productruimte van maten en de Euclidische ruimte ($\mathcal{M}(\mathbb{R}^d) \times \mathbb{R}^d$), is de toestandsruimte oneindig-dimensionaal. Directe simulatie is onmogelijk omdat men het volledige maat $O_t$ niet kan opslaan of updaten bij elke tijdstap. Bestaande methoden worstelen om dit oneindig-dimensionale object te benaderen terwijl de fout onder controle wordt gehouden.

2. Methodologie: Cylindrische Projecties

De auteurs stellen een numerieke benaderingsmethode voor genaamd Cylindrische Projecties om het oneindig-dimensionale occupatiemaat te reduceren tot een eindig-dimensionaal systeem.

A. Cylindrische Structuur

De methode berust op het identificeren van een eindig maat $o$ met een rij van zijn integralen tegenover een aftelbare familie van testfuncties $(f_k)_{k \in \mathbb{N}}$.

• Cylindrische Norm: Een norm wordt gedefinieerd op de ruimte van maten met behulp van een scheidende familie van functies $(g_j)$:
  $$\|o\| := \left( \sum_{j=0}^\infty |o(g_j)|^2 \right)^{1/2}$$
• Projectie: Voor een truncatieniveau $K$ wordt het occupatiemaat $O_t$ geprojecteerd op een eindig-dimensionale vector $Z_t^K \in \mathbb{R}^{N_K+1}$:
  $$Z_t^K := (O_t(f_0), \dots, O_t(f_{N_K}))$$
  waarbij $(f_k^K)$ een eenheidspartitie is ondergeschikt aan een overdekking van de toestandsruimte.

B. Het Benaderde Systeem

De oorspronkelijke OSDE wordt benaderd door een eindig-dimensionaal systeem $(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}$$
Hier worden de coëfficiënten via een reconstructie-operator $o^K(z) = \sum z_k \delta_{x_k}$ terug "gelift" naar de maatr uimte. Dit transformeert het probleem in een standaard eindig-dimensionale SDE die kan worden gesimuleerd met standaard schema's zoals Euler–Maruyama.

3. Belangrijkste Bijdragen

1. Theoretisch Kader: Introductie van cylindrische projecties als een rigoureus hulpmiddel om pad-afhankelijke dynamiek gedreven door occupatiematen te benaderen.
2. Sterke Convergentiesnelheden: De auteurs bewijzen dat de benadering sterk convergeert naar het ware proces. Specifiek, voor een vaste tijdhorizon $T$:
   $$\mathbb{E}\left[ \sup_{0 \le t \le T} |X_t - X_t^K|^2 \right]^{1/2} \le \frac{C}{K}$$
   Dit vestigt een lineaire convergentiesnelheid met betrekking tot het truncatieniveau $K$.
3. Zwakke Foutanalyse: Het artikel breidt de analyse uit naar zwakke convergentie (verwachtingen van functionalen), waarbij wordt aangetoond dat voor Lipschitz-functionalen $\phi$:
   $$|\mathbb{E}[\phi(O_T, X_T)] - \mathbb{E}[\phi(o^K(Z_T^K), X_T^K)]| \le \frac{C}{K}$$
   Dit is cruciaal voor het prijzen van financiële derivaten.
4. Omgaan met Ongelimiteerde Coëfficiënten: De auteurs behandelen gevallen waarin de coëfficiënten lineair kunnen groeien (Lipschitz-voorwaarde) door lokalisatietechnieken (stopstoptijden) te gebruiken en convergentiesnelheden af te leiden die afhankelijk zijn van logaritmische factoren als de groei ongelimiteerd is.

4. Numerieke Resultaten en Voorbeelden

De methode wordt gevalideerd via Euler–Maruyama-simulaties op drie verschillende voorbeelden:

• Cranston–Le Jan Zelf-aantrekkende Diffusie:
  • Een proces dat wordt aangetrokken door zijn eigen lopend gemiddelde.
  • Er bestaat een expliciete sterke oplossing, wat directe foutvergelijking mogelijk maakt.
  • Resultaat: De numerieke fout neemt lineair af met $K$, wat de theoretische $O(1/K)$-snelheid bevestigt.
• Raimond's Zelf-interagerende Diffusie:
  • Een meervoudig dimensionale diffusie waarbij de drift afhankelijk is van de richting van de posities in het verleden.
  • Omdat er geen expliciete oplossing bestaat, vergelijkt de methode de projectie met een referentie met hoge resolutie ($K=75$).
  • Resultaat: Convergentie wordt waargenomen, beperkt alleen door bemonsteringsruis en discretisatiefout.
• Lokale Occupatie Volatiliteit (LOV) Model:
  • Een financieel model waarbij de volatiliteit afhankelijk is van de exponentiële tijd-occupatiestroom.
  • Resultaat: De methode simuleert succesvol pad-afhankelijke volatiliteit, wat de toepasbaarheid op complexe financiële derivaten aantoont.

5. Betekenis en Toepassingen

• Computationele Haalbaarheid: Het artikel biedt de eerste praktische, computationeel hanteerbare methode voor het simuleren van algemene occupied diffusions, die eerder als onberekenbaar werden beschouwd vanwege de oneindige dimensionaliteit.
• Prijzen van Financiële Derivaten: De zwakke foutanalyse ondersteunt direct het gebruik van Monte Carlo-methoden voor het prijzen van pad-afhankelijke opties (bijvoorbeeld Aziatische opties) in modellen met pad-afhankelijke volatiliteit (LOV).
• Generaliseerbaarheid: De aanpak van cylindrische projecties is veelzijdig en kan worden toegepast op andere proceswaarden-processen, waaronder McKean-Vlasov-type dynamiek en systemen met vertraagde interacties.
• Efficiëntie: In tegenstelling tot algemene pad-afhankelijke SDE's waarbij de toestandsdimensie groeit met de tijdstappen, behoudt de cylindrische projectie een vaste eindige dimensie $K$, wat aanzienlijke opslag- en computationele voordelen biedt.

Conclusie

Tissot-Daguette en Zhang slaan succesvol de brug tussen het theoretische bestaan van occupied diffusions en hun praktische simulatie. Door de introductie van cylindrische projecties reduceren ze een oneindig-dimensionaal probleem tot een eindig-dimensionaal probleem met rigoureuze foutgrenzen ($O(1/K)$), waardoor de simulatie van complexe zelf-interagerende systemen en het prijzen van geavanceerde financiële derivaten mogelijk wordt.