transformée de fourier sous matlab

[heruhur]

Salut à tous

voilà, j'ai un probleme en progammation (la question est dans le nom du sujet).
Si quelqu'un saurait comment on écrit une transformée de fourier sous matlab.

Merci d'avance
-----

[invite208c538e]

tape:
help fft
ça devrait répondre à ta question

[heruhur]

merci je vais essayer

[invite208c538e]

si c'est pour afficher le spectre d'un signal, utilise ensuite fftshift pour recentrer ton graphe.
help fftshift

[A voir en vidéo sur Futura]

[invite53b1262f]

Bonjour à tous! j'ai un problème, peut être avez vous des idées.

Alors voila j'ai deux images sources, que je décompose chacune suivant leur composante Rouge, Vert et Bleu.
Ensuite je fais la différence entre composantes, en gros :

IR=I2R-I1R
IV=I2V-I1V
IB=I2B-I1B

qui forment l'image de la différence IDIFF (fichier joint)

En gros, disons que je fais fftshift(fft2(X))) de ces trois images IR, IV et IB ;
ensuite abs(ifft2(X)) des trois spectres obtenus ;
enfin la reconstruction d'une image IOUT avec les trois résultats (fichier joint)

Quand je compare IDIFF à IOUT je ne trouve pas les mêmes images. Le passage dans la représentation fréquentielle à changé quelque chose, mais je ne sais pas quoi. IOUT à l'air de faire un mix des deux images source, c'est bizarre.


Avez vous une idée de ce qui cause dette différence? est-ce une conséquence irréversible de la FFT2? Est-ce que ca à un rapport avec les énergies?

PS :

max(max(IDIFF)) =>
ans(:,:,1) = 0.6588 (rouge)
ans(:,:,2) = 0.6627 (vert)
ans(:,:,3) = 0.6039 (bleu)

max(max(IOUT)) =>
ans(:,:,1) = 0.7137 (rouge)
ans(:,:,2) = 0.7020 (vert)
ans(:,:,3) = 0.6235 (bleu)

[invite53b1262f]

rebonjour à tous, j'ai trouvé mon problème. En fait les images sources ont, pour chacune de leur composante RGB, un "range" d'intensité de lumière allant de 0 à 1.

Si je fais la différence de 2 images, ce "range" va donc d'une valeur <0 à une valeur <1, hors la FFT demande des valeurs positives seulement. ou bien, sons résultat n'est que positif, a vrai dire je sais plus

le résultat, cest que pour que ca marche il me suffit de faire +1 à chaque composante; ensuite je peut faire ifft2(fft2(X)) ou bien ce que je veux, il me suffit de faire -1 à l'image finale.

je sais pas si vous avez compris mon histoire, en tout cas moi je comprends et ca résoud mon problème, mais si ca peut aider des gens alors tant mieux!

[invite9b8a9ef4]

bonjour
je commence à travailler avec transformée de fourier sous matlab.
ça m'aiderai bien un exemple simple de l'oscillateur de Duffing de la forme dx/dt)2+ 2ε w dx/dt + w[EXP]2 x + α x3 = - dx0(t)/dt
comment visualiser sa reponse en amplitude sous matlab en utilisant la transformée de fourrier
merci

[invite9b8a9ef4]

bjr

comment peut on visualiser la réponse d'un système masse ressort amorti défini soit par son équation différentielle ou sa forme matricielle
dx/dt= Ax+Bf
merci bien