% ---------------------------------------------------------------------------------------------------
% 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);