clear;??
close?all;
clc;
Sensor=8;??????????????????????????????????%陣元數(shù)
bw1=1e6;????????????????????????????????????%?帶寬
bw2=1e6;?
bw3=1e6;?
f0=1e8;?
T1=1e-2;????????????????????????????????????%時(shí)間
T2=1e-4;????
T3=1e-6;????
L=512;??????????????????????????????????????%快拍數(shù)
dt1=T1/L;?????????????????????????????????%時(shí)間步長
dt2=T2/L;??
dt3=T3/L;??
gc?=?3e+8;?????????????????????????????????%速度
snr?=?2;???????????????????????????????????%信噪比
degrad=pi/180;?
p=3;???????????????????????????????????????%信源數(shù)
angle=10;????????????
angle2=20;??????????????????????????????????%DOA方向
angle3=-40;???
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?

t1=（0:L-1）*dt1;???????????????????????????%時(shí)間
t2=（0:L-1）*dt2;???
t3=（0:L-1）*dt3;???
sig1=exp（j*pi*bw1*t1.*t1/T1+j*2*pi*f0*t1）;???????????????????%信號(hào)
sig2=exp（j*pi*bw2*t2.*t2/T2+j*2*pi*f0*t2）;???
sig3=exp（j*pi*bw3*t3.*t3/T3+j*2*pi*f0*t3）;???
sig1=awgn（sig1snr）;?
sig2=awgn（sig2snr）;?
sig3=awgn（sig3snr）;?
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?
?
F=bw1/（T1）*t1;????????????????????????????
f1=F（64）;?
f2=F（2*64）;?
f3=F（3*64）;?
f4=F（4*64）;?
f5=F（5*64）;?
f6=F（6*64）;?
f7=F（7*64）;?
f8=F（8*64）;?
F=[f1f2f3f4f5f6f7f8];?
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?
?
?
%???????????????
Lambda=gc/f8;???????????????????????????????%?波長
d=0.5*Lambda;???????????????????????????????%陣元間距
tmp2=d*[0:Sensor-1?]‘;?????????????
?
v11=sig1（1:64）;?
v12=sig1（65:128）;?
v13=sig1（129:192）;?
v14=sig1（193:256）;?
v15=sig1（257:320）;?
v16=sig1（321:384）;?
v17=sig1（385:448）;?
v18=sig1（449:512）;?
v21=sig2（1:64）;?
v22=sig2（65:128）;?
v23=sig2（129:192）;?
v24=sig2（193:256）;?
v25=sig2（257:320）;?
v26=sig2（321:384）;?
v27=sig2（385:448）;?
v28=sig2（449:512）;?
v31=sig3（1:64）;?
v32=sig3（65:128）;?
v33=sig3（129:192）;?
v34=sig3（193:256）;?
v35=sig3（257:320）;?
v36=sig3（321:384）;?
v37=sig3（385:448）;?
v38=sig3（449:512）;?
?
s1=[（v11）;（v21）;v31];?
s2=[（v12）;（v22）;v32];?
s3=[（v13）;（v23）;v33];?
s4=[（v14）;（v24）;v34];?
s5=[（v15）;（v25）;v35];?
s6=[（v16）;（v26）;v36];?
s7=[（v17）;（v27）;v37];?
s8=[（v18）;（v28）;v38];?
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?
?
%導(dǎo)向矢量
A1=[exp（-j*2*pi*f1/gc*tmp2*sin（angle*degrad））exp（-j*2*pi*f1/gc*tmp2*sin（angle2*degrad））exp（-j*2*pi*f1/gc*tmp2*sin（angle3*degrad））];??
A2=[exp（-j*2*pi*f2/gc*tmp2*sin（angle*degrad））exp（-j*2*pi*f2/gc*tmp2*sin（angle2*degrad））exp（-j*2*pi*f2/gc*tmp2*sin（angle3*degrad））];??
A3=[exp（-j*2*pi*f3/gc*tmp2*sin（angle*degrad））exp（-j*2*pi*f3/gc*tmp2*sin（angle2*degrad））exp（-j*2*pi*f3/gc*tmp2*sin（angle3*degrad））];??
A4=[exp（-j*2*pi*f4/gc*tmp2*sin（angle*degrad））exp（-j*2*pi*f4/gc*tmp2*sin（angle2*degrad））exp（-j*2*pi*f4/gc*tmp2*sin（angle3*degrad））];??
A5=[exp（-j*2*pi*f5/gc*tmp2*sin（angle*degrad））exp（-j*2*pi*f5/gc*tmp2*sin（angle2*degrad））exp（-j*2*pi*f5/gc*tmp2*sin（angle3*degrad））];??
A6=[exp（-j*2*pi*f6/gc*tmp2*sin（angle*degrad））exp（-j*2*pi*f6/gc*tmp2*sin（angle2*degrad））exp