版塊導(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小時熱門版塊排行榜

返回列表

cb123123

新蟲 (小有名氣)

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

[求助] 用相位變換法（PHAT）計(jì)算時延

采集到四個通道音頻以后在MATLAB中用相位變換法（ＰＨＡＴ）處理得之間的時延，用軟件可以看出四個通道的波形有卻別，MATLAB中的實(shí)驗(yàn)結(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ù)此樓

» 猜你喜歡

請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有32人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有261人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有78人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復(fù)
面向全球申請人的中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)

1樓 2014-07-10 08:53:13

已閱回復(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 的主題更新