matlab宽带doa,请教 宽带DOA经典算法CSM

最近作宽带，用经典的CSM算法测向，其中聚焦矩阵是用信号在空间法得到的 跑出来的结果有些问题 请各位帮忙指点一下阿！十分感谢！

clc;

clear;

r=0.5;                % space of array

sensor=8;             % number of sensor

s_s=[r*(0:(sensor-1)/sensor);zeros(1,sensor)];

c=340;                % velocity of sound in air

t=10;                % 观测时间

fl=200;fh=300;        % 频带范围

f0=(fl+fh)/2;         % 中心频率

fs=1000;              % 采样频率

snap=fs*t;            % 采样点数

K=100;                 % 频域块拍

snap_s=snap/K;        % 频域点数

snaps=1000;           % fft点数

%----随机高斯信号通过带通滤波器--------

snr=10^(20/10);        % 信噪比

doa_azimuth=[20 0]*pi/180;                       % DOA

nu=length(doa_azimuth);                           % 信号个数

p=[cos(doa_azimuth);sin(doa_azimuth)];