matlab运行后没结果怎么办？

冷凌川

A=[1,5,-3;5,-2,1;2,1,-5];
b=[2;4;-11];
x0=ones(size(b));
D=diag(diag(A));
U=triu(A,1);
L=tril(A,-1);
G=-(D+L)\U;
d1=(D+L)\b;
x=G*x0+d1;
n=1;
while norm(x-x0)>=1e-4
x0=x;
x=G*x0+d1;
n=n+1;
end
x
n

x =

  -Inf
   NaN
  -Inf


n =

   315

ChaChing

回复 1 # 冷凌川 的帖子

检查下算法是否正确? 发散!?

meiyongyuandeze

我看了你的程序，应该是用Gauss-Seidel 迭代方法求线性方程组的解的问题。程序编写有个错的地方：G=-(D+L)\U，应该是G=(D+L)\U。
但即使修改了还是出现发散的结果，原因在于该线性方程组的迭代矩阵G并不满足迭代收敛的充分必要条件即谱半径小于1. 你的这个问题的谱半径为9.5958>1,因此得不到收敛解。
你可以事先根据下面的两个条件来判断使用迭代法是不是收敛：
1.看系数矩阵是不是严格对角占优
2. 若系数矩阵为正定的对称阵则可采用Gauss-Seidel 迭代方法。
你可以看下，你的系数矩阵A不是正定的。
希望对你有帮助！