版塊導航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

返回列表

cb123123

新蟲 (小有名氣)

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

[求助] 相位變換法計算時延

采集到四個通道的音頻信號后利用MATLAB進行處理，使用相位變換法（PHAT）計算時延，可是得出來的結(jié)果卻都是0，程序如下求高手指點
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');

回復此樓

» 猜你喜歡

求調(diào)劑已經(jīng)有3人回復
265求調(diào)劑已經(jīng)有4人回復
085700資源與環(huán)境308求調(diào)劑已經(jīng)有6人回復
一志愿吉林大學材料學碩321求調(diào)劑已經(jīng)有12人回復
286分人工智能專業(yè)請求調(diào)劑愿意跨考！已經(jīng)有3人回復
329求調(diào)劑已經(jīng)有5人回復
申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復已經(jīng)有4人回復
材料學碩318求調(diào)劑已經(jīng)有5人回復
一志愿中國海洋大學，生物學，301分，求調(diào)劑已經(jīng)有6人回復
081700化工學碩調(diào)劑已經(jīng)有3人回復

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

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

智能機器人

Robot (super robot)

我們都愛小木蟲

找到一些相關的精華帖子，希望有用哦~

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

點擊這里搜索更多相關資源

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

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