clc
clear
format?long
N=30;
A=zeros（NN）;
for?i=1:N-1
???A（ii）=3;
???A（ii+1）=-1;
???A（i+1i）=-1;
end
A（NN）=3;
b=ones（N1）;
b（1）=2;
b（N）=2;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%雅克比迭代
D=diag（diag（A））;
L=tril（A-1）;
U=triu（A1）;
BJ=inv（D）*（L+U）;
fJ=inv（D）*b;

x=ones（N1）;
newx=x;
while?true
???newx=BJ*x+fJ;?
???error=abs（norm（newx1）-norm（x1））;%用范數(shù)作為終止條件
???if?error<1e-6
???????break
???end
???x=newx;
end
%雅克比迭代結(jié)果
fprintf（‘1.雅克比迭代結(jié)果為:‘）?
Jx=round（newx6）;
vpa（Jx6）

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%高斯-賽德爾迭代
BG=inv（（D-L））*U;
fG=inv（（D-L））*b;
while?true
???newx=BG*x+fG;?
???error=abs（norm（newx1）-norm（x1））;%用范數(shù)作為終止條件
???if?error<1e-6
???????break
???end
???x=newx;
end
%高斯-賽德爾迭代結(jié)果
fprintf（‘2.高斯-賽德爾迭代結(jié)果:‘）?
Gx=round（newx6）;
vpa（Gx6）
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????51152??2019-12-17?11:07??2-7\2-7文檔.docx
?????文件?????????969??2019-12-17?10:54??2-7\m2_7.m
?????目錄???????????0??2019-12-26?09:22??2-7\