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