Java迭代方程求方程近似解_线性方程组的迭代求解java.docx

线性方程组的迭代求解java.docx

线性方程组的迭代求解线性方程组的迭代求解摘要迭代法是一种逐次逼近方法，在使用迭代法解方程组时，其系数矩阵在计算过程中始终不变。它利用计算机运算速度快、适合做重复性操作的特点，让计算机对一组指令(或一定步骤)进行重复执行。迭代法具有循环的计算方法，方法简单，适宜解大型稀疏矩阵方程组本文总结了解线性方程组的三个迭代法，Jacobi迭代法，Gauss-Seidel迭代法，SOR迭代法，并且介绍了软件JAVA在这方面的应用。关键词：Jacobi迭代法；Gauss-Seidel迭代法；SOR迭代法；计算线性方程组的迭代求解SOLUTIONOFLINEAREQUATIONSOFITERATIONWITHTHEEXPERIMENTALABSTRACTIterationisakindoftosolvequestionsbystep-by-stepapproximation.Whenwearegettingthesolutionoflinearequationsbyusingiteration,thecoefficientmatrixisalwaysstayingthesameincomputationprocess.Computercouldoperatefastlysothatitissuitableforoperatingagainandagain.Iterationiseasytooperatetosolvethelargematrixequationsbyusingacalculatecalledcirculation.Thissummaryunderstandingoflinearequationsthreekindofiteration,Jacobiiteration,Gauss-Seideliteration,successiveoverrelaxation,andintroducemodernsoftwareJAVAinthisrespect.Keywords:Jacobiiteration;Gauss-Seideliteration;SuccessiveOverRelaxation;calculating线性方程组的迭代求解目录1迭代法概述……………………………………………….………………………….11.1迭代法定义……………………………………………………………………………11.2迭代法基本原理………………………………………………………………………12迭代法解线性方程组…………………………………………………………………….12.1雅克比(Jacobi)迭代法……………………………………….………………………12.2高斯—赛德尔(Gauss-Seidel)迭代法………………………………………….42.3超松弛(SOR)迭代法…………………………………………………………….…73总结……………………………………………………………………………….9参考文献………………………………………………………………………….……10附录…………………………………………………………………………….……………11线性方程组的迭代求解第1页共19页1迭代法概述迭代法也称辗转法，是一种逐次逼近方法，在使用迭代法解方程组时，其系数矩阵在计算过程中始终不变。它利用计算机运算速度快、适合做重复性操作的特点，让计算机对一组指令(或一定步骤)进行重复执行，在每次执行这组指令(或步骤)时，都从变量的原值推出它的一个新值。迭代法具有循环的计算方法，方法简单，适宜解大型稀疏矩阵方程组，在用计算机计算时只需存储A的非零元素(或可按一定公式形成系数，这样A就不需要存储)[1]。1.1迭代法定义(1)对于给定的方程组，用式子xBf(1-1)(1)(0)2(1)()kkxfBf逐步代入求近似解的方法称为迭代法(或称为一阶定常迭代法，这里与B和k无关)(2)如果存在(记作)，称此迭代法收敛，显然就是方程组的解，否()limkxxx则称此迭代法发散。1.2迭代法基本定理设有方程组，对于任意初始向量及任意f,解此方程组的迭代法(即xBf(0)x)收敛的充要条件是.(1)()kkxf()1B2迭代法解线性方程组2.1雅克比(Jacobi)迭代法线性方程组的迭代求解第2页共19页2.1.1Jacobi迭代法的定义设有方程组(),n1ijiaxb1,2in记作(2-1)AA为非奇异阵且。将A分裂为,其中a0(1,2)ijnDLU，，12nDa21312,100nnaL。123121,00naUa将式(2-1)第i()个方程用去除再移项，得到等价方程组,2inija(),(2-2)1()niijjxbax1,i简单记作，其中,0Bf110()BIDALU1fDb对方程组(2-2)应用迭代法，得到(2-1)的迭代公式(2-3)(0)()0()12()()1,,Tnkkiiijjxxba线性方程组的迭代求解第3页共19页其中为第k次迭代向量，设已经算出，由式(2-3)可计(0)()0()12,,Tnxx()kx算下一次迭代向量。(1),2;1,kin显然迭代公式(2-3)的矩阵形式为(2-4)(0)1()+fkkixB)初始向量其中称为Jacobi方法迭代矩阵。02.1.2JAVA程序实现Jacobi迭代法编写java程序用Jacobi迭代法解如下方程组：例1：1235+40 x实验结果如下图所示(JAVA程序设计详见附录源程序1)：线性方程组的迭代求解第4页共19页2.2Gauss-Seidel迭代法线性方程组的迭代求解第5页共19页2.2.1高斯—赛德尔(Gauss-Seidel)迭代法的定义雅克比迭代法的优点是公式简单，迭代矩阵容易计算。在每一步迭代时，用的全部分量求出的全部分量，因此称为同步迭代法，计算时需保留两个)(kx)1(kx近似解和。)(k)1(但在雅克比迭代过程中，对已经计算出的信息未能充分利用，即在计算第i个分量时，已经计算出的最新分量没有被利用。从直观上看，)1(kix)1()1(,kikx在收敛的前提下，这些新的分量应比旧的分量更好，更精)(1)1(,ikx)(1)(1,kikx确一些。因此，如果每计算出一个新的分量便立即用它取代对应的旧分量进行迭代，可能收敛的速度更快，并且只需要储存一个近似解向量即可。据此思想可构造高斯—赛德尔(Gauss-Seidel)迭代法，其迭代公式为(i=1,2,…,n)(2-5))(1(1)()(ikjijnijkjikibxaxax也可以写成矩阵形式SGkSkfxBx)()1(仍将系数矩阵A分解为ULD则方程组变为b)(得(2-6)x将最新分量代替为旧分量，得bULDxkkk)(