Jacobi迭代法和Gauss-Seidel迭代法Matlab程序

整理文档很辛苦,赏杯茶钱您下走!

免费阅读已结束,点击下载阅读编辑剩下 ...

阅读已结束,您可以下载文档离线阅读编辑

资源描述

.分别Jacobi迭代法和Gauss-Seidel迭代法求解线性方程组12346212425027,2085113270xxxx迭代法计算停止的条件为:6)()1(3110maxkjkjjxx.解(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

1 / 4
下载文档,编辑使用

©2015-2020 m.777doc.com 三七文档.

备案号:鲁ICP备2024069028号-1 客服联系 QQ:2149211541

×
保存成功