Jacobi迭代求解线性方程组
对于线性方程组Ax=b,A是nxn矩阵,b是一个nx1向量。若A非奇异,则有A的如下分解
A=D-L-U
其中,D,L,U都是nxn矩阵,A=(aij)nxn,D=diag(a11,…,ann),-L和U分别是A的除了对角线外的上下三角部分。那么原方程变为
(D-L-U)x=b,因为A可逆,所以D可逆,所以有
Dx=b+(L+U)x,
进而
x=D-1b+D-1(L+U)x.
那么就有迭代格式
xk+1=D-1b+D-1(L+U)xk.
令B=D-1(L+U),通过分析B的相应性质,可知道收敛性(有相应定理的)
Jacobi迭代求解线性方程组
最新推荐文章于 2024-05-25 12:48:51 发布