entropy_min.m

化学计量学中利用信息上进行分辨的程序！用到了pso优化算法！可以得到很真实的光谱！

function entropy=entropy_min(T)

load('X.mat');
[U,S,V]=svd(X);
V=V(:,1:2);
%%%intial T--->transformation%%%%%%%%
%[m,n]=size(V);
%T=randn(1,m);
x=T*V;
c=A*x'*inv(a*a');
%%%the derivative spectrum%%%%%%%%%%%%
 G1=orient(x);
 [m1,n1]=size(x);
 %%%%%%penalty constraint%%%%%%%%%%%%%
 for i=1:n1
     if x(i)<0
         F1=x*x';
     else 
         F1=0;
     end
 end
 [m2,n2]=size(c);
 for j=1:n2
     if c(j)<0
         F2=c*c';
     else
         F2=0;
 end
 
 lamda_1=0.001;lamda_2=0.01;
 if F1<lamda_1
     gam_a=0;
 elseif F1<lamda_2&F1>=lamda_1
         gam=10;
     else
         gam=10000;
     end
 end
 gam_c=1000;
 
alfa=5;
if max(x)>alaf
    gam_max=10000;
else 
    gam_max=0;
end

 %%%%the optimum function%%%%%%%%%%%%%%%%
 G=G1+gam_a*F1+gam_c*F2+gam_max;