第2章 解方程组(二).pdfVIP

§2 解线性方程组的迭代法 问题1：需要求解一系列方程： AA xx bb ，，AA xx bb ，，，，AA xx bb ，， 1 1 2 2 m m 其中A , A , …, A 之间和b ,b , …,b 之间相 1 2 m 1 2 m 差很小，那么怎么求解能节省计算量？ 问题2：大规模稀释矩阵计算，例如求解Ax b , 其中 ⎡⎡ 33 −11 1/1/ 22⎤⎤ ⎢ ⎥ ⎢−1 ⎥ ⎢⎢ 33 −−11 1/1/ 22 ⎥⎥ ⎢ ⎥ A ⎢ −1 3 −1 ⎥ ⎢⎢ −11 33 −11 ⎥⎥ ⎢ ⎥ ⎢ 1/ 2 −1 3 ⎥ ⎢⎢ −1 ⎥⎥ ⎢ ⎥ ⎢⎣1/ 2 −1 3 ⎦⎥ §2.1 Jacobi迭代法 设设AA ((aij ))n×n 为为非非奇奇异异矩矩阵阵且且aii =≠≠0(0(ii 11, 22,..., n)), 现需要求解Ax b。 令D表示A的主对角部分，L表示A的下三角部分，U 表示 AA的的上上三三角角部部分分，则则AA LL =+=+DD ++UU ,因因此此 Ax b =⇔ D +L +U x b =⇔Dx b =− L +U x ( ) ( ) ⇔x D −1 (b =−(L +U )x )，这就是矩阵的不动点方程形式。 Jacobi方法： x (0)(0) 初初始始向向量量， xx (k +1) D −1 (b(b −((L +UU))xx(k ) ),), kk 0,0, 1,, 2,, 。。 (k ) * 向量序列x 收敛时，一定收敛到方程Ax