[音源位置] MATLABアレイマニホールドマトリックスに基づく信号表示[Matlabソースコード549を含む]

1.はじめに

1アレイマニホールドAには、アレイの幾何学的構造、アレイ要素のモード、アレイ要素間の結合、周波数、およびその他の影響が含まれます。
解像度：アレイの方向探知では、特定の方向のソースの解像度は、その方向に近いアレイの方向ベクトルの変化率に直接関係します。方向ベクトルが急激に変化する方向付近では、光源角度の変化に伴いアレイスナップショットデータが大きく変化し、対応する解像度も高くなります。
2均一な線形配列の場合：
D（θ）-cosθは
、信号の解像度が0方向で最も高く、60方向の解像度が半分になっていることを意味します。したがって、一般的な線形の方向探知範囲は配列は-60-60です。

3アレイの共分散行列：1。M
> K、つまり、アレイ要素の数Mは、アレイシステムが受信する可能性のある共周波数空間信号の数よりも多くなります。
2.異なる信号方向θの場合、信号方向ベクトルa（θ）は線形独立です。
3.アレイ内のノイズは、ガウス分布特性を持つランダムプロセスです。
4.空間ソース信号ベクトルS（t）の共分散行列Rsは対角非特異行列であり、空間ソース信号がインコヒーレントであることを示しています。また、空間ソース信号が出力とインコヒーレントである必要があります。配列要素。

第二に、ソースコード

clear all;
close all;
clc
%****************************************************%
M=2;                                %阵元数目
theta0=30;                          %声源到达方向
p=1;                                 %声源个数
lamda=1.6;                        %波长
d=lamda/2;                       %阵元间距
a=[0:1:M-1];
a_theta0=exp(j*pi*2*d/lamda*sin(theta0*pi/180)*a);%第一阵列流型矢量
%****************************************************%
[s,fs]=wavread('C9_3_y.wav');   
s=s/max(abs(s));
sound_length = 6400;
%装载房间冲激响应
load h.mat;            
%麦克风1的信号
s1=conv(s,h1);
s1=s1(1:sound_length);
%麦克风2的信号
s2=conv(s,h2);
s2=s2(1:sound_length);
%****************************************************%
s1=s1*a_theta0(1);                  %叠加角度信息
s2=s2*a_theta0(2);
%****************************************************%
%添加白噪声
%SNR=20dB
s1_20db = awgn(s1,20,'measured','db');
s2_20db = awgn(s2,20,'measured','db');
%SNR=10dB
s1_10db = awgn(s1,10,'measured','db');
s2_10db = awgn(s2,10,'measured','db');
%SNR=0dB
s1_0db = awgn(s1,0,'measured','db');
s2_0db = awgn(s2,0,'measured','db');
%SNR=-5dB
s1_m5db = awgn(s1,-5,'measured','db');
s2_m5db = awgn(s2,-5,'measured','db');
figure(1);
%显示原始输入语音信号波形
subplot(5,2,1),plot(s1,'k'),title('麦克风1（原始）'),xlabel('采样点'),ylabel('幅度');
subplot(5,2,2),plot(abs(s2),'k'),title('麦克风2 （原始）'),xlabel('采样点'),ylabel('幅度');
subplot(5,2,3),plot(s1_20db,'k'),title('麦克风1（SNR=20dB）'),xlabel('采样点'),ylabel('幅度');
subplot(5,2,4),plot(abs(s2_20db),'k'),title('麦克风2（SNR=20dB）'),xlabel('采样点'),ylabel('幅度');
subplot(5,2,5),plot(s1_10db,'k'),title('麦克风1（SNR=10dB）'),xlabel('采样点'),ylabel('幅度');
subplot(5,2,6),plot(abs(s2_10db),'k'),title('麦克风2（SNR=10dB）'),xlabel('采样点'),ylabel('幅度');
subplot(5,2,7),plot(s1_0db,'k'),title('麦克风1（SNR=0dB）'),xlabel('采样点'),ylabel('幅度');

3、実行中の結果

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

四、備考