线性方程组迭代法x.pptVIP

思路;§4.1 三种基本的迭代法 4.1.1．雅克比（Jacobi）迭代法 （以三阶方程 组为例）;假设;任选一初值向量：;对于n阶方程组;n阶方程的 Jacobi 迭代格式:;若用矩阵来表示雅可比迭代，则如下:;A x=b,  ( D – L – U) x = bDx = (L+U) x+b;4.1.2 高斯—赛德尔迭代法 (Gauss-Seidel);这种迭代称为 高斯—赛德尔(Gauss-Seidel)迭代法（G-S）;n阶方程的G-S迭代格式:;用矩阵表示：;例1：分别Jacobi及G-S迭代法解下列线性方程组。初值均取(0，0，0)T;%用雅可比法解例1 a=[9,-1,-1;-1,8,0;-1,0,9]; D=-(a-triu(a)-tril(a)); % D=diag(diag(a)); L=-(tril(a)-D); U=-(triu(a)-D); xo=[0;0;0];bo=[7;7;8]; ep=0.0001;dx=1;k=0; while dxep k=k+1; x=D\(L+U)*xo+D\bo; dx=abs(norm(x)-norm(xo)); xo=x; end k,x;%用雅可比法解例1 a=[9,-1,-1;-1,8,0;-1,0,9]; bo=[7;7;8]; ep=