模拟退火算法Matlab源程序

MCM战备历程3（模拟退火算法Matlab源程序）For glory

2007-02-03 11:20:04| 分类： 数学建模 | 标签：学习 |字号 订阅

%模拟退火算法程序

T_max=input('please input the start temprature'); T_min=input('please input the end temprature');

iter_max=input('please input the most interp steps on the fit temp'); s_max=input('please input the most steady steps ont the fit temp'); T=T_max; load d:\\address.txt;

order1=randperm(size(address,1))';%生成初始解。 plot(address(order1,1),address(order1,2),'*r-') totaldis1=distance(address,order1); while T>=T_min iter_num=1; s_num=1; plot(T,totaldis1) hold on

while iter_num

totaldis2=distance(address,order2); R=rand;

DeltaDis=totaldis2-totaldis1; if DeltaDis<0; order1=order2; totaldis1=totaldis2;

elseif (exp((totaldis1-totaldis2)/(T))>R) order1=order2; totaldis1=totaldis2; else s_num=s_num+1;

end

iter_num=iter_num+1; end T=T*0.99; end order1 totaldis1 figure(2)

plot(address(order1,1),address(order1,2),'*r-')

function y=distance(address,order) nmb=size(address,1); y=0;

for i=1:nmb-1

y=y+sqrt((address(order(i+1),1)-address(order(i),1))^2+(address(order(i+1),2)-address(order(i),2))^2); end

y=y+sqrt((address(order(i+1),1)-address(order(1),1))^2+(address(order(i+1),2)-address(order(1),2))^2);

function y=exhgpath(order) while 1

b=size(order,1); r=unidrnd(b,1,2); if r(1)-r(2)~=0 break end end

b=order(r(2)); order(r(2))=order(r(1)); order(r(1))=b;

（模拟退火算法Matlab源程序）

学术讨论 2008-04-01 18:05:56 阅读1322 评论6 字号：大中小 订阅

（模拟退火算法Matlab源程序）For glory

%模拟退火算法程序

T_max=input('please input the start temprature'); T_min=input('please input the end temprature');

iter_max=input('please input the most interp steps on the fit temp'); s_max=input('please input the most steady steps ont the fit temp'); T=T_max; load d:\\address.txt;

order1=randperm(size(address,1))';%生成初始解。 plot(address(order1,1),address(order1,2),'*r-') totaldis1=distance(address,order1); while T>=T_min iter_num=1; s_num=1; plot(T,totaldis1) hold on

while iter_num

totaldis2=distance(address,order2); R=rand;

DeltaDis=totaldis2-totaldis1; if DeltaDis<0; order1=order2; totaldis1=totaldis2;

elseif (exp((totaldis1-totaldis2)/(T))>R) order1=order2; totaldis1=totaldis2; else s_num=s_num+1; end

iter_num=iter_num+1; end T=T*0.99; end order1 totaldis1 figure(2)

plot(address(order1,1),address(order1,2),'*r-')