Trying to run Gauss Siedel method, I have no idea where I'm going wrong in this code. I get an error of "Index exceeds matrix dimensions" Not sure how? Just need this to run Thanks

Question

Thomas Susi 2016-10-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/305945-trying-to-run-gauss-siedel-method-i-have-no-idea-where-i-m-going-wrong-in-this-code-i-get-an-erro

回答： Walter Roberson 2016-10-5

% Solve system from HW using Gauss Seidel
% System of Equations Simplified Below
% 2x1-3x2=-160 ---> x1=(3x2-160)/2
% x1+2x2-3x3 -----> x2=(3x3-x1+40)/2
% 2x3-x2=80 ------> x3=(x2+80)/2
clear;clc;
i=1;
x2(i)=0;x3(i)=0;
error_x1(i)=9999;
while error_x1(i) >= .01;
        x1(i+1)=(3*x2(i)-160+0*x3(i))/2;
        x2(i+1)=(3*x3(i)-x1(i+1)+40)/2;
        x3(i+1)=(x2(i+1)+80+0*x1(i+1))/2;
          errorx1(i+1)=abs((x1(i+1)-x1(i))/x1(i+1))*100;
          errorx2(i+1)=abs((x2(i+1)-x2(i))/x2(i+1))*100;
          errorx3(i+1)=abs((x3(i+1)-x3(i))/x3(i+1))*100;
          i=i+1;
  end
disp('                  x1               error(%)');
disp([x1',error_x1']);
disp('                  x2               error(%)');
disp([x2',error_x2']);
disp('                  x3               error(%)');
disp([x3',error_x3']);

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Marc Jakobi 2016-10-5

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/305945-trying-to-run-gauss-siedel-method-i-have-no-idea-where-i-m-going-wrong-in-this-code-i-get-an-erro#answer_237625

在 MATLAB Online 中打开

You have to pre-initialize x1, x2, x3, errorx1, errorx2 and errorx3. The index exceeds the matrix dimensions, because you index x1(i+1), which is x1(2) in the first iteration of your loop. Since x1 is initialized with a length of 1, the index 2 exceeds the dimension 1.

One solution would be

clear;clc;
i=1;
errorx1=9999;
errorx2 = 9999;
errorx3 = 9999;
x1 = 0;
x2 = 0;
x3 = 0;
while errorx1(i) >= .01
        x1 = [x1; (3*x2(i)-160+0*x3(i))/2];
        x2 = [x2; (3*x3(i)-x1(i+1)+40)/2];
        x3 = [x3; (x2(i+1)+80+0*x1(i+1))/2];
          errorx1 = [errorx1; abs((x1(i+1)-x1(i))/x1(i+1))*100];
          errorx2 = [errorx2; abs((x2(i+1)-x2(i))/x2(i+1))*100];
          errorx3 = [errorx3; abs((x3(i+1)-x3(i))/x3(i+1))*100];
          i=i+1;
end
disp('                  x1               error(%)');
disp([x1',errorx1']);
disp('                  x2               error(%)');
disp([x2',errorx2']);
disp('                  x3               error(%)');
disp([x3',errorx3']);