Discussion :

[Benoît1989]

Bonjour,

J'aimerais savoir comment exactement représenter les résultats d'une fft. En effet, je trouve l'aide de MATLAB peu descriptive à ce sujet.

Voici mon code :

Code :

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

1
2
3
4
5
6
7
8
9
10
 
n=length(a);
y=abs(fft(a,n))/n;
fe=1/(time(2)-time(1));
f=fe*(0:n-1)/n;
 
figure
plot(time,a)
figure
plot(f,y)

où a sont les ordonnées du signal à transformer et time (le temps) représente les abscisses du signal à transformer.

J'aimerais savoir comment déterminer la fréquence d'échantillonnage (fe dans mon code). En effet, on aura pas du tout la même échelle de fréquence (plot(f,y)) en fonction du fe choisi!

Merci et bonne journée!

[Invité]

Bonjour,

En effet, je trouve l'aide de MATLAB peu descriptive à ce sujet.

Je ne suis pas vraiment d'accord avec toi As-tu regardé les exemples de la documentation?

J'aimerais savoir comment déterminer la fréquence d'échantillonnage (fe dans mon code).

N'est-ce pas ce que tu fais avec fe=1/(time(2)-time(1));?
Pour plus de précision, je te conseillerais de prendre peut-être plus de points.

[Benoît1989]

Je me suis peut être mal exprimé : je trouve l'aide MATLAB au sujet de la fft assez bien fournie (d'ailleurs, mon code se base dessus), excepté en ce qui concerne la fréquence d'échantillonnage (fe dans mon exemple, Fs dans l'exemple de MATLAB).

En premier lieu, j'avais donné une fréquence fe de manière arbitraire, et, bien entendu, les résultats présentés (diagramme transformée/fréquence) dépendaient de cela.

Et effectivement, j'ai déjà donné une valeur de fe... mais je ne suis pas certain de mon coup, tout simplement... (point de vue dimensions (s^-1), c'est correct, mais...).


En conclusion, j'aimerais simplement savoir si le choix de fe comme

Code :

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

1
2
 
fe=1/(time(2)-time(1))

est correct. J'ai également une autre question : lorsque tu écris

Pour plus de précision, je te conseillerais de prendre peut-être plus de points.

qu'entends tu par là? Les mesures que j'ai en entrée : a (qui est une accélération en fait) et time (le temps) sont des signaux discrets issus de mesure en laboratoire (en soufflerie, pour être plus précis) et je ne vois pas ce que la fréquence d'échantillonnage peut apporter comme précision supplémentaire à ma transformée.

Merci et bonne soirée!

[Invité]

Eh bien je peux comprendre que MATLAB ne s'étende pas trop sur je sujet vu que c'est une donnée que tu es censé avoir, étant donné que c'est toi-même qui fait les mesures
Oui ce que tu fais est tout à fait correct.
Pour ce qui est de ma remarque, je te renvoie à tes TP de physique lorsque tu prenais plusieurs intervalles sur ton oscilloscope pour calculer la période de ton signal.