zkp.m

来自「该文件主要包括了matlab的一些主要应用例程!对学习matlab很有帮助」· M 代码 · 共 42 行

function [X,P]=zkp(w,c,M,t0,tf)
[m,n]=size(w);

L=100*n;
t=t0;

clear m;
x=zeros(1,n);
xmax=x;
fmax=0;

while t>tf
   for k=1:L
      xr=change(x);
      gx=g_0_1(w,x);
      gxr=g_0_1(w,xr);
      if gxr<=M
         fr=f_0_1(xr,c);
         f=f_0_1(x,c);
         df=fr-f;
         if df>0
            x=xr;
            if fr>fmax
               fmax=fr;
               xmax=xr;
            end
         else
            p=rand;
            if p<exp(df/t)
               x=xr;
            end
         
         end
      end
   end
   t=0.87*t;
end
P=fmax;
X=xmax;