MATLAB样例之雅克比迭代法.docx

_要求：下面分别使用雅克比迭代法和高斯-赛德尔迭代法求一个方程组的近似解用的线性方程组是按实验要求给的：7*x1+x2+2*x3=10x1+8*x2+2*x3=82*x1+2*x2+9*x3=6雅克比迭代法的matlab代码：（老师写的）A=7,1,2;1,8,2;2,2,9;b=10;8;6;if(any(diag(A)=0)error(error,pause)endeps=input(误差限eps=);N=input(迭代次数N=);D=diag(diag(A);B=inv(D)*(D-A);f=inv(D)*b;K=0;x0=zeros(size(b);while1x1=B*x0+fK=K+1;fprintf(第-次迭代的近似解为,K)disp(x1);ifnorm(x1-x0,inf)Nfprintf(迭代超限)endx0=x1;end高斯-赛德尔迭代法matlab代码:（自己改的）A=7,1,2;1,8,2;2,2,9;b=10;8;6;if(all(diag(A)=0)error(error,pause)endeps=input(误差限eps=);N=input(迭代次数N=);D=diag(diag(A);B=inv(D)*(D-A);f=inv(D)*b;K=0;x0=zeros(size(b);x00=x0;while1x11=B*x0+f;x00(1,1)=x11(1,1);x12=B*x00+f;x00(2,1)=x12(2,1);x13=B*x00+f;x00(3,1)=x13(3,1);x1=x00K=K+1;fprintf(第-次迭代的近似解为,K)disp(x1);ifnorm(x1-x0,inf)Nfprintf(迭代超限)endx0=x1;endWelcome ToDownload !欢迎您的下载，资料仅供参考！精品资料