1:n,1))/a(i,i);enddisp(k)disp(x')end>>jacobi
1:n,1))/a(i,i);enddisp(k)disp(x')end>>jacobi
x(i)=(b(i)-l(i,1:i-1)*x0(1:i-1,1)-u(i,i+1:n)*x0(i+
x(i)=(b(i)-l(i,1:i-1)*x0(1:i-1,1)-u(i,i+1:n)*x0(i+
,1);while norm(x0-x,inf)>e,k=k+1;x0=x;for i=1:n
,1);while norm(x0-x,inf)>e,k=k+1;x0=x;for i=1:n
=length(b);x=x0;x0=x+2*e;k=0;l=tril(a,-1);u=triu(a
=length(b);x=x0;x0=x+2*e;k=0;l=tril(a,-1);u=triu(a
tion x=jacobi(a,b,x0,e)% 用途：用分量形式的sor迭代解线性方程组ax=bn
tion x=jacobi(a,b,x0,e)% yong tu yong fen liang xing shi de sor die dai jie xian xing fang cheng zu ax=bn
代次数上限');end&g网通传奇sft;> sor(a,b,1.45,[1 11]',1e-6)func
dai ci shu shang xian ');end&g wang tong chuan qi sft;> sor(a,b,1.45,[1 11]',1e-6)func
i)+omega*x1(i);enddisp(x')endif k==n,warning('已达到迭
i)+omega*x1(i);enddisp(x')endif k==n,warning(' yi da dao die
-u(i,网通传奇sfi+1:n)*x0(i+1:n,1))/a(i,i);x(i)=(1-omega)*x0(
-u(i, wang tong chuan qi sfi+1:n)*x0(i+1:n,1))/a(i,i);x(i)=(1-omega)*x0(
k+1x0=x;for i=1:nx1(i)=(b(i)-l(i,1:i-1)*x(1:i-1,1)
k+1x0=x;for i=1:nx1(i)=(b(i)-l(i,1:i-1)*x(1:i-1,1)
-1);u=triu(a,1);whilenorm(x0-x,inf)>e&k<n,k=
-1);u=triu(a,1);whilenorm(x0-x,inf)>e&k<n,k=
f nargin<6,n=500;endx=x0;x0=x+2*e;k=0;网通传奇sfl=tril(a,
f nargin<6,n=500;endx=x0;x0=x+2*e;k=0; wang tong chuan qi sfl=tril(a,
ega,x0,e,n)% 用途：用分量形式的sor迭代解线性方程组ax=bn=length(b);i
ega,x0,e,n)% yong tu yong fen liang xing shi de sor die dai jie xian xing fang cheng zu ax=bn=length(b);i
=gaussseidel(a,b,[1 11]网通传奇sf',1e-6)functionx=sor(a,b,om
=gaussseidel(a,b,[1 11] wang tong chuan qi sf',1e-6)functionx=sor(a,b,om
p(x')endif k==n,warning('已达到迭代次数上限');end>> x
p(x')endif k==n,warning(' yi da dao die dai ci shu shang xian ');end>> x
f)>e&k<n,k=k+1;x0=x;x=-a2*(a-a网通传奇sf1)*x0+a2*b;dis
f)>e&k<n,k=k+1;x0=x;x=-a2*(a-a wang tong chuan qi sf1)*x0+a2*b;dis
=x+2*e;k=0;a1=tril(a);a2=inv(a1);whilenorm(x0-x,in
=x+2*e;k=0;a1=tril(a);a2=inv(a1);whilenorm(x0-x,in
-seidel迭代解线性方程组ax=bif nargin<5,n=500;e网通传奇sfndx=x0;x0
-seidel die dai jie xian xing fang cheng zu ax=bif nargin<5,n=500;e wang tong chuan qi sfndx=x0;x0
function x=gaussseidel(a,b,x0,e,n)% 用途：用向量形式的gauss
function x=gaussseidel(a,b,x0,e,n)% yong tu yong xiang liang xing shi de gauss