zzf.m

#include "stdio.h" #include "math.h" #include "malloc.h" #include "string.h" #define m 6 #defin

minz=3000000;
for r=1:15
A=0.01/3;
B=0.05/3;
C=0.1/3;
e=500;
q=0;
a=[0.075,0.225,0.075,0.075,1.125,12,0.5625];
b=[0.125,0.375,0.125,0.125,1.875,20,0.935];
x1=ones(1,7);
for i=1:7
     x1(1,i)=unifrnd(a(i),b(i),1,1);
end
D1=[B,C,C,C,C,C,B];
while 1
[z1,x1,D1]=YWSSS(x1,D1);
A=0.01/3;
B=0.05/3;
C=0.1/3;
DJ=[A,B,C];
D=ones(108,7);
D(:,1)=B*ones(108,1);
D(:,5)=C*ones(108,1);
for i=2:3
    for j=1:3
        for k=1:3
            for l=1:3
                for m=1:2
                    D((i-2)*2*3^3+(j-1)*3*3*2+(k-1)*3*2+(l-1)*2+m,2)=DJ(i);
                    D((i-2)*2*3^3+(j-1)*3*3*2+(k-1)*3*2+(l-1)*2+m,3)=DJ(j);
                    D((i-2)*2*3^3+(j-1)*3*3*2+(k-1)*3*2+(l-1)*2+m,4)=DJ(k);
                    D((i-2)*2*3^3+(j-1)*3*3*2+(k-1)*3*2+(l-1)*2+m,6)=DJ(l);
                    D((i-2)*2*3^3+(j-1)*3*3*2+(k-1)*3*2+(l-1)*2+m,7)=DJ(m);
                end;end;end;end;end
for i=1:108
    zz(i)=Czfee(x1,D(i,:));
end
z2=min(zz);
for i=1:108
    if z2==zz(i)
        x2=x1;
        D2=D(i,:);
        break;
    end
end
if z1-z2<e
    zzz=z2;
    xxx=x2;
    DDD=D2;
    break;
end
z1=z2;
x1=x2;
D1=D2;
q=q+1;
end
if minz>z1
    minz=z1;
    minx=x1;
    minD=D1;
end
minz
r
end