Optimal Calibration of the Endpoint-corrected Hilbert Transform

Résumé Technique : Calibrage Optimal de la Transformée de Hilbert avec Correction de l'Extrémité (Endpoint-corrected)

Énoncé du Problème
L'estimation précise et à faible latence de la phase instantanée des oscillations est cruciale pour la détection et l'actionnement en boucle fermée, particulièrement dans les systèmes de neurostimulation à verrouillage de phase et les systèmes de contrôle en temps réel. L'approche standard consiste à construire un signal analytique via la Transformée de Hilbert (TH) sur une fenêtre glissante. Cependant, les implémentations en temps discret utilisant la Transformée de Fourier Discrète (TFD) supposent une périodicité au sein de la fenêtre. Lorsque le signal n'est pas parfaitement périodique à travers les limites de la fenêtre, cette hypothèse crée des discontinuités qui se manifestent par un phénomène de Gibbs (ringing), provoquant des distorsions significatives d'amplitude et de phase aux bords de la fenêtre. Puisque les systèmes en temps réel reposent sur l'échantillon le plus récent (l'extrémité ou endpoint), ces artefacts de bord corrompent directement l'estimation critique de la phase.

Bien que la Transformée de Hilbert avec correction de l'extrémité (ecHT) ait été introduite pour atténuer ces problèmes en appliquant un filtre à bande étroite causal au spectre analytique, ses distorsions systématiques à l'extrémité ont manqué d'une analyse théorique fondée sur des principes. Les évaluations précédentes étaient largement empiriques, laissant la nature du biais résiduel et la stratégie de correction optimale indéfinies.

Méthodologie
Cette étude fournit une caractérisation analytique de l'ecHT basée sur la TFD. Les auteurs dérivent analytiquement l'opérateur d'extrémité pour un signal monocromatique à bande étroite, décomposant la sortie de l'ecHT en deux composantes distinctes :

1. Un gain complexe déterministe ($G_+$) : Un terme de fréquence positive qui induit un biais d'amplitude et de phase calibrable.
2. Un terme de fuite résiduelle ($G_-$) : Un terme de fréquence négative qui échappe au masque analytique et au filtre, créant une ondulation dépendante de la phase qui varie avec la phase instantanée.

Sur la base de cette décomposition, l'article dérive un calibrage scalaire optimal en erreur quadratique moyenne (MSE). Les auteurs prouvent que pour toute implémentation de l'ecHT, il existe un scalaire complexe unique $C_{opt}$ qui minimise l'erreur entre la sortie calibrée et le signal analytique idéal. Ce scalaire est dérivé analytiquement comme suit :
$$C_{opt} = \frac{E[\hat{Z}^* Z]}{E[|\hat{Z}|^2]}$$
où $Z$ est le véritable point d'extrémité et $\hat{Z}$ est l'estimation par l'ecHT. Géométriquement, cela représente une projection orthogonale (un filtre de Wiener scalaire) qui extrait l'information linéairement récupérable de manière maximale. Les auteurs démontrent que $C_{opt}$ dépend uniquement du gain de fréquence positive et de la fuite de fréquence négative, lesquels peuvent être pré-calculés pour toute longueur de fenêtre, largeur de bande de filtre, ordre et fréquence centrale sans recours à un ajustement de données (data-driven fitting).

Contributions Clés

1. Décomposition Analytique : L'article établit que les erreurs de l'ecHT consistent en un biais systématique correctible (déterminé par $G_+$) et un plancher de variance irréductible (déterminé par $G_-$).
2. Calibrage Optimal (cecHT) : Les auteurs introduisent la « ecHT calibrée » (cecHT), qui applique le scalaire dérivé $C_{opt}$ pour éliminer le biais de phase systématique. Cela permet d'obtenir une erreur de phase moyenne quasi nulle dans des conditions stables.
3. Règles de Conception : L'étude fournit des limites explicites pour les erreurs de phase et d'amplitude à l'extrémité et relie les paramètres de conception pratiques (longueur de fenêtre, largeur de bande/ordre du filtre et décalage de la fréquence centrale) au biais résiduel via un délai de groupe à l'extrémité.
4. Efficacité Computationnelle : Le calibrage ajoute un surcoût négligeable (une seule multiplication complexe) au pipeline existant de l'ecHT basé sur la FFT en $O(L \log L)$, préservant ainsi la compatibilité avec le temps réel.

Résultats
L'efficacité de la cecHT a été validée par des simulations et des données réelles :

• Simulations : Dans les simulations de signaux monocromatiques, la cecHT a systématiquement réduit l'erreur de phase moyenne d'environ 8,53° (ecHT non calibrée) à 0,40°, et la réduction de l'erreur maximale de 11,70° à 1,23°. Le calibrage a rendu l'algorithme virtuellement indépendant des choix de paramètres (bande passante, ordre du filtre, dérive de fréquence), permettant de régler les paramètres pour d'autres contraintes (ex: robustesse au bruit) sans sacrifier la précision de la phase.
• Ondes Alpha EEG : Appliquée au jeu de données des ondes Alpha de l'EEG, la cecHT a réduit l'erreur de phase moyenne de 4,9° à 0,0° ($p < 10^{-6}$) tout en préservant la valeur de verrouillage de phase (PLV) et en réduisant significativement l'indice de retard de phase (PLI), indiquant l'élimination d'une asymétrie de phase.
• Tremblement Essentiel : Dans l'analyse des données de tremblement, le calibrage a réduit l'erreur de phase moyenne de -15,9° à 5,2°. En combinant cela avec des mises à jour périodiques de la fréquence centrale ($f_0$) pour suivre la dérive spectrale, l'erreur moyenne a été réduite à 0,0°, démontrant que les erreurs résiduelles étaient principalement dues à la dérive plutôt qu'intrinsèques à la méthode de calibrage.

Signification et Revendications
L'article affirme que la cecHT offre un « remplacement direct et pratique » pour les pipelines ecHT existants, fournissant un estimateur statistiquement robuste avec des règles de conception claires. En séparant le biais correctible de la fuite intrinsèque, le cadre offre à la fois un outil pratique pour une amélioration immédiate et un langage diagnostique pour comprendre les limites de l'estimation de phase.

Les auteurs soulignent que si le calibrage élimine le biais systématique encodé dans la réponse de phase du filtre, il ne peut pas supprimer le plancher d'erreur irréductible causé par la fuite spectrale et la non-stationnarité. Par conséquent, la méthode permet d'obtenir des estimations sans biais dans des conditions stables, le jitter restant reflétant soit les limitations de l'estimateur, soit l'instabilité réelle du signal. L'étude positionne la cecHT comme une solution intermédiaire entre les filtres causaux simples et les modèles adaptatifs complexes, adaptée au traitement de signaux neuraux en temps réel, sous réserve que la stabilité de la fréquence et la robustesse aux artefacts soient gérées.