%非相干信號子空間（ISM）算法
clc
clear?all
close?all
L=8;???????????????????????????????????????%陣元數
K=33;????????????????????????????????????%快拍數
ts=0.0125;?????????????????????????????????%時域采樣間隔
f0=100;????????????????????????????????????%入射信號中心頻率?
f1=80;?????????????????????????????????????%入射信號最低頻率
f2=120;????????????????????????????????????%入射信號最高頻率
c=1500;????????????????????????????????????%聲速
lambda=c/f0;???????????????????????????????%波長
d=lambda/2;????????????????????????????????%陣元間距
b=pi/180;??
theat1=-20*b;???????????????????????????????%入射信號波束角1
theat2=20*b;????????????????????????????????%入射信號波束角2
n=ts:ts:K*ts;????????????????????????????
theat=[theat1?theat2]‘;

s1=chirp（n801120）;?%余弦調頻信號（t0f0t1f1），信號的頻率隨時間的增長而變化
sa=fft（s12048）;
%?figure（2）;
%???????specgram（s12561E3256250）;?%?Display?the?spectrogramsa=fft（s12048）;??????????????????????????????%進行FFT變換
s2=chirp（n+0.12580