关于一维MUSIC的仿真问题

请教各位，我做了个MUSIC的仿真，可是结果不对劲，没有估计出方向，不知道错在哪，程序和图如下。
% 一维music
clear;
close al
%----产生数据-------------
signal_a=sqrt(10^(snr/10));
t=[0:1:L-1)];
f1 = 0.21;S1=sin(2*pi*f1*t);
f2 = 0.33;S2=sin(2*pi*f2*t);
S=[S1;S2];
A=exp(i*2*pi*ddd*(0:M-1)'*sin(signal_doa*pi/180));
X=A*S*signal_a;
R=X*X'/L;
[U D]=eig(R);
US=U(:,1:N);
UN=U(:,N+1:M);
ddoa=-90:0.1:90;
steer=exp(i*2*pi*ddd*(0:M-1)'*sin(ddoa*pi/180));
for ij=1:length(ddoa)
aa=steer(:,ij);
p1(ij)=1/(aa'*UN*UN'*aa);
end
%----画图得到角度-------------
figure;box on;grid on;
plot(ddoa,20*log10(abs(p1)));
legend('MUSIC')