gde.m

高频地波雷达中双排线形天线阵列的盖氏圆盘法估计信号源数

function m=GDE(t,k)

rsub=t(1:(k-1),1:(k-1));

% [v,d]=eig(rsub); 
% diad=abs(diag(d));
% [diad,indd]=sort(diad);
% diad=flipud(diad);
% indd=flipud(indd);
% v=v(:,indd);%7阶的特征向量按降序排列构成的酉矩阵
[u,d,v]=svd(rsub);
diad=diag(d);

dat=v;
dat(k,k)=1;%dat是酉变换矩阵T

rt=dat'*t*dat;
getrt=abs(rt(1:k-1,k));%getrt是前7行的半径
 
getradii=zeros(k-1,1);%求ri一撇的值
for i=1:(k-1)
    getradii(i,1)=getrt(i)/(sqrt(diad(i))*sqrt(abs(rt(k,k))));
end 
 
%%%%%%%%%%%%% %NEDGE的判别式
dm=1;
m=0;
sum1=0;
for l=1:k-1
    sum2=0;
    for i=l:k-1
        sum2=sum2+getradii(i);
    end
    sum1=sum1+sum2; 
end
for n=1:k-2
    sum3=0;  
    for l=1:k-1
        sum3=sum3+getradii(n);
    end
    sum4=0;
    sum4=sum1-sum3;
    edge(n)=getradii(n)-sum4*dm/(k-1);
    if edge(n)<0
        m=n-1;
        break
    end
end