版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

cb123123

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 392
帖子: 142
在線: 26.3小時(shí)
蟲號: 3041743
注冊: 2014-03-12
專業(yè): 導(dǎo)航、制導(dǎo)與傳感技術(shù)

[求助] 相位變換法計(jì)算時(shí)延

采集到四個(gè)通道的音頻信號后利用MATLAB進(jìn)行處理，使用相位變換法（PHAT）計(jì)算時(shí)延，可是得出來的結(jié)果卻都是0，程序如下求高手指點(diǎn)
a=wavread('KinectAudio 211');
b=wavread('KinectAudio 212');
c=wavread('KinectAudio 213');
d=wavread('KinectAudio 214');
l1=length(a);
e=1;
l2=e*l1;
signal=zeros(l2,1);
signal1=zeros(l2,1);
signal2=zeros(l2,1);
signal3=zeros(l2,1);
for m=1:e;
for n=1:l1;
signal(n+(m-1)*l1)=a(n);
signal1(n+(m-1)*l1)=b(n);
signal2(n+(m-1)*l1)=c(n);
signal3(n+(m-1)*l1)=d(n);
end
end
signallength=length(signal);
p=1;
for n=1:signallength;
if 2.^n<(signallength*2);
p=p+1;
end
end
assumesignallength=2.^p;
signal4=zeros(assumesignallength,1);
signal5=zeros(assumesignallength,1);
signal6=zeros(assumesignallength,1);
signal7=zeros(assumesignallength,1);
for n=1:signallength;
signal4(n)=signal(n);
signal5(n)=signal1(n);
signal6(n)=signal2(n);
signal7(n)=signal3(n);
end
cccorrelation=xcorr(signal4,signal5);
cccorrelation1=xcorr(signal4,signal6);
cccorrelation2=xcorr(signal4,signal7);
gcc=zeros((assumesignallength*2-1),1);
gcc1=zeros((assumesignallength*2-1),1);
gcc2=zeros((assumesignallength*2-1),1);
phatfilter=zeros((assumesignallength*2-1),1);
phatfilter1=zeros((assumesignallength*2-1),1);
phatfilter2=zeros((assumesignallength*2-1),1);
crossspectrum=fft(cccorrelation);
crossspectrum1=fft(cccorrelation1);
crossspectrum2=fft(cccorrelation2);
for n=1

assumesignallength*2-1);
phatfilter(n)=abs(crossspectrum(n));
phatfilter1(n)=abs(crossspectrum1(n));
phatfilter2(n)=abs(crossspectrum2(n));
gcc(n)=crossspectrum(n)/phatfilter(n);
gcc1(n)=crossspectrum1(n)/phatfilter1(n);
gcc2(n)=crossspectrum2(n)/phatfilter2(n);
end
gcccorrelation=ifft(gcc);
gcccorrelation1=ifft(gcc1);
gcccorrelation2=ifft(gcc2);
for n=1

assumesignallength*2-1);
gcccorrelation(n)=abs(gcccorrelation(n));
gcccorrelation1(n)=abs(gcccorrelation1(n));
gcccorrelation2(n)=abs(gcccorrelation2(n));
end
[gccmaximum,gcctime]=max(gcccorrelation);
[gccmaximum1,gcctime1]=max(gcccorrelation1);
[gccmaximum2,gcctime2]=max(gcccorrelation2);
gccestimation=abs(assumesignallength-gcctime);
gccestimation1=abs(assumesignallength-gcctime1);
gccestimation2=abs(assumesignallength-gcctime2);
gccestimation,gccestimation1,gccestimation2
lag=zeros((assumesignallength*2-1),1);
for n=1

assumesignallength*2-1);
lag(n)=assumesignallength-n;
end
subplot(3,1,1);
plot(lag,gcccorrelation,'b')
axis([-inf inf -inf inf]);
legend('MIC12');
subplot(3,1,2);
plot(lag,gcccorrelation1,'r')
axis([-inf inf -inf inf]);
ylabel('cross-correlation');
legend('MIC13');
subplot(3,1,3);
plot(lag,gcccorrelation2,'g')
axis([-inf inf -inf inf]);
legend('MIC14');
xlabel('time lag');

回復(fù)此樓

» 猜你喜歡

085600材料與化工求調(diào)劑已經(jīng)有16人回復(fù)
317求調(diào)劑已經(jīng)有6人回復(fù)
304求調(diào)劑已經(jīng)有4人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有7人回復(fù)
265求調(diào)劑已經(jīng)有5人回復(fù)
279分求調(diào)劑一志愿211 已經(jīng)有3人回復(fù)
求調(diào)劑已經(jīng)有3人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有12人回復(fù)
286分人工智能專業(yè)請求調(diào)劑愿意跨考！已經(jīng)有3人回復(fù)
329求調(diào)劑已經(jīng)有5人回復(fù)

1樓 2014-07-10 08:47:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

智能機(jī)器人

Robot (super robot)

我們都愛小木蟲

找到一些相關(guān)的精華帖子，希望有用哦~

零相位濾波器、傅里葉變換及逆變換濾波器已經(jīng)有4人回復(fù)
【求助】在matlab中如何用希爾伯特變換定義相位已經(jīng)有22人回復(fù)

點(diǎn)擊這里搜索更多相關(guān)資源

科研從小木蟲開始，人人為我，我為人人

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 cb123123 的主題更新