互相關(guān)延時估計加權(quán)函數(shù)性能分析(共12頁)_第1頁
互相關(guān)延時估計加權(quán)函數(shù)性能分析(共12頁)_第2頁
互相關(guān)延時估計加權(quán)函數(shù)性能分析(共12頁)_第3頁
互相關(guān)延時估計加權(quán)函數(shù)性能分析(共12頁)_第4頁
互相關(guān)延時估計加權(quán)函數(shù)性能分析(共12頁)_第5頁
已閱讀5頁,還剩10頁未讀, 繼續(xù)免費閱讀

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認領(lǐng)

文檔簡介

1、互相關(guān)延時估計加權(quán)函數(shù)(hnsh)性能分析廣義互相關(guān)函數(shù)法是通過首先求出倆信號之間的互功率譜,然后在頻域內(nèi)給予一定的加權(quán),以此對信號和噪音進行白化處理,從而增強信號中信噪比較高的頻率成分,抑制(yzh)噪聲的影響,最后再反變換到時域,得到兩信號之間的互相關(guān)函數(shù),即: (1)其中(qzhng)是廣義互相關(guān)加權(quán)函數(shù)。廣義互相關(guān)加權(quán)函數(shù)的選擇主要基于倆個方面:噪聲和反射情況。根據(jù)不同的情況選擇加權(quán)函數(shù),其目的就是使具有比較尖銳的峰值。峰值處就是倆個傳感器之間的時延。由于來自同一聲源的信號存在一定的相關(guān)性,通過計算不同麥克風(fēng)所接受到的信號之間的相關(guān)函數(shù),就可以估計出TDOA值。然而在實際環(huán)境中,由于噪

2、聲和混響的影響,相關(guān)函數(shù)的最大峰會被弱化,有時還會出現(xiàn)多個峰值,這些都造成了實際峰值的檢測困難。此時就通過加權(quán)的方法來銳化峰值,通常我們通過時間、精度來確定算法的合理性。廣義互相關(guān)函數(shù)模擬clear all; clc; close all;N=1024;%長度Fs=500;%采樣頻率n=0:N-1;t=n/Fs; %時間序列a1=5; %信號幅度a2=5;d=2; %延遲點數(shù)x1=a1*cos(2*pi*10*n/Fs); %信號1x1=x1+randn(size(x1); %加噪聲x2=a2*cos(2*pi*10*(n+d)/Fs); %信號2x2=x2+randn(size(x2);su

3、bplot(211);plot(t,x1,r);axis(-0.2 1.5 -6 6);hold on;plot(t,x2,:);axis(-0.2 1.5 -6 6);legend(x1信號, x2信號);xlabel(時間/s);ylabel(x1(t) x2(t);title(原始信號);grid on;hold off%互相關(guān)函數(shù)X1=fft(x1,2*N-1);X2=fft(x2,2*N-1);Sxy=X1.*conj(X2);Cxy=fftshift(ifft(Sxy);%Cxy=fftshift(real(ifft(Sxy);subplot(212);t1=(0:2*N-2)/F

4、s; %注意plot(t1,Cxy,b);title(互相關(guān)函數(shù));xlabel(時間/s);ylabel(Rx1x2(t);grid onmax,location=max(Cxy);%求出最大值max,及最大值所在的位置(第幾行)location;%d=location-N/2-1 %算出延遲了幾個點d=location-NDelay=d/Fs %求得時間延遲運行程序得到(d do)的結(jié)果是:d = 2Delay = 0.0040可以看出,通過(tnggu)互相關(guān)函數(shù)的求解d=2,delay=0.0040,這和我們給出的信號的時延d/Fs=0.0040是一致的。這表明互相關(guān)函數(shù)可以給出信號的

5、時延估計。PHAT-GCC模擬(mn)clear all; clc; close all;N=1024; %長度(chngd)Fs=500; %采樣(ci yn)頻率n=0:N-1;t=n/Fs; %時間(shjin)序列a1=5; %信號幅度a2=5;d=9; %延遲點數(shù)x1=a1*cos(2*pi*10*n/Fs); %信號1x1=x1+randn(size(x1); %加噪聲%x1=x1.*hamming(max(size(x1);%加窗x2=a2*cos(2*pi*10*(n+d)/Fs); %信號2x2=x2+randn(size(x2);%x2=x2.*hamming(max(si

6、ze(x2);%加窗subplot(211);plot(t,x1,r);axis(-0.2 2 -6 6);hold on;plot(t,x2,:);axis(-0.2 2 -6 6);legend(x1信號, x2信號);xlabel(時間/s);ylabel(x1(t) x2(t);title(原始信號);grid on;hold off%互相關(guān)函數(shù)X1=fft(x1,2*N-1);X2=fft(x2,2*N-1);Sxy=X1.*conj(X2);%Cxy=fftshift(ifft(Sxy);Cxy=fftshift(ifft(Sxy./abs(Sxy);subplot(212);t1

7、=(-N+1:N-1)/Fs;plot(t1,Cxy,b);title(Rx1x2);xlabel(t/s);ylabel(Rx1x2(t);grid onmax,location=max(Cxy);%d=location-N/2-1 d=location-NDelay=d/Fs %求得時間延遲運行程序得到的結(jié)果是:d = 1 Delay = 0.0020我們可以看見結(jié)果是d=1,delay=0.0020,而實例中給出的時延為d/fs=0.016,這并不表示PHAT-GCC算法是錯誤的,只是(zhsh)因為,我們在信號中加入了均值為0,方差為1的高斯白噪音,所以才會導(dǎo)致了誤差的存在。ROTH-

8、GCC模擬(mn)clear;N=1024;%信號(xnho)長度fs=500;%采樣頻率n=0:N-1;t=n/fs;%時間序列a1=5;%信號幅度a2=5;%信號幅度d=2;%延遲點數(shù)x1=a1*sin(2*pi*10*n/fs);x2=a2*sin(2*pi*10*(n+d)/fs);%x2=awgn(x1./4,-3); %噪聲強度大于信號%x2=x2 .* hamming(N);x1=x1+randn(size(x1); %加入噪聲x2=x2+randn(size(x2);S1=fft(x1,2*N-1);S2=fft(x2,2*N-1);S12 = S1.* conj(S2);S1

9、1 = S1.* conj(S1);R1 =real(fftshift(ifft(S12./abs(S11);ts=(-N+1:N-1)/fs;plot(ts,R1);xlabel(時間(shjin)/s);ylabel(R1(t);title(互相(h xing)關(guān)函數(shù));max,location=max(R1);%d=location-N/2-1 d=location-NDelay=d/fs 運行(ynxng)程序的結(jié)果為:d = 4Delay =0.0080SCOT-GCC模擬clear;N=1024;%信號(xnho)長度fs=1000;%采樣(ci yn)頻率n=0:N-1;t=n/

10、fs;%時間(shjin)序列ts = 1/fs * (-N + (1 : 2*N - 1); %互相關(guān)時間序列a1=5;%信號幅度a2=5;%信號幅度d=26;%延遲點數(shù)x1=a1*sin(2*pi*10*t)+1.9*sin(2*pi*18*t)+2.8*sin(2*pi*55*t);x2=a2*sin(2*pi*10*(n+d)/fs)+1.9*sin(2*pi*18*(n+d)/fs)+2.8*sin(2*pi*55*(n+d)/fs);%x2=awgn(x1./4,-3); %噪聲強度大于信號%x2=x2 .* hamming(N);x=awgn(x1,20); %加入噪聲y=awg

11、n(x2,0.001);S1=fft(x,2*N-1);S2=fft(y,2*N-1);X = S1.* conj(S2);X11 = S1.* conj(S1);X22 = S2.* conj(S2);Y=sqrt(X11.*X22);R1 =real(fftshift(ifft(X./Y);plot(ts,R1);xlabel(時間/s);ylabel(R1(t);title(ifft計算結(jié)果)max,location=max(R1);%d=location-N/2-1 d=location-NDelay=d/fs 運行程序的結(jié)果是:d = 8Delay =0.0080相同信噪比不同(b

12、tn)算法的比較clear all; clc; close all;N=1024; %長度(chngd)Fs=500; %采樣(ci yn)頻率n=0:N-1;t=n/Fs; %時間序列a1=30; %信號幅度a2=30;d=9; %延遲點數(shù)x1=a1*cos(2*pi*10*n/Fs); %信號1x1=awgn(x1,20); %加噪聲%x1=x1.*hamming(max(size(x1);%加窗x2=a2*cos(2*pi*10*(n+d)/Fs); %信號2x2=awgn(x2,20);%x2=x2.*hamming(max(size(x2);%加窗subplot(511);plot(

13、t,x1,r);axis(-0.2 2 -40 40);hold on;plot(t,x2,:);axis(-0.2 2 -40 40);legend(x1信號, x2信號);xlabel(時間(shjin)/s);ylabel(x1(t) x2(t);title(原始(yunsh)信號);grid on;%互相(h xing)關(guān)函數(shù)ticX1=fft(x1,2*N-1);X2=fft(x2,2*N-1);Sxy=X1.*conj(X2);%Cxy=fftshift(ifft(Sxy);%GCCCxy=fftshift(ifft(Sxy);subplot(512);t1=(-N+1:N-1)/

14、Fs;plot(t1,Cxy,b);title(GCC);xlabel(t/s);ylabel(Cxy);grid on;max1,location1=max(Cxy);%d=location-N/2-1 d1=location1-NDelay1=d1/Fs %求得時間延遲toc%phat-gccticX1=fft(x1,2*N-1);X2=fft(x2,2*N-1);Sxy=X1.*conj(X2);Pxy=fftshift(ifft(Sxy./abs(Sxy);subplot(513);t1=(-N+1:N-1)/Fs;plot(t1,Pxy,b);title(phat-gcc);xlab

15、el(t/s);ylabel(Pxy);grid on;max2,location2=max(Pxy);%d=location-N/2-1 d2=location2-NDelay2=d2/Fs %求得時間延遲toc%rhat-gccticX1=fft(x1,2*N-1);X2=fft(x2,2*N-1);Sxy=X1.*conj(X2);S11 = X1.* conj(X1);Rxy=fftshift(ifft(Sxy./abs(S11);subplot(514);t1=(-N+1:N-1)/Fs;plot(t1,Rxy,b);title(phat-gcc);xlabel(t/s);ylabe

16、l(Rxy);grid on;max3,location3=max(Rxy);%d=location-N/2-1 d3=location3-NDelay3=d3/Fs %求得時間延遲toc%scot-gccticX1=fft(x1,2*N-1);X2=fft(x2,2*N-1);Sxy=X1.*conj(X2);S11 = X1.* conj(X1);S22 = X2.* conj(X2);Y=sqrt(S11.*S22);SCxy=fftshift(ifft(Sxy./Y);subplot(515);t1=(-N+1:N-1)/Fs;plot(t1,SCxy,b);title(scot-gc

17、c);xlabel(t/s);ylabel(SCxy);grid on;max4,location4=max(SCxy);%d=location-N/2-1 d4=location4-NDelay4=d4/Fs %求得時間延遲tocSNR=0時SNR=10時:SNR=20時SNR=50時分析(fnx):從運行(ynxng)結(jié)果上來看,在時間上,基本(jbn)互相關(guān)、PHAT加權(quán)、ROTH加權(quán)和SCOT加權(quán)四種算法的運行時間基本相同;但是從峰度的銳化來說,這四種方式的時延估計的準確性隨著信噪比的降低而惡化,互相關(guān)函數(shù)峰值的尖銳程度隨信噪比的降低而降低。對于SCOT加權(quán)來說,隨著信噪比的降低,性能急劇下降?;净ハ嚓P(guān)函數(shù)和RHOT加權(quán)雖然有一定的抗噪能力,但隨著信噪比的降低,其波動程度明顯加強,特別是對外圍的噪聲、反射和有限觀測數(shù)據(jù)很敏感,會造成峰值不明顯;對于PHAT加權(quán),

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論