matlab运行结果显示不出,程序运行结果不对，自己查不出来错，请各位大侠帮忙...

下面是写的一段UCA ROOT-MUSIC程序，但是运行结果不对，不是设置的[20,50]。同样是本程序的噪声子空间，用谱峰搜索的MUSIC就是对的，但是用ROOT-MUSIC就不对。应该就是root-music那段程序有错，但是自己查也查不出来错误，请各位大侠帮帮忙，非常感谢！

sig_num = 1000;%%信号的快拍数

snr=10;

rou=0;

x1=randint(sig_num,1,4);%产生sig_sum*1的矩阵，数值在0-3之间

x2=randint(sig_num,1,4);

x3=randint(sig_num,1,4);

x222=rou*x1+sqrt(1-(abs(rou))^2)*x3;%产生相关的两路信号

x2=round(x222);

for i=1:sig_num

if x2(i)>=4

x2(i)=3;

elseif x2(i)<=0

x2(i)=0;

end

end

y1=dmod(x1,1,0.1,4,'psk',4);

y1=y1(1:sig_num);

y2=dmod(x2,1,0.1,4,'psk',4);

y2=y2(1:sig_num);

x=[y1';y2'];

s=x;

q=2;

p=12;%智能天线阵元数

fc=10.7*10^6;

<