function [Xt,t,MahalaD]=Mahalanob(X,e)
% X(m x n)為光譜矩陣(row is samples,column is spectrum)
% e(設定設置閾值范圍)馬氏距離標準差的權重系數
[m,n]=size(X);
Xm=mean(X);
CenterX=X-ones(m,1)*Xm;
M=CenterX'*CenterX./(m-1);
Minv=inv(M);
% (計算校正集樣品的馬氏距離)
for i=1:m
    MahalaD(i)=(X(i,:)-Xm)*Minv*(X(i,:)-Xm)';
end
MahalaDm=mean(MahalaD);
Dstd=std(MahalaD);
p=1;
q=1;
for i=1:m
    if MahalaD(i)<(MahalaDm+e*Dstd)
        Xt(p,:)=X(i,:);
        p=p+1;
    else
        t(q)=i;
        q=q+1;
    end
end