第6章解线性方程组的迭代法

所以，序列收敛
Gk 0
与初值的选取无关
数 学 系 University of Science and Technology of China
DEPARTMENT OF MATHEMATICS
定义6.1：（收敛矩阵） 定理：
Gk 0
矩阵G为收敛矩阵，当且仅当G的谱半径<1
G k பைடு நூலகம் (G) 1
数 学 系 University of Science and Technology of China
DEPARTMENT OF MATHEMATICS
第6章 解线性方程组的迭代法
直接法得到的解是理论上准确的，但是我们可以看得出，它们的计算量都是n3 数量级，存储量为n2量级，这在n比较小的时候还比较合适（n<400），但是对于现 在的很多实际问题，往往要我们求解很大的n的矩阵，而且这些矩阵往往是系数矩阵 就是这些矩阵含有大量的0元素。对于这类的矩阵，在用直接法时就会耗费大量的时 间和存储单元。因此我们有必要引入一类新的方法：迭代法。 迭代法具有的特点是速度快。与非线性方程的迭代方法一样，需要我们构造一 个等价的方程，从而构造一个收敛序列，序列的极限值就是方程组的根
Jacobi迭代算法 1、输入系数矩阵A和向量b，和误差控制eps 2、x1={0,0,…..,0} , x2={1,1,…..,1} //赋初值 3、while( ||A*x2-b||>eps) { x1=x2; for(i=0;i<n;i++) { x2[i]=0; for(j=0;j<i;j++) { x2[i] += A[i][j]*x1[j] } for(j=i+1;j<n;j++) { x2[i] += A[i][j]*x1[j] } x2[i]=-(x2[i]-b[i])/A[i][i] } } 4、输出解x2
则，我们可以构造序列 若
x( k 1) G x( k ) g
x ( k ) x * x* G x * g Ax* b
同时：
x( k 1) x* Gx( k ) Gx* G( x( k ) x*) G k 1 ( x(0) x*)
数 学 系 University of Science and Technology of China
DEPARTMENT OF MATHEMATICS
( k 1) x1 x2 ( k 1) ( k 1) xn
格式很简单：
1 (k ) (k ) (a12 x2 a1n xn b1 ) a11 1 (k ) (k ) (k ) (a21 x1 a23 x3 a1n xn b2 ) a22 1 (k ) (k ) (an1 x1 an n 1 xn 1 bn ) ann
x1 x2 xn 1 (a12 x2 a1n xn b1 ) a11 1 (a21 x1 a23 x3 a1n xn b2 ) a22 1 (an1 x1 an n 1 xn 1 bn ) ann
DEPARTMENT OF MATHEMATICS
• 收敛条件 迭代格式收敛的充要条件是G的谱半径<1。对于Jacobi迭代，我们有一些保证收敛 的充分条件
定理：若A满足下列条件之一，则Jacobi迭代收敛。
a1n 0 a12 0 U 0 an 1n 0 0
数 学 系 University of Science and Technology of China
DEPARTMENT OF MATHEMATICS
易知，Jacobi迭代有
数 学 系 University of Science and Technology of China
DEPARTMENT OF MATHEMATICS
• 迭代矩阵
记
A D L U
0 a11 D 0 ann
0 0 a21 0 L 0 a a 0 nn 1 n1
xi
( k 1)
n 1 i 1 (k ) (k ) ( aij x j aij x j bi ) aii j 1 j i 1
数 学 系 University of Science and Technology of China
DEPARTMENT OF MATHEMATICS
知，若有某种范数
由
(G) G
G
p
1
则，迭代收敛
数 学 系 University of Science and Technology of China
DEPARTMENT OF MATHEMATICS
6.1 Jacobi迭代
a11 x1 a1n xn b1 a x a x b nn n n n1 1
(D L U ) x b Dx ( L U ) x b
x D1 ( L U ) x D1b G D1 ( L U ) I D1 A , g D1b
数 学 系 University of Science and Technology of China
数 学 系 University of Science and Technology of China
DEPARTMENT OF MATHEMATICS
对方程组 如：令
Ax b
AM N
做等价变换 ，则
x Gx g
Ax b (M N ) x b Mx b Nx x M 1 Nx M 1b