雅克比迭代法的迭代矩阵J是Bx+f中的B,那高斯-赛德尔迭代法的迭代矩阵G,好像和J是一样的?是怎么算出来的?

发布网友 发布时间：2022-04-24 03:55

我来回答

共1个回答

热心网友 时间：2023-10-26 08:39

迭代法写成x<-Bx+f之后迭代矩阵当然是B
不过问题在于不同的迭代法产生的B和f是不同的

在Jacobi迭代中A=D-L-U，Ax=b <=> Dx=(L+U)x+b <=> x=D^{-1}(L+U)x+D^{-1}b
所以B=D^{-1}(L+U)

在Gauss-Seidel迭代中同样A=D-L-U，但是Ax=b <=> (D-L)x=Ux+b <=> x=(D-L)^{-1}Ux+(D-L)^{-1}b
所以B=(D-L)^{-1}U

热心网友 时间：2023-10-26 08:39

迭代法写成x<-Bx+f之后迭代矩阵当然是B
不过问题在于不同的迭代法产生的B和f是不同的

在Jacobi迭代中A=D-L-U，Ax=b <=> Dx=(L+U)x+b <=> x=D^{-1}(L+U)x+D^{-1}b
所以B=D^{-1}(L+U)

在Gauss-Seidel迭代中同样A=D-L-U，但是Ax=b <=> (D-L)x=Ux+b <=> x=(D-L)^{-1}Ux+(D-L)^{-1}b
所以B=(D-L)^{-1}U

热心网友 时间：2023-10-26 08:39

迭代法写成x<-Bx+f之后迭代矩阵当然是B
不过问题在于不同的迭代法产生的B和f是不同的

在Jacobi迭代中A=D-L-U，Ax=b <=> Dx=(L+U)x+b <=> x=D^{-1}(L+U)x+D^{-1}b
所以B=D^{-1}(L+U)

在Gauss-Seidel迭代中同样A=D-L-U，但是Ax=b <=> (D-L)x=Ux+b <=> x=(D-L)^{-1}Ux+(D-L)^{-1}b
所以B=(D-L)^{-1}U

热心网友 时间：2023-10-26 08:39

迭代法写成x<-Bx+f之后迭代矩阵当然是B
不过问题在于不同的迭代法产生的B和f是不同的

在Jacobi迭代中A=D-L-U，Ax=b <=> Dx=(L+U)x+b <=> x=D^{-1}(L+U)x+D^{-1}b
所以B=D^{-1}(L+U)

在Gauss-Seidel迭代中同样A=D-L-U，但是Ax=b <=> (D-L)x=Ux+b <=> x=(D-L)^{-1}Ux+(D-L)^{-1}b
所以B=(D-L)^{-1}U