clear?all;close?all;clc;
sensor_number=8;
source_number=2;
snapshort_number=250;
source_theta=[?-3050];
source_fei=[20-60];
snr=10;
A=zeros（sensor_numbersource_number）;
k=0:sensor_number-1;
?for?i=1:source_number
???A（:i）=exp（-j*pi*cos（2*pi*k‘/sensor_number-source_theta（i）*pi/180）*sin（source_fei（i）*pi/180））;
?end

length=500;
M1=2;
M2=4;
S1=randi（[0?M1-1]length1）;
S2=randi（[0?M2-1]length1）;
ps1=pskmod（S1M1）;
ps2=pskmod（S2M2）;
t=151:150+snapshort_number;
S=[ps1（t）?ps2（t）]‘;
s=A*S;
noise=randn（sensor_numbersnapshort_number）+j*randn（sensor_numbersnapshort_number）;
ss=s+noise;
%ss=awgn（ssnr）;

Rs=ss*ss‘/snaps