解線性方程组的迭代法_计算方法大作业.docxVIP

解线性方程组的迭代法1.1 方法概述对于线性方程组Ax=b，其中A为非奇异矩阵，当A为大型稀疏矩阵时，考虑用迭代法求解上述方程组，其基本思想是求不动点,即构造一个向量系列，使其收敛至某个极限，则就是要求的方程组的准确解。这儿主要介绍Jacobi迭代法和Gauss-Seidel迭代法的算法理论及数值实验。1）Jacobi迭代法算法概述Jacobi迭代法推导过程为将方程组：在假设下，改写成如果引用系数矩阵, 及向量 ，，，方程组（1）和（2）分别可写为：及，这样得到了迭代格式用迭代解方程组时，就可任意取初值带入迭代可知式，然后求。但是，比较大的时候，写方程组和是很麻烦的，如果直接由，能直接得到，就是矩阵与向量的运算了，那么如何得到，呢实际上，如果引进非奇异对角矩阵将分解成：要求的解，实质上就有而奇异的，所以存在，从而有我们在这里不妨令就得到迭代格式： 。用向量的分量来表示为：其中为初始向量.由此看出,雅可比迭代法公式简单,每迭代一次只需计算一次矩阵和向量的乘法.在电算时需要两组存储单元,以存放及.2）Gauss-Seidel迭代法算法概述由雅可比迭代公式可知,在迭代的每一步计算过程中是用的全部分量来计算的所有分量,显然在计算第i个分量时,已经计算出的最新分量没有被利用，从直观上看,最新计算出的分量可能比旧的分量要好些.因此，对这些最新计算出来的第次近似的分量加以利用,就得到