clear?all;close?all;clc;



snapshot_number=500;?????????????%快拍數(shù)
Ts=0.001;????????????????????????%采樣周期
SNR=40;??????????????????????????%信噪比


snr=15;%信噪比
tangle1=10;??%三個(gè)信號(hào)的發(fā)射角度
tangle2=35;?%三個(gè)信號(hào)的入射角度
tangle3=60;

c1=kron（exp（-1i*pi*（0:9）*sin（tangle1*pi/180））‘exp（-1i*pi*（0:9）*sin（tangle1*pi/180））‘）;??%3個(gè)信號(hào)導(dǎo)向矢量?????
c2=kron（exp（-1i*pi*（0:9）*sin（tangle2*pi/180））‘exp（-1i*pi*（0:9）*sin（tangle2*pi/180））‘）;???????
c3=kron（exp（-1i*pi*（0:9）*sin（tangle2*pi/180））‘exp（-1i*pi*（0:9）*sin（tangle2*pi/180））‘）;

A=[c1c2?c3];?
for?r=1:50

s=5*（2*round（rand（31*snapshot_number））-1）;???%仿真信號(hào)
z1=A*s;
sigPower=mean（sum（abs（z1）.^2）/length（c1））;
noisePower=sigPower*10.^（-SNR/10）;
x=z1+randn（length（c1）snapshot_number）*sqrt（noisePower）;