Discussion :

[sw7gm]

Bonjour à tous !

Je suis nouveau sur le forum (mais je me débrouille assez bien sur Matlab).
Alors voilà le but de mes manips est de déterminer le retard entre un signal d'émission et 4 signaux de réception.
J'utilise l'intercorrelation pour mes mesures mais en fait je ne suis pas sûr que ce soit vraiment ce que je cherche.

J'ai fait un schéma en PJ : si en rouge est mon signal d'émission, et en bleu & vert réception; quel(s) seraient les meilleurs outils selon vous pour arriver à déterminer les 2 retards (rouge - bleu et rouge - vert).

Une autre question : est-ce qu'il faut mesure selon vous "pic à pic" ou "pied à pied" ?

Merci par avance pour votre aide !

[FLB]

Salut,
pour affiner la mesure de l'instant de réception je pense que tu peux t'en sortir avec un produit de convolution

Code :

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

doc conv;

Ensuite si tu as plusieurs pics, il faut que tu cherches les différents maximums locaux (tu as surement une information de la séparation minimale entre 2 maximums consécutifs?).
Je penche pour du pic à pic, avec la convolution ce sera d'autant plus marqué.
La complexité du filtrage (convolution, etc) dépend beaucoup de la qualité de ton signal : si tu as vraiment les courbes que tu nous montres, il suffit de chercher 4 différents maximums locaux.

[sw7gm]

Bonjour FLB,
Helas non les signaux que j'obtiens sont plus compliqué que les gaussiennes représentées, pour te montrer j'obtiens ca :

Le traitement du signal est assez ardu (sinus avec enveloppe gaussienne comme tu peux le voir). Et ma question sur le pic/pic ou pied/pied vient du fait que pour moi je ne sais pas si le capteur "sent" passer l'onde, quand le signal devient different de 0 ou quand il "sent" le maximum d'énergie...
Dans un premier cas il n'y a aucun moyen de détecter calculer le temps de retard, dans le second je m'en sors avec l'intercorrélation.

Pour l'émission le signal est en fait une moitié de sinus-gaussienne (coupée en 2 centrée en 0).

Ok tu parles de conv, j'ai utilisé xcorr jusqu'ici c'est plus précis tu ne penses pas ?

[FLB]

Ton problème est pas simple!
Tu peux nous montrer ce que donne le résultat de la corrélation?
Je pense que la convolution t'aidera plus que la corrélation, mais faut voir à l'usage...

[sw7gm]

Voici le résultat des 4 correlations
Comme mon picoscope ne possède que 4 voies j'ai d'abord enregistré les signaux de réception lorsque ceux ci sont excités par l'emetteur.

Puis dans un second temps j'ai enregistré le signal de l'émetteur.
J'ai ensuite fais l'intercorrelation entre l'émission et les 4 réceptions (cf image) grâce à Matlab.

[FLB]

Si tu prends le max de chaque courbe et l'instant associé, tu n'as pas l'information que tu veux?

Pour peaufiner tu peux effectuer une démodulation en amplitude. Matlab a une toolbox qui possède l'outil amdemod. Si tu n'as pas accès à la toolbox ça peut se coder sans grand souci technique.

[sw7gm]

Salut FLB,
L'ennui (c'est pas mal une question de physique) c'est que je ne sais pas si conceptuellement il convient de prendre le maximum de courbe à courbe ou bien dès que les courbes sont différentes de 0 à la fin ou au début (comme montré sur cette image ).

Tu en penses quoi ? C'est vraiment sur le concept de l'onde (du paquet d'onde même pour être plus précis) que je me pose des questions...

Penses-tu toujours que le produit de convolution soit la solution la plus adaptée pour ce problème ? Car au final convolution et intercorrelation c'est presque la même chose non ? Je ne vois pas quelle sera la valeur ajoutée du produit de convolution...

Sinon (désolé ) je ne vois pas à quoi va servir la démodulation en amplitude... et non je viens de vérifier sur mon matlab je n'ai pas accès à cette option..

[FLB]

Salut,
de mon expérience (mais je ne fais plus de traitement du signal depuis quelques années) il vaut mieux chercher les pics, tu seras plus précis, d'autant plus que le bruit de mesure te compliquera la tache pour trouver le pied de courbe.
La grosse différence entre la convolution et l'intercorrelation est que le premier est une fonction (vecteur), alors que le second est un scalaire. En fait ce que tu cherches c'est trouver pour quel offset de temps tu dois décaler le signal émis pour avoir un maximum d'intercorrelation avec le signal recu. Ca revient à faire un produit de convolution.

La démodulation en amplitude te permet de séparer la porteuse (sinus haute fréquence) de ta gaussienne (qui détient ton information). Mais tu n'en a pas forcement besoin car sur ta figure de ton post d'hier tu peux facilement chercher le maximum, qui correspond quasiment au maximum qu'aurait la gaussienne sans le sinus : tu n'as pas besoin d'interpoler si la précision te suffit.

[sw7gm]

Okay très bien merci FLB, je vais continuer sur cette voie du coup.
Serais-je plus précis :
- si je prends la valeur absolue des 2 signaux avant de les intercorreler
- si je prends la valeur absolue du signal d'intercorrelation ?

[FLB]

Je n'ai pas de réponse définitive à t'apporter sur ces points, une analyse qualitative présente des pours (utiliser la valeur absolu des signaux revient un peu à comparer les puissances des signaux) et des contres (tu perds de l'information). Je ne peux que te conseiller de tester les différentes possibilités et comparer quantitativement les résultats.
Si tu trouves des réponses je suis intéressé

[sw7gm]

OK FLB merci.
Je vais faire mes manipulations ainsi que le post-traitement et je reviendrai vers toi (cela risque tout de même de prendre un peu de temps !).
Je te tiens au courant !

[soft001]

Personnellement je te conseille d'utiliser la fonction xcross
Par exemple si j'ai deux signaux Sign1 et Sign2, pour calculer le retard :

Code :

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

1
2
3
C = xcorr(Sign1,Sign2,'coeff');
[~, I] = max(abs(C));
Retard=I- max(numel(Sign1),numel(Sign2));

Retard : représente le nombre de point et non pas la composante du temps par exemple (l'axe des x)