Matlab求解117×117线性方程组遇奇异矩阵报错，如何解决？

1. 先明确矩阵的本质问题 #

• 计算矩阵秩：执行rank(A)，若结果小于117，说明矩阵确实奇异，不存在唯一解；若秩等于117但cond(A)远大于1e10，说明矩阵是病态矩阵（接近奇异），单纯提高精度无法解决核心问题。
• 排查矩阵有效性：用any(isnan(A(:)))或any(isinf(A(:)))检查矩阵中是否混入了NaN/Inf元素，这类元素会直接导致求解失败。

2. 奇异矩阵的处理方案（无唯一解场景） #

若矩阵秩不足：

• 求最小范数解：使用lsqminnorm(A,b)，该函数专门针对奇异/病态方程组，返回稳定性最优的最小范数解，不会输出NaN。
• 构造全解空间：结合伪逆与零空间基，示例代码：

  x_special = pinv(A)*b; % 特解
  null_space = null(A); % 零空间基向量组
  % 所有解可表示为：x = x_special + null_space*c，c为任意系数向量
  

3. 病态矩阵的优化方法（秩满但条件数极高场景） #

• 正则化修正：加入小正则项平衡稳定性与精度，示例：

  lambda = 1e-8; % 可根据条件数调整，范围1e-10~1e-6
  x = (A'*A + lambda*eye(size(A)))\(A'*b);
  

• 稳定分解求解：用QR分解替代直接求解，QR分解对病态矩阵的稳定性更强：

  [Q,R] = qr(A);
  x = R\(Q'*b);
  

• 删除冗余方程：手动检查矩阵行/列，删除重复或近似线性相关的行，降低矩阵的病态程度。

4. 纠正digits(64)的使用误区 #

digits(64)仅对符号计算生效，常规数值求解（\、linsolve）不会读取该设置。若要结合高精度符号计算，需将矩阵转为符号类型：

A_sym = sym(A);
b_sym = sym(b);
x_sym = linsolve(A_sym, b_sym);
x = double(x_sym);

但需注意：若矩阵本质奇异，符号计算也会提示无解，此时仍需回到前面的解空间构造方法。

5. 排查建模根源错误 #

• 检查方程组构建：确认物理模型、边界条件、参数输入是否正确，是否存在矛盾约束或重复方程。
• 回溯矩阵生成过程：检查循环赋值、索引操作是否正确，避免因代码逻辑错误导致矩阵结构异常。

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