ISM算法的matlab代码

%非相干信号子空间（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