ywsss.m

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

function [minz,minx,minD]=YWSSS(x,D)
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];
dh=[0.001,-0.001];
minx2=ones(2,7);
minx1=ones(1,7);
        xx1=x;
        z=Czfee(xx1,D);
        minz=3000000;
        minz1=z;
        for k=1:7
            for l=1:2
                z1=z;
                xx2=xx1;
                for j=1:1
                   u=1;
                   z2=z1;
                   while z2<=z1
                       z1=z2;
                       xx1(k)=u*dh(l)+xx1(k);
                       if xx1(k)>b(k)|xx1(k)<a(k)
                           break;
                       end
                       z2=Czfee(xx1,D);
                       u=u*2;
                   end
                   if z2>z1
                       u=u/2;
                   end
                   xx1(k)=xx1(k)-dh(l)*u;
                end
                z2=z1;
                while z2<=z1
                    z1=z2;    
                    xx1(k)=xx1(k)+dh(l);
                    if xx1(k)>b(k)|xx1(k)<a(k)
                        break;
                    end
                    z2=Czfee(xx1,D);
                end
                xx1(k)=xx1(k)-dh(l);
                minz2(l)=z1;
                minx2(l,:)=xx1;
                xx1=xx2;
            end
            minz1=min(minz2);
            if minz1==minz2(1)
                n=1;
            else
                n=2;
            end
            minx1=minx2(n,:);
            xx1=minx1;
            if minz>minz1
                minx=minx1;
                minz=minz1;
            end
        end
minD=D;