04.08.2013 Views

La Transformée de Fourier - LPSC

La Transformée de Fourier - LPSC

La Transformée de Fourier - LPSC

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Les <strong>Transformée</strong>s du<br />

Traitement du Signal<br />

Cours signaux et systèmes<br />

M1 physique<br />

<strong>Transformée</strong>s <strong>de</strong> <strong>Fourier</strong><br />

Cours signaux et systèmes<br />

M1 physique<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

<strong>Transformée</strong> <strong>de</strong> <strong>Fourier</strong><br />

Définition<br />

Échantillonnage et périodisation<br />

Signaux <strong>de</strong> durée limitée et signaux périodiques<br />

Signaux échantillonnés <strong>de</strong> durée limitée<br />

Signaux discrets<br />

<strong>Transformée</strong> <strong>de</strong> <strong>La</strong>place<br />

Définition<br />

Relation avec la transformée <strong>de</strong> <strong>Fourier</strong><br />

<strong>Transformée</strong> en Z<br />

Les transformées du TS<br />

Définition<br />

Relation avec la transformée <strong>de</strong> <strong>Fourier</strong><br />

Relation avec la transformée <strong>de</strong> <strong>La</strong>place<br />

<strong>La</strong> <strong>Transformée</strong> <strong>de</strong> <strong>Fourier</strong><br />

Pour les signaux à énergie finie<br />

vérifiant :<br />

<strong>La</strong> transformée <strong>de</strong> <strong>Fourier</strong> donne le spectre en fréquence<br />

<strong>de</strong> x(t) :<br />

Le spectre X( f ) est une fonction à valeur complexe :<br />

2<br />

4


Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

<strong>La</strong> <strong>Transformée</strong> <strong>de</strong> <strong>Fourier</strong><br />

Le spectre X(f) est une fonction à valeur complexe :<br />

Forme polaire du spectre X(f) :<br />

Linéarité :<br />

Parité :<br />

Quelques propriétés <strong>de</strong> la TF<br />

x(t)<br />

Réelle et paire<br />

Réelle et impaire<br />

Imaginaire et paire<br />

Imaginaire et impaire<br />

Complexe et paire<br />

Complexe et impaire<br />

Réelle quelconque<br />

Imaginaire quelconque<br />

Partie réelle paire et partie<br />

imaginaire impaire<br />

Partie réelle impaire et<br />

partie imaginaire paire<br />

X( f )<br />

Réelle et paire<br />

Imaginaire et impaire<br />

Imaginaire et pair<br />

Réelle et paire<br />

Complexe et paire<br />

Complexe et impaire<br />

Partie réelle paire et partie<br />

imaginaire impaire<br />

Partie imaginaire paire et<br />

partie réelle impaire<br />

Réelle<br />

Imaginaire<br />

5<br />

7<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

<strong>La</strong> <strong>Transformée</strong> <strong>de</strong> <strong>Fourier</strong> Inverse<br />

TF -1<br />

Inversibilité :<br />

Similitu<strong>de</strong> :<br />

x(t) peut donc être interprété comme<br />

une somme d’harmoniques pures<br />

Quelques propriétés <strong>de</strong> la TF<br />

6<br />

8


Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Translation :<br />

Quelques propriétés <strong>de</strong> la TF<br />

x<br />

-T/2 T/2<br />

x<br />

T<br />

t<br />

t<br />

ReX<br />

ReX<br />

f f<br />

Quelques propriétés <strong>de</strong> la TF<br />

I<strong>de</strong>ntité <strong>de</strong> Parseval :<br />

f<br />

ImX<br />

ImX<br />

Densité Spectrale d’Énergie<br />

f<br />

9<br />

11<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Quelques propriétés <strong>de</strong> la TF<br />

Translation en<br />

fréquence<br />

ou modulation :<br />

x(t)<br />

t<br />

x(t) cos(2πf p)<br />

t<br />

|X|<br />

|X mod |<br />

TF et Convolution<br />

Convolution continue<br />

TF<br />

Convolution discrète<br />

TF<br />

f<br />

f<br />

10<br />

12


Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

A<br />

Quelques signaux et leur TF<br />

Séries <strong>de</strong> <strong>Fourier</strong> complexe<br />

13<br />

15<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Tout signal périodique (T) <strong>de</strong> puissance finie peut être décomposé<br />

en une somme <strong>de</strong> sinus et <strong>de</strong> cosinus.<br />

Continu / Fondamental / harmoniques<br />

Signal créneau :<br />

Séries <strong>de</strong> <strong>Fourier</strong><br />

Temporel : Fréquentiel :<br />

A n =0<br />

1(4/π)<br />

1+ 3 (4/3π)<br />

1+ 3+5 (4/5π)<br />

1+ 3+ 5 + 7 (4/7π)<br />

T=5τ Regraduons l ’axe <strong>de</strong>s n<br />

en fréquence ...<br />

14<br />

16


Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

<strong>Transformée</strong> <strong>de</strong> <strong>Fourier</strong> <strong>de</strong>s signaux<br />

continus périodiques<br />

Signal périodique ⇔ Spectre discret<br />

Du périodique au non-périodique<br />

17<br />

19<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Et la puissance d’un signal périodique ?<br />

I<strong>de</strong>ntité <strong>de</strong> Parseval<br />

Densité Spectrale <strong>de</strong> Puissance<br />

<strong>Transformée</strong> <strong>de</strong> <strong>Fourier</strong><br />

<strong>de</strong>s signaux continus<br />

18<br />

20


Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Module / Argument<br />

Parties réelle & imaginaire<br />

amplitu<strong>de</strong><br />

amplitu<strong>de</strong><br />

1<br />

0.5<br />

0<br />

-0.5<br />

1<br />

0.8<br />

0.6<br />

0.4<br />

0.2<br />

Représentation <strong>de</strong> la TF<br />

|X(ω)|<br />

Arg(X(ω))<br />

représentation <strong>de</strong>s signaux en temps<br />

-1<br />

0 1 2 3 4 5 6 7 8 9 10<br />

temps (sec)<br />

représentation <strong>de</strong>s signaux en fréquence<br />

0<br />

0 1 2 3 4 5 6 7 8 9 10<br />

fréquence (Hz)<br />

21<br />

23<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

amplitu<strong>de</strong><br />

amplitu<strong>de</strong><br />

1<br />

0.5<br />

0<br />

-0.5<br />

représentation <strong>de</strong>s signaux en temps<br />

-1<br />

0 1 2 3 4 5 6 7 8 9 10<br />

temps (sec)<br />

1<br />

0.8<br />

0.6<br />

0.4<br />

0.2<br />

amplitu<strong>de</strong><br />

amplitu<strong>de</strong><br />

représentation <strong>de</strong>s signaux en fréquence<br />

0<br />

0 1 2 3 4 5 6 7 8 9 10<br />

fréquence (Hz)<br />

1<br />

0.5<br />

0<br />

-0.5<br />

représentation <strong>de</strong>s signaux en temps<br />

-1<br />

0 1 2 3 4 5 6 7 8 9 10<br />

temps (sec)<br />

1<br />

0.8<br />

0.6<br />

0.4<br />

0.2<br />

représentation <strong>de</strong>s signaux en fréquence<br />

0<br />

0 5 10 15<br />

fréquence (Hz)<br />

20 25 30<br />

22<br />

24


Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

0.4<br />

0.2<br />

0<br />

amplitu<strong>de</strong><br />

amplitu<strong>de</strong><br />

1<br />

0.5<br />

représentation <strong>de</strong>s signaux en temps<br />

0<br />

-5 -4 -3 -2 -1 0 1 2 3 4 5<br />

temps (µsec)<br />

1<br />

0.8<br />

0.6<br />

0.4<br />

0.2<br />

représentation <strong>de</strong>s signaux en fréquence<br />

0<br />

0 1 2 3 4 5 6 7 8 9 10<br />

fréquence (MHz)<br />

[e]<br />

-0.2<br />

0 50<br />

t (ms)<br />

100<br />

(dB)<br />

50<br />

40<br />

30<br />

20<br />

10<br />

DFT<br />

PSD [e]<br />

0<br />

0 1 2 3 4 5<br />

f (kHz)<br />

25<br />

27<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

amplitu<strong>de</strong><br />

Énergie (dB)<br />

2<br />

1<br />

0<br />

-5 -4 -3 -2 -1 0 1 2 3 4 5<br />

temps (µsec)<br />

10<br />

0<br />

-10<br />

-20<br />

-30<br />

représentation <strong>de</strong>s signaux en temps<br />

représentation <strong>de</strong>s signaux en fréquence<br />

0 1 2 3 4 5 6 7 8 9 10<br />

fréquence (MHz)<br />

et l'énergie d'un signal ?<br />

I<strong>de</strong>ntité <strong>de</strong> Parseval :<br />

Densité Spectrale d’Énergie<br />

26<br />

28


Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

<strong>Transformée</strong> <strong>de</strong> <strong>Fourier</strong> & Systèmes<br />

Un SLTI va être caractérisé par sa réponse impulsionnelle h(t)<br />

<strong>La</strong> transformée <strong>de</strong> <strong>Fourier</strong> <strong>de</strong> h(t) donne la réponse en fréquence du système H(f)<br />

TF TF TF<br />

L’inverse est aussi vrai<br />

Fréquence d’échantillonnage Fe =1/Te <strong>La</strong> transformée <strong>de</strong> <strong>Fourier</strong> est discrète et donne un spectre périodique<br />

Le spectre est représenté <strong>de</strong> 0 à Fe ou <strong>de</strong> -Fe /2 à Fe /2<br />

Remarque : pour <strong>de</strong>s signaux discrets, on pose T e = 1 = F e, (et n = t)<br />

0<br />

x e (t)<br />

<strong>Transformée</strong> <strong>de</strong> <strong>Fourier</strong><br />

<strong>de</strong>s signaux échantillonnés<br />

NTe<br />

t<br />

TF<br />

0<br />

|X(f)|<br />

1<br />

Te<br />

f<br />

29<br />

31<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Ban<strong>de</strong> passante et largeur <strong>de</strong> ban<strong>de</strong><br />

Ban<strong>de</strong> passante<br />

Caractérise un système<br />

Module <strong>de</strong> la rép. en fréquence<br />

Définie à −3 dB soit Gainmax /√2<br />

ou Pmax /2.<br />

<strong>La</strong>rgeur <strong>de</strong> ban<strong>de</strong><br />

Caractérise un signal<br />

Densité Spectrale<br />

Espace <strong>de</strong>s fréquences utiles !<br />

TF <strong>de</strong>s signaux échantillonnés périodiques<br />

Fréquence d’échantillonnage F e = 1/T e ,<br />

Pério<strong>de</strong> du signal NT e , Fréquence du signal F = 1/NT e<br />

<strong>La</strong> transformée <strong>de</strong> <strong>Fourier</strong> est discrète et donne un spectre périodique et discret<br />

Le spectre est constitué <strong>de</strong> N raies, il est représenté <strong>de</strong> 0 à Fe ou <strong>de</strong> -Fe /2 à Fe /2<br />

Remarque : pour <strong>de</strong>s signaux discrets, on pose T e = 1 (et n = t)<br />

x e (t)<br />

0 NT e 2NTe<br />

<strong>Transformée</strong><br />

<strong>de</strong> <strong>Fourier</strong><br />

X e (f)<br />

0 1<br />

NTe 1<br />

T e<br />

f<br />

30<br />

32


Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Échantillonnage idéal...<br />

⇒ <strong>Transformée</strong> <strong>de</strong> <strong>Fourier</strong>...<br />

... périodisation en fréquence.<br />

Du continu au discret, on<br />

échantillonne<br />

Échantillonnage temporel ⇔ Périodisation en fréquence<br />

Échantillonnage en fréquence ⇔ Périodisation temporelle<br />

x e (t)<br />

TF en Temps Discret<br />

(TFTD)<br />

0 NT t 0<br />

1<br />

f<br />

Périodisation<br />

T<br />

x (t) Te<br />

<strong>Transformée</strong><br />

<strong>de</strong> <strong>Fourier</strong> Discrète<br />

(TFD)<br />

X (f) e<br />

Echantillonnage<br />

en fréquence<br />

0 NT 2NT<br />

0<br />

X(f)<br />

1<br />

NT<br />

1<br />

T<br />

f<br />

33<br />

35<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

x(t)<br />

Périodisation<br />

<strong>Transformée</strong> <strong>de</strong><br />

<strong>Fourier</strong><br />

(TF)<br />

0 T t f<br />

x T (t)<br />

0 T 2T<br />

Série <strong>de</strong><br />

<strong>de</strong> <strong>Fourier</strong><br />

0<br />

X(f)<br />

X e (f)<br />

Et le numérique ?<br />

Un signal numérique est fini (N points).<br />

Pour faire sa TF, on le périodise implicitement<br />

On rajoute éventuellement <strong>de</strong>s zéros (Nz),<br />

pour avoir une TF sur (N+Nz) points.<br />

1<br />

T<br />

2<br />

T<br />

Echantillonnage<br />

en fréquence<br />

f<br />

34<br />

36


Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

TF complexe / réelle<br />

Transformation <strong>de</strong> <strong>Fourier</strong> (TF)<br />

Analyse :<br />

X<br />

TF complexe TF réelle<br />

Synthèse :<br />

Synthèse :<br />

+∞<br />

j ωt x()<br />

t = ∫ X(<br />

ω)<br />

e dω<br />

−∞<br />

+∞<br />

x()<br />

t = ∫ Re X(<br />

ω)<br />

cos( ωt)<br />

0<br />

− ImX(<br />

ω)<br />

sin( ωt)<br />

dt<br />

+∞<br />

1<br />

− jωt<br />

ω<br />

π ∫ e dt<br />

2 −∞<br />

( ) = x()<br />

t<br />

Domaine en temps :<br />

x(t) est complexe, continu et apériodique<br />

t va <strong>de</strong> -∞ à+∞<br />

Domaine fréquentiel :<br />

X(ω) est complexe, continue et apériodique<br />

ω va <strong>de</strong> -∞ à+∞<br />

Fréquences positives <strong>de</strong> ω > 0<br />

Fréquences négatives <strong>de</strong> ω < 0<br />

X<br />

1<br />

( ) ∑ +∞<br />

ω =<br />

2π<br />

n=<br />

−∞<br />

Analyse :<br />

X<br />

+∞<br />

1<br />

− jωt<br />

ω e dt<br />

( ) = ∫ x()<br />

t<br />

π<br />

−∞<br />

Domaine en temps :<br />

x(t) est réel, continu et apériodique<br />

t va <strong>de</strong> -∞ à+∞<br />

Domaine fréquentiel :<br />

X(ω) est complexe, continue et apériodique<br />

ω va <strong>de</strong> 0 à +∞<br />

TFDT complexe / réelle<br />

Transformation <strong>de</strong> <strong>Fourier</strong> à Temps Discret (TFDT)<br />

Synthèse :<br />

Analyse :<br />

TFDT complexe TFDT réelle<br />

x[<br />

n]<br />

=<br />

2π<br />

∫ X(<br />

ω)<br />

e<br />

0<br />

jωn dω<br />

x[<br />

n]<br />

e<br />

− jωn<br />

Domaine en temps :<br />

x[n] est complexe, discret et apériodique<br />

n va <strong>de</strong> -∞ à+∞<br />

Domaine fréquentiel :<br />

X(ω) est complexe, continue et périodique<br />

ω parcours une seule pério<strong>de</strong> <strong>de</strong> 0 à 2π<br />

Fréquences positives <strong>de</strong> ω = 0 à ω = π<br />

Fréquences négatives <strong>de</strong> ω = π à ω = 2π<br />

Synthèse :<br />

Analyse :<br />

x[<br />

n]<br />

X<br />

π<br />

= ∫<br />

0<br />

Re X<br />

− Im X<br />

1<br />

( ) ∑ +∞<br />

ω =<br />

π n=<br />

−∞<br />

( ω)<br />

cos( ωn)<br />

( ω)<br />

sin( ωn)<br />

dω<br />

x[<br />

n]<br />

e<br />

− jωn<br />

Domaine en temps :<br />

x[n] est réel, discret et apériodique<br />

n va <strong>de</strong> -∞ à+∞<br />

Domaine fréquentiel :<br />

X(ω) est complexe, continue et périodique<br />

ω parcours une <strong>de</strong>mi pério<strong>de</strong> <strong>de</strong> 0 à π<br />

37<br />

39<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Série <strong>de</strong> <strong>Fourier</strong> complexe / réelle<br />

Séries <strong>de</strong> <strong>Fourier</strong> ou transformation en fréquences réduites<br />

Synthèse :<br />

x<br />

Analyse :<br />

() ∑ +∞<br />

t =<br />

k=<br />

−∞<br />

1<br />

X[<br />

k]<br />

=<br />

T<br />

complexe réelle<br />

T<br />

X[<br />

k]<br />

e<br />

∫ x()<br />

t e<br />

0<br />

j2πkt<br />

/ T<br />

− j2πkt<br />

/ T<br />

Domaine en temps :<br />

x(t) est complexe, continu et périodique<br />

t parcours une pério<strong>de</strong> <strong>de</strong> 0 à T<br />

Domaine fréquentiel :<br />

X[k] est complexe, discret et apériodique<br />

k va <strong>de</strong> -∞ à+∞<br />

Fréquences positives <strong>de</strong> k > 0<br />

Fréquences négatives <strong>de</strong> k < 0<br />

1<br />

X[<br />

k]<br />

=<br />

N<br />

∑ − N 1<br />

n=<br />

0<br />

dt<br />

Synthèse :<br />

x()<br />

t = ∑ Re X[<br />

k]<br />

cos(<br />

2πkt<br />

/ T)<br />

k 0<br />

− ImX[<br />

k]<br />

sin(<br />

2πkt<br />

/ T)<br />

+∞<br />

=<br />

Analyse :<br />

2<br />

X[<br />

k]<br />

=<br />

T<br />

T<br />

∫ x()<br />

t e<br />

0<br />

− j2πkt<br />

/ T<br />

Domaine en temps :<br />

x(t) est réel, continu et périodique<br />

t parcours une pério<strong>de</strong> <strong>de</strong> 0 à T<br />

Domaine fréquentiel :<br />

X[k] est complexe, discret et apériodique<br />

k va <strong>de</strong> 0 à +∞<br />

NB: Avant synthèse, il faut diviser la valeur <strong>de</strong> ReX[0]<br />

par <strong>de</strong>ux.<br />

TFD complexe / réelle<br />

Transformation <strong>de</strong> <strong>Fourier</strong> Discrète (TFD)<br />

Synthèse :<br />

x[<br />

n]<br />

=<br />

Analyse :<br />

TFD complexe TFD réelle<br />

∑ − N 1<br />

k=<br />

0<br />

X[<br />

k]<br />

e<br />

j2πkn<br />

/ N<br />

x[<br />

n]<br />

e<br />

− j2πkn<br />

/ N<br />

Domaine en temps :<br />

x[n] est complexe, discret et périodique<br />

n parcours une pério<strong>de</strong>, <strong>de</strong> 0 à N-1<br />

Domaine fréquentiel :<br />

X[k] est complexe, discret et périodique<br />

k parcours une pério<strong>de</strong> <strong>de</strong> 0 à N-1<br />

Fréquences positives <strong>de</strong> k = 0 à k=N/2<br />

Fréquences négatives <strong>de</strong> N/2 à N-1<br />

Synthèse :<br />

Analyse :<br />

x[<br />

n]<br />

N / 2<br />

= ∑<br />

k=<br />

0<br />

2<br />

X[<br />

k]<br />

=<br />

N<br />

dt<br />

Re X[<br />

k]<br />

cos( 2πkn<br />

/ N)<br />

− Im X[<br />

k]<br />

sin( 2πkn<br />

/ N)<br />

∑ − N 1<br />

n=<br />

0<br />

x[<br />

n]<br />

e<br />

− j2πkn<br />

/ N<br />

Domaine en temps :<br />

x[n] est réel, discret et périodique<br />

n parcours une pério<strong>de</strong>, <strong>de</strong> 0 à N-1<br />

Domaine fréquentiel :<br />

X[k] est complexe, discret et périodique<br />

k parcours une <strong>de</strong>mi pério<strong>de</strong> <strong>de</strong> 0 à N/2<br />

NB: Avant synthèse, il faut diviser les valeurs <strong>de</strong> ReX[0]<br />

et ReX[N/2] par <strong>de</strong>ux.<br />

38<br />

40


Signaux et systèmes - M1R Physique UJF F.Montanet<br />

x[<br />

i]<br />

avec:<br />

Re X[<br />

k]<br />

Re X[<br />

k]<br />

= pour k ≠ 0 et<br />

N / 2<br />

Re X[<br />

0]<br />

Re X[<br />

0]<br />

=<br />

N<br />

Re X[<br />

N / 2]<br />

Re X[<br />

N / 2]<br />

=<br />

N<br />

et :<br />

Im X[<br />

k]<br />

=<br />

+<br />

N / 2<br />

∑<br />

k=<br />

0<br />

N / 2<br />

∑<br />

k=<br />

0<br />

Re X[<br />

k]<br />

cos( 2πki<br />

/ N)<br />

Im X[<br />

k]<br />

sin( 2πki<br />

/ N)<br />

ImX[<br />

k]<br />

= −<br />

N / 2<br />

TFD réelle<br />

TFD réelle : TFD réelle<br />

inverse :<br />

k ≠<br />

N / 2<br />

Re X[<br />

k]<br />

=<br />

N−1<br />

∑<br />

i=<br />

0<br />

Im X[<br />

k]<br />

= −<br />

x[<br />

i]<br />

cos( 2πki<br />

/ N)<br />

N−1<br />

∑<br />

i=<br />

0<br />

x[<br />

i]<br />

sin( 2πki<br />

/ N)<br />

• Asymétrique<br />

• Cas particulier <strong>de</strong> k=0 et k=N/2<br />

• Mathématiquement peu pratique<br />

Transformation <strong>de</strong> <strong>La</strong>place<br />

Cours signaux et systèmes<br />

M1 physique<br />

41<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

TFD complexe :<br />

X[<br />

k]<br />

=<br />

=<br />

=<br />

−<br />

1<br />

N<br />

N−1<br />

∑<br />

k=<br />

0<br />

N−1<br />

∑<br />

k=<br />

0<br />

N−1<br />

∑<br />

k=<br />

0<br />

N−1<br />

∑<br />

n=<br />

0<br />

x[<br />

n]<br />

e<br />

TFD complexe<br />

− j2πkn<br />

/ N<br />

TFD complexe inverse :<br />

x[<br />

n]<br />

X[<br />

k]<br />

e<br />

Re X[<br />

k]<br />

j2πkn<br />

/ N<br />

Im X[<br />

k]<br />

=<br />

1<br />

N<br />

N−1<br />

∑<br />

n=<br />

0<br />

x[<br />

n]<br />

( cos(<br />

2πkn<br />

/ N)<br />

+ jsin(<br />

2πkn<br />

/ N)<br />

)<br />

( sin(<br />

2πkn<br />

/ N)<br />

− jcos(<br />

2πkn<br />

/ N)<br />

)<br />

( cos(<br />

2πkn<br />

/ N)<br />

− jsin(<br />

2πkn<br />

/ N)<br />

)<br />

Transformation <strong>de</strong> <strong>La</strong>place<br />

<strong>La</strong> définition mathématique rigoureuse<br />

<strong>de</strong> la transformée <strong>de</strong> <strong>La</strong>place est la suivante:<br />

42<br />

44


Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Transformation <strong>de</strong> <strong>La</strong>place<br />

Analyse non plus seulement en terme <strong>de</strong> composantes fréquentielles,<br />

mais en termes <strong>de</strong> sinusoï<strong>de</strong>s et d’exponentielles<br />

TF:<br />

TL:<br />

Le plan complexe ou domaine s<br />

Exponentielles<br />

décroissantes<br />

jω<br />

s<br />

Exponentielles<br />

croissantes<br />

σ<br />

45<br />

47<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Le plan complexe ou domaine s<br />

jω<br />

s<br />

σ<br />

NB: on trouve souvent la variable complexe notée p au lieu <strong>de</strong> s<br />

dans les ouvrages francophones<br />

46<br />

Le plan complexe ou domaine s<br />

jω<br />

s<br />

Fréquences<br />

positives<br />

σ<br />

Fréquences<br />

négatives<br />

48


Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Le plan complexe ou domaine s<br />

A<br />

A’<br />

Domaine s Signaux associés<br />

B<br />

B’<br />

C<br />

C’<br />

Transformation <strong>de</strong> <strong>La</strong>place<br />

Domaine temporel<br />

Domaine Fréquentiel Domaine s<br />

TF TL<br />

49<br />

51<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

1 ère étape:<br />

Signal dans le domaine temporel x(t)<br />

2 ème étape:<br />

Multiplication par e -σt<br />

pour tout σ∈]−∞, ∞[<br />

3 ème étape:<br />

Calculer la TF <strong>de</strong><br />

chaque signal pondéré<br />

par une exponentielle<br />

4 ème étape:<br />

Construire la TL dans le plan s<br />

constituée d’une infinité <strong>de</strong><br />

spectres en fréquence juxtaposés<br />

le long <strong>de</strong> l’axe σ<br />

Analyse <strong>de</strong> <strong>La</strong>place<br />

Relations entre <strong>La</strong>place et <strong>Fourier</strong><br />

-a<br />

σ<br />

ϕ<br />

s=jω<br />

50<br />

52


Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

2) Impulsion <strong>de</strong> Dirac<br />

aire sous la courbe = 1<br />

Zéros et Pôles<br />

Signaux Tests<br />

53<br />

55<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

1) Échelon unitaire :<br />

3) Rampe at<br />

Signaux Tests<br />

C’est également la fonction <strong>de</strong> transfert<br />

d’un intégrateur pur.<br />

Signaux Tests<br />

54<br />

56


Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

4) Exponentielle décroissante<br />

Signaux Tests<br />

Propriétés générales<br />

57<br />

59<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

<strong>Transformée</strong>s Classiques<br />

Propriétés générales<br />

58<br />

60


Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Théorèmes <strong>de</strong> la valeur initiale<br />

et <strong>de</strong> la valeur finale<br />

Théorème <strong>de</strong> la valeur finale :<br />

Théorème <strong>de</strong> la valeur initiale :<br />

Propriétés <strong>de</strong>s systèmes linéaires<br />

Superposition<br />

décomposition et synthèse<br />

décomposition synthèse<br />

TF(réponse impulsionnelle) = fonction <strong>de</strong><br />

transfert<br />

61<br />

63<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Transformation <strong>de</strong> <strong>La</strong>place<br />

& Systèmes Linéaires<br />

Un système linéaire va être caractérisé par sa réponse impulsionnelle h(t)<br />

<strong>La</strong> transformée <strong>de</strong> <strong>Fourier</strong> <strong>de</strong> h(t) donne la réponse en fréquence du système H(f)<br />

On peut écrire (1) :<br />

TF TF TF<br />

Réponse impusionnelle<br />

Convolution<br />

L ’inverse est aussi vrai<br />

Soit h(t) la réponse impulsionnelle du système invariant en<br />

temps :<br />

Par linéarité, (1) <strong>de</strong>vient :<br />

62<br />

64


Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Applications aux fonctions <strong>de</strong> transfert<br />

<strong>de</strong>s système linéaires<br />

Entrée x(t) Sortie y(t)<br />

H(s)<br />

Système régit par une équation différentielle<br />

<strong>de</strong> forme générale :<br />

Vrai pour tout système linéaire invariant en temps.<br />

De plus, pour un système physique<br />

le <strong>de</strong>gré du dénominateur est supérieur<br />

à celui du numérateur : m 6 n<br />

Fonction <strong>de</strong> transfert<br />

<strong>La</strong> fonction <strong>de</strong> transfert, définie par la réponse impulsionnelle :<br />

donc :<br />

et la réponse temporelle du système s'obtient<br />

donc par :<br />

65<br />

67<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Fonction <strong>de</strong> transfert<br />

Cas <strong>de</strong> systèmes causals:<br />

x(t) = 0 et y(t) = 0 et toutes les dérivées <strong>de</strong> x et <strong>de</strong> y pour t < 0<br />

En prenant la TL <strong>de</strong> l'équation différentielle<br />

on obtient:<br />

Fonction <strong>de</strong> transfert<br />

Cas <strong>de</strong> conditions initiales non-nulles :<br />

on a <strong>de</strong>s termes supplémentaires dus aux limites en 0+<br />

dans les dérivations :<br />

conditions initiales<br />

même dénominateur que H(s)<br />

P(s) est un polynôme<br />

en s qui dépend <strong>de</strong>s<br />

conditions initiales<br />

66<br />

68


Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Applications aux fonctions <strong>de</strong> transfert<br />

Equations différentielles linéaires<br />

Considérons le système :<br />

Entrée e(t) Sortie s(t)<br />

Processus<br />

dont le processus est régit par les équations différentielles suivantes:<br />

d y<br />

+ 5<br />

dt<br />

dy<br />

+ 6y = e(t)<br />

dt<br />

Conditions initiales : dy<br />

2<br />

2<br />

(0) = 2 et y(0) = 2<br />

dt<br />

A t = 0, on applique une entrée constante e(t) = 6. Que vaut y(t) ?<br />

Application aux fonctions <strong>de</strong><br />

transfert<br />

2<br />

2<br />

2s + 12s + 6 2s + 12s + 6<br />

Y(s) =<br />

=<br />

2<br />

s(s + 5s + 6) s(s + 3)(s + 2)<br />

Décomposition en éléments simples :<br />

2<br />

2s + 12s + 6 a<br />

Y(s) =<br />

=<br />

s(s + 3)(s + 2) s<br />

+<br />

b c<br />

+<br />

s + 3 s + 2<br />

Pour trouver a, on multiplie les 2 membres <strong>de</strong> l’égalité par s et on fait s = 0<br />

Pour trouver b, on multiplie les 2 membres <strong>de</strong> l’égalité par s + 3 et on fait s = -3<br />

Pour trouver a, on multiplie les 2 membres <strong>de</strong> l’égalité par s + 2 et on fait s = -2<br />

On trouve alors :<br />

soit :<br />

1 4 5<br />

Y(s) = - +<br />

s s + 3 s + 2<br />

y<br />

() t<br />

= 1−<br />

4e<br />

− 3t<br />

+<br />

5e<br />

−2t<br />

69<br />

71<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Application aux fonctions <strong>de</strong><br />

transfert<br />

Appliquons les transformées <strong>de</strong> <strong>La</strong>place à :<br />

2<br />

d y dy<br />

+ 5 + 6y<br />

= e(t)<br />

2<br />

dt dt<br />

On obtient :<br />

s Y s 2s<br />

2 5 sY s 2 6Y<br />

s<br />

2<br />

− − + − + =<br />

() ( () ) ( ) E(<br />

s)<br />

Conditions initiales :<br />

e t = 6u(<br />

t)<br />

⇒ E s = 6<br />

On en tire Y(s):<br />

Y(s)=<br />

() ( ) s<br />

2s<br />

2<br />

s(s<br />

+ 12s<br />

+ 6<br />

+ 5s<br />

+ 6)<br />

2<br />

=<br />

2s<br />

2<br />

+ 12s<br />

+ 6<br />

s(s+ 3)(s+<br />

2)<br />

Applications aux systèmes électriques<br />

Soit u(t) la tension aux borne d'un composant, U(s) sa<br />

TL, i(t) le courant qui le traverse et I(s) sa TL:<br />

Résistance :<br />

Capacité:<br />

Inductance:<br />

70<br />

72


Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Applications aux systèmes électriques<br />

R1<br />

i(t)<br />

x(t) R2<br />

y(t)<br />

C<br />

i out =0<br />

Zéros et Pôles<br />

73<br />

75<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Zéros et Pôles<br />

Zéros Pôles<br />

Zéros et Pôles<br />

74<br />

76


Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Analyse <strong>de</strong> filtres dans le domaine s<br />

Circuit Sallen-Key<br />

Pôles en σ + jω<br />

Transformation en z<br />

Cours signaux et systèmes<br />

M1 physique<br />

77<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Analyse <strong>de</strong> filtres dans le domaine s<br />

Relations entre <strong>La</strong>place et TZ<br />

<strong>Transformée</strong> <strong>de</strong> <strong>La</strong>place <strong>de</strong> x[kT], signal échantillonné :<br />

78<br />

80


Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Transformation en Z<br />

<strong>La</strong> transformation en z reprend le formalisme et les propriétés <strong>de</strong> la<br />

transformation <strong>de</strong> <strong>La</strong>place et les applique aux signaux discrets.<br />

<strong>Transformée</strong> en Z = <strong>Transformée</strong> <strong>de</strong> <strong>La</strong>place<br />

<strong>de</strong> la fonction échantillonnée avec z = e Ts ou T est la pério<strong>de</strong> d’échantillonnage.<br />

Transformation <strong>de</strong> <strong>La</strong>place Transformation en Z<br />

Exemple :<br />

e(t)<br />

1<br />

1/2<br />

t<br />

Le plan ou domaine z<br />

81<br />

83<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Définition<br />

Quelques propriétés<br />

Linéarité<br />

Décalage temporel :<br />

Convolution :<br />

Multiplication par série<br />

exponentielle :<br />

<strong>Transformée</strong> en Z<br />

Somme <strong>de</strong> série... donc problèmes <strong>de</strong> convergence !<br />

Propriétés <strong>de</strong> la transformées en Z<br />

Retard :<br />

Cette propriété est particulièrement importante. Elle donne<br />

tout son intérêt à la TZ. En effet, la multiplication par z -1<br />

correspond à la translation en temps d'un échantillon au<br />

suivant :<br />

Avance :<br />

conditions initiales<br />

82<br />

84


Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Propriétés <strong>de</strong> la transformées en Z<br />

Conditions initiales et finales :<br />

Convolution :<br />

Exponentielle :<br />

( )<br />

Initiales : lim X z = lim x[ n]<br />

z→0 n→0 n→∞ z→1 [ ( ) ( ) ]<br />

Finales : lim xn [ ] = lim z−1X z<br />

Soit<br />

alors<br />

n<br />

∑<br />

[ ] = [ − ] ⋅ [ ]<br />

yn hn k xk<br />

k = 0<br />

( ) = ( ) ⋅ ( )<br />

Y z H z X z<br />

<strong>Transformée</strong>s en Z<br />

Sinus et Cosinus : TZ[<br />

cos( at)<br />

u(<br />

t)<br />

]<br />

Signaux test élémentaires<br />

(échelon unitaire, rampe,<br />

paraboles, etc ...) :<br />

TZ<br />

−at<br />

z<br />

[ e u()<br />

t ] = aT<br />

TZ −<br />

z − e<br />

1<br />

2!<br />

1<br />

n!<br />

[ sin( at)<br />

u(<br />

t)<br />

]<br />

f<br />

t<br />

t<br />

( z − cos( aT)<br />

)<br />

z<br />

= 2<br />

z − 2cos(<br />

aT)<br />

z + 1<br />

zsin(<br />

aT)<br />

= 2<br />

z − 2cos(<br />

aT)<br />

z + 1<br />

() t L[<br />

f () t ] Z[<br />

f ( t)<br />

]<br />

u(<br />

t)<br />

tu(<br />

t)<br />

u(<br />

t)<br />

2<br />

u(<br />

t)<br />

n<br />

<br />

1<br />

s<br />

1<br />

2<br />

s<br />

1<br />

3<br />

s<br />

<br />

1<br />

n+<br />

s<br />

1<br />

z<br />

z −1<br />

Tz<br />

2 ( z −1)<br />

N2<br />

() z<br />

2 ( z −1)<br />

<br />

Nn<br />

() z<br />

n ( z −1)<br />

si ces limites existent<br />

() z , N () z , N<br />

() z<br />

avec N2<br />

3<br />

polynômes en z<br />

n<br />

85<br />

87<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Y(<br />

z)<br />

+ a Y(<br />

z)<br />

z<br />

1<br />

<strong>Transformée</strong>s en Z<br />

Échelon unitaire :<br />

Rampe unitaire :<br />

à démontrer<br />

−1<br />

+ ... + a<br />

P<br />

( )<br />

z<br />

TZ ⎡<br />

⎣u() t ⎤<br />

⎦ =<br />

z −1<br />

* −Ts −2Ts −3Ts<br />

U s = 1+<br />

e + e + e + <br />

−T<br />

comme e < 1, la série converge<br />

⎛ 1 ⎞ 1 z<br />

TZ ⎣<br />

⎡u() t ⎦<br />

⎤ = ⎜ ⎟ = =<br />

⎝1−e ⎠ 1−z<br />

z −1<br />

Y(<br />

z)<br />

z<br />

()<br />

TZ ⎡<br />

⎣tut⎤ ⎦ =<br />

−M<br />

= b<br />

0<br />

−Ts −1<br />

Ts<br />

z= e<br />

Tz<br />

( ) 2<br />

z −1<br />

Systèmes différentiels et TZ<br />

<strong>La</strong> sortie y[n] et x[n] d’un système physique linéaire<br />

échantillonné sont reliés par une équation aux différences :<br />

y( n)<br />

+ a1y(<br />

n −1)<br />

+ ... + a My(<br />

n − M)<br />

= b0x(<br />

n)<br />

+ b1x(<br />

n −1)<br />

+ ... + b Nx(<br />

n − N)<br />

y[<br />

n]<br />

=<br />

M<br />

∑<br />

k=<br />

0<br />

a<br />

k<br />

x[n]<br />

Système<br />

y[n]<br />

x[<br />

n<br />

− k]<br />

+<br />

N<br />

∑<br />

k=<br />

1<br />

b<br />

k<br />

y[<br />

n<br />

− k]<br />

X(<br />

z)<br />

+ b X(<br />

z)<br />

z<br />

1<br />

Causal : N ≥ M<br />

la sortie au temps nT dépend <strong>de</strong>s valeurs <strong>de</strong>s échantillons connu à ce<br />

moment c.a.d x[n], x[n-1], x[n-2]… et y[n-1], y[n-2], … et <strong>de</strong>s<br />

caractéristiques du systèmes définies par les a k et b k (système linéaire<br />

causal)<br />

!<br />

−1<br />

+ ... + b<br />

Q<br />

X(<br />

z)<br />

z<br />

−N<br />

86<br />

88


Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Fonction <strong>de</strong> transfert et Pôles<br />

dans le domaine Z<br />

En prenant la TZ et en utilisant les propriétés <strong>de</strong> linéarité et <strong>de</strong> translation en<br />

temps, on obtient :<br />

⎛ M ⎞ ⎛ M ⎞<br />

−k −k<br />

Y( z) ⎜1− bkz ⎟= X ( z) ⎜ akz ⎟<br />

⎜ ∑ ⎟ ⎜∑ ⎟<br />

⎝ k= 1 ⎠ ⎝ k=<br />

0 ⎠<br />

On défini la fonction <strong>de</strong> transfert H(z) qui ne dépend que <strong>de</strong>s a k et b k et qui<br />

résume les propriétés du systèmes :<br />

Y( z) = H( z) X ( z)<br />

Donc :<br />

( )<br />

( )<br />

Y z<br />

H( z)<br />

= =<br />

X z<br />

1−<br />

M<br />

−k<br />

∑ az k<br />

k = 0<br />

N<br />

−k<br />

∑bz<br />

k<br />

k = 1<br />

Pôles dans le domaine Z<br />

H ( z)<br />

=<br />

avec :<br />

z = j /<br />

z<br />

2<br />

1<br />

= − j /<br />

( z − z1)(<br />

z − z2<br />

)<br />

( z − p )( z − p )<br />

LC<br />

LC<br />

1<br />

2<br />

− R +<br />

p1<br />

=<br />

2<br />

R − 4L<br />

/ C<br />

2L<br />

− R −<br />

p2<br />

=<br />

2<br />

R − 4L<br />

/ C<br />

2L<br />

89<br />

91<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Fonction <strong>de</strong> transfert et Pôles<br />

dans le domaine Z<br />

Fonction <strong>de</strong> transfert en z, un rapport <strong>de</strong> polynômes en z avec M ≤ N pour un<br />

système physique.<br />

M<br />

−k<br />

( ) ∑ az k<br />

Y z k = 0<br />

H( z)<br />

= =<br />

N<br />

X ( z)<br />

−k<br />

1−∑<br />

bz k<br />

qui peut se mettre également sous la forme k = 1:<br />

M<br />

∏(<br />

z−zi) P− M i=<br />

1<br />

H( z) = a0z N<br />

∏(<br />

z−pj) j=<br />

1<br />

Les zi sont les zéros et les pj sont les pôles du système.<br />

Les zéros et les pôles décrivent le système à un facteur d’amplitu<strong>de</strong> près a0 (gain<br />

statique) et à un retard près zN-M avec N-M ≤ 0.<br />

En particulier, on appellera<br />

• filtres MA (à Moyenne Adaptée) les filtres n’ayant que <strong>de</strong>s zéros (et <strong>de</strong>s pôles en z = 0)<br />

• filtres AR (Auto Régressifs) les filtres n’ayant que <strong>de</strong>s pôles (et <strong>de</strong>s zéros en z =0)<br />

Relations entre TZ et <strong>Fourier</strong><br />

z = exp(j2pf) ⇒ on restreint z au cercle unité<br />

X() z = x[]exp( k −j2 π f k) = X( f )<br />

z = 1<br />

+∞<br />

∑<br />

k=−∞<br />

On retrouve la transformée<br />

<strong>de</strong> <strong>Fourier</strong> discrète du signal<br />

x[k], et sa périodicité<br />

f=1/2<br />

f=1/4<br />

-1<br />

Im(z)<br />

1<br />

f croissante<br />

f=0<br />

Re(z)<br />

f=1<br />

90<br />

92


Signaux et systèmes - M1R Physique UJF F.Montanet<br />

k<br />

s = r + j2π(<br />

f + )<br />

T<br />

k = −∞,...,<br />

0,...,<br />

+∞<br />

2πF e =2π/T<br />

0<br />

Im(s)=ω<br />

, k entier<br />

Re(s)=r<br />

6πFe<br />

4πFe<br />

2πFe<br />

0<br />

-2πFe<br />

-4πFe<br />

-6πFe<br />

Im(s)<br />

Plan <strong>de</strong> <strong>La</strong>place<br />

Plan <strong>de</strong> <strong>La</strong>place Plan <strong>de</strong>s Z<br />

0<br />

Im(z)<br />

1<br />

Re(s)<br />

f=0<br />

f=1<br />

Re(z)<br />

Im(z)<br />

Plan <strong>de</strong>s Z<br />

1<br />

f=0<br />

f=1<br />

Re(z)<br />

93<br />

Signaux et systèmes - M1R Physique UJF F.Montanet<br />

Interprétation géométrique <strong>de</strong> la TZ<br />

X() z<br />

-a<br />

1<br />

=<br />

z+ a<br />

ρ<br />

0<br />

ϕ<br />

Plan <strong>de</strong>s Z<br />

|a|

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!