[音源定位] matlab一般化相互相関に基づく音源定位[Matlabソースコード548を含む]

1.はじめに

最近、スペクトル監視システムのTDOA測位精度が最適化されています。TDOAで最も重要な問題は、信号の時間差の抽出にすぎません。国内外の文献を読んで、彼らは満場一致で「一般化相互相関（GCC）アルゴリズム」を推奨しています。
文献を確認し、MATLABでシミュレートされた正弦波シミュレーションを使用して一般化された相関アルゴリズムをテストしましたが、期待した効果が得られず、戸惑いました。インターネットでGCCのコメントを閲覧したところ、国内外の多くの人が期待したシミュレーション結果が得られないと不満を漏らしました。最後に、一般化された相関アルゴリズムの式では、計算結果に影響を与える主なパラメータは位相情報であり、シミュレートされた変調信号でテストする必要があると思います。Matlabを使用してFM / FSK / PSKをシミュレートした後、効果は非常に良好です。このアルゴリズムをC ++プログラムに実装するのは、作業負荷の問題です。参考までに、Matlabシミュレーションのソースコードをここに示します。Matlabのバージョンが異なれば、異なる変調関数を使用する必要があることに注意してください。
Generalized Cross-Correlationを使用して、信号が2つのセンサーに到達するまでの時間遅延を推定できます。たとえば、マイクアレイを使用して音源の位置を特定するアプリケーションの背景。この問題は、次のモデルで説明できます。
ここに画像の説明を挿入します

第二に、ソースコード

clear all;
close all;
clc
load s.mat;
s1=s1_10db(:)';             %变为行向量
s2=s2_10db(:)';             %变为行向量
wnd=512;
inc=256;

[delay]=GCC_Method('standard',s1,s2,wnd,inc);
subplot(411)
plot(delay-wnd,'*')
ylim([0,12])
title('标准GCC')
xlabel('帧数')
ylabel('延时/点')

[delay]=GCC_Method('phat',s1,s2,wnd,inc);
subplot(412)
plot(delay-wnd,'*')
ylim([0,12])
title('Phat-GCC')
xlabel('帧数')
ylabel('延时/点')

[delay]=GCC_Method('scot',s1,s2,wnd,inc);
subplot(413)
plot(delay-wnd,'*')
ylim([0,12])
title('Scot-GCC')
xlabel('帧数')
ylabel('延时/点')

[delay]=GCC_Method('ml',s1,s2,wnd,inc);
 subplot(414)
plot(delay-wnd,'*')
ylim([0,12])
unction frameout=enframe(x,win,inc)

nx=length(x(:));            % 取数据长度
nwin=length(win);           % 取窗长
if (nwin == 1)              % 判断窗长是否为1，若为1，即表示没有设窗函数
   len = win;               % 是，帧长=win
else
   len = nwin;              % 否，帧长=窗长
end
if (nargin < 3)             % 如果只有两个参数，设帧inc=帧长
   inc = len;
end
nf = fix((nx-len+inc)/inc); % 计算帧数
frameout=zeros(nf,len);            % 初始化
indf= inc*(0:(nf-1)).';     % 设置每帧在x中的位移量位置
inds = (1:len);             % 每帧数据对应1:len
frameout(:) = x(indf(:,ones(1,len))+inds(ones(nf,1),:));   % 对数据分帧

3、実行中の結果

ここに画像の説明を挿入します

四、備考