迭代法解线性方程组数值分析实验报告.docVIP

修I劝械京呼協 数学与计算科学学院 《数值分析》课程设计 题 目: 专 业: 学 号: 姓 名: 指导教师: 迭代法解线性方程组 信息与计算科学 1309302-24 成 绩： 二零一六年六月二十日 —、前言：（目的和意义） 实验目的 掌握用迭代法求解线性方程组的基本思想和步骤。 了解雅可比迭代法，高斯-赛德尔法和松弛法在求解方程组过程中的优缺点。 实验意义 迭代法是用某种极限过程去逐步逼近线性方程组精确解的方法，它是解高阶稀疏方 程组的重要方法。迭代法的基本思想是用逐次逼近的方法求解线性方程组。比较雅可比迭代 法，高斯-赛德尔迭代方法和松弛法，举例子说明每种方法的试用范I韦I和优缺点并进行比较。 二、数学原理： 设有方程组 …①…② …① …② …③ 将其转化为等价的，便于迭代的形式 x= Bx^ f （这种转化总能实现，如令B = I-AJ = h\ 并由此构造迭代公式 严）=加）+ / 式中B称为迭代矩阵，f称为迭代向量。对任意的初始向量x（0）,由式③可求得 向量序列｛兀⑹斤，若则/就是方程①或方程②的解。此时迭代公 式②是收敛的，否则称为发散的。构造的迭代公式③是否收敛，取决于迭代矩阵 B的性 1 ?雅可比迭代法基本原理 …① …① X aijxi 二 bj (i = 1,2,3,…,川) 矩阵形式为Ax = bf设系数矩阵A为非奇异矩阵，且山H0,（/二1,2,3,???,力 从式①中笫i个方程中解出x,得其等价形式 ] n 无=—0-工…② 取初始向量+°)=(聲),雋°),???,对°)),对式②应用迭代法，可建立相应的迭代公式: 1 “宀訂―,盘屛*) 1 “ 宀訂―,盘屛*) …③ 也可记为矩阵形式:严)二矿+ Fj…④若将系数矩阵A 也可记为矩阵形式: 严)二矿+ Fj …④ 若将系数矩阵A分解为A二D-L-U, A=D—L—U= a\n a2n 式屮0… 0 ■ 0… 00a” 式屮 0 … 0 ■ 0 … 0 0 a” ? ? ? ? ? ? ~a2l ° ■ ■ ■ ■ ■ ■ ???0 ? ? ? ? ? ? ? ? ? ? ? ? 0 ? ? ? ° %. __ an\ … 一 ann-\ an\ 0 — 6f|2 0 勺2 -a\n ~an-\n 0 %1 0丿 ro an 0 则方程Ax 则方程Ax二b变为 (D-L-U)x = b Dx = (L + U)x^b = D-\D-A)x-^Dib = (/-￡_1A)x+D_1Z? 于是式中④中的 Bj=I — DTA,fj=D-b° 式③和式④分别称为雅克比迭代法的分量形式和矩阵形式，分量形式用于编程计算，矩阵型 式用于讨论迭代法的收敛性。 2.高斯一赛德尔迭代法 高斯一赛德尔（Gauss-Seidel）迭代法，其迭代公式为 也可以写成矩阵形式 仍将系数矩阵A分解为则方程组变为得将最新分量代替为旧分量,A=D-L-U 仍将系数矩阵A分解为 则方程组变为 得 将最新分量代替为旧分量, A=D-L-U (D- L-U)x = b Dx = Lx+ Ux + b 得 D严）严）+S⑹+b (D-L)x(k+l} =Ux(k)+b 于是有严)=(D-L)~[Ux(k) + (D-LYlb 于是有 Bg_s=(D — u fas =—L「b 超松弛迭代法 设已知第k次迭代向量兀⑷，及第k+1次迭代向量的前i-1个分量￡Z）,（j=l,2,???iT）, 现在研究如何求向Mxu+,）的第i个分量#3）。 首先，有高斯一赛徳尔迭代法求出一个值，记为 再将第k次迭代向量的第i个分量X,）与列⑷）进行加权平均, 得兀:⑷），即： 君7=(1-劲护+应严) = h)+秋#Z-聆) 于是的SOR迭代公式 TOC \o 1-5 \h \z zn T xj1+l)=兀丫)■! (? -￡吗兀严-￡坷兀丫)) (i=l, 2,???n) …① aii 戶 1 7=1 ^ 或 X.) /-I ” 才+】)=(1 _血)才)+ — (/?. ■-工知对⑷—￡知兀畀)(i=h 2, ???n) …② aii j=\ j=i+\ 当e=i吋，式①即为高斯一赛徳尔迭代法； 当06yl时，式①称为低松弛方法，当某些方程组用高斯一赛德尔迭代法不收敛时, 可以用低松弛方法获得收敛； 当。1时，式①称为超松弛方法，可以用来提髙收敛速度。 将式②写成矩阵的形式，得： DX a+1) =(1- a))DX 也 + a)(b + LX a+1) + UX ⑷) 即 (D — coL)x(k+}) = [(1-(6)D + coU]x(k) + ebb 于是得SOR迭代的矩阵表示 x{k+l) = Ba)X(k) fco 式屮 Bco={D-泌尸[(1 一効D + qU] fa)=