MATLAB编写Jacobi迭代法报Index exceeds array bounds错误如何解决

该报错核心原因是代码执行时访问了超出数组实际长度的索引位置，针对Jacobi迭代的常见代码实现，触发该报错的常见原因及修复方案如下：

• 维度不匹配：首先确认系数矩阵A为n阶方阵，右端向量b、初始迭代向量x0的长度均为n，三者维度只要有1位偏差就会触发越界，可通过size(A)、length(b)、length(x0)输出维度信息确认。
• 循环下标越界：MATLAB数组下标从1开始计数，若遍历未知数的循环变量范围写成了0:n-1，或者超出1:n的区间，会直接触发越界，需将循环范围修正为1:size(A,1)。
• 求和逻辑边界错误：Jacobi迭代更新公式中需要排除对角线元素求和，若你手动编写了j的循环做求和，需确认i=1时不会访问x(0)、i=n时不会访问x(n+1)，推荐直接使用向量化求和写法，MATLAB对空数组求和默认返回0不会触发错误：

  s = A(i,1:i-1)*x(1:i-1) + A(i,i+1:end)*x(i+1:end);
  

• 迭代变量未正确初始化：若迭代过程中用到的x_new这类变量未提前初始化为长度n的向量，动态扩展失败也会触发越界，需在迭代开始前统一初始化所有迭代变量的维度。

你可以参考以下标准Jacobi迭代实现逐行对比修正：

function [x, iter_num] = jacobi_solve(A, b, tol, max_iter)
    n = size(A, 1);
    % 前置维度校验，避免维度不匹配问题
    if size(A, 2) ~= n
        error('系数矩阵A必须为方阵');
    end
    if length(b) ~= n
        error('右端向量b长度需等于矩阵A的阶数');
    end
    % 初始化迭代变量
    x = zeros(n, 1);
    x_new = x;
    iter_num = 0;
    
    while iter_num < max_iter
        for i = 1:n
            sum_term = A(i, [1:i-1, i+1:n]) * x([1:i-1, i+1:n]);
            x_new(i) = (b(i) - sum_term) / A(i, i);
        end
        % 收敛判断
        if norm(x_new - x, inf) < tol
            x = x_new;
            return;
        end
        x = x_new;
        iter_num = iter_num + 1;
    end
    warning('达到最大迭代次数未收敛');
end

若仍无法定位问题，可在报错行前添加调试语句输出对应数组的维度和当前索引值，即可快速定位故障点。

内容的提问来源于stack exchange，提问作者Mimi