fft sous matlab

[invite1178bc9c]

Bonjour,
Je ne connais pas matlab et je doit faire un petit programme qui me permette de:
* créer un sinus
* afficher mon sinus dans le domaine temporel
* faire sa fft pour l'afficher dans le domaine fréquenciel
Pour les 2 premiers points j'y suis arrivé par contre je n'arrive pas a faire le passage dans le domaine fréquenciel
Quand je regarde ce qu'il m'affiche cela ne correspond pas du tout a la fréquence voulu.
J'utilise le mot clef y=fft(signal_en temporel);
quelqu'un peut-il m'aider
-----

[inviteb9f49292]

Tu es en numérique, donc au même titre qu'en temporel où ton échantillon n correspond à l'instant Tn = (n - 1) / Fech en secondes, tu as un échantillonnage fréquenciel. De plus, méfies toi la fft de MATLAB te renvoies les fréquencesb positives en premier, puis les fréquences négatives.

[invite1178bc9c]

ok. mais du coup je doit taper quoi, pour que cela ressemble vraiment a mon signal?
Car a part des raies qui s'affichent mais qui ne correspondent pas du tout a ma fréquence mais plutôt au borne de ma fenêtre et un message d'erreur me disant qu'il y a un pb avec la partie complexe de mon signal.

[inviteb9f49292]

ok. mais du coup je doit taper quoi, pour que cela ressemble vraiment a mon signal?

Je pense que la réponse est quelque part dans ton cours... pour t'aider, quelles sont les fréquences min et max que t'affiche la transformée... Quelle est la longueur (le nb de points) de ta transformée?

mais plutôt au borne de ma fenêtre

Je t'ai expliqué que les fréquences négatives se retrouvent après les fréquences positives, un "help fft" devrait te donner toutes les indications supplémentaires.

et un message d'erreur me disant qu'il y a un pb avec la partie complexe de mon signal.

Je te rappelle qu'une transformée de Fourier d'un signal réel est complexe... donc utiliser "plot" directement sur le résultat de fft n'est pas approprié... Mais c'est une supposition puisque tu ne dis pas quelle opération te donne une erreur, ni quel est le message exact.

[A voir en vidéo sur Futura]

[invite1178bc9c]

sur mon signal je fait (de mémoire) car pas sur ce pc
Tn = (n - 1) * Fech pour le domaine fréquentiel
y=fft(signal);
plot(y);
J'ai quelque chose comme cela
Désolé mais se sont mets premières lignes de programmation en matlab. du coup même avec l'aide de matlab j'ai un peu du mal

[invite3c59dc04]

Bonjour Isabelle38,

Voilà : tu fait un nouveau fichier .m avec ça :

fe=100000;
m=300;
f=[0:m-1]*fe/m;
t=[0: (m-1)]/fe;
%Création sigal s
s=5*sin(2*1000*pi*t);
plot(t,s);
figure;
%fft
h=abs(fft(s));
stem(f(1: (m/2)),h(1: (m/2))/(m/2));

Save and Run et tu verras deux figures : l'une le signal de 1000 Hz au point de vue temporel et l'autre le signal dans le domaine fréquentiel à 1000 Hz.

@+

[inviteb9f49292]

sur mon signal je fait (de mémoire) car pas sur ce pc

Pour t'amuser chez toi, installe 'OCTAVE' c'est un clône gratuit de MATLAB

[invite1178bc9c]

merci de vos réponses
je vais essayer au plus vite.
mais pas facile l'apprentissage sur le tas

[Jack]

Pour t'amuser chez toi, installe 'OCTAVE' c'est un clône gratuit de MATLAB

Ou bien scilab:
http://www.scilab.org/

A+

[inviteb9f49292]

Ou bien scilab:

Le problème de scilab c'est qu'il n'est pas franchement compatible avec MATLAB (à ma connaissance en tout cas). OCTAVE en revanche c'est compatible à plus de 99%, tu peux utiliser des scripts de l'un sur l'autre sans problème.

[invite0ab6bca5]

Bonsoir,

Désormais en retraite, je n'utilise plus Matlab mais Scilab ; donc voici, traduit en Scilab, le programme proposé par le précédent intervenant, et à la suite un autre ; en effet, je trouve que pour bien représenter le spectre avec fft, il vaut mieux donner UNE SEULE période du signal.
L'illustration du spectre est avec des fréquences négatives et positives.

Bonne continuation

 //pour forum futura sciences
fe=100000;
m=300;
f=[0:m-1]*fe/m;
t=[0: (m-1)]/fe;
//Création sigal s
s=5*sin(2*1000*%pi*t);
xset("window",1), xbasc(1)
plot(t,s);
figure;
//fft
h=abs(fft(s));
xset("window",2),xbasc(2)
plot(f(1: (m/2)),h(1: (m/2))/(m/2));
//
//autre proposition
dt=1/64 ; 
N=64 ; temps=[-N:N-1]*dt ;
signal=3*sin(2*%pi*temps/(2*N*dt));
xset("window",3),xbasc(3)
plot(temps,signal,'b'),xgrid()
df=1/(2*N*dt) ;
frequence=[-N:N-1]*df;
spectre=fft(signal,-1);
xset("window",4),xbasc(4)
specsor=fftshift(spectre);
plot(frequence,real(specsor),'k',frequence,imag(specsor),'r')