%*********ZIBAI*********
%Jacobi迭代求解

clear?all;
format?long;

%設(shè)置方程參數(shù)
A?=?[10-1-2;-110-2;-1-15];
b?=?[7.2;8.3;4.2];
x0?=?[1.1;1.2;1.3];
%設(shè)置精度，最大迭代次數(shù)，
eps?=?1.0e-8;
M?=?200;

%Jacobi迭代求解
D?=?diag（diag（A））%求A的對(duì)角矩陣
L?=?-tril（A-1）%求A的下三角矩陣
U?=?-triu（A1）%求A的上三角矩陣
B?=?D\（L+U）;
f?=?D\b;
x00?=?[0;0;0];
x?=?B*x00+f;
n?=?1;%迭代次數(shù)
while?norm（（x-x0）?‘inf‘）>=eps
????xk?=?x;
????x?=?B*xk+f;
????n?=?n+1;
????if（n>=M）
????????disp（‘Warning:迭代次數(shù)太多可能不收斂!‘）
????????break;
????end
end

fprintf?（‘%d次迭代后，相鄰無窮范數(shù)誤差小于：1.0e-8\n‘?n）
fprintf?（‘迭代解：‘）;x

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-10-12?20:12??Jacobi迭代\
?????文件?????????658??2018-10-12?20:11??Jacobi迭代\Jacobi.m