宽带信号的DOA估计.m

clf;
clear all;
N=1024;%快拍数
N1=6;%Rx矩阵大小
j=sqrt（-1）;
n=0:N-1;
k=0.5*2*pi;%chirp 信号调频斜率 k
M=4; %阵元数
q=1;%信号源数目
p=1;
sita=30*pi/180;%角度->弧度
d=900;        %间距
c=1500;       %速度
f0=2000*2*pi; %信号中心频率）
 

 for m=1:M
tao（m）=（m-1）*d*sin（sita）/c;%tao的表达式
 end
 %for n=1:1:N 
     s1=exp（j*（（1/2）*k*[0:N]^2+f0*[0:N]））;
    x1=awgn（s140‘measured‘）;
    for m=1:M
    s2（m）=exp（-j*（k*[0:N]+f0）*tao（m）+j*pi*k*tao（m）.^2）*exp（j*（（1/2）*k*[0:N].^2+f0*[0:N]））;
    x2（m）=awgn（s240‘measured‘）;
    s=conj（x2）.*x1;%共轭乘积法构造正弦信号
    end
% for m=1:M-1
  %