MATLABで線形システムを解くJacobiのメソッド

Question

MATLABでこれをどのようにコード化しますか？

これは私が試したものですが、正しく動作しません。

function x = my_jacobi(A,b, tot_it) 
%Inputs: 
%A: Matrix 
%b: Vector 
%tot_it: Number of iterations 
%Output: 
%:x The solution after tot_it iterations 
    n = length(A); 
    x = zeros(n,1); 
    for k = 1:tot_it 
     for j = 1:n 
     for i = 1:n 
      if (j ~= i) 
       x(i) = -((A(i,j)/A(i,i)) * x(j) + (b(i)/A(i,i))); 

      else 
       continue; 
      end 
      end 
     end 
    end 
end 

Answer 1

j各iにわたる和のイテレータであるので、あなたは彼らの順序を変更する必要があります。また、数式には合計があり、コードには何も追加していないので、考慮する必要があります。あなたが省略している最後の事は、式の右辺にそれが必要なので、前の状態であるxを保存しなければならないということです。次のようなものを試してください：

function x = my_jacobi(A,b, tot_it) 
    %Inputs: 
    %A: Matrix 
    %b: Vector 
    %tot_it: Number of iterations 
    %Output: 
    %:x The solution after tot_it iterations 
    n = length(A); 
    x = zeros(n,1); 
    s = 0; %Auxiliar var to store the sum. 
    xold = x 
    for k = 1:tot_it 
    for i = 1:n 
     for j = 1:n 
     if (j ~= i) 
      s = s + (A(i,j)/A(i,i)) * xold(j); 
     else 
      continue; 
     end 
     end 
     x(i) = -s + b(i)/A(i,i); 
     s = 0; 
    end 
    xold = x; 
    end 
end