宽带信号的波达方向估计.m

clear;  
close all;
clc;
Sensor=8;                                  %阵元数
bw1=1e6;                                    % 带宽
bw2=1e6; 
bw3=1e6; 
f0=1e8; 
T1=1e-2;                                    %时间
T2=1e-4;    
T3=1e-6;    
L=512;                                      %快拍数
dt1=T1/L;                                 %时间步长
dt2=T2/L;  
dt3=T3/L;  
gc = 3e+8;                                 %速度
snr = 2;                                   %信噪比
degrad=pi/180; 
p=3;                                       %信源数
angle=10;            
angle2=20;                                  %DOA方向
angle3=-40;   
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

t1=（0:L-1）*dt1;                           %时间
t2=（0:L-1）*dt2;   
t3=（0:L-1）*dt3;   
sig1=exp（j*pi*bw1*t1.*t1/T1+j*2*pi*f0*t1）;                   %信号
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]; 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
 
%导向矢量
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