Transformation de Hilbert

En mathématiques et en traitement du signal, la transformation de Hilbert, ici notée ${\mathcal {H}}$ , d'une fonction de la variable réelle est une transformation linéaire qui permet d'étendre un signal réel dans le domaine complexe, de sorte qu'il vérifie les équations de Cauchy-Riemann.

La transformation de Hilbert tient son nom en honneur du mathématicien David Hilbert, mais fut principalement développée par le mathématicien anglais G. H. Hardy^[1].

Définition[modifier | modifier le code]

Soit $s$ une fonction définie sur $\mathbb {R}$ à valeur dans $\mathbb {R}$ , on appelle transformée de Hilbert la fonction ${\widehat {s}}$ définie par:

{\widehat {s}}(t)={\mathcal {H}}\{s\}(t)={\text{vp}}\left\{(h*s)(t)\right\}={\text{vp}}\left\{\int _{-\infty }^{\infty }s(\tau )h(t-\tau )~{\rm {d}}\tau \right\}={\frac {1}{\pi }}{\text{vp}}\left\{\int _{-\infty }^{\infty }{\frac {s(\tau )}{t-\tau }}~{\rm {d}}\tau \right\}

où ${\mathcal {H}}$ est la transformation de Hilbert et où

h(t)={\frac {1}{\pi t}}

et

{\text{vp}}\left\{\int _{-\infty }^{\infty }s(\tau )h(t-\tau )~{\rm {d}}\tau \right\}=\lim _{\epsilon \to 0}\left\{\int _{-\infty }^{t-\epsilon }s(\tau )h(t-\tau )~{\rm {d}}\tau +\int _{t+\epsilon }^{+\infty }s(\tau )h(t-\tau )~{\rm {d}}\tau \right\}.

$vp$ étant l'abréviation de valeur principale de Cauchy.

On peut montrer que pour tout réel p > 1, ${\mathcal {H}}$ est un opérateur borné de l'espace $L p$ (ℝ) dans lui-même.

Réponse fréquentielle[modifier | modifier le code]

Il s'ensuit que la transformée de Hilbert d'un signal peut être calculée dans le domaine fréquentiel en remarquant que la transformée de Fourier de la fonction $h$ , notée $H$ est :

H(\omega )={\mathcal {F}}\{h\}(\omega )=-{\rm {i}}\operatorname {sgn}(\omega ),

où

${\mathcal {F}}\{h\}$ désigne la transformée de Fourier de $h$ ,
$i$ désigne l'unité imaginaire (parfois notée $j$ )
$\omega \,$ est la fréquence angulaire, et
$\operatorname {sgn}(\omega )={\begin{cases}\ \ 1,&{\mbox{si }}\omega >0,\\\ \ 0,&{\mbox{si }}\omega =0,\\-1,&{\mbox{si }}\omega <0,\end{cases}}$ qui est souvent appelée fonction de signe.

Ainsi, comme la transformée de Fourier d'une convolution est égale au produit des transformées de Fourier de ses constituants, on a:

{\mathcal {F}}\{{\widehat {s}}\}(\omega )={\mathcal {F}}\{h*s\}(\omega )=H(\omega ){\mathcal {F}}\{s\}(\omega )=-{\rm {i}}\operatorname {sgn}(\omega ){\mathcal {F}}\{s\}(\omega ).

La transformation de Hilbert a pour effet de tourner de +90° la composante de fréquence négative de $s (t)$ et de −90° la composante de fréquence positive. Autrement dit, la transformée de Hilbert tourne de 90° le plan de projection du signal analytique.

Transformée de Hilbert inverse[modifier | modifier le code]

On peut remarquer que $H^{2}(\omega )=1^{-}$ (excepté en $\omega =0$ , où la fonction signe est nulle, ce qui implique $H^{2}(0)=0^{-}$ ). Donc la transformée de Hilbert inverse apparait clairement :

s(t)=(h*h*h*{\widehat {s}})(t)={\mathcal {H}}^{3}\{{\widehat {s}}\}(t)={\mathcal {H}}^{-1}\{{\widehat {s}}\}(t).

Exemples de transformées de Hilbert[modifier | modifier le code]

Remarque. Certains auteurs, par exemple Bracewell 1986, utilisent notre $-{\mathcal {H}}$ comme définition de la transformation présentée dans la suite. La conséquence est qu'il faut prendre l'opposé de la colonne de droite du tableau.

Signal $s(t)\,$	transformée de Hilbert ${\mathcal {H}}\{s\}(t)$
$\sin(t)\,$	$-\cos(t)\,$
$\cos(t)\,$	$\sin(t)\,$
$1 \over t^{2}+1$	$t \over t^{2}+1$
Sinus cardinal $\sin(t) \over t$	$1-\cos(t) \over t$
fonction porte $\sqcap (t)$	${1 \over \pi }\ln \left\|{t+{1 \over 2} \over t-{1 \over 2}}\right\|$
impulsion de Dirac $\delta (t)\,$	${1 \over \pi t}:=\delta (jt)$

Modèle à bande étroite[modifier | modifier le code]

De nombreux signaux peuvent être modélisés par le produit d'un signal harmonique à support borné, $s m (t)$ , et d'une « porteuse » sinusoïdale : $s(t)=s_{m}(t)\cos(\omega t+\varphi ).$

Lorsque $s m (t)$ n'a pas de composante fréquentielle au-delà de la fréquence de la porteuse, ${\tfrac {\omega }{2\pi }}$ Hz, alors : ${\widehat {s}}(t)=s_{m}(t)\sin(\omega t+\varphi ).$

Donc, la transformation de Hilbert peut être simplement vue comme un circuit qui produit un déphasage de 90° de la fréquence de la porteuse. De plus : $(\omega t+\varphi )_{\mathrm {mod} \,2\pi }=\arctan \left({{\widehat {s}}(t) \over s(t)}\right)$ d'où on peut reconstruire la porteuse. Puis le message peut être extrait de $s (t)$ par une démodulation cohérente.

En biologie et particulièrement pour l'analyse des signaux électrencéphalographiques (EEG), la démodulation complexe^[2], l'analyse par transformée de Hilbert^[3]^,^[4], l'analyse en bande étroite^[5], l'analyse de modulation d'amplitude et de phase^[6]^,^[7], ont été appliquées en électroencéphalographie quantitative, en neurosciences et en imagerie cérébrale par cartographie EEG^[8].

Représentation analytique[modifier | modifier le code]

Une représentation analytique d'un signal est définie ainsi à l'aide de la transformée de Hilbert :

s_{a}(t)=s(t)+{\rm {i}}{\widehat {s}}(t).

Par exemple pour le modèle à bandes étroites, la représentation analytique est :

$s_{a}(t)$	$=s_{m}(t)\cos(\omega t+\varphi )+{\rm {i}}s_{m}(t)\sin(\omega t+\varphi )$
	$=s_{m}(t)\left[\cos(\omega t+\varphi )+{\rm {i}}\sin(\omega t+\varphi )\right]$
	$=s_{m}(t){\rm {e}}^{{\rm {i}}(\omega t+\varphi )}$ (par la formule d'Euler)

Cette opération complexe hétérodyne enlève les composantes de fréquence négative de $s m (t)$ . Dans ce cas, la partie imaginaire du résultat est la transformée de Hilbert de la partie réelle. C'est donc une manière indirecte de former une transformée de Hilbert.

Alors que la représentation analytique d'un signal n'est pas nécessairement analytique, il existe un lien avec les fonctions analytiques, qui est en fait la façon dont la transformation de Hilbert est apparue historiquement. L'idée est la suivante. Commençons avec une fonction

f:\mathbb {R} \to \mathbb {R} .

On peut l'étendre à une fonction harmonique sur le demi-plan de Poincaré (demi-plan supérieur du plan complexe), à l'aide d'une convolution avec le noyau de Poisson. Toute fonction harmonique est la partie réelle d'une fonction analytique. Maintenant on considère la partie imaginaire d'une fonction analytique, en particulier sa valeur à la frontière. Il apparaît que les valeurs à la frontière sont ${\mathcal {H}}(f)$ . Il s'ensuit que la fonction analytique peut être décrite à l'aide de l'intégrale de Poisson de $f+{\rm {i}}{\mathcal {H}}(f)$ .

Considérations pratiques[modifier | modifier le code]

La fonction $h$ avec $h (t)= 1 ⁄ (π t)$ est un filtre non causal et donc ne peut pas être faite telle quelle, si $s$ est un signal dépendant du temps. Si $s$ est une fonction à variable non temporelle, par exemple des variables spatiales, la non-causalité ne doit pas être un problème. Le filtre est de plus à support non borné, ce qui peut être un problème dans certaines applications. Un autre problème peut apparaître du fait du comportement à fréquence nulle, ce qui peut être évité en s'assurant que $s$ ne contient pas de composante continue.

Une implémentation pratique implique dans de nombreux cas qu'un filtre à support fini, qui peut de plus être rendu causal grâce à un délai raisonnable, est utilisé pour approximer une simulation informatique. Cette approximation peut aussi impliquer que seule une plage de fréquence spécifique est sujette au changement de phase dû à la transformation de Hilbert.

Transformée de Hilbert discrète[modifier | modifier le code]

Soit $s[n]$ un signal discret et ${\hat {s}}[n]$ sa transformation de Hilbert discrète. On a

{\hat {s}}[n]=s[n]*h[n]

où $h[n]$ représente un noyau de Hilbert discret tel que

h[n]={\begin{cases}0,&n{\text{ pair}}\\{\frac {2}{\pi n}}&n{\text{ impair}}\end{cases}}

.

Références[modifier | modifier le code]

(en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Hilbert transform » (voir la liste des auteurs).
(en) Ronald Bracewell, The Fourier Transform and Its Applications, McGraw-Hill, 1986, 2^e éd..
(en) A. Bruce Carlson, Paul B. Crilly et Janet C. Rutledge, Communication Systems, McGraw-Hill, 2002, 4^e éd..

Voir aussi[modifier | modifier le code]

Articles connexes[modifier | modifier le code]

Liens externes[modifier | modifier le code]

(en) Eric W. Weisstein, « Hilbert Transform », sur MathWorld
(en) Analytic Signals and Hilbert Transform Filters, par Julius O. Smith III, université Stanford
(en) The Hilbert transform, par Mathias Johansson (mémoire de maîtrise, université de Växjö)

Notes et références[modifier | modifier le code]

↑ (en) King, Frederick W., Hilbert Transforms. Volume 1., Cambridge University Press, 2009, 744 p. (ISBN 978-0-511-72145-8, 0511721455 et 9780521887625, OCLC 776965734, lire en ligne), p. 3
↑ Walter D.O. The method of complex demodulation. In D.O. Walter and M.A.B. Brazier M.A.B., eds. Advances in EEG Analysis. Electroenceph. Clin. Neurophysiol., 1968, Suppl.27.
↑ [1]. Pierre Etevenon. Étude méthodologique de l'électroencéphalographie quantitative. Application à quelques exemples. Thèse de doctorat d'état ès sciences naturelles, 29 novembre 1977, Université Pierre et Marie Curie, Paris VI, 303 p., Copédith, Paris, 1978. Archive ouverte HAL Université Pierre et Marie Curie.
↑ John S. Barlow. The ELECTROENCEPHALOGRAM. Its Patterns and Origins. MIT Press, 1993, (ISBN 978-0262023542).
↑ B. Ribemont, P.Etevenon, F. Giannella Une méthode d’approche en électroencéphalographie : l’analyse de Fourier en bande étroite. Septième Colloque sur le Traitement du Signal et ses Applications. GRETSI. 101/1, Nice, 28 mai-2 juin 1979.
↑ P. Etevenon, F. Giannella, F. Abarnou. Modèle d'EEG par modulation radioélectrique d'amplitude et de phase. Conséquences sur l'aspect morphologique du signal. Rev. EEG Neurophysiol. clin. 10, 1, 69 - 80, 1980.
↑ Etevenon, P., Lebrun, N., Clochon, P., Perchey, G., Eustache, F., Baron, J.C. High Temporal Resolution Dynamic Mapping of Instantaneous EEG Amplitude Modulation after Tone-burst Auditory Stimulation. Brain Topography, 12, 2, 129-137, 2000.
↑ H. Witte et al. Using discrete Hilbert transformation to realize a general methodological basis for dynamic EEG mapping. A methodological investigation. Automedica, 13:1-13, jan. 1990.

[1] (en) King, Frederick W., Hilbert Transforms. Volume 1., Cambridge University Press, 2009, 744 p. (ISBN 978-0-511-72145-8, 0511721455 et 9780521887625, OCLC 776965734, lire en ligne), p. 3

[2] Walter D.O. The method of complex demodulation. In D.O. Walter and M.A.B. Brazier M.A.B., eds. Advances in EEG Analysis. Electroenceph. Clin. Neurophysiol., 1968, Suppl.27.

[3] [1]. Pierre Etevenon. Étude méthodologique de l'électroencéphalographie quantitative. Application à quelques exemples. Thèse de doctorat d'état ès sciences naturelles, 29 novembre 1977, Université Pierre et Marie Curie, Paris VI, 303 p., Copédith, Paris, 1978. Archive ouverte HAL Université Pierre et Marie Curie.

[4] John S. Barlow. The ELECTROENCEPHALOGRAM. Its Patterns and Origins. MIT Press, 1993, (ISBN 978-0262023542).

[5] B. Ribemont, P.Etevenon, F. Giannella Une méthode d’approche en électroencéphalographie : l’analyse de Fourier en bande étroite. Septième Colloque sur le Traitement du Signal et ses Applications. GRETSI. 101/1, Nice, 28 mai-2 juin 1979.

[6] P. Etevenon, F. Giannella, F. Abarnou. Modèle d'EEG par modulation radioélectrique d'amplitude et de phase. Conséquences sur l'aspect morphologique du signal. Rev. EEG Neurophysiol. clin. 10, 1, 69 - 80, 1980.

[7] Etevenon, P., Lebrun, N., Clochon, P., Perchey, G., Eustache, F., Baron, J.C. High Temporal Resolution Dynamic Mapping of Instantaneous EEG Amplitude Modulation after Tone-burst Auditory Stimulation. Brain Topography, 12, 2, 129-137, 2000.

[8] H. Witte et al. Using discrete Hilbert transformation to realize a general methodological basis for dynamic EEG mapping. A methodological investigation. Automedica, 13:1-13, jan. 1990.

[1]

[2]

[3]

[4]

[5]

[6]

[7]

[8]