Discussion :

[carterslap]

Bonjour à tous,

je viens vers vous car j'ai un problème. Avant toute chose, je précise que je suis totalement débutant sur Matlab, pas de notions du tout sur le comment ça fonctionne (quelques trucs quand même, mais considérez que vous avez affaire à un neuneu.)

J'ai un appareil hydrographique qui me sort une analyse spectrale de houle réalisé à partir d'un enregistrement brutes à 2hz contenant 1024 données.

J'ai les données brutes, j'ai le spectre, mais je n'ai pas le déroulement du calcul et j'aimerai la retrouver pour pouvoir l'appliquer à d'autre jeu de données.

Je vous joins les données en question: brutes.xlsx et spectre.xlsx

Ce que je sais faire sur Matalab dans ce problème: importer les données brutes et créer une variable, faire une FFT, faire un graphique. Mais je ne sais pas comment paramétrer la FFT pour pouvoir retomber sur celle que l'instrument me sort.

J'espère être assez clair

D'avance merci pour votre aide!

[vbbarent]

Bonjour,
Voici le pas de fréquence vous devriez calculer la fft sur 128 points (df=fréquence d'échantillonnage(2 Hz ici)/nbre de points de la fft). Ensuite on ne considère que les fréquences positives.

Code :

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

1
2
3
4
5
6
7
8
 
fe           = 2;
nfft         = 128;
freq        = (0:nfft-1)*fe/nfft-fe/2;
psd        = (1/(fe*nfft))*abs(fftshift(fft(houle(:,2),nfft)).^2);
spectre   = 4*psd(freq>0); % on tient compte des fréquences négatives
freq        = freq(freq>0); % frequence positive, on laisse la composante continue
figure, plot(freq, spectre);

Notez par contre l'utilisation d'une fenêtre d'opodisation peut faire varier légèrement l'amplitude.

@++

[carterslap]

Bonjour Vbparrent,

tout d'abord merci pour ta réponse!

je reçois un message d'erreur quand je teste ton code --> "Attempted to access houle(:,2); index out of bounds because size(houle)=[1024,1]."

c'est quoi le soucis? Car je créé bien une variable houle à partir des données brutes.?

[carterslap]

ok c'est bon j'ai trouvé, j'importais une variable et non une matrice.

Je reviendrais peut être pour des précisions

encore merci!

[carterslap]

Re,

quand je compare le spectre que mon instrument me donne et celui que tu as fait il est vraiment différent, est-ce qu'il n'y a pas un moyen d'avoir le même ?