transformée de Fourier sans fft

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

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

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??

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.

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:

---

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
http://img4.hostingpics.net/pics/237496fft.jpg

Merci de me répondre

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.

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

Code:

---

1 2 3 4 5

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

---

Bonjour,

En reprenant les valeurs de votre exemple :

Code:

---

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:

---

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.