CHAPITRE II
LA TRANSFORMEE DE FOURIER DISCRETE
I. DISCRÉTISATION DE LA FRÉQUENCE
II. PROPRIÉTÉS DE W
III. INFLUENCE DE LA DISCRÉTISATION
IV. PROPRIÉTÉS DE LA TFD
V. CONVOLUTION LINÉAIRE DE SÉQUENCES FINIES
VI. CONVOLUTION PAR TFD
1
Transformée de Fourier d'un signal discret :
La transformée de Fourier d'un signal discret est définie par l'expression suivante :
on reconnait alors l'expression de la décomposition en série de Fourier d'une fonction
périodique de période 1 (ou Fe).
Il sera donc suffisant d'étudier celle-ci sur une période, soit par exemple -≤ f ≤
ou 0 ≤ f < 1.
D'un point de vue pratique, il est impossible d'évaluer sur calculateur la transformée de
Fourier d'un signal de durée finie en raison de l'aspect continu de la variable
fréquentielle f. Il faut donc discrétiser l'axe des fréquences, et évaluer le résultat en ces
points (il faut échantillonner la transformée de Fourier). Le résultat obtenu sera appelé
transformée de Fourier discrète.
I. DISCRETISATION DE LA FREQUENCE
Soit un signal x(k), k=0,...,M-1, dont nous souhaitons évaluer numériquement la
transformée de Fourier pour des fréquences fnéquiréparties (fn= n f). Nous allons
choisir
f = , et pour -≤ f < nous obtenons : -≤ n ≤ - 1
(nous aurions pu choisir 0 ≤ f ≤ 1 et donc 0 ≤ n ≤ N - 1).
Evaluons x(k) par la méthode des rectangles : x(k) = X(f) e2jπkf df devient
x (k) = X(nf) e2jπn
p
f= en utilisant la notation X(n)≜X(nf)
on pose W= e2jπkn/N
II. PROPRIETES DE W
W = e2jπ/N
N
W= e2jnkπ/N
2
W= W. W
W= 1
W= Wpour l entier
W= - 1
W= -W
Orthogonalité : =
III. INFLUENCE DE LA DISCRETISATION
x (k) = e2jπkn/N
p
= e2jπkn/N
= e-2jπ(l-k) n/N
Appelons S : e-2jπ(l-k) n/N)
S=
S=
xp(k) =
Le résultat de ce calcul est le signal x(k) rendu périodique de période N. Si l'on
souhaite que les différentes périodes ne se recouvrent pas, il faut que N soit
supérieur à M-1. La valeur la plus faible pour N est donc N=M.
Retenons simplement :
La TFD n'est utilisable que pour des signaux à durée finie.
Le nombre d'échantillons de la TFD doit être ≥ au nombre d'échantillons du
signal.
3
On appellera TFD d'un signal de durée finie (k = 0,..., N-1)
X(n) = e-2jπkn/N
Ceci doit être évalué pour n=0,…,N-1 (périodicité)
x(k) = e2jπkn/N
Ceci doit être évalué aussi pour k=0,…,N-1 (périodicité)
IV. PROPRIETES DE LA TFD
- Linéarité : soit x1(k) 0 ≤ k ≤ N1-1
x2(k) 0 ≤ k ≤ N2-1
N = max(N1, N2)
TFD (a1 x1(k) + a2 x2(k)) = a1 TFD(x1(k)) + a2 TFD (x2(k))
X1(n) X2(n)
il faut concaténer des zéros (ajouter des zéros) à la fin du signal le plus court.
En général, pour la FFT, on prend N = puissance de 2.
- décalage cyclique
x(k) X(n)
x(k-k0) X(n) e-2jπ
x(k) . e2jπ X(n - n0)
V. CONVOLUTION LINEAIRE DE SEQUENCES FINIES
y(k) =
x(l) ≠ 0 pour 0 ≤l ≤ Nx - 1
g(l) ≠ 0 0 ≤ l ≤ Ng - 1
x(k - l) non nul pour 0 ≤ k - l ≤ Nx - 1
4
l ≤ k ≤ Nx + l - 1
Donc y(k) non nul pour 0 ≤ k ≤ Nx + (Ng - 1) - 1
pour k [0, Nx + Ng - 2] c'est-à-dire
y(k) possède Nx + Ng - 1 échantillons non nuls.
* Si l'on veut faire le calcul par TFD (puis TFDI) :
- on doit prendre le même nombre de points N pour les deux signaux
(Concaténation de zéros) jusqu'à Nx + Ng - 1 échantillons.
- N > Nx + N g - 1
- sinon le calcul n'est pas correct.