matlab宽带doa,低频宽带信号的空间重采样实现doa估计

我想对宽带信号实现DOA估计，应用fft空间重采样方法聚焦，MUSIC算法估计，

问题： 只能识别-90°到0°的信号，不能识别0°到90°的信号！！！

跪求大神指点！！！！

%%

clear all;

close all;

clc;

c=340;                  %波速

M=16;                    %阵元数

M_M=[(0:(M-1))];

f=[80:120];            %信号频率

SNR=8;                  %信噪比

N=800;                  %信号长度

fs=800;                  %信号采样频率

l=c/f(end);            %波长

d=l/2;                  %阵元间距

%%

doa_azimuth=[-20 50]*pi/180;                         % D方向

p=[sin(doa_azimuth)];                               % 方向向量

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

[b,a]=butter(5,[f(1) f(end)]/fs*2);                 % 带通滤波器

x=awgn((randn(D,N)+j*