Bonjour !
J'ai un petit problème...
Sous MatLab, quand on écrit :
x = 1:1000
plot(abs(fft( cos(x/1000 *2*pi)))) soit une période
on obtient un pic fréquentiel d'ordre 2
Maintenant :
plot(abs(fft( cos(x/200 *2*pi)))) soit 5 périodes
on obtient un pic d'ordre 6
Maintenant :
plot(abs(fft( cos(x/100 *2*pi)))) soit 10 périodes
on obtient un pic d'ordre 11, impair, je comprends pas
Y'a surement un truc qui m'échape...
Si vous avez une idée, n'hésitez pas.
Merci !
2006-09-16 11:27:39 · 6 réponses · demandé par jojo 1 dans Sciences et mathématiques ➔ Mathématiques
ça y'est j'ai trouvé. abs(fft(...)) donne la valeur d'amplitude de la composante fréquentielle. j'aurais du comparer abs(fft(...)) et imag(fft(...)) pour pouvoir différencier les sin des cos.
2006-09-17 12:23:01 · update #1
tu dois mal lire la sortie car y a un (+1) sur l'ordre de ton harmonique, c'est pas une coïncidence...
2006-09-17 04:10:22 · answer #1 · answered by Ramis V 7 · 0⤊ 0⤋
Pour le calcul des coeffs de Fourrier, le nombre de valeur servant à leurs déterminations doit être une puissance de 2 alors un graphe entre 0 et 1000 ne comporte pas assez de valeurs donc un risque d'erreur.
Essaye entre 0 et 1023 par exemple avec un pas de 1
2006-09-16 20:25:39 · answer #2 · answered by Anonymous · 1⤊ 0⤋
Ben oui, pour la FFT, la fréquence 0 (c'est à dire le continu) correspond au rang 1. Les index vont de 0 à N-1 pas de 1 à N.
2006-09-17 10:13:56 · answer #3 · answered by Champoleon 5 · 0⤊ 0⤋
Une possibilité (?) : Pour plot(abs(fft( cos(x/500 *2*pi)))) , as-tu deux ou trois pics ? Vu la périodicité, le décompte ne prend t'il pas deux fois le pic "début période, fin période" ?
C'est juste une idée : Je n'ai pas utilisé Matlab depuis deux ans
2006-09-16 13:17:33 · answer #4 · answered by Obelix 7 · 0⤊ 0⤋
Peut-être est-ce en rapport avec l'unité pour ton intervalle au sujet de x : puisque c'est du "1:1000" et pas du "0:1000" c'est ptêt ça qui bloque ? À mon avis il faudrait déjà revoir l'aide par "help fft", ça doit y être détaillé.
Bon, je ne sais pas trop, vu que je n'ai pas touché à Matlab, mais juste un peu à Scilab (et encore, pile depuis deux semaines) : une solution pour toi serait d'aller voir sur le forum assez fréquenté de développez.com, qui est à la page http://www.developpez.net/forums/forumdisplay.php?f=148
Bon clavier !
2006-09-16 12:23:25 · answer #5 · answered by myahoo 3 · 0⤊ 0⤋
trop compliquer pour moi
2006-09-16 11:42:26 · answer #6 · answered by Rouky 4 · 0⤊ 1⤋