宽带信号DOA估计算法

%初始条件
clear all; 
Sensor=8;                                  %天线阵元数
bw1=1e7;                                    % 信号带宽
bw2=1e7;
bw3=1e7;
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;                                  %信号源角度
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+2*pi*f0*t1）;                   %产生线性调频信号
% 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））;                   %产生线性调频信号
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）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%线性调频信号频率
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））];