btypegreyincidence.m

基于灰色系统方法的离散多目标优化设计程序

function f2=BtypeGreyIncidence(f0,f1)  
%%%%%%%%%%%%%B TYPE degree of grey incidence%%%%%%%%%%%%%%%%%%%%%%
global GreyEntropy;
    [n,n1]=size(f0);
    d0=0;
    for i=1:n
       d0=d0+abs(f0(i,1)-f1(i,1)); 
    end
    d1=0;
    for i=1:n-1
       d1=d1+abs(f1(i+1,1)-f0(i+1,1)-f1(i,1)+f0(i,1)); 
    end
    d2=0;
    for i=2:n-1
       d2=d2+0.5*abs((f1(i+1,1)-f0(i+1,1))-2*(f1(i,1)+f0(i,1))+(f1(i-1,1)-f0(i-1,1))); 
    end
    %%%%%%%%%%%%%
    f2=1/(1+d0/n+d1/(n-1)+d2/(n-2));
    %%%%%%%%%%%%%%%%%%灰熵计算(改进并取负）%%%%%%%%%%%%%%%5
    ff2=(2+f2)*log(2+f2);
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 %%%%%%%%%%%%%%%%%%优度计算  %%%%%%%%%%%%%%%%%%%%%
    ff3=1/(1+((1-f2)/f2))^2;
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if GreyEntropy==1;
    f2=ff2;% 使用灰熵计算;
elseif GreyEntropy==2;
    f2=ff3;
else
    f2=f2;%灰关联计算。
end
 %%%%%%2007.8.18%%%%%%%%%%%%%%%%%%%%%