Discussion :

[driou]

Bonjour,

J'ai une courbe de pluviométrie et une courbe d'énergie sismiques. Les deux sont représentées en cumulé avec un point toute les heures pour l'énergie et un point toutes les 30min pour la pluvio.

J'aimerais calculer un coefficient de corrélation pour voir si il existe une relation entre mes deux jeux de données.

Comment pourrais-je procéder sous Matlab?

Merci pour votre aide.

[christophe_halgand]

Une fois que tu nous auras donné ton algorithme, on pourra t'aider dans la programmation sous matlab mais ce n'est pas à nous de faire tes statistiques...

Christophe

[Jean Dumoncel]

Salut,

tu peux faire une recherche :
lookfor correlation
pour voir ce qui est déjà implémenté dans matlab.

[driou]

J'ai bien trouvé les fonction xcorr et corrcoef mais j'avoue ne pas bien comprendre leur fonctionnement. Elles me renvoient une matrice alors que j'aimerais simplement un "coefficient de similarité" entre mes deux courbes.

Pour etre plus clair j'ai deux jeux de données (les vitesse de vent et l'enerie simique) toutes deux dépendent d'un vecteur temps différent tel que:

Code :

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

1
2
3
4
5
6
7
8
9
wind_cum = cumsum(wind_speed);
E_sism_cum = cumsum(E_sism);
time_sism % est le vecteur temps associé à l'énergie sismique (pas d'échantilonage de 1h)
time_wind % tps associé aux vitesses de vent (pas d'echantillonage de 30min)
 
plot(time_sism,E_sism_cum)
hold on;
plot(time_wind,wind_cum)
set(gca,'Xlim',[time_wind(1) time_wind(end)])

j'ai donc tenté

Code :

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

[X,Y] = xcorr(wind_cum,E_sism_cum);

J'obtiens bien une matrice mais je ne comprends pas sa signification

j'ai égalment tenté

Code :

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

corrcoef(wind_cum,interp1(time_sism,E_sism_cum,time_wind);

dans ce cas j'obtient NaN,(je suis obligé d'interpoler pour avoir le même nombre de valeur dans mes 2 vecteurs).

J'aimerais simplement trouver un coefficient de comparaisons me permettant de montrer que le vent et résponsable de l'augmentation de l'énergie simique.


Merci

[Jerome Briot]

Tu pourrais nous montrer ces deux courbes ?

[driou]

Voici les courbes.
En fait il y en a trois AVR et AVM sont les deux courbes d'énergies de deux stations sismologique et la courbe en vert correspond aux vitesses de vent. (toutes les courbes sont normalisées)

Après quelques tests, si j'ai bien compris, xcorr permet de tracer un coefficient de correlation en fonction du décalage (en l'occurence temporel) entre mes deux courbes.
tel que

Code :

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

[Ycorel, Xdecalage]=xcorr(wind_cum,E_sism_cum);

ceci est donc interessant pour moi si je veux connaitre le décalage temporel entre l'action du vent et la réponse des stations mais ceci ne me donne pas vraiment d'indication sur la similarité des courbes.

Corrigez moi si je me trompe

[driou]

et voici la courbe de correlation

Code :

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

1
2
[Ycorel, Xdecalage]=xcorr(wind_cum,E_sism_cum);
plot(Xdecalage,Ycorel)

[abidineb]

Bonjour

Juste une information, le coefficient de corrélation ne tient pas compte des unités, donc c'est pas la peine d'écrire, les unités (heure et min)

Cordialement