求解线性方程组的方法课题.ppt

计 算 方 法;§8.0 引 言; 在自然科学和工程技术中许多问题的解决转化为解线性方程组，而这些方程组的系数矩阵大致分为两种，一种是低阶稠密矩阵，一种是高阶稀疏矩阵。; ;例1 用消去法求解方程组;设有方程组;将原方程组记为;第二步，设经过k-1次消元后的同解方程组为;重复上述过程，可以将方程组化为等价的简单方程组;乘除运算量;定理2 约化的主元素;定理2 约化的主元素;function x=gauss(A,b) n=length(b); x=zeros(n,1); for i=1:n-1 %消元过程 for k=i+1:n for j=i+1:n A(k,j)=A(k,j)+A(i,j)*(-A(k,i)/A(i,i)); end b(k)=b(k)+b(i)*(-A(k,i)/A(i,i)); A(k,i)=0; end end disp(A) disp(b) pause;A=[1 1 1;0 4 -1;2 -2 1]; b=[6 5 1]; x=gauss(A,b);§8.2 雅可比Jacobi迭代法 (AX=b); 与解f (x)=0 的不动点迭代相类似，将AX=b改写为X=BX+f 的形式，建立雅可比方法的迭代格式：