Chapitre III
Matériels et Méthodes
Page | 55
environ 4h et 24h, avec respectivement, les fréquences de 0.000069 et 0.000011 Hz. Quant à
la gamme fréquentielle des autres ondes dues à l’élévation soudaine du niveau de la mer, nous
avons retenu celles qui vont de 5 min à 3h, avec respectivement les fréquences de 0.0033 et
0.00009 Hz.
Dans le but de séparer ou extraire les différentes ondes contenues dans le signal brut du
niveau de la mer, nous allons réaliser un filtre avec Butterworth sur MATLAB pour le signal
enregistré avec un pas de mesure de 1 min (0.016Hz).
Le code utilisé pour le filtre pass-band (bande-passante) est le suivant :
load inct.txt;
fs=0.016;
x=inct(:,1);
t=inct(:,2);
plot(t,x) ;
grid
Fc1=0.0000069;
Fc2=0.000011;
[b,a]=butter(2,[Fc1 Fc2]/(fs/2));
f=filter(b,a,x);
hold on
plot(t,f,'r') ;
Où :
Fc1, Fc2 : représentant les fréquences de la bande passante pour la période entre 4h et 24h.
[b,a]=butter(2,[Fc1 Fc2]/(fs/2)): filtre bande passante de type Butterworth.
Matériels et Méthodes
Page | 55
environ 4h et 24h, avec respectivement, les fréquences de 0.000069 et 0.000011 Hz. Quant à
la gamme fréquentielle des autres ondes dues à l’élévation soudaine du niveau de la mer, nous
avons retenu celles qui vont de 5 min à 3h, avec respectivement les fréquences de 0.0033 et
0.00009 Hz.
Dans le but de séparer ou extraire les différentes ondes contenues dans le signal brut du
niveau de la mer, nous allons réaliser un filtre avec Butterworth sur MATLAB pour le signal
enregistré avec un pas de mesure de 1 min (0.016Hz).
Le code utilisé pour le filtre pass-band (bande-passante) est le suivant :
load inct.txt;
fs=0.016;
x=inct(:,1);
t=inct(:,2);
plot(t,x) ;
grid
Fc1=0.0000069;
Fc2=0.000011;
[b,a]=butter(2,[Fc1 Fc2]/(fs/2));
f=filter(b,a,x);
hold on
plot(t,f,'r') ;
Où :
Fc1, Fc2 : représentant les fréquences de la bande passante pour la période entre 4h et 24h.
[b,a]=butter(2,[Fc1 Fc2]/(fs/2)): filtre bande passante de type Butterworth.
