基于矢量水听器阵列的MUSIC算法.m

clear all
close all
derad = pi/180;        
radeg = 180/pi;
twpi = 2*pi;
kelm = 8;               % 阵列数量
dd = 0.5;               % 阵元间距 
d=0:dd:（kelm-1）*dd;     % 阵元位置
iwave = 3;              % 信源数
theta = [10 30 60];     % 波达方向
snr = 10;               %  信噪比）
n = 500;                %  采样点数
l=[111];
u=[l;cos（theta*derad）;sin（theta*derad）];   %  加矢量水听器



Ap1=exp（-j*twpi*d.‘*sin（theta（1）*derad））;%%%% 声压阵列流行
Ap2=exp（-j*twpi*d.‘*sin（theta（2）*derad））;
Ap3=exp（-j*twpi*d.‘*sin（theta（3）*derad））;


A1=kron（Ap1u（:1））;
A2=kron（Ap2u（:2））;
A3=kron（Ap3u（:3））; %  直积得出矢量水听器阵列流行
A=[A1A2A3];
S=randn（iwaven）;             %  构造随机信号
X=A*S;
X1=awgn（Xsnr‘measured‘）;    %  水听器接收
Rxx=X1*X1‘/n;                 %  协方差