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;?????????????????%??協方差