Jacobi iteration what the wrong of my code please help

2 次查看（过去 30 天）

ayah zaghlol 2021-5-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/827915-jacobi-iteration-what-the-wrong-of-my-code-please-help

评论： Walter Roberson 2021-5-13

% jacobi Method
 clc
 clear
 format short
 A=[2.08 -1 0 0;-1 2.08 -1 0;0 -1 2.08 -1;0 0 -1 2.08];
 b=[200.0 0.8 0.8 40.8];
 l=tril(-A,-1);
 u=triu(-A,1);
 d=diag(diag(A));
 tj=inv(d)*(l+u);
 cj=inv(d)*b;
 xj=[1;1;1];
 N=5;
 for i=1:N
   xj=tj*xj+cj;
 [i xj'] 
end

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

请先登录，再进行评论。

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

采纳的回答

Walter Roberson 2021-5-12

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/827915-jacobi-iteration-what-the-wrong-of-my-code-please-help#answer_697755

size(d) is 4 x 4, so size(inv(d)) is 4 x 4.

size(b) is 1 x 4.

You cannot use matrix multiplication * operator between a 4 x 4 and a 1 x 4. The second dimension of the first operand (4 here) must match the first dimension of the second operand (1 here)