Discussion :

[Invité]

Bonsoir à tous,


je suis nouveau ici, je suis étudiant et j'ai besoin d'utiliser Matlab pour un projet.

Je rencontre une difficulté, et j'aimerai savoir si vous pourriez m'aider :

J'ai un fichier de données que je charge sous Matlab, il s'agit d'un simple signal temporel.
Je souhaite détecter la présence d'oscillations régulières dans ce signal, j'ai donc besoin de faire une transformée de Fourier.

J'ai trouvé des algorithmes permettant de tracer une transformée de Fourier d'un signal, mais je ne souhaite pas un résultat graphique :
Je souhaiterai savoir si il serait possible, au lieu d'obtenir une courbe, d'avoir les résultats dans un tableau ?

Ce tableau aurait pour indice la fréquence et contiendrait l'amplitude du signal de la transformée.

En effet je souhaite pouvoir traiter ces données pour détecter des oscillations, et un tableau de valeurs me serait très utile pour pouvoir analyser les données dans des boucles...

Je sais pas si vous comprenez ce que je veux dire, mais si jamais vous avez des solutions ou des algos je suis preneur !!

Merci beaucoup d'avance

[Matersss]

Il me semble que la fonction fft fait ce dont tu as besoin...

Matersss

[Mac52]

Salut,

je travaille également sur les transformées de Fourier, le code que j'utilise généralement :

Code :

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

1
2
3
4
5
% Soit At la trace temporelle à transformer et Aw sa transformée de Fourier
Aw = fft(fftshift(fft(At))).*dt     % "*dt" pour la conservation d'énergie
 
% Transformée de Fourier inverse pour revenir à At
At = ifft(ifftshift(ifft(Aw))).*dw

At et Aw étant des vecteurs tu as directement ton tableau.
Pour une étude précise je te conseille tout de même une étape visuelle intermédiaire pour régler la résolution fréquentielle de Aw.
Si ta plage temporelle est trop large, ta résolution temporelle sera diminuée mais ta résolution fréquentielle sera grande (et réciproquement).
Cela est à considérer si tu as des événements ponctuels subtiles à traiter.

[Invité]

Merci beaucoup à vous deux. J'ai encore quelques réglages à faire pour obtenir ce que je veux du coup je passerai peut être vous demander quelques conseils, mais je pense tenir le bon bout. Bonne continuation !

[Mac52]

Je tiens à m'excuser j'ai écris ma réponse un peu vite.
J'ai inversé certaines fonctions !

la vraie réponse :

Code :

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

1
2
3
4
5
% Transformée de Fourier
Aw = ifftshift(ifft(ifftshift(At)));
 
% Transformée de Fourier inverse
At = fftshift(fft(fftshift(Aw)));

Ca m'apprendra à répondre trop vite