Changements récents - Rechercher:

edit SideBar

Introduction À Lanalyse Temps-échelle Et Temps-fréquence

Page issue de la présentation de Bruno Torresani (LATP, Université de Provence, Marseille) au groupe de travail Signal MEEG, le 12 janvier 2006. http://www.cmi.univ-mrs.fr/~torresan/

Introduction

Le principe des représentations temps-fréquences (ou temps-échelle) appliquées au traitement du signal est de donner des points de vue "alternatifs" sur le signal étudié.

T : x(t)\to X(\lambda)\ ,\quad\lambda\in\Lambda

On distingue deux types de problématiques:

  • Analyse, détection, estimation: on cherche \`a savoir ce que le signal

contient, éventuellement caractériser la présence ou l'absence de certaines composantes.

  • Traitement (séparation de composantes, débruitage,...): dans ce cas,

il est nécessaire de pouvoir inverser la transformation

Illlustration: mélange d'un son de flute et d'un chant d'oiseau

Séparation des deux composantes:

Il existe deux types de représentations temps-fréquence:

  • Représentations bilinéaires: Fonction de Wigner
W_x(b,\nu) = \int \overline{x}(b-\tau/2) x(b+\tau/2) e^{-2i\pi\nu\tau}\,d\tau

interprétable comme spectre local pour les signaux non-stationnaires, mais difficile à interpréter pour des signaux un tant soit peu complexes. De plus, pas de transformation inverse.

  • Représentations "atomiques": nécessitent de choisir un outil d'analyse

(une fenêtre) g:

G_x(b,\nu) = \int x(t)\overline{g}(t-b) e^{-2i\pi\nu t}\, dt

Il existe une (infinité de) transformation(s) inverse(s).

L'alternative temps-échelle

Les représentations temps-fréquence s'interprêtent comme "spectre variant dans le temps". Dans certains cas, il est plus naturel d'introduire une variable d'échelle plutôt que de fréquence.

  • Fonctions de Wigner affines: bilinéaires,
  • Représentations atomiques: ondelettes
T_x(b,a) = \frac1{\sqrt{a}}\, \int x(t) \overline{\psi}\left(\frac{t-b}a\right)\,dt
  • Variantes: entre temps-fréquence et temps-échelle: un zoo complet de transformations possibles

La transformation de Gabor

On considère une fenêtre g; à tout signal x on associe sa transformée de Fourier à court terme Gx:

G_x(b,\nu) = \int x(t)\overline{g}(t-b) e^{-2i\pi\nu t}\, dt

Il existe une (infinité de) transformation(s) inverse(s); par exemple, pour une fenêtre (dite de synthèse) h, on peut écrire

x(t) = \int\int G_x(b,\nu) e^{2i\pi\nu t} h(t-b)\,dbd\nu\

c'est à dire exprimer x comme combinaison linéaire de Gaborettes

h_{(b,\nu)}(t) = e^{2i\pi\nu t} h(t-b)\

Le rôle de la fenêtre

La fenêtre joue un rôle d'outil d'analyse; deux fenêtres différentes peuvent conduire à des lectures très différentes d'un même signal.

Illustration d'analyse d'un son de clarinette avec une fenêtre étroite:

Avec fenêtre large:

Un autre exemple, sur la phrase "how are you":

Version discrète

Pour des raisons pratiques, il est nécessaire de discrétiser la transformation de Fourier à fenêtre: on obtient la transformation de Gabor

G_x(m,n) = \int x(t) \overline{g}(t-mb_0) e^{-2i\pi n\nu_0 t}\, dt\

b_0 et \nu_0 étant des constantes. Leur choix conditionne l'inversibilité de la transformation x\to G_x. Pour une fenêtre fixée, il existe une valeur critique \tau <1 telle que pour b_0\nu_0 \le\tau, on puisse écrire

x(t) = \sum_{m,n} G_x(m,n) \tilde g(t-mb_0) e^{2i\pi n\nu_0 t}\ ,

\tilde g est une fenêtre duale.

Il existe une infinité de fenêtres duales, dont certaines sont calculables explicitement.

Le compromis de la redondance

Plus le paramètre de redondance b_0\nu_0 est petit:

  • plus l'image de la transformée Gx du signal x est précise
  • plus la fenêtre de synthèse \tilde g peut être prise proche de la

fenêtre d'analyse.

Mais

  • plus la transformée est redondante (les coefficients G_x(m,n)

sont corrélés): les modifications de la transformée temps-fréquence sont plus difficiles à contrôler.

  • plus le temps de calcul est élevé.

Exemple de la fenêtre Gaussienne, avec deux valeurs de b_0\nu_0:

La transformation en ondelettes

Les ondelettes jouent un rôle assez similaire à celui des "Gaborettes": elles fournissent des représentations alternatives des signaux.

Une ondelette est une fonction \psi qui oscille suffisamment pour être d'intégrale nulle. Les ondelettes associées sont des copies translatées et dilatées de \psi

\psi_{(b,a)}(t) = \frac1{\sqrt{a}}\,\psi\!\left(\frac{t-b}a\right)\ .

La transformation en ondelettes (continue) associe à tout signal x la fonction de deux variables

T_x(b,a) = \frac1{\sqrt{a}}\,\int x(t)\overline{\psi}\!\left(\frac{t-b}a\right)\,dt\ .

Pour des signaux suffisamment oscillants, la variable d'échelle joue un rôle comparable à l'inverse de la variable de fréquence de la transformation de Gabor.

Ondelettes: discrétisation

  • Discrétisation "naturelle": discrétisation dyadique:
\psi_{jk}(t) = s^{-j/2}\,\psi\!\left(s^{-j}t -k\tau\right)\ ,\quad a_j = s^j,\ b_{jk} = k\tau s^{-j}

Pour des \psi,s,\tau bien choisis, ces ondelettes peuvent former des bases orthonormées de l'espace des signaux: aucune redondance.

  • Alternative "populaire": Littlewood-Paley: dyadique en échelle,

uniforme en temps: représentations nécessairement redondantes, mais invariantes par translation.

\psi_j^k(t) = s^{-j/2}\,\psi\!\left(s^{-j}(t -k\tau)\right)\ ,\quad a_j = s^j,\ b_{k} = k\tau

Ondelettes vs Gaborettes

Les Gaborettes ont une taille constante (de même que leur transformée de Fourier: résolution temps-fréquence uniforme).

Les ondelettes sont de forme constante, et de taille variable:

  • Bonne résolution temporelle aux petites échelles: grande précision pour

détecter et caractériser des év\`enements localisés.

  • Bonne résolution fréquentielle aux grandes échelles.

A quoi ça sert ?

  • Représentations alternatives des signaux, à partir desquelles certaines caractéristiques peuvent être mises en évidence, et estimées.
  • Représentations parcimonieuses: les signaux "intéressants" sont souvent caractérisés par un petit nombre de coefficients de Gabor ou d'ondelettes. Application naturelle à la compression, et au débruitage (par seuillage).
  • Modélisation de la variabilité, par transformations.

Illustration: Gabor

Illustration: Ondelettes

Illustration: Ondelettes et bruit

Éditer - Historique - Imprimer - Changements récents - Rechercher
Page mise à jour le 31 janvier 2007 à 12h11