calculer la profondeur d'un local a l'aide de la correlation

Version imprimable

06/12/2006, 16h01
mariogarcia

calculer la profondeur d'un local a l'aide de la correlation

Bonjour à tous,
le problème est le suivant, à l'aide de matlab, je dois calculer la profondeur d'un local type garage, je lance un son on m'a dit un 'chirp', et avec le microphone j'imagine je dois capturer le son de retour et calculer ainsi a quel distance je me trouve du bout du local, un peu comme un sonar qui calcule a quel distance se trouve un obstacle, mais je ne suis pas dans l'eau, je suis dans l'air, et je n'ai pas de capteur d'ultrasons, je lance un son bien audible et que la carte son d'un portable puisse gérer.
on m'a dit d'utiliser la corrélation et matlab a fin de pouvoir réaliser cela, on m'a donné aussi la fonction xcorr2, mais je ne vois pas comment convertir, et calculer les distances.
si vous pouvez me donner un peu d'aide, je serai très reconnaissant.

comment produire le son qui va être envoyé, il faudrai avoir un microphone directionnel? comment peut on fabriquer un? ou connaissez vous de bonnes marques pour ce type de microphones?

merci d'avance.
06/12/2006, 16h13
Matthieu Brucher

Tu prends ton signal original, tu cherches le maximum de la correlation entre lui et ce que tu obtiens, tu utilises la fréquence d'échantillonage du signal d'origine pour calculer le temps, puis tu as la distance.
07/12/2006, 00h26
mariogarcia

comment ça chercher le maximum de la corrélation? je ne comprends pas beaucoup ce qui est la corrélation ou la "cross-correlation", surtout pour l'implémenter en matlab.
en tout cas merci,
07/12/2006, 08h05
Matthieu Brucher

Tu n'as pas à l'implémenter, Mtlab le fournit, tu as toi-même donné les fonctions à utiliser.
07/12/2006, 08h27
progfou

Et vu la seconde question posée, tu ne sais pas ce qu'est la corrélation...
Tu as eu des outils pour faire ton système (ton prof ?), un petit :

Code:

doc xcorr2

dans Matlab devrait te fournir une bonne explication, pour le reste, le nombre de sites expliquant cela est suffisamment grand pour que je n'ai pas à mettre un lien. :google2: et tant qu'à faire :recherch: car le forum a déjà plusieurs sujets touchant à la corrélation.
17/12/2006, 21h08
mariogarcia

J ai envoye un signal et j en recois un en retour; j obtiens un delai entre les 2 signaux en echatillon mais je ne sais pas comment obtenir la distance avec le frequence d echantillonage et ce delai.
si vous pouvez m aider je vous serait reconnaissant
18/12/2006, 09h08
PRomu@ld

Un fois que tu as le délai, tu utilises tes formules de lycée :

Code:

v = d/t

Tu connais v, tu connais t, pour calculer d ça ne devrait pas être bien dur ... (attention tout de même, le son est parti , et il est revenu ;) )

que veux tu dire par le son est parti et puis est revenu, j'ai du faire qqch d'illogique dans mon raisonnement car j'ai des resultats pas logiques du tout.

avec matlab j'ai fait le suivant:

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
 
clear all;
close all;
 
%analog output of chirp
 
ao = analogoutput('winsound');
chan = addchannel(ao,1);
duration = 3.4547;
set(ao,'SampleRate',11025)
set(ao,'TriggerType','Manual')
ActualRate = get(ao,'SampleRate');
% len = ActualRate*duration;
x = wavread('sonarping.wav');
putdata(ao,x)
%start(ao)
%trigger(ao)
%waittilstop(ao,5)
 
%capture using a micro and sound card of analog output object
 
AI = analoginput('winsound');
chan_in = addchannel(AI,1);
duration = 3.4547; %acquisition
set(AI,'SampleRate',11025)
ActualRate = get(AI,'SampleRate');
set(AI,'SamplesPerTrigger',duration*ActualRate)
set(AI,'TriggerType','Manual')
% blocksize = get(AI,'SamplesPerTrigger');
Fs = ActualRate;
%start(AI)
%trigger(AI)
 
start([AI ao])
trigger([AI ao])
 
y = getdata(AI);
y = y-mean(y);
plot(x);
%plot(y);
autocorr = xcorr(x);
%autocorr = fftshift(autocorr);
crosscorr = xcorr(x,y);
%crosscorr = fftshift(crosscorr);
[value_auto,position_auto]=max(autocorr);
[value_cross,position_cross]=max(crosscorr);
delay=position_cross-position_auto

est ce que j'ai calculé le delai correctement?
ainsi quels sont vos conseils pour realiser l'expérience, (tests)?
merci d'avance

18/12/2006, 20h40
Matthieu Brucher

Tu ne sais même pas ce que tu fais... Tu ne lis même pas ce qu'on t'indique. Tu n'as pas un problème de dimension dans ce que tu nous donnes ?
18/12/2006, 21h07
progfou

C'est ce que j'avais pressenti dans mon premier post.
Je ne sais pas quel est ton niveau, mais évite de dire que tu as plus que bac+2 :aie:
18/12/2006, 22h01
mariogarcia

déjà; je ne sais pas ce qu'est 'bac +2' je suis en Belgique.

ce que je fais,
ce que je fais c'est d'envoyer un son par la carte son, le récupérer avec un microphone.
je calcule la auto corrélation du signal d'origine, puis la inter corrélation entre le signal d'origine que j'envoie et le signal que je capture avec le microphone. et je compare les deux (c'est cela la corrélation, non?)
vous voyez la fonction 'max' avec celle la, j'essaye de calculer le maximum de la corrélation et je compare pour obtenir le délai.
Je dois faire ceci sur une matière non vue au cours,
je vous prie d'excuser mon ignorance, je ne fais pas des études d'ingénieur, je reconnais que j'ai des problèmes d'incompréhension, mais bon dois je laisser tomber?
18/12/2006, 22h31
Matthieu Brucher

Tant que tu ne te décides pas à suivre tous les conseils qu'on te donne, ça ne servira pas à grand chose. La question que tu poses, on y a répondu à plusieurs reprises.