Jacobi迭代法分量形式matlab,有无穷解的线性方程组的迭代法_田学全

第16卷　第2期塔　里　木　农　垦　大　学　学　报Vol.16No.2

2004年6月JournalofTarimUniversityofAgriculturalReclamationJun.2004

①　　文章编号:1009-0568(2004)02-0055-02

有无穷解的线性方程组的迭代法

田学全1　朱世建2

(1　塔里木农垦大学文理学院,新疆阿拉尔　843300)

(2　农十师北屯181团中学,新疆阿勒泰　836001)

随着计算机技术的发展,计算机的存储量日益增大,计算机速度也迅速提高,直接法(如GAUSS消去法,平方根法等)在计算机上可以求解的线性方程组的规模也越来越大,但直接法大多数均需要对系数矩阵A进行分解,因而一般不能保持A的稀疏性,而实际应用中,特别是偏微分方程的数值求解时,常常遇到大型稀疏线性方程组的求解问题,迭代法是能保持线性方程组稀疏性的有效算法,它也是数值

代数中的一种常用的非常重要的方法。本文运用雅可比(Jacobi)迭代法,高斯-塞得尔(Gauss-Seidel)迭代法探讨有无穷解的线性方程组的迭代法。

1　方程组Ax=b的等价方程组:

设方程组Ax=b有无穷解,并且不妨令化简后的方程组有如下形式:

a11x1+a12x2+…+a1mxm=b1-a1,m+1xm+1-…-a1,nxna21x1+a22x2+…+a2mxm=b2-a2,m+1xm+1-…-a2,nxn……

am1x1+am2x2+…+ammxm=bm-am,m+1xm+1-…-am,nxn

其中(n>m),aii≠0,A1=(aij)mxm是非奇异方阵。

n

n

令g=(b1-

j=m+1

<