Chapitre II
Matériels et méthodes
55
Nous définissons x avec :
X =Alcudia ( :,1) ;
X : représente la hauteur du niveau de la mer
Alcudia: Le nom du fichier dans le fichier
1 : la colonne que prend x
On utilise le point-virgule pour terminer l’instruction
Par la suite, on va définir la méthode utilisée dans ce travail qui est Pwelch pour
estimer la densité spectrale de puissance.
Langage de programmation Matlab :
La méthode de Welch est représentée sur Matlab par la fonction pwelch. Elle s’écrit
comme suit :
load Alcudia.txt
x=Alcudia(:,1);
fs=0.016;
n = length(x)
[Pxx,F] = pwelch(x,500,250,n,fs)
loglog(F,Pxx)
Où :
 Alcudia : Le titre du fichier contenant les données à utiliser dans les calculs.
 X : Correspond aux variations du paramètre à étudier en fonction du temps et qui
est représenté dans le fichier par la 1ère colonne du tableau des données.
 Fs (Sampling frequency) : La fréquence d’échantillonnage en (Hz).
 N : nombre de point FFT
 [Pxx,F]:le spectre de pwelch
 Pwelch : type de spectre
 500 : la valeur de la fenêtre ou window
 250 : la valeur de recouvrement ou de Noverlap