Transformées de Fourier avec fftw3

Bonjour, je vais devoir utiliser pas mal la bibli fftw, et j'ai déjà des problèmes dès les premiers tests.... (oui ça commence bien)

Alors voilà, j'essaye juste de faire une fft sur un cos, la pas de pb j'ai bien un spectre nul partout sauf en deux pics (c'est bien ça que je suis censé avoir??)

Après, c'est la fft inverse qui me pose souci, elle ne me redonne pas le signal d'origine.
(en fait j'ai l'amplitude du signal multipliée par N)
Je me demande si ce n'est pas une question de taille N dans:

p = fftw_plan_dft_1d(N, in, out, FFTW_FORWARD, FFTW_ESTIMATE)
(je ne suis pas sûr de devoir utiliser le même N lors de la fft inverse que pour la fft forward, mais je ne vois pas quelle taille utiliser alors...)

Merci de votre aide...

Salut

Citation:

---

Envoyé par BackInBlack

Après, c'est la fft inverse qui me pose souci, elle ne me redonne pas le signal d'origine.
(en fait j'ai l'amplitude du signal multipliée par N)

---

C'est normale. Entre la fft et la fft inverse, un facteur N (nombre d'échantillon) apparait. A toi de choisir sur qui tu veut l'appliquer.
Tu peut aussi appliquer sqrt(N) sur les deux.

Oui je sais ça, mais je pensais que fftw bacward le prenait en compte....
(je vais le faire à la main :D)

Sinon, j'ai une autre petite question:
je voudrais filter un signal venu d'un accéléromètre, j'ai donc 3 données à traiter (X, Y et Z), est-ce que dans fftw, une fft_3d est équivalente à 3 fft_1d? Et y a t-il une différence au niveau temps de calcul?

Merci

bon... ya pas beaucoup de réponses. Du coup j'ai trouvé tout seul!
Alors je donne la réponse au cas ou un jour, qqun se pose un jour la question aussi et qu'il n'ai pas de réponse non plus:

fftw_3d est equivalent à 3 fftdw_1d

Citation:

---

Envoyé par BackInBlack

fftw_3d est equivalent à 3 fftdw_1d

---

?
pour une matrice de taille M*N*O
fftw_3d ne correspond pas à M*N*O fftdw_1d?

Bonjour, svp j'ai besoin de faire une fft sur un cos
est ce que vous pouvez m'aider ?

merci d'avance.:)

Salut.
As tu vue ce tuto?
http://humbert-florent.developpez.co...aitement/fftw/