數(shù)值線性代數(shù)課程設(shè)計(jì)資料

1、數(shù)值線性代數(shù)課程設(shè)計(jì)報(bào)告（2014-2015 第二學(xué)期）姓名:王美玲學(xué)號(hào):081310104任課教師:楊熙南京航空航天大學(xué)2015年6月18Jacobi 迭代法)Gauss-Seidel 迭代法)SOR 迭代法求解線性方程組的數(shù)值效果比較摘要：Jacobi迭代法，Gauss-Seide迭代法，SOR迭代法是三種經(jīng)典的用于求解線性方程組的迭代方法，本文主要對(duì)這三種方法的數(shù)值逼近效果進(jìn)行比較。關(guān)鍵詞：Jacobi迭代法；Gauss-Seidel迭代法；SOR迭代法；線性方程組線性方程組的求解方法可歸納為直接法和迭代法。 迭代法中有三種最為經(jīng)典的迭代方法，就是Jacobi迭代法，Gauss-Seid

2、el迭代法和SOR迭代法。然而三種方法的收斂性，近似解的逼近效果有不同。本文將對(duì)三種方法求解線性方程組的迭代效果做相應(yīng)探討。設(shè)有線性方程組Ax=b,A為非奇異矩陣，求x的近似解，三種迭代方法如下。1. Jacobi迭代法算法：(1)選取初始點(diǎn)x(0)化迭代次數(shù)k=0。,精度要求ep,最大迭代次數(shù)N,初始(3)相對(duì)誤差err小于等于精度要求ep時(shí)，輸出(k+1)作為方程的近似解(4)x(k)=x(k+i),k=k+1,轉(zhuǎn)步驟(2)。2. Gauss-Seidel迭代法算法：(1)選取初始點(diǎn)x(0),精度要求ep,最大迭代次數(shù)N,初x始化迭代次數(shù)k=0。(2)由Gauss-Seide迭代法計(jì)算公式

3、計(jì)算點(diǎn)x(k+1)。Gau-Seidel迭代法計(jì)算公式為：7。-值7(初始向期=b工a吟川工的甘施j=1j=a-l(13)(3)相對(duì)誤差err小于等于精度要求ep時(shí)，輸出丫(k+1)作為方x程的近似解。(4)x(k)=x(k+1),k=k+1,轉(zhuǎn)步驟(2)。3. SOR迭代法算法：(2)由Jacobi迭代法計(jì)算公式計(jì)算點(diǎn)Jacobi迭代法的計(jì)算公式為：=/v網(wǎng)TX(Xin)ni2ayxj)/an;j百(k+1)x(L2)(1)選取初始點(diǎn)x(0),精度要求ep,最大迭代次數(shù)N,初始化迭代次數(shù)k=0。(2)由SOR迭代法計(jì)算公式計(jì)算點(diǎn)x(k+1)。(14)(3)相對(duì)誤差err小于等于精度要求ep時(shí)

4、，輸出(k+i)作為方x程的近似解。(4)x(k)=x(k+i),k=k+1,轉(zhuǎn)步驟(2)。上述三種經(jīng)典迭代法收斂的充分必要條件是迭代矩陣譜半徑小于1。譜半徑不易求解，而在一定條件下，通過(guò)系數(shù)矩陣A的性質(zhì)可判斷迭代法的收斂性。定理1:若系數(shù)矩陣A是嚴(yán)格對(duì)角占優(yōu)或不可約對(duì)角占優(yōu)，則Jacobi迭代法和Gauss-Seidel迭代法均收斂。定理2：(1)SOR迭代法收斂的必要條件是0w2;(2)若系數(shù)矩陣A嚴(yán)格對(duì)角占優(yōu)或不可約對(duì)角占優(yōu)且0wGauss-SeidelSORb=第29卷第6期,2008年11月2徐樹(shù)方.數(shù)值線i代數(shù).北京:北京大學(xué)出版社,1995.3馬昌鳳.現(xiàn)代數(shù)值分析.北京:國(guó)防工業(yè)

5、出版社.2013.4劉春鳳，米翠蘭.實(shí)用數(shù)值分析教程.北京冶金工業(yè)出版社.2006附錄：源代碼I.Jacob迭代法functionx,k=myjacobi(A,b,x,ep,N)%雅可比迭代法解線性方程組Ax=b%A為系數(shù)矩陣，b為右端向量，x為初始向量(默認(rèn)為零向量)，ep為精度%N為最大迭代次數(shù)(默認(rèn)最大值為500),x為近似解，k為迭代次數(shù)clearkx;ifnargin5N=500;Endifnargin4ep=1e-6;endifnargin3x=zeros(size(b);endD=diag(diag(A);fork=1:Nx=D(D-A)*x+b);err=norm(b-A*x)

6、/norm(b);iferrepbreakendendend2.Gauss-Seidellt代法functionx,k=myseidel(A,b,x,ep,N)%高斯-賽德?tīng)柕ń饩€性方程組Ax=b%A為系數(shù)矩陣，b為右端向量，x為初始向量(默認(rèn)為零向量)，ep為精度%N為最大迭代次數(shù)(默認(rèn)最大值為500),x為近似解，k為迭代次數(shù)clearkx;ifnargin5N=500;endifnargin4ep=1e-6;endifnargin3x=zeros(size(b);endD=diag(diag(A);L=D-tril(A);U=D-triu(A);fork=1:Nx=(D-L)(U*x+b);err=norm(b-A*x)/norm(b);iferrepbreak;endend3.SORt代法functionx,k=mysordd(A,b,w,x,ep,N)%超松弛迭代法(SOR)解線性方程組Ax=b%w為松弛因子%A為系數(shù)矩陣，b為右端向量，x為初始向量(默認(rèn)為零向量)，ep為精度%N為最大迭代次數(shù)(默認(rèn)最大值為500),x為近似解，k為迭代次數(shù)clearkx;ifnargin6N=500;endifnargin5ep=1e-6;endifnargin4x=zeros(size(b);endifnargin3w=0.8;endD=diag(dia