Discussion :

[NulleenMatlab]

Bonjour ,

Etant donné que je n'arrivais pas à écrire le programme permettant de faire la transformée de Fourier d'une fonction, un ami m'a passé le m-file suivant:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
 % compute horizontal vector (1,n) fourier from -xm to xm with zoom effect on
% frequency, window range=n*dx, dx pixel size, xm=(n-1)*dx/2 ; zoom range
% from 1 to 100
 
function F=fur1(v,xm,zoom)
n=length(v);
fxm=(n-1)^2/4/n/xm/zoom;
x=linspace(-xm,xm,n)
fx=linspace(-fxm,fxm,n)
dx=2*xm/(n-1)
F=dx*v*exp(-2*pi*i*x'*fx);

C'est censé faire la TF 1D de la fonction v, le problème c'est que je ne comprends ni comment s'en servir, ni à quel moment ce programme fait l'intégrale. Et apparemment il y a aussi un paramètre qui permet aussi de faire un zoom ??

L'ami en question est parti en vacances pr 2 semaines, je ne peux donc pas lui demander!!

Quelqu'un pourrait-il m'aider à comprendre?

Merci.

[Dam2227]

Est-ce que tu as regardé du côté de la fonction fft? (dans l'aide de Matlab)

[NulleenMatlab]

Oui mais l'objectif est que je code moi même la transformée de Fourier pas que j'utilise celle de Matlab!!

[StitchP]

Bonjour,

à mon avis partir sur un code déjà à moitié fait n'est pas une bonne idée
surtout si tu n'as pas compris le principe de la transformée de Fourier...

Essaie à la rigueur de regarder le code de la fonction fft de matlab...

Bon courage !

[robotron]

Essaie à la rigueur de regarder le code de la fonction fft de matlab...

Je ne suis pas sûr que ce soit la meilleure idée sachant que l'algorithme qui se cache derrière la FFT de matlab est l'algorithme de Cooley & Tukey qui est assez complexe.

[robotron]

C'est censé faire la TF 1D de la fonction v, le problème c'est que je ne comprends ni comment s'en servir, ni à quel moment ce programme fait l'intégrale. Et apparemment il y a aussi un paramètre qui permet aussi de faire un zoom ??

L'ami en question est parti en vacances pr 2 semaines, je ne peux donc pas lui demander!!

Quelqu'un pourrait-il m'aider à comprendre?

Merci.

Sinon pour répondre à l'auteur du msg :
Je ne vois pas trop pourquoi tu parles d'intégrale. Tu travailles sur un ordinateur, ton signal est donc discret, tu dois donc utiliser la transformée de Fourier discrète (TFD) qui elle n'utilise pas d'intégration mais une somme d'exponentielle.

Ce qu'on t'as donné est peut etre pas mal mais ce bout de programme n'est pas très clair, rien qu'à en voir le nom de la fonction. En gros pour t'en servir, tu dois entrer à la main le signal dont tu souhaites calculer la TFD, c'est à dire 2 matrices, l'une (par ex. t) pour les abscisses (le temps) et l'autre (par ex. x) pour les ordonnées, ensuite tu fais simplement

Code :

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

TFD=fur1(x,t,1)

Et le résultat devrait être dans la matrice TFD. J'ai testé le code rapidement, il ne fonctionne pas mais là j'ai pas trop le temps. Déjà tu n'as pas besoin de la fonction de zoom...

[NulleenMatlab]

Merci de votre aide a tous, mais j'ai a peu près compris , je me débrouillerai...