dengtypegreyincidence.m

基于粒子群优化与灰色系统理论的多目标优化程序

function f2=DengTypeGreyIncidence(f0,f1,factor)
% Deng type degree of  Grey Incidence 
%f0=[1000000;1000000;1000000;1000000;1000000;1000000;1000000];
%f1=[30622.0927; ;1477.74233;1214.19412;59811.2107;211.272102;604.756576;1485.3764];
%factor=[1/7;1/7;1/7;1/7;1/7;1/7;1/7];
global GreyEntropy;
[n,n1]=size(f0);
f00=f0/f0(1);
f10=f1/f1(1);
for i=1:n,
    da(i)=f00(i)-f10(i);
end;
da1=min(abs(da));da2=max(abs(da));kk=0.5;
for i=1:n, 
    fe(i)=(da1+kk*da2)/(abs(da(i))+kk*da2);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%
f2=0;
for i=1:n 
    f2=f2+fe(i)*factor(i);
end
%%%%%%%%%%%%%%%%%%熵权关联计算%%%%%%%%%%%%%%%效果？？？
ff=0;
for i=1:n
    ff=ff+fe(i);
end
ff1=fe/ff;
H1=0;
for i=1:n
    H(i)=-(1+ff1(i)).*(1+ff1(i));
    H1=H1+H(i);
end
for i=1:n
    w(i)=(1-H(i))/(n-H1);
end
ff2=0;w1=0;
for i=1:n 
    ff2=ff2+fe(i)*w(i);
    w1=w1+w(i);
end
%w,w1
%%%%%%%%%%%%%%%%%%优度计算  %%%%%%%%%%%%%%%%%%%%%
    ff3=1/(1+((1-f2)/f2))^2;
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if GreyEntropy==1;
    f2=ff2;% 使用灰熵计算;
elseif GreyEntropy==2;
    f2=ff3;
else
    f2=f2;%灰关联计算。
end
%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%2007.818%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%