LSQR算法是Paige和Sanders在1982年提出的最小二乘法和线性系统的方法,适用于解决大型稀疏线性方程组,并且,再相比较其他方法,LSQR算法的求解量较小,能够使用简化来计算或求解大型稀疏线性方程组,其有效性是毋庸置疑的。
LSQR算法的优点是占有用的存储空间少,运算速度较快,特别适用于求解系数为大型稀疏矩阵的方程组,和其他的迭代方程组相比较下,在求解奇异或者病态的问题时,LSQR能够显示其更快的收敛性以及更好地可接受的结果。文献综述
LSQR算法是一种特别适用于求解大型稀疏矩阵的迭代算法,在数学上就相当于共轭梯度法。
第二章 线性方程组
线性方程组即线性系统的求解问题是一个探究很久的数学问题,早在中国古代的《九章算术》中,就已经详尽的描写了线性系统的消元法。而到了19世纪初,西方也有Gauss消去法。但是求解多未知数的线性方程组则是在20世纪中叶,计算机问世以后才得以求解出来。
线性方程组的求解的数值方法有很多种,但是一般可分为直接法和迭代法两种。直接法是在没有舍入误差的情况下,经过了有限次的运算来求得方程组的精确的解法。为此,直接法又可称为精确法。而迭代法则是采用逐次逼近的解法,就是从一个初始向量出发,按照一定量的计算,构造出一个向量的无穷序列,其极限是其方程组的精确解,若是只经过有限次运算是得不到其精确解的,只能获得近似解。
2。1线性方程组的直接解法
高斯消去法是很早以前提出的一种求解线性系统的方法,其最早出现在公元前250年前,数学家就有研究此方法的应用,可见其悠久。高斯消去法的基础是以行初等变化为主,而后通过改进和变化等一系列的进化,变成了计算机算法。其中高斯消去法又有两种分支的出现,即高斯选主元消去法以及三角分解法。
高斯消去法是目前求解中小规模线性方程组(即阶数不能太高,例如:小于等于1000)常用的方法,一般是用在稠密的系数矩阵(矩阵的绝大多数元素为非零)且没有任何特殊结构的线性方程组。若系数矩阵存在特殊形式,为了尽可能减少计算量和储存量,就需要采取特别的方法来求解。
高斯消去法的核心思想是采用初等行变换,改变其系数矩阵,并且化为较为简单的矩阵或者线性系统,也就是将线性方程组进行一些变换来使其简单化,例如,左变换等,化为上三角矩阵,接着转变为简单的线性方程组,然后再进行将公式回代求解,最近进行求解。
高斯消去法的基本步骤:先进行消元计算,然后再回代求解。来.自^优+尔-论,文:网www.youerw.com +QQ752018766-
算法公式:
消元公式:有,如果,取消元因子,
对公式进行回代:如果,
2。2 线性方程组的迭代法
迭代法是按照规则构造一个向量序列,使它的极限向量是方程组的精确解。一个方法的有效性要看得到的近似解的精确程度,对其付出的代价如何。一般以运算量和存储量要求为标准。在这个标准下,直接法在很多方程组中会比迭代法有效,但是对于大型的稀疏方程组,迭代法会更加适用。
2。2。1 雅克比迭代法
1、雅可比迭代法
在阶线性方程组中,是非奇异性的矩阵,。在线性方程组中第个方程可以呈现为可以把它改变成为和它一样的等价值的线性方程组
取任意的初始向量,在上式的计算中我们可以得出一个迭代序列,若是要满足,(为精度),终止计算的话,那么就是满足了的近似解。 用LSQR迭代法求解线性方程组(2):http://www.youerw.com/shuxue/lunwen_112968.html