Discussion :

[johnnybgood1602]

Bonjour à tous,

Je me trouve face à un problème comme vous vous en doutez, je vous expose d'abord le contexte :

Je dispose de deux signaux sinusoïdaux enregistrés par des géophones. Le premier à été placé près d'une source vibratoire et l'autre à une distance donnée.





Ce n'est pas moi qui ai écrit le début du programme, vous remarquerez qu'il manque les titres des axes, en abscisse ce sont les indices des point (que je n'ai pas encore repassé en temps ^^) et en ordonnée l'amplitude en volt.

L'objectif est d’éliminer l'effet de la source sur la réponse du second géophone, on m'a pour cela conseillé d'utiliser la fonction xcorr. Je l'ai donc appliqué à mes deux signaux ...

Code :

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

1
2
[Ycorel, Xdecalage]=xcorr(WaveformData_V(:,1),WaveformData_Source_VLT(:,1));
plot(Xdecalage,Ycorel);

n'étant pas calé en traitement du signal et ne sachant pas ce que la fonction xcorr fait réellement je ne sais pas si c'est juste et quoi en déduire !

Que dois-je faire pour remonter au signal du second géophone sans l'influence de la source ?

Je vous remercie par avance pour votre aide !

[jyber]

Regardée à cette échelle, l'inter-corrélation nous dit juste qu'il s'agit de deux signaux d'enveloppes rectangulaires. Pour le visualiser, imaginez que le second rectangle glisse sur l'axe des temps. Lorsqu'il est à gauche du premier, le produit des deux est nul. Le produit augmente au fur-et-à-mesure que les deux rectangles se superposent. Il atteint un maximum lorsque les deux coïncident, puis décroît pour atteindre à nouveau zéro lorsqu'il est à droite. D'où ce triangle.

Il faudrait expliciter "éliminer l'effet de la source" et regarder ces signaux de plus près (zoom). Comment la propagation modifie le signal? Que cherchez-vous?

[johnnybgood1602]

Bonjour,

Ce que je cherche a avoir c'est le retard du second géophone par rapport à la source pour connaitre la vitesse de propagation de l'onde.

Merci !

[jyber]

Ok, mais quel est l'indice sur lequel on s'appuie pour mesurer ce retard? Le début des signaux? Il est zéro pour les deux. Il n'y a pas de retard, ce que confirme la corrélation qui "pic" à zéro.

Si la source est en régime permanent (fréquence constante) on peut mesurer le déphasage, qui se traduit en retard mais avec une ambiguïté de 2pi (un décalage de 0.5*longueur d'onde, 1.5*longueur d'onde, 2.5*longueur d'onde donneront le même déphasage).

Si la source est modulée, on peut corréler après démodulation mais n'étant pas géologue je préfère m'abstenir de faire des suppositions hasardeuses.

Votre programme MATLAB fait bien ce qu'on attend de lui mais ce n'est pas ce que vous cherchez.