Jacobi迭代法和Gauss-Seidel迭代法Matlab程序(共4页).docx

《Jacobi迭代法和Gauss-Seidel迭代法Matlab程序(共4页).docx》由会员分享，可在线阅读，更多相关《Jacobi迭代法和Gauss-Seidel迭代法Matlab程序(共4页).docx（4页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上分别Jacobi迭代法和Gauss-Seidel迭代法求解线性方程组迭代法计算停止的条件为： 解（1）：采用Jacobi迭代法时，Matlab计算程序为：clearclci=1;a=6 2 1 -2;2 5 0 -2;-2 0 8 5;1 3 2 7;d=diag(diag(a);l=d-tril(a);u=d-triu(a);d0=inv(d);b=4;7;-1;0;x0=zeros(4,1);B=d0*(l+u);f=d0*b;x=B*x0+f;while norm(x-x0,inf)=1e-6 x0=x; x=B*x0+f; i=i+1;endxi计算结果为：x

2、 =（0.0520；1.1509；0.2446；-0.5706） 迭代次数i =15专心-专注-专业采用Gauss-Seidel迭代法计算时，Matlab计算程序为：clearclci=1;a=6 2 1 -2;2 5 0 -2;-2 0 8 5;1 3 2 7;d=diag(diag(a);l=d-tril(a);u=d-triu(a);b=4;7;-1;0;x0=zeros(4,1);B=inv(d-l)*u;f=inv(d-l)*b;x=B*x0+f;while norm(x-x0,inf)=1e-6 x0=x; x=B*x0+f; i=i+1;endxi计算结果为：x =（0.0520；1.1509；0.2446；-0.5706） 迭代次数i =8