论文部分内容阅读
在数值计算领域中,对线性方程组进行求解是十分活跃的研究课题,一般可分为两种情况,一是方程是良态的,即方程的解存在唯一且稳定,可采用常规的算法,例如直接解法和迭代算法。另一种为方程是超定的或不定的,也就是病态线性方程组,寻找病态线性方程组有效求解方法一直以来是数值代数领域研究的热点。本文主要以病态线性方程组的求解为重点,研究基于Krylov子空间的投影算法解决此类问题的可行性和有效性,并且将研究的成果运用于一类热传导方程反问题的求解当中。本文所开展的主要研究工作如下: (1)对广义极小残余(GMRES)方法、RRGMRES方法、限制值域的广义极小残量(RRGMRES)算法、LSQR方法、共轭梯度法等进行了详细的描述。通过分析表明基于Krylov子空间的投影算法能够有效求解病态线性方程组,避免了所求数值解的不稳定。 (2)热传导方程反问题是一种高度不适定问题,将其化为积分方程,然后采用数值积分离散,最终转化为病态的线性方程组,采用本文方法对病态的线性方程进行求解,数值模拟结果表明了这些方法是可行和有效的。 (3)本文尝试将基于krylov子空间的投影算法应用到热传导方程反问题中,不仅可以解出线性方程组的数值解,并且数值模拟效果不错。