dmatrix.m

宽带信号源测向的一种新方法,TOPS方法的源程序

% ---------------------------------------------------------------------------------------------------
% TOPS法子程序，由接受数据得出噪声子空间
%对多次快拍用music方法，形参Xt为矩阵，每一列为一次快拍
% ---------------------------------------------------------------------------------------------------

function [NoisespaceM,SignalspaceM]=Dmatrix(Xt,Nsensor,Nsignal)

CovM=CorGenM(Xt,Xt);
[EigvectorM,DiagM]=eig(CovM);
Temp1=diag(abs(DiagM));%diag是取对角元素为列矢量

[Temp2,in]=sort(Temp1);
Temp3=EigvectorM(:,in);%将特征向量按照in的顺序排列，即排列后前面的列向量对应小特征值，后面的对应大特征值，那样从前向后开始取为噪声空间就可以了。
NoisespaceM=Temp3(:,1:Nsensor-Nsignal);
SignalspaceM=Temp3(:,Nsensor-Nsignal+1:Nsensor);

% %NoisespaceM=Temp3(:,1)
% theta=-90:0.1:90;
% 
% for m=1:length(theta)
%    aa=steer(Nsensor,Radius,mid_f,theta(m),fre);
%   % s(m)=norm(a)/norm(NoisespaceM'*a);
%   s(m)=1/((NoisespaceM'*aa)'*(NoisespaceM'*aa));
% end
% 
% s=s/max(s);