宽带非相干信号DOA估计

clc
clear all
close all
M=12;%阵元数
N=200;%快拍数
ts=0.01;%时域采样间隔
f0=100;%入射信号中心频率
f1=80;%入射信号最低频率
f2=120;%入射信号最高频率
c=1500;%声速
lambda=c/f0;%波长
d=lambda/2;%阵元间距
SNR=15;
b=pi/180;
theat1=30*b;%入射信号波束角1
theat2=20*b;%入射信号波束角2
n=ts:ts:N*ts;
theat=[theat1 theat2]‘;
s1=chirp（n801120）;%生成线性调频信号1；
sa=fft（s12048）;%进行fft变换
figure
specgram（s12561E3256250）%频谱图
s2=chirp（n+0.100801120）;%生成线性调频信号2；
sb=fft（s22048）;%进行fft变换
%%%%%%%%%%%%%ISM算法
P=1:2;
a=zeros（M2）;
sump=zeros（1181）;
for i=1:N
    f=80+（i-1）*1.0;
    s=[sa（i） sb（i）]‘;
    for m=1:M
         a（mP）=exp（-j*2*pi*f*d/c*si