模拟退火算法-MATLAB

%% SA 模拟退火: 求解函数y = 11*sin(x) + 7*cos(5*x)在[-3,3]内的最大值(动画演示）

tic

clear; clc

%% 绘制函数的图形

x = -3:0.1:3;

y = 11*sin(x) + 7*cos(5*x);

figure

plot(x,y,'b-')

hold on % 不关闭图形，继续在上面画图

%% 参数初始化

narvs = 1; % 变量个数

T0 = 100; % 初始温度

T = T0; % 迭代中温度会发生改变，第一次迭代时温度就是T0

maxgen = 200; % 最大迭代次数

Lk = 100; % 每个温度下的迭代次数

alfa = 0.95; % 温度衰减系数

x_lb = -3; % x的下界

x_ub = 3; % x的上界

%% 随机生成一个初始解

x0 = zeros(1,narvs);

for i = 1: narvs

x0(i) = x_lb(i) + (x_ub(i)-x_lb(i))*rand(1);

end

y0 = Obj_fun1(x0); % 计算当前解的函数值

h = scatter(x0,y0,'*r'); % scatter是绘制二维散点图的函数（这里返回h是为了得到图形的句柄