重心法选址matlab代码

当在 MATLAB 中使用重心法进行选址时，可以采用以下代码实现：

function [x_opt, y_opt] = gravity_method(locations, masses)
    % 初始化变量
    n = size(locations, 1);
    x = locations(:, 1);
    y = locations(:, 2);
    
    % 定义迭代次数和阈值
    max_iter = 1000;
    eps = 1e-6;
    
    % 初始化重心位置
    x_center = sum(x) / n;
    y_center = sum(y) / n;
    
    % 迭代计算重心位置
    for iter = 1:max_iter
        % 计算每个位置与重心的距离的平方
        d_square = (x - x_center).^2 + (y - y_center).^2;
        
        % 根据质量和距离更新重心位置
        x_center_new = sum(masses .* x) / sum(masses);
        y_center_new = sum(masses .* y) / sum(masses);
        
        % 判断是否收敛
        if abs(x_center_new - x_center) < eps && abs(y_center_new - y_center) < eps
            break;
        end
        
        % 更新重心位置
        x_center = x_center_new;
        y_center = y_center_new;
    end
    
    % 返回最优重心位置
    x_opt = x_center;
    y_opt = y_center;
end

在使用该函数时，`locations` 参数是一个 n 行 2 列的矩阵，每一行表示一个候选位置的坐标；`masses` 参数是一个 n 行 1 列的矩阵，表示每个位置的质量。函数会返回选址结果的 x 坐标和 y 坐标。