Jacobi method not working

1 次查看（过去 30 天）
Joe Smith
Joe Smith 2016-7-9
评论： Azzi Abdelmalek 2016-7-9
Not sure what's wrong with my code. It does the first iteration correctly, then goes off in the wrong direction.
if true
% function x = Playing(A,b)
n = length(b);
x = zeros(1,n)';
y = x;
for k = 1:20
for i = 1:n
for j = 1:i-1,i+1:n
y(i) = y(i) + A(i,j)*x(j)
end
y(i) = (b(i) - y(i))/A(i,i)
end
end
end
end

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

回答（1 个）

Azzi Abdelmalek
Azzi Abdelmalek 2016-7-9
What is this: for j = 1:i-1,i+1:n ?
Maybe you mean for j = [1:i-1,i+1:n]
  2 个评论
Joe Smith
Joe Smith 2016-7-9
I've tried both, it didn't make a difference.
Azzi Abdelmalek
Azzi Abdelmalek 2016-7-9
for ii= 1:n
s=0
for jj = setdiff(1:n,ii)
s = s + A(ii,jj)*y(jj);
end
y(ii) = (b(ii) - s)/A(ii,ii);
end
and don't forget the initial guess y=ones(n,1) for example

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Logical 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by