本文共 1409 字，大约阅读时间需要 4 分钟。

模拟退火算法求最大值：

clear,close all;clc;y=@(x)10*(x-2*(sin(x)).^2).^2 %目标函数%……初始参数表……%t1=100; %初始温度t0=1; %终止温度w=0.9 %衰减系数(经验在0.9-0.99之间，w越大搜索次数越多)xmax=13,xmin=-13;%搜索范围在[-13，13]之间xcurrent=0 %当前解xnew=0 %新解xbest=0%最优解%……初始参数表……%[hp,hn,hb]=show(y,xmax,xmin,xcurrent,xnew,xbest);%作图while t1>t0    for i=1:100%降一次温搜索100次        step=0.3;%搜索步长        xnew=xcurrent+step*xmax*(rand-0.5);%扰动(rand-0.5)是(-0.5,0.5)之间的随机数        if xnew>xmax||xnew
   
    y(xcurrent)            xcurrent=xnew;            %记忆化过程        if y(xnew)>y(xbest)            xbest=xnew;        end        %Metroplis准则        elseif exp(y(xnew)-y(xcurrent)./t1)>rand %(y(xnew)-y(xcurrent))<0            xcurrent=xnew;%以一定的概率接受新解        else            xcurrent=xcurrent;%否则不接受        end        %%更新图像        set(hn,'xdata',xnew,'ydata',y(xnew));        pause(0.1); %0.1秒暂停一下        set(hp,'xdata',xcurrent,'ydata',y(xcurrent));        set(hb,'xdata',xbest,'ydata',y(xbest));            end    t1=t1*w;%降温    xlabel(sprintf('最大值=%.1f',y(xbest)),'fontsize',13);%x轴显示动态的最大值，即从局部最优解变化到全局最优解的过程    title(sprintf('T=%.1f',t1),'fontsize',13);%标题显示温度的变化    endfunction [hp,hn,hb]=show(y,xmax,xmin,xcurrent,xnew,xbest)x=xmin:0.1:xmax;plot(x,y(x),'-k');hold onhp=plot(xcurrent,y(xcurrent),'o','markersize',10,'markerfacecolor','r');hn=plot(xnew,y(xnew),'o','markersize',10,'markerfacecolor','y');hb=plot(xbest,y(xbest),'o','markersize',10,'markerfacecolor','b')end
   

转载地址：http://nsgwi.baihongyu.com/