mSAA程序及流程图.docxVIP

 PAGE \* MERGEFORMAT 5 带记忆的模拟退火算法流程图 开始 输入参数（T,Ts,Lm,Km）和初始解（X1） 调用目标函数objf(X)计算f(X1) 用xp和fp分别记忆X1和f(X1) 开始迭代计数L=1 开始mapkob链计数K=1 利用解变换函数Generate(X)和X1生成X2 调用目标函数objf(X)计算f(X1) 用0-1均布随机矩阵rand生成r 计算X2与X1的函数差df=f(X2)-f(X1) df=0 yes no exp((-df)/T) r X1=X2 f(X1)=f(X2) yes no 比较fp和f(X1),若f(X1)小，用xp和fp分别对X1和f(X1) 予以记忆 比较fp和f(X1),若f(X1)小，用xp和fp分别对X1和f(X1) 予以记忆 L=L+1 K=K+1 KKm L Lm no yes T=0.99*T TTs no yes 结束 总程序 %******************带记忆式模拟退火算法总程序***************** clear; clc; T=input(please input initial temperature T:); %输入初始温度 Ts=input(please input stop temperature Ts:); %输入终止温度 Km=input(please input available maximum iterative number Km:); %输入允许的最大的次数 Lm=input(please input available maximum Mapkob chain number Lm:); %输入允许的最大Mapkob链的长度 X1=input(please input initial solution X1:); %输入初始解 f1=objf(X1); %调用目标函数子程序计算初始解的函数值 xp=X1; %记忆初始解 fp=f1; %记忆初始解的函数值 while T=Ts %开始退火 K=1; %迭代计数开始 L=1; %链长计数开始 while KKm LLm X2=Generate(X1); %调用发生器函数子程序生成X2 f2=objf(X2); r=rand; f=f2-f1; if f0 X1=X2; %接受新解 f1=f2; T=T*0.99; %退火降温 if fpf1 xp=X1; %记忆新解 fp=f1; end elseif (exp(f1-f2)/(T)r) X1=X2; f1=f2; T=0.9