帮我用matlab写一个从22个坐标中运用禁忌搜索算法选取符合规定的坐标进行重心法计算的坐标。其步骤如下： 1、要求22个点用重心法确定的选址地点为禁忌搜索算法的初始点x。 2、并且已经有三个固定的点A、B、C 3、判断22个点分别到点x的距离是不是比这个点到ABC三点都近。 4、如果22中的点到x最近，那么将站点留下用于下次重心法选址，否则加入禁忌表。并在下次循环中也用于判断。 5、保留下的站点用重心法确定新的站点x2 6、禁忌搜索的目标函数是各站点到A、B、C三个原有站点和禁忌搜索算法确定的站点距离最短值得和。

好的，这里给您提供一下 Matlab 的代码实现。代码如下：

% 确定初始点 x
x = centroid([x1, y1; x2, y2; x3, y3; ... ; x22, y22]);

% 确定固定点 A、B、C
A = [xA, yA];
B = [xB, yB];
C = [xC, yC];

% 禁忌表的初始值为空
tabuList = [];

% 禁忌搜索的循环次数
numIterations = 100;

% 禁忌搜索算法的主循环
for i = 1:numIterations
    
    % 判断每个站点是否符合规定
    for j = 1:22
        
        % 如果该站点已经在禁忌表中，则跳过
        if ismember(j, tabuList)
            continue;
        end
        
        % 判断该站点是否比 x 到 ABC 更近
        if norm([x(1), x(2)] - [x1, y1]) <= norm([x(1), x(2)] - A) && ...
           norm([x(1), x(2)] - [x2, y2]) <= norm([x(1), x(2)] - A) && ...
           norm([x(1), x(2)] - [x3, y3]) <= norm([x(1), x(2)] - A)
            % 如果是，则更新 x 为该站点
            x = centroid([x; [x(j, 1), x(j, 2)]]);
        else
            % 如果不是，则将该站点加入禁忌表
            tabuList = [tabuList, j];
        end
    end
    
    % 计算各站点到 ABC 和 x 的距离
    distances = zeros(1, 22);
    for j = 1:22
        distances(j) = norm([x(1), x(2)] - [x(j, 1), x(j, 2)]) + ...
                       norm(A - [x(j, 1), x(j, 2)]) + ...
                       norm(B - [x(j, 1), x(j, 2)]) + ...
                       norm(C - [x(j, 1), x(j, 2)]);
    end
    
    % 找到距离最短的站点，并将其从禁忌表中移除
    [minDistance, minIndex] = min(distances);
    tabuList = setdiff(tabuList, minIndex);
    
    % 如果禁忌表已经包含了所有站点，则退出循环
    if length(tabuList) == 22
        break;
    end
end

% 计算各站点到 ABC 和 x 的距离
distances = zeros(1, 22);
for j = 1:22
    distances(j) = norm([x(1), x(2)] - [x(j, 1), x(j, 2)]) + ...
                   norm(A - [x(j, 1), x(j, 2)]) + ...
                   norm(B - [x(j, 1), x(j, 2)]) + ...
                   norm(C - [x(j, 1), x(j, 2)]);
end

% 输出各站点到 ABC 和 x 的距离最短值得和
totalDistance = sum(distances);
disp(['The total distance is ', num2str(totalDistance)]);

需要注意的是，此代码中的 `centroid` 函数用于计算给定点集的重心坐标。您需要自己实现这个函数。