在数值代数中考虑线性方程组 Ax=b A x = b ,其中 A∈Rn×n A ∈ R n × n , x∈Rn x ∈ R n 。一般有两种思路:
- 当矩阵 A A 为低阶稠密矩阵时,考虑使用Gauss消元思想进行求解;
- 当 为大型稀疏矩阵时,考虑使用迭代法进行求解。
本文对两类基本的迭代方法进行一些归纳讨论。
一、迭代法的基本概念
迭代法的基本思想是将 Ax=b A x = b 转化为便于迭代的等价方程组: x=Bx+c x = B x + c ,进而构造迭代格式
xk+1=Bxk+c(1) (1) x k + 1 = B x k + c
对于给定的初始向量 x0 x 0 ,由( 1 1 )产生的序列 { xk} { x k } 满足 limk→∞xk=x∗ lim k → ∞ x k = x ∗ 且 Ax∗=b A x ∗ = b ,称迭代格式( 1 1 ) 收敛。
一般情形下迭代格式( 1 1 )收敛的必要条件为:
1. 矩阵