Résumé technique : Équations de noyaux Signature et Schwinger-Dyson comme équations différentielles rugueuses à deux paramètres

1. Énoncé du problème

L'article aborde la formulation mathématique et l'approximation numérique de deux équations de noyaux spécifiques — le noyau Signature et le noyau de Schwinger-Dyson — lorsqu'ils sont pilotés par des signaux hautement irréguliers qui sortent du régime classique à variation bornée.

Bien que ces noyaux aient été précédemment bien compris pour des trajectoires à variation bornée (où l'intégration classique s'applique), de nombreux signaux d'intérêt dans les applications (par exemple, le mouvement brownien, les solutions d'équations différentielles stochastiques) sont trop irréguliers pour une intégration standard à deux variables. Les auteurs visent à :

1. Développer un cadre de trajectoires rugueuses rigoureux pour les équations différentielles à deux paramètres sur des domaines rectangulaires et simpliciaux.
2. Interpréter les équations de noyaux Signature et de Schwinger-Dyson comme des instances d'équations différentielles rugueuses à deux paramètres (2P-RDE).
3. Établir le bon positionnement (existence, unicité et stabilité) de ces équations dans des espaces de trajectoires rugueuses conjointement contrôlées.
4. Dériver des schémas numériques pour l'équation de Schwinger-Dyson rugueuse et analyser leur complexité computationnelle.

2. Méthodologie

2.1 Cadre théorique : Trajectoires rugueuses conjointement contrôlées

La théorie est construite sur l'espace des trajectoires rugueuses conjointement contrôlées, étendant les travaux de [GH19] et [CP23].

• Définitions : Une trajectoire à deux paramètres $Z$ est $(X, Y)$-contrôlée si ses incréments dans le premier paramètre sont contrôlés par la trajectoire rugueuse $X$ et dans le second par $Y$, avec des termes de reste spécifiques satisfaisant des conditions de régularité conjointe.
• Stabilité sous composition : Une étape technique critique consiste à prouver que la classe des trajectoires contrôlées à deux paramètres est stable sous composition avec des fonctions lisses. Cela permet aux termes non linéaires dans les équations de noyaux (par exemple, les produits du noyau avec lui-même) de rester dans l'espace contrôlé, permettant ainsi la définition de 2P-RDE non linéaires.
• Intégration rugueuse : Les auteurs utilisent un cadre robuste pour l'intégration rugueuse sur des rectangles $[a, b] \times [c, d]$, définie comme la limite de sommes discrètes sur des partitions rectangulaires.

2.2 Extension aux domaines simpliciaux

Une innovation méthodologique clé est l'extension de l'intégration rugueuse aux simplices bidimensionnels (domaines triangulaires), ce qui est requis pour l'équation de Schwinger-Dyson.

• Symétrisation : Les auteurs définissent une trajectoire « symétrisée » $\hat{Z}$ à travers la diagonale d'un domaine carré.
• Théorème de type Fubini : Ils prouvent que l'intégrale sur un simplexe peut être récupérée à partir de l'intégrale sur un rectangle via cette symétrisation. Plus précisément, l'intégrale sur le simplexe est définie comme la moitié de l'intégrale de la trajectoire symétrisée sur le carré. Cela comble le fossé entre la théorie d'intégration rectangulaire de [CP23] et les exigences simpliciales de l'équation de Schwinger-Dyson.

2.3 Formulation des équations de noyaux

L'article formule les noyaux comme des points fixes de 2P-RDE :

• Noyau Signature : Défini sur un domaine rectangulaire $D_{t,v} = [a, t] \times [c, v]$. L'équation est linéaire en le noyau $K$ :
  $$K(t, v) = 1 + \int_{[a,t]\times[c,v]} K(s, u) \langle dX_s, dY_u \rangle$$
  Il est démontré qu'il s'agit d'une 2P-RDE où l'intégrande est une fonction lisse de la trajectoire contrôlée.
• Noyau de Schwinger-Dyson : Défini sur un domaine simplicial $D_{t,v} = \Delta_{[t,v]}$. L'équation est quadratique :
  $$K(s, t) = 1 - \int_{\Delta_{[s,t]}} K(s, u) K(u, r) \langle dX_u, dX_r \rangle$$
  Les auteurs étendent cela des trajectoires à variation bornée aux trajectoires rugueuses géométriques de régularité arbitrairement faible ($p \ge 1$).

2.2 Régime lisse et reformulation intégrale-différentielle

Lorsque les signaux directeurs sont des trajectoires rugueuses lisses (admettant une dérivée temporelle le long de la diagonale), les auteurs relient les 2P-RDE aux EDP classiques et aux équations intégrales-différentielles :

• Pour le noyau Signature, la 2P-RDE retrouve l'EDP linéaire dérivée dans [LLS24].
• Pour le noyau de Schwinger-Dyson, les auteurs dérivent un système d'équations intégrales-différentielles. Cela est réalisé en exploitant des identités combinatoires impliquant des appariements partiels non croisés (liés à la probabilité libre) pour relier la composante de trace du noyau à ses composantes tensorielles d'ordre supérieur.

2.5 Schéma numérique

Un schéma numérique est proposé pour l'équation de Schwinger-Dyson rugueuse, basé sur le système intégral-différentiel dérivé et les identités combinatoires.

• Algorithme : La méthode utilise la programmation dynamique sur une grille, résolvant un système linéaire à chaque étape. Elle exploite la structure des partitions non croisées pour calculer efficacement les termes d'ordre supérieur.
• Analyse de complexité : Les auteurs fournissent des estimations rigoureuses pour la complexité temporelle et mémoire, montrant une dépendance exponentielle au niveau de troncature de l'algèbre tensorielle et une dépendance polynomiale à la dimension de la grille.

3. Contributions clés

1. Cadre unifié de trajectoires rugueuses : L'article établit une théorie générale pour les 2P-RDE sur les domaines rectangulaires et simpliciaux, unifiant le traitement des noyaux Signature et de Schwinger-Dyson sous l'égide des trajectoires rugueuses conjointement contrôlées.
2. Intégration rugueuse simpliciale : Il introduit une définition rigoureuse de l'intégration rugueuse sur des simplices bidimensionnels à faible régularité, étendant les résultats précédents limités aux approximations lisses ou à des régimes de variation spécifiques.
3. Bon positionnement et stabilité : Les auteurs prouvent l'existence, l'unicité et la dépendance continue par rapport aux trajectoires rugueuses directrices pour les deux équations de noyaux. Cela fournit une fondation théorique solide pour l'approximation de ces noyaux avec des données irrégulières.
4. Reformulation combinatoire : Pour le noyau de Schwinger-Dyson, l'article dérive une nouvelle reformulation intégrale-différentielle dans le régime lisse en utilisant des identités combinatoires, offrant une description analytique alternative.
5. Implémentation et analyse numériques : Un algorithme numérique concret est présenté avec une analyse détaillée de la complexité. L'article valide les estimations de complexité théoriques et démontre les performances de la méthode à travers différents régimes de régularité.

4. Résultats

• Théoriques : Le noyau Signature est prouvé être l'unique solution d'une 2P-RDE linéaire. Le noyau de Schwinger-Dyson est prouvé être l'unique solution d'une 2P-RDE quadratique sur un simplexe. Les deux solutions sont montrées comme étant stables par rapport aux perturbations des trajectoires rugueuses directrices.
• Régime lisse : Dans le cas de trajectoires rugueuses lisses, la 2P-RDE pour le noyau Signature se réduit à l'EDP connue. Le noyau de Schwinger-Dyson est montré satisfaisant un système spécifique d'équations intégrales-différentielles dérivé des propriétés combinatoires des composantes tensorielles du noyau.
• Expériences numériques :
  • Les auteurs ont testé le schéma proposé (désigné $SDK_\kappa$) contre une méthode CDE randomisée ($SDK_r$) de [CT26] en utilisant un mouvement brownien fractionnaire avec des indices de Hurst variables ($H$).
  • Constats : À mesure que la régularité du signal directeur diminue (H plus faible), l'écart (Erreur Absolue Moyenne et Écart Type) entre les différents schémas numériques augmente.
  • Les résultats confirment que les schémas avec des ordres d'approximation plus élevés (κ plus élevé) maintiennent un meilleur accord entre eux lorsque la régularité diminue, ce qui est cohérent avec l'attente théorique selon laquelle une régularité plus faible nécessite des approximations d'ordre supérieur.
  • Les expériences de temps d'exécution ont confirmé l'échelle exponentielle avec la taille de la grille et la taille des blocs, et l'échelle polynomiale avec la dimension, correspondant aux estimations de complexité théorique.

5. Signification et revendications

L'article revendique de fournir un cadre robuste de trajectoires rugueuses permettant l'interprétation, l'analyse et l'approximation numérique des identités de noyaux (Signature et Schwinger-Dyson) pilotées par des signaux qui ne sont pas à variation bornée.

• Signification théorique : En encadrant ces noyaux comme des 2P-RDE, les auteurs exploitent les puissants résultats de stabilité et de bon positionnement de la théorie des trajectoires rugueuses. Cela est particulièrement pertinent pour les applications impliquant des processus stochastiques (comme le mouvement brownien) où l'intégration classique échoue.
• Signification pratique : Le développement d'un schéma numérique pour l'équation de Schwinger-Dyson rugueuse, accompagné d'une analyse de complexité, permet le calcul de ces noyaux dans des contextes précédemment inaccessibles aux méthodes standard. Les expériences suggèrent que les méthodes proposées sont cohérentes et que les exigences de régularité théoriques sont nécessaires pour une approximation précise.
• Limites : Les auteurs notent que certaines comparaisons numériques impliquent des signaux directeurs en dehors de la justification théorique de certains schémas (par exemple, l'utilisation de méthodes à variation bornée sur des trajectoires rugueuses), traitant ces comparaisons comme des références heuristiques plutôt que comme des preuves de convergence rigoureuses.

En résumé, ce travail comble le fossé entre la structure algébrique des trajectoires rugueuses et les propriétés analytiques des équations de noyaux, fournissant à la fois une fondation théorique pour les données irrégulières et un outil pratique pour leur calcul.