clear

clc

a=rand(8,21)

lamuda=4

for i=1:7

b(:,i)=(a(:,3*i-1)+(a(:,3*i)+a(:,3*i-2)*lamuda)/(lamuda+1))/2;

end

b

w=[0.1415,0.1698,0.1604,0.1321,0.1792,0.0943,0.1226]

for i=1:7

   c(:,i)=b(:,i).*w(i);

end

c

S_Max=max(c);

S_Min=min(c);

   for j=1:7

     tt1(:,j)=  min(abs(c(:,j)-S_Max(j)));

     tt2(:,j)=  max(abs(c(:,j)-S_Max(j)));

   end

 MMMin=min(tt1)  

 MMMax=max(tt2)

 rou=0.6

 for i=1:8

     for j=1:7

         kesi1(i,j)=(MMMin+rou*MMMax)/(abs(S_Max(j)-c(i,j))+rou*MMMax);

     end

 end

 kesi1

for i=1:8

   fa1(i)=-kesi1(i,:)*(log(kesi1(i,:)))';

end

fana=fa1./(max(fa1))

   for j=1:7

     tt3(:,j)=  min(abs(c(:,j)-S_Min(j)));

     tt4(:,j)=  max(abs(c(:,j)-S_Min(j)));

   end

 MMMin=min(tt3)  

 MMMax=max(tt4)

 rou=0.6

 for i=1:8

     for j=1:7

         kesi2(i,j)=(MMMin+rou*MMMax)/(abs(S_Min(j)-c(i,j))+rou*MMMax);

     end

 end

 kesi2

for i=1:8

   fa2(i)=-kesi2(i,:)*(log(kesi2(i,:)))';

end

fan2=fa2./(max(fa2))

fanZ=fana+fan2

for i=1:8

   p(i)=fana(i)/fanZ(i);

end

p