雅可比(jacobi)迭代法 matlab實現
阿新 • • 發佈:2018-12-15
clc clear n = input('請輸入矩陣階數:\n'); for i = 1:n fprintf('請輸入矩陣第%d行\n',i); A(i,:) = input(''); end A B(:,1) = input('請輸入B向量:\n'); B for i = 1:n x(i) = 0; x2(i) = 0; end for i =1:4 for j = 1:n for k = 1:n if j ~= k x(j) = x(j)-A(j,k)*x2(k); end end x(j)=x(j) + B(j); x(j)=x(j)/A(j,j); end x2=x; for j = 1:n x(j)=0; end fprintf('第%d次迭代結果:\n',i); x2 end