%初始條件
clear?all;?
Sensor=8;??????????????????????????????????%天線陣元數(shù)
bw1=1e7;????????????????????????????????????%?信號帶寬
bw2=1e7;
bw3=1e7;
f0=1e8;
T1=1e-2;????????????????????????????????????%信號脈沖寬度
T2=1e-4;???
T3=1e-6;???
L=512;??????????????????????????????????%采樣點數(shù)
dt1=T1/L;?????????????????????????????????%采樣間隔
dt2=T2/L;?
dt3=T3/L;?
gc?=?3e+8;?????????????????????????????????%光速
snr?=?2;???????????????????????????????????%信躁比
degrad=pi/180;
p=3;???????????????????????????????????????%?信號源個數(shù)?
angle=10;???????????
angle2=20;??????????????????????????????????%信號源角度
angle3=-40;??
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%構(gòu)造采樣序列
t1=（0:L-1）*dt1;???????????????????????????%時域采樣點
t2=（0:L-1）*dt2;??
t3=（0:L-1）*dt3;??

%?sig1=exp（j*pi*bw1*t1.*t1/T1+2*pi*f0*t1）;???????????????????%產(chǎn)生線性調(diào)頻信號
%?sig2=exp（j*pi*bw2*t2.*t2/T2+2*pi*f0*t2）;??
%?sig3=exp（j*pi*bw3*t3.*t3/T3+2*pi*f0*t3）;??
%上面三行代碼是原程序的，有點小問題

sig1=exp（j*（pi*bw1*t1.*t1/T1+2*pi*f0*t1））;???????????????????%產(chǎn)生線性調(diào)頻信號
sig2=exp（j*（pi*bw2*t2.*t2/T2+2*pi*f0*t2））;??
sig3=exp（j*（pi*bw3*t3.*t3/T3+2*pi*f0*t3））;??
sig1=awgn（sig1snr）;
sig2=awgn（sig2snr）;
sig3=awgn（sig3snr）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%線性調(diào)頻信號頻率
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];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%構(gòu)造采樣數(shù)據(jù)矩陣
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];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%構(gòu)造方向矩陣
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））];