您好,欢迎访问三七文档
当前位置:首页 > 电子/通信 > 综合/其它 > Jacobi迭代法和Gauss-Seidel迭代法Matlab程序
.分别Jacobi迭代法和Gauss-Seidel迭代法求解线性方程组12346212425027,2085113270xxxx迭代法计算停止的条件为:6)()1(3110maxkjkjjxx.解(1):采用Jacobi迭代法时,Matlab计算程序为:clearclci=1;a=[621-2;250-2;-2085;1327];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;whilenorm(x-x0,inf)=1e-6x0=x;x=B*x0+f;i=i+1;endxi计算结果为:x=(;;;)迭代次数i=15采用Gauss-Seidel迭代法计算时,Matlab计算程序为:clearclci=1;a=[621-2;250-2;-2085;1327];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;whilenorm(x-x0,inf)=1e-6x0=x;x=B*x0+f;i=i+1;endxi计算结果为:x=(0.0520;1.1509;0.2446;-0.5706)迭代次数i=8
本文标题:Jacobi迭代法和Gauss-Seidel迭代法Matlab程序
链接地址:https://www.777doc.com/doc-4872037 .html