simulated.asv

模拟退火matlab源码

function [x0,y0,k]=simulated
% aim0=forward(x0,y0);
T0=10;
K=4000;
k=0;
esp=forward(3,0);
syms x0 y0 x1 y1;
while k<=K&aim0>esp
%       x0=12.0*rand-6.0;
%        y0=12.0*rand-6.0; 
% deltax=0.01*rand;
% deltay=0.01*rand;
% x1=x0+deltax;
% y1=y0+deltay;
%   model_xy=12.0*rand(1,2)-6.0;
% x1=x0*rand;
% y1=y0*rand;
x1=12*(rand(1)-0.5);
y1=12*(rand(1)-0.5);
% x0=12.0*rand(1)-6.0;
% y0=12.0*rand(1)-6.0;
l=0;
 while  l<10%(-6<=x1<=6)&(-6<=y1<=6)&  %x1<max&x1>min&y1<max&y1>min
% OBJFUN=100;
%        aim0=forward(x0,y0);
%         x0=unifrnd(-6,6);
%         y0=unifrnd(-6,6);
deltax=0.01*rand(1);
deltay=0.01*rand(1);
%         model_XY=model_xy+model_change;
x1=x0+deltax;
y1=y0+deltay;
aim0=forward(x0,y0);
aim1=forward(x1,y1);
delta_aim=aim1-aim0;
p_delta_aim=exp(-delta_aim/T0);
if(delta_aim<0)
    x0=x1;
    y0=y1;
elseif p_delta_aim>rand
    x0=x1;
    y0=y1;
end
l=l+1;
% aim0=aim1;
end
% obj_fun=forward(x0,y0);
% if(obj_fun<=OBJFUN)
%             xmin=x0;
%             ymin=y0;
%         end
%         OBJFUN=obj_fun;
%         l=l+1;
%        x0=xmin;
%        y0=ymin;
%     end
% aim0=aim1;
T0=T0*0.99;
if(abs(aim0-esp)<=10.^(-4))
            break;
end
        
k=k+1;
end