模拟退火算法matlab代码

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

% 模拟退火算法
% f: 目标函数
% x0: 初始解
% T0: 初始温度
% alpha: 温度衰减因子
% L: 每个温度下的迭代次数
% T_min: 最小温度
% delta: 邻域半径
function [x, fx] = SA(f, x0, T0, alpha, L, T_min, delta)
    x = x0;
    fx = f(x0);
    T = T0;
    while T > T_min
        for i = 1:L
            % 在邻域内随机生成新解
            x_new = x + delta * (2*rand(size(x)) - 1);
            fx_new = f(x_new);
            % 计算接受概率
            delta_fx = fx_new - fx;
            p = exp(-delta_fx / T);
            % 根据接受概率决定是否接受新解
            if delta_fx < 0 || rand() < p
                x = x_new;
                fx = fx_new;
            end
        end
        % 降温
        T = alpha * T;
    end
end

其中，`f` 表示目标函数，`x0` 表示初始解，`T0` 表示初始温度，`alpha` 表示温度衰减因子，`L` 表示每个温度下的迭代次数，`T_min` 表示最小温度，`delta` 表示邻域半径。函数返回最优解 `x` 和最小值 `fx`。