main.m

本文以著名的旅行商问题(TSP)为例说明如何利用MATLAB语言实现模拟退火算法。

%TSP问题模拟退火算法

global T_cur;   %当前温度
global T_ini;   %初始温度
global T_end;   %终止温度
global city_n;  %城市数量
global Markov_len;  %内层循环马式链长度
global solution;
global solution_everbest;   %历代最优解路径
global adapt_everbest;  %历代最优解值
global adapt_cur;       %当前解值
global adapt_last;      %上一代解值
global solution_temp;   %复制临时解集

adapt_everbest=inf;
adapt_last=inf;
solution(30)=0;
T_ini=50;
T_end=1;
T_cur=T_ini;
city_n=30;
Markov_len=10;

initial;    %初始化解、温度等参数

while T_cur>T_end   %未退火到终止温度
    for i=1:Markov_len  
        pause(0.01);
        create;     %状态产生函数
        adapting;   %计算新解适应值
        keepbest;
        accept;     %状态接收函数
        paint;      %绘制回路图形
    end
    my_ada_file=fopen('value.dat','a');
    fprintf(my_ada_file,',%f',adapt_cur);
    fclose(my_ada_file);
    annealing;  %退温函数
end

result;