第三章 线性方程组的迭代解法.pptVIP

第1页，共23页，星期日，2025年，2月5日任取一个向量x(0)作为x的近似解，用迭代公式则有x*=Bx*+f，即x*是(3.2)的解，当然x*也就是Ax=b的解.1如何构造迭代公式x(k+1)=Bx(k)+f.这样的构造形式不止一种，它们各对应一种迭代法.2迭代法产生的向量序列{x(k)}的收敛条件是什么，收敛速度如何.结束*从以上的讨论中，可以看出，迭代法的关键有：x(k+1)=Bx(k)+f,(k=0,1,2,?)(3.3)产生一个向量序列{x(k)}，若第2页，共23页，星期日，2025年，2月5日先看一个算例：按下式进行迭代(k=0,1,2,?)结束*3.2几种常用的迭代法公式例13.2.1Jacobi迭代法从以上三个方程中分别解出x1,x2,x3第3页，共23页，星期日，2025年，2月5日任取一初始向量，例如x(0)=(0,0,0)T，得到迭代序列{x(k)}(k=0,1,2,?)，列表如下表3-1。容易验证，原方程组的精确解为x=(1,2,3)T，从上面的计算可看出，{x(k)}收敛于精确解.一般说来，对方程组：设aii≠0(i=1,2,?,n)，从第i个方程解出xi，得等价的方程组：k012345678x100.30000.80000.91800.97160.98040.99620.99850.9998x201.50001.76001.92601.97001.98971.99611.99861.9998x302.00002.66002.95402.95402.98232.99382.99772.9997结束*第4页，共23页，星期日，2025年，2月5日迭代公式为：i=1,2,?,n(3.5)(3.6)这种迭代形式称为Jacobi迭代法(雅可比迭代法)，也称简单迭代法.Jacobi迭代法的矩阵迭代形式:(推导)结束*其中第5页，共23页，星期日，2025年，2月5日3.2.2Gauss-Seidel迭代法在Jacobi迭代法的迭代形式(3.6)中，可以看出，在计算时，要使用.但此时已计算出来.看来此时可提前使用代替，一般地，计算(n≥i≥2)时，使用代替(ip≥1)，这样可能收敛会快一些，这就形成一种新的迭代法，称为Gauss-Seidel迭代法。例2用Gauss-Seidel迭代法计算例1并作比较.(k=0,1,2,?)结束*解迭代公式为：第6页，共23页，星期日，2025年，2月5日用它计算得到的序列{x(k)}列表如表3-2：可见它对这一方程组比Jacobi迭代法收敛快一些。Gauss-Seidel迭代法的公式如下：(3.9)Gauss-Seidel迭代法的矩阵迭代形式：（推导）k0123456x100.30000.88040.98430.99780.99971.0000x201.50001.94481.99221.99891.99982.0000x302.68402.95392.99382.99912.99993.0000结束*其中第7页，共23页，星期日，2025年，2月5日SOR迭代法也称为逐次超松弛法，它可看成Gauss-Seidel迭代法的加速，Gauss-Seidel迭代法是SOR迭代法的特例.改写为结束*3.2.3SOR迭代法将Gauss-Seidel迭代法的(3.9)式记则第8页，共23页，星期日，2025年，2月5日为加快收敛，在增量前加一个因子,得称此公式为SOR法(逐次超松弛法).从(3.13)式不难推