phasedelayestimated.txt

本例程是针对不同通道之间采集的声音信号之间波形的相位时延差进行估计

figure;
load A150.TXT;
time=A150(:,1);
A1=A150(:,2);
A2=A150(:,3);
A3=A150(:,4);
%A4=A174(:,5);
len=length(time)
fs=10000;

%A1=filter(b,a,A1);A2=filter(b,a,A2);A3=filter(b,a,A3);A4=filter(b,a,A4);
%time=time((8300*2.5):(8700*2.5));A1=A1((8300*2.5):(8700*2.5));A2=A2(1:22000);A3=A3(1:22000);A4=A4(1:22000);
subplot(421);plot(time,A1);
y=fft(A1,len);
y(1)=0;
z=abs(y(1:len./2));
w=(0:1:(len/2-1))./len.*fs;
subplot(422);plot(w,z);
hold on;
subplot(423);plot(time,A2);
y=fft(A2,len);
y(1)=0;
z=abs(y(1:len./2));
w=(0:1:(len/2-1))./len.*fs;
subplot(424);plot(w,z);
hold on;
subplot(425);plot(time,A3);
y=fft(A3,len);
y(1)=0;
z=abs(y(1:len./2));
w=(0:1:(len/2-1))./len.*fs;
subplot(426);plot(w,z);
hold on;
subplot(427);plot(time,A4);
y=fft(A4,len);
y(1)=0;
z=abs(y(1:len./2));
w=(0:1:(len/2-1))./len.*fs;
subplot(428);plot(w,z);

kb=8000;ko=8010; fb=0;fo=500;
subplot(421);axis([kb ko -200 200]);axis('auto y');
subplot(423);axis([kb ko -200 200]);axis('auto y');
subplot(425);axis([kb ko -200 200]);axis('auto y');
subplot(427);axis([kb ko -200 200]);axis('auto y');
subplot(422);axis([fb fo -200 200]);axis('auto y');
subplot(424);axis([fb fo -200 200]);axis('auto y');
subplot(426);axis([fb fo -200 200]);axis('auto y');
subplot(428);axis([fb fo -200 200]);axis('auto y');


% sound output from speaker
ao=analogoutput('winsound');   
addchannel(ao,1:2);
putdata(ao,[A3 A3]);
start(ao);

%calculate the time delay between each two sensor;
fs=10000;
lag=100;
kb=4.8*fs; ko=6.0*fs;
t=time(kb:ko);
x=A1(kb:ko);
y=A2(kb:ko);
z=A3(kb:ko);
%w=A4(kb:ko);

figure;
subplot(411);plot(t,x); axis([kb/2.5 ko/2.5 -200 200]);axis('auto y');
subplot(412);plot(t,y); axis([kb/2.5 ko/2.5 -200 200]);axis('auto y');
subplot(413);plot(t,z); axis([kb/2.5 ko/2.5 -200 200]);axis('auto y');
%subplot(414);plot(t,w); axis([kb/2.5 ko/2.5 -200 200]);axis('auto y');

figure;
[c,lags]=xcorr(x,y,lag,'unbiased');
[b,a]=max(c),c1=(a-lag-1)/fs,
subplot(221);
plot(lags/fs,c);
xlabel('t');
ylabel('x xcorr y');
grid;

[c,lags]=xcorr(x,z,lag,'unbiased');
[b,a]=max(c),c2=(a-lag-1)/fs,
subplot(222);
plot(lags/fs,c);
xlabel('t');
ylabel('y xcorr z');
grid;

%[c,lags]=xcorr(z,w,lag,'unbiased');
[c,lags]=xcorr(y,z,lag,'unbiased');
[b,a]=max(c),c3=(a-lag-1)/fs,
subplot(223);
plot(lags/fs,c);
xlabel('t');
ylabel('z xcorr w');
grid;

[c,lags]=xcorr(y,w,lag,'unbiased');
[b,a]=max(c),c4=(a-lag-1)/fs,
subplot(224);
plot(lags/fs,c);
xlabel('t');
ylabel('x xcorr w');
grid;