clf;
clear?all;
N=1024;%快拍數
N1=6;%Rx矩陣大小
j=sqrt（-1）;
n=0:N-1;
k=0.5*2*pi;%chirp?信號調頻斜率?k
M=4;?%陣元數
q=1;%信號源數目
p=1;
sita=30*pi/180;%角度->弧度
d=900;????????%間距
c=1500;???????%速度
f0=2000*2*pi;?%信號中心頻率）
?

?for?m=1:M
tao（m）=（m-1）*d*sin（sita）/c;%tao的表達式
?end
?%for?n=1:1:N?
?????s1=exp（j*（（1/2）*k*[0:N]^2+f0*[0:N]））;
????x1=awgn（s140‘measured‘）;
????for?m=1:M
????s2（m）=exp（-j*（k*[0:N]+f0）*tao（m）+j*pi*k*tao（m）.^2）*exp（j*（（1/2）*k*[0:N].^2+f0*[0:N]））;
????x2（m）=awgn（s240‘measured‘）;
????s=conj（x2）.*x1;%共軛乘積法構造正弦信號
????end
%?for?m=1:M-1
??%