Summary: La Décomposition Modale Empirique (ou EMD pour Empirical Mode Decomposition) est une méthode de décomposition adaptative, non paramétrique et locale de signaux non stationnaires. On en décrit le principe ainsi que la structure de l'algorithme associé, et on donne un exemple sur des données réelles.
Principe — La technique de Décomposition Modale Empirique (ou EMD pour “Empirical Mode Decomposition”) a été introduite par N.E. Huang et al., “The empirical mode decomposition and Hilbert spectrum for nonlinear and non-stationary time series analysis," Proc. Roy. Soc. London A, Vol. 454, pp. 903–995, 1998), avec pour objectif de décomposer tout signal en une somme de composantes oscillantes extraites directement de celui-ci de manière adaptative. Ces composantes (ou IMF pour “Intrinsic Mode Functions”) s'interprètent comme des formes d'ondes non stationnaires (i.e., modulées en amplitude et en fréquence) pouvant être éventuellement associées à des oscillations non linéaires.
EMD vs. ondelettes
— L'EMD partage avec la transformée en ondelettes (TO) l'idée d'effectuer une analyse récursive dans laquelle un signal
![]() |
On aboutit ainsi à une décomposition de la forme :
où, tout comme pour la TO, l'approximation à l'“échelle”
En contraste notable avec la TO, les “échelles” mises en jeu dans l'EMD ne correspondent pas à une grille fixée a priori mais sont définies à partir du signal de manière locale et adaptative.
Algorithme — Plus précisément, la décomposition (Équation 2) qui effectue la dichotomie entre oscillation rapide et oscillation lente est obtenue en itérant un opérateur non linéaire agissant à une échelle de temps variable définie par la distance entre extrema locaux consécutifs. Cet opérateur, dit de tamisage (ou “sifting”) est décrit par la procédure suivante :
Si on itère
D'un point de vue pratique :
Il peut anisi être décrit par le pseudo-code suivant:
![]() |
Un code Matlab est disponible à l'URL http://perso.ens-lyon.fr/patrick.flandrin/emd.html .
Un exemple réel — La Figure 3 présente un résultat de la décomposition par EMD des données mensuelles de variabilité de la température moyenne globale à la surface de la Terre durant la période 1850–2007, mesurée (en o C) comme écart par rapport à la moyenne durant la période 1961–1990. Ces données sont disponibles à l'URL http://www.cru.uea.ac.uk/cru/data/temperature et on pourra en trouver une analyse plus complète de même nature dans Z. Wu, N.E. Huang, S.R. Long and C.-K. Peng, “On the trend, detrending, and variability of nonlinear and nonstationary time series,” Proc. Nat. Acad. Sc., Vol. 104, No. 38, pp. 14889–14894, 2007. En identifiant et regroupant les différents modes obtenus, la méthode permet de donner une représentation synthétique sous la forme d'une tendance très basse fréquence à laquelle se superposent des oscillations représentatives de cycles permettant, de façon hiérarchique, de construire des approximations successives du signal complet de la résolution la plus grossière à la plus fine, celles-ci étant sélectionnées de manière adaptative à partir des seules données.
![]() |