<page>
    <title>Traitement numérique du signal</title>
    <id>111193</id>
    <revision>
      <id>31793682</id>
      <timestamp>2008-07-22T15:47:57Z</timestamp>
      <contributor>
        <username>Sylenius</username>
        <id>65080</id>
      </contributor>
      <minor />
      <comment>Annulation des modifications 31792422 par [[Special:Contributions/Gmp.sti|Gmp.sti]] ([[User talk:Gmp.sti|Discuter]])</comment>
      <text xml:space="preserve">Le '''traitement numérique du signal''' étudie les techniques de [[Traitement du signal|traitement]] ([[filtrage]], [[compression]], etc), d'analyse et d'interprétation des signaux numérisés. À la différence du traitement des signaux [[analogique]]s qui est réalisé par des dispositifs en [[électronique]] analogique,  le traitement des signaux [[numérique]]s est réalisé par des machines numériques (des [[ordinateur]]s ou des circuits dédiés). 
&lt;br&gt;Ces machines numériques donnent accès à des algorithmes puissants, tel le calcul de la [[transformée de Fourier]]. Les signaux naturels n’étant généralement pas numériques, il faut les numériser par un [[convertisseur analogique-numérique]].

== Avantages du traitement numérique ==

En comparaison du traitement du signal analogique, le traitement numérique présente un certain nombre d'avantages :
* Immunité au bruit : un signal codé numériquement n'est pas bruité lors de son traitement. On peut donc réaliser des calculs de précision arbitraire et mettre en cascade un grand nombre d'opérations.

* Souplesse : un traitement numérique est facilement ajustable ou re-paramétrable en cours de fonctionnement. Le traitement peut même s'adapter de lui même à la situation (évolution du signal d'entrée au cours du temps). L’aptitude des machines numériques à enchaîner des séquences d’opérations facilite également la réalisation d’algorithmes de traitement complexes. Un exemple important est la [[transformée de Fourier]].

* Mémorisation : la facilité de mettre un signal en mémoire permet de réaliser des retards facilement et donc une grande variété de filtres ou des opérations de [[corrélation]]. On a également accès par ce biais à des traitements itératifs, fonctionnant par raffinements successifs, ce qui ouvre considérablement les possibilités.

Ces avantages sont tempérés par quelques limites :
* La numérisation en elle-même peut dégrader le signal, du fait d'une [[quantification]] ou d'un [[échantillonnage (signal)|échantillonnage]] insuffisants. De plus, les calculs numériques ne sont pas exempts d'erreurs, du fait des arrondis. Les calculs effectués en virgule fixe sont particulièrement vulnérables à ces problèmes.

* Le traitement numérique est en général plus lent, plus consommateur et plus gourmant en ressources matérielles qu'une approche analogique. Cependant, les progrès de la microélectronique  réduisent de plus en plus l'impact de ces défauts et autorisent le développement du traitement numérique temps réel sur un vaste domaine d'applications, s'étendant jusque à des fréquences élévés (&gt;100 MHz).

== Spécificités analytiques du traitement numérique ==

Alors qu’un signal analogique est représentable [[Mathématiques|mathématiquement]] par une [[fonction continue]] d’une variable continue (par exemple, modélisation d'une [[tension électrique]] variable au cours du [[Le temps en physique|temps]]), un signal numérique est une [[Topologie discrète|suite de nombres]]. Il faut donc des outils mathématiques différents pour les manipuler. Citons :
* la [[transformée en Z]] — l’équivalent discret de la [[transformée de Laplace]] en domaine continu. Un cas particulier en est la [[transformée de Fourier]] d’un signal discret ;
* la [[transformée de Fourier discrète]](TFD) — elle a de nombreuses applications en analyse spectrale et en filtrage, du fait de l’existence d’algorithmes de calcul rapide regroupés sous le vocable de [[transformée de Fourier rapide]].

== Conversion analogique-numérique ==

Puisque les signaux naturels sont presque tous des signaux continus et que les circuits numériques ne manipulent que des données discrètes, il faut d'abord ''transformer'' ces signaux avant de pouvoir leur appliquer un traitement numérique. Cette transformation s'appelle la ''discrétisation'' ou ''numérisation'' ; elle est réalisée par un [[convertisseur analogique-numérique]]. Cette opération comprend deux actions :
* l’[[échantillonnage (signal)|échantillonnage]], qui transforme le signal à support continu en un signal à support discret (une suite de valeurs) ;
* la [[Quantification (signal)|quantification]], qui remplace les valeurs continues par des [[mathématiques discrètes|valeurs discrètes]].

== Théorème de Shannon-Nyquist  ==
{{Article détaillé|Théorème d'échantillonnage de Nyquist-Shannon}}
L'échantillonnage transforme un signal continu &lt;math&gt;S(t)\,&lt;/math&gt; en un signal discret &lt;math&gt;S_i\,&lt;/math&gt; composé de plusieurs mesures du signal continu, relevées à des instants successifs séparés par un pas temporel constant. Ainsi : &lt;math&gt;S_i = S(i \cdot T_e )\,&lt;/math&gt; où &lt;math&gt;T_e\,&lt;/math&gt; est la ''période d’échantillonnage''. &lt;math&gt;F_e =1/T_e\,&lt;/math&gt; est la fréquence d’échantillonnage.

D’après le [[théorème d'échantillonnage de Nyquist-Shannon]], un signal dont les composantes fréquentielles sont inférieures à &lt;math&gt;F_e/2\,&lt;/math&gt; est échantillonné de façon réversible : il est possible de reconstituer le signal continu à partir des valeurs discrètes. La fonction interpolante utilisée à cette fin est le [[sinus cardinal]] (noté sinc) :

&lt;math&gt;S(t)=\sum_i S_i \cdot \operatorname{sinc}\left( \pi (\frac{t}{T_e}-i) \right)&lt;/math&gt;

De façon à éviter les artéfacts de « repliements de spectre », il est nécessaire d'éliminer les composantes fréquentielles au delà de &lt;math&gt;F_e/2\,&lt;/math&gt;. Cela est réalisé par un filtre (analogique), appelé [[filtre anticrènelage]], possédant en général une coupure raide (d'ordre élevé).

== Filtrage linéaire ==

Un filtre calcule un signal discret &lt;math&gt;y(n)\,&lt;/math&gt; à partir d’un signal discret &lt;math&gt;x(n)\,&lt;/math&gt;. Chaque échantillon de &lt;math&gt;y(n)\,&lt;/math&gt; est une combinaison linéaire d’échantillons de &lt;math&gt;x(n)\,&lt;/math&gt;.

Un filtre est dit causal si l’échantillon &lt;math&gt;y(n)\,&lt;/math&gt; ne dépend exclusivement que des valeurs de &lt;math&gt;x(i)\,&lt;/math&gt; pour &lt;math&gt;i&lt;n\,&lt;/math&gt;, c’est-à-dire des valeurs antérieures. C'est le cas de tous les systèmes réels, qui ne peuvent pas anticiper les valeurs futures. Toutefois, grâce à la mémorisation des signaux rendue possible par la numérisation, on peut réaliser des réponses impulsionnelles d'apparence non-causale.

Un filtre linéaire et invariant est caractérisé par sa réponse impulsionnelle, ''ie''. la réponse à une entrée ne comportant qu’une valeur non nulle. La nature de cette réponse caractérise le filtre qui peut être à réponse impulsionnelle finie (RIF) ou infinie (RII). La sortie d’un filtre est la [[convolution]] de l’entrée par la réponse impulsionnelle. Cette convolution donne un algorithme de calcul d’un filtre RIF. Cette convolution peut être également calculée par [[transformée de Fourier rapide]]. Pour réaliser des filtres de type
RII, on utilise des techniques récursives (la sortie du filtre dépend des échantillons de sortie précédents).

La [[transformée en Z]] de la réponse impulsionnelle est la fonction de transfert du filtre. Sa [[factorisation]] permet de décrire un filtrage par une équation aux différences discrètes. La transformée de Fourier — c'est-à-dire la transformée en Z pour &lt;math&gt;z = e^{\,j2\pi\lambda}\,&lt;/math&gt; — de la réponse impulsionnelle est la « [[réponse en fréquence]] » du filtre : elle permet de passer dans un domaine spectral. &lt;math&gt;\lambda\,&lt;/math&gt; est la fréquence réduite : &lt;math&gt;\lambda  = F/Fe\,&lt;/math&gt;.

Le calcul des coefficients d’un filtre pour obtenir une réponse en fréquence spécifiée est appelé « synthèse du filtre. »

== Spectre unilatéral ==
Tout signal périodique de période &lt;math&gt;T_0=1/F_0&lt;/math&gt; se décompose en [[série de Fourier]] :
&lt;math&gt;S(t)=A_0 + \sum_{n=1}^\infty A_n   cos(2\pi\cdot n\cdot F_0\cdot t) +  \sum_{n=1}^\infty  B_n sin(2\pi\cdot n\cdot F_0\cdot t )&lt;/math&gt;
Ce qui conduit directement au spectre discret du signal.

== Analyse spectrale ==

La [[transformée de Fourier discrète]] (TFD) est utilisée pour  calculer numériquement le [[Spectre (physique)|spectre]] d'un signal, c’est-à-dire sa représentation fréquentielle. L’utilisation de cette TFD nécessite cependant quelques précautions. D’une part, parce qu’elle n’est applicable que sur un signal de durée limité : le signal doit donc souvent être tronqué, ce qui entraîne l’apparition d’ondulations parasites sur le spectre — elles peuvent être atténuées par les techniques d'[[apodisation]]. D’autre part, parce que le spectre obtenu est lui-aussi échantillonné, ce qui rend son interprétation plus difficile. Une [[interpolation (mathématiques)|interpolation]] de ce spectre peut être souhaitable.

En pratique, on utilise pour réaliser la transformée de Fourier discrète un [[Algorithmique|algorithme]] connu sous le nom de [[transformée de Fourier rapide]] (FFT pour ''Fast Fourier Transform'' en [[anglais]]). Il consiste à ordonner astucieusement les calculs de manière à réduire le nombre total de multiplications nécessaires. La complexité de l'algorithme passe alors de &lt;math&gt;N^2&lt;/math&gt;
à &lt;math&gt;N.\log_2(N)&lt;/math&gt; (N étant le nombre d'éléments du tableau).

== Mise en œuvre ==

Le traitement de signaux déjà numérisés peut se faire sur des ordinateurs d’usage général. Dans les applications de traitement en temps réel, des microprocesseurs spécialisés sont utilisés ([[Digital Signal Processor|DSP]]). Pour plus de rapidité, des fonctions spécifiques sont réalisées directement sous la forme de circuits intégrés ([[Application Specific Integrated Circuit|ASIC]]) ou implantées sur des composants en logique programmable ([[FPGA]]).

== Signal causal ==

&lt;math&gt;s(t)=0&lt;/math&gt; pour &lt;math&gt;t&lt;0&lt;/math&gt;

Un signal est dit causal si ce signal est nul quand t&lt;0 

Par opposition, un signal non-causal est défini sur l'ensemble des réels : le signal existe même pour t&lt;0.

== Bibliographie ==
*elliot price traitement numerique du signal pour debutant
* Maurice Bellanger, ''Traitement numérique du signal'', [[éditions Masson|Masson]], [[Paris]], 2006 {{ISBN|2100501623}}
* Murat Kunt, ''Techniques modernes de traitement numérique des signaux'', [[Presses Polytechniques et Universitaires Romandes]], [[Lausanne]], [[1991]] {{ISBN|2-88074-297-2}}
* A.W.M. Van Den Enden, N.A.M. Verhoeckx, ''Traitement numérique du signal : une introduction'', [[éditions Masson|Masson]], [[Paris]], [[1992]]  {{ISBN|2100073753}}
* René Boite, Henri Leich, ''Les Filtres numériques — Analyse et synthèse des filtres unidimensionnels'', Collection [[CNET]]-[[ENST]], [[éditions Masson|Masson]] {{ISBN|2225818843}}

== Voir aussi ==
=== Articles connexes ===
* [[Distribution de Dirac]]
* [[Sinus cardinal]]

=== Liens externes ===
* [http://www.e-dsp.com Des projets de traitement du signal, des tutoriaux et des livres pour s'instruire]

[[Catégorie:Traitement numérique du signal|*]]

[[ar:معالجة الإشارة الرقمية]]
[[da:Digital signalbehandling]]
[[de:Digitale Signalverarbeitung]]
[[en:Digital signal processing]]
[[es:Procesamiento digital de señales]]
[[fi:Digitaalinen signaalinkäsittely]]
[[he:עיבוד אותות דיגיטלי]]
[[hu:Digitális jelfeldolgozás]]
[[it:Elaborazione numerica dei segnali]]
[[ja:デジタル信号処理]]
[[mn:Тоон дохио боловсруулалт]]
[[ms:Pemprosesan isyarat digital]]
[[nn:Digital signalhandsaming]]
[[no:Digital signalbehandling]]
[[pl:Cyfrowe przetwarzanie sygnałów]]
[[ru:Цифровая обработка сигналов]]
[[sl:Digitalno procesiranje signalov]]
[[sv:DSP-processor]]
[[th:การประมวลผลสัญญาณดิจิทัล]]
[[tr:Sayısal sinyal işleme]]
[[zh:数字信号处理]]</text>
    </revision>
  </page>