32Jacobi迭代法和Gauss-Seidel迭代法.pptVIP

3.2 Jacobi迭代法和Gauss-Seidel迭代法 本节主要内容 Jacobi迭代法 Gauss-Seidel迭代法 3.2.1 Jacobi迭代法 【解】 k k 0 0 0 0 3 0.9950 0.9850 0.9900 1 0.9000 0.7000 0.8000 4 0.9985 0.9975 0.9970 2 0.9700 0.9500 0.9400 5 0.9998 0.9992 0.9995 卡尔.雅可比 其中A是n阶非奇异矩阵. 且其主对角元素 从而得到Jacobi迭代法的分量形式： 下面推导Jacobi迭代法的矩阵形式： 把系数矩阵A分解成三部分： 任取向量 ，则Jacobi迭代法可写成如下的矩阵形式： 称矩阵J 为Jacobi迭代法的迭代矩阵. 3.2.2 算法与程序 算法3.1 Jacobi迭代法 说明：为简单起见，假定系数矩阵A非奇异，且 ，且假设Jacobi迭代法收敛. 步骤1 输入系数矩阵A，右端向量b，以及初始向量 算法3.1的 Matlab 程序 %Jacobi.m function x = Jacobi (A, b, x0, eps, N) % 功能：用Jacobi迭代法解n 阶线性方程组 Ax=b n=length(b); x=o