%---雅克比迭代法-----
%---Jacobi?iteration?method
clear;clc;
%?A=[10-1-2;-110-2;-1-15];
%?b=[728342]‘;
A=[?28-3000;
????-338-100-5;
????-10025-150;
????00-15450;
????0-50030];
b=[100000]‘;
N=length（b）;
fprintf（‘庫函數計算結果：‘）;
x=inv（A）*b???%庫函數計算結果
B=zeros（NN）;g=zeros（N1）;
x=zeros（N1）;%迭代初始值
eps=0.001;%相鄰解的距離小于該數時，結束迭代
%?for?i=1:N
%?????B（i:）=-A（i:）/A（ii）;
%?????B（ii）=0;???????%迭代矩陣
%?????g（i）=b（i）/A（ii）;
%?end
D=diag（diag（A））;
E=-tril（A-1）;%下三角
F=-triu（A1）;%上三角
B=inv（D）*（E+F）;g=inv（D）*b;
for?k=1:100?%最大迭代次數為100
????fprintf（‘第%d次迭代：‘k）;
????y=B*x+g;
????fprintf（‘\n與上次計算結果的距離（2范數）:%f?\n‘norm（x-y）^2）;
????if?norm（x-y）????????break;
????end
????x=y
end
x

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????813??2020-06-10?22:23??程序\Jacobi.m

?????文件????????636??2020-06-10?22:34??程序\SOR.m

?????文件????????606??2020-06-10?22:49??程序\SOR1.m

?????目錄??????????0??2020-06-10?23:04??程序

-----------?---------??----------?-----??----

?????????????????2055????????????????????4