Discussion :

[angelmn]

bonjour,
j'essais de faire la transformée de Fourier d'un signal sauf qu'il n'est pas forcément en puissance de 2 donc je ne peux pas utiliser de fft donc je voudrais savoir comment je pourrais faire une transformée de fourier discrète sans fft et qui s'adapte à un signal qui n'est pas en puissance de deux
Merci de votre aide

[Invité]

Bonjour,

La fonction fft accepte tout aussi bien des signaux de taille en puissance de 2 ou non, il n'y a aucune problème
Tu peux consulter sa documentation, et plus particulièrement sa partie algorithmes

[angelmn]

Merci pour votre réponse, en fait je croyais qu'il fallait la faire sur un nombre de points qui soit en 2^n c'est ce qu'on m'a dit du moins , et en plus j'ai testé avec la fft j'ai eu des résultats pas très exacts. je fesais ma fft sur 295 points ce n'est pas de là que viens le problème vous pensez??

[Invité]

C'est en général le cas pour de nombreux algorithmes effectuant le calcul de la FFT, mais MATLAB lève cette obligation en utilisant comme tu as pu le voir différents algorithmes selon le nombre d'échantillons du signal.

Si tu veux affiner ton résultat, tu peux effectuer du zéro padding en spécifiant le deuxième argument "n" de la fonction.

[angelmn]

je vous remercie énormément pour vos réponses, j'ai l'impression que j'utilise mal la fft par exemple quand j'essais de faire la fft d'un signal sinusoïdal de fréquence 100 hz

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
 
>> t=0:0.001:1;
>> x=sin(2*pi*t*100);
>> plot(abs(fft(x)));

je suis sensé obtenir un pic en 100 mais le résultat que j'obtiens n'est pas très exact; j'ai un pic en 101 cette inexactitude pourrais me fausser mes calculs donc je me demande à quoi elle est due


Merci de me répondre

[Invité]

Attention à bien mettre des fréquences en abscisses et non des échantillons

Sinon, tu auras toujours une précision de fs/N (fs : fréquence d'échantillonnage ; N : nombre d'échantillons)
Je t'invite à lire ce passage de Wikipédia.

[angelmn]

je vous remercie, mais si j'ai bien compris je dois modifier l'axe des abscisses comme ceci??

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
 
f=abs(fft(x));
N=length(f);
w=-fs/2:fs/N:fs/2;
plot(w,f)

[Krillz]

Bonjour,

En reprenant les valeurs de votre exemple :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
 
f0  = 100;
fs  = 1000; 
 
t   = 0:1/fs:1;
x   = sin(2*pi*f0*t);
y   = abs(fft(x)); 
f   = fs*linspace(0,1,length(x));
 
plot(f,y);

A noter qu'il suffit de tracer jusqu'à fs/2 :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
 
f0   = 100;
fs   = 1000; 
 
t    = 0:1/fs:1;
x    = sin(2*pi*f0*t);
 
nfft = 2^nextpow2(length(x)); % prochaine puissance de 2 >= à la durée de votre signal
y    = abs(fft(x,nfft)); 
f    = fs/2*linspace(0,1,nfft/2+1);
 
plot(f,y(1:nfft/2+1))

Cordialement,

K.