模拟退火matlab代码

以下是一个简单的模拟退火算法的 MATLAB 代码示例：

% 目标函数
f = @(x) sin(x) + sin(10/3 * x);

% 初始解和温度
x0 = 0;
T0 = 100;

% 终止温度和迭代次数
Tend = 1e-8;
maxiter = 1e4;

% 退火参数
alpha = 0.99;

% 初始化
x = x0;
T = T0;
fval = f(x);
iter = 0;

% 迭代
while T > Tend && iter < maxiter
    % 随机生成新解
    xnew = x + randn() * T;
    
    % 计算新解的目标函数值
    fnew = f(xnew);
    
    % 判断是否接受新解
    deltaf = fnew - fval;
    if deltaf < 0 || exp(-deltaf / T) > rand()
        x = xnew;
        fval = fnew;
    end
    
    % 降低温度
    T = alpha * T;
    
    % 更新迭代次数
    iter = iter + 1;
end

% 输出结果
fprintf('最优解为 %.4f，目标函数值为 %.4f\n', x, fval);

其中，`f` 是目标函数，`x0` 是初始解，`T0` 是初始温度，`Tend` 是终止温度，`maxiter` 是最大迭代次数，`alpha` 是退火参数。在迭代过程中，使用随机生成的新解，根据 Metropolis 准则判断是否接受新解，降低温度，更新迭代次数。最终输出最优解和对应的目标函数值。