gusucode.com > 十大算法matlab程序说明 > 十大算法matlab程序说明/遗传退火法/Anneal.m
function [X2,fX2] = Anneal(x,fx,FitnessFcn,B,stepSize,A,LB,UB,R,L,t) GenomeLength = length(x); Lbnd = B(1,:); Ubnd = B(2,:); span = max([sqrt(eps)*x;min([x-Lbnd;Ubnd-x])]); scale = stepSize * span; X2=[];fX2=[]; r = 0; l = 0; while r<R && l<L x2 = x + scale .* randn(1,GenomeLength); y2 = (A*x2')'; if all(LB'<y2 & y2<UB'); l = l+1; fx2 = feval(FitnessFcn,x2); delta_f = fx2 - fx; if delta_f < 0 r = r+1; x = x2; fx = fx2; span = min([x-Lbnd;Ubnd-x]); scale = stepSize * span; X2 = [X2;x2]; fX2 = [fX2;fx2]; elseif exp(-delta_f/t) > rand r = r+1; x = x2; fx = fx2; span = min([x-Lbnd;Ubnd-x]); scale = stepSize * span; X2 = [X2;x2]; fX2 = [fX2;fx2]; end end end