数值分析雅克比高斯迭代法MATLAB程序

1.雅可比迭代法：function[x,n]=jaccbi(A,b,x0,eps,t)ifnargin==3;eps=1e-6;m=200;elseifnargin3error('输入的数有误');return;elseifnargin==5m=t;endD=diag(diag(A));L=-tril(A,-1);U=-triu(A,1);B=D\(L+U);f=D\b;x=B*x0+f;n=1;whilenorm(x-x0)=epsx0=x;x=B*x0+f;n=n+1;if(n=m)disp('可能不收敛');return;end;end2.高斯—赛德尔迭代法：function[x,n]=gsdddy(A,b,x0,eps,t)ifnargin==3;eps=1e-6;m=200;elseifnargin3error('输入有误');return;elseifnargin==5m=t;endD=diag(diag(A));L=-tril(A,-1);U=-triu(A,1);B=(D-L)\U;f=(D-L)\b;x=B*x0+f;n=1;whilenorm(x-x0)=epsx0=x;x=B*x0+f;n=n+1;if(n=m)disp('迭代次数过多，可能不收敛');return;end;endPS:有兴趣小朋友试一下，结果与x0无关