clear;
clc;
Sensor=8;?????????????????????????????????????%天線陣元數(shù)
bw1=1e7;??????????????????????????????????????%?信號(hào)帶寬
bw2=1e6;
bw3=51e6;
f1=1e8;
f2=1e7;
f3=1e7;
T1=1e-2;??????????????????????????????????????%信號(hào)脈沖寬度??
L=1024;???????????????????????????????????????%采樣點(diǎn)數(shù)
dt1=T1/L;?????????????????????????????????????%采樣間隔
gc=3e8;???????????????????????????????????????%光速
snr=-500;???????????????????????????????????????%信躁比
degrad=pi/180;
p=3;??????????????????????????????????????????%信號(hào)源個(gè)數(shù)?
angle=10;???????????
angle2=20;????????????????????????????????????%信號(hào)源角度
angle3=-40;??
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%構(gòu)造采樣序列
t1=（0:L-1）*dt1;???????????????????????????????%時(shí)域采樣點(diǎn)
sig1=exp（j*pi*bw1*t1.*t1/T1+j*2*pi*f1*t1）;????%產(chǎn)生寬帶信號(hào)
sig2=exp（j*pi*bw2*t1.*t1/T1+j*2*pi*f2*t1）;??
sig3=exp（j*pi*bw3*t1.*t1/T1+j*2*pi*f3*t1）;??
sig1=awgn（sig1snr）;
sig2=awgn（sig2snr）;??????????????????????????%加白噪音
sig3=awgn（sig3snr）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%寬帶信號(hào)頻率
F=f1+bw1/（T1）*t1;???????????????????????????
f1=F（128）;
f2=F（2*128）;
f3=F（3*128）;
f4=F（4*128）;
f5=F（5*128）;
f6=F（6*128）;
f7=F（7*128）;
f8=F（8*128）;
F=[f1f2f3f4f5f6f7f8];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%構(gòu)造采樣數(shù)據(jù)矩陣
Lambda=gc/f8;?????????????????????????????????%波長(zhǎng)
d=0.5*Lambda;?????????????????????????????????%陣元間距離
tmp2=d*[0:Sensor-1?]‘;????????????

v11=sig1（1:128）;
v12=sig1（129:256）;
v13=sig1（257:384）;
v14=sig1（385:512）;
v15=sig1（513:640）;
v16=sig1（641:768）;
v17=sig1（769:896）;
v18=sig1（897:1024）;
v21=sig2（1:128）;
v22=sig2（129:256）;
v23=sig2（257:384）;
v24=sig2（385:512）;
v25=sig2（513:640）;
v26=sig2（641:768）;
v27=sig2（769:896）;
v28=sig2（897:1024）;
v31=sig3（1:128）;
v32=sig3（129:256）;
v33=sig3（257:384）;
v34=sig3（385:512）;
v35=sig3（513:640）;
v36=sig3（641:768）;
v37=sig3（769:896）;
v38=sig3（897:1024）;

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））];?
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（-j*2*pi*f6/gc*tmp2*sin（angle3*degrad））];?
A7=[exp（-j*2*pi*f7/g