newstate1.m

来自「从国外专业网站下载的非常不错的模拟退火算法」· M 代码 · 共 11 行

function [A,X,Y]=NEWSTATE1(Aceng,Xceng,Yceng,EE)
Emax=48.118;
if round(rand(1))==1
    A=Aceng+(5-Aceng)*(1-rand(1).^((EE/Emax).^2));%(1-rand(1).^((EE/Emax).^2)当rand(1)固定时为EE的增函数，而rand(1)越小此值越大
    X=Xceng+round((20-Xceng)*(1-rand(1).^((EE/Emax).^2)));
    Y=Yceng+round((20-Yceng)*(1-rand(1).^((EE/Emax).^2)));
else 
    A=Aceng-Aceng*(1-rand(1).^((EE/Emax).^2));
    X=Xceng-round(Xceng*(1-rand(1).^((EE/Emax).^2)));
    Y=Yceng-round(Yceng*(1-rand(1).^((EE/Emax).^2)));
end