论文部分内容阅读
算法的准确性和稳定性是数值代数中重要而基础的研究分支.随着科学的进步,计算机技术的发展,数值计算实际问题的规模不断扩大,产生的矩阵越来越大,那么解线性方程组就变得更加困难,需要提出一些切实可行的算法来解决这样的问题.在提出方法之后就会面临一个一般性的问题,即这个方法可不可行,当轻微的扰动引起后,会不会出现数值解离真实解差距比较大的问题.所以对算法的稳定性和准确性进行研究是相当必要的.本文研究块三对角矩阵,一般非奇异矩阵和非对称严格t-对角占优矩阵的相关算法的稳定性和准确性;探讨几类增长因子;分析极分解的更新和秩r更新的准确性;分别研究鞍点问题和加权线性最小二乘问题的准确性和稳定性.
研究块三对角矩阵相关算法的准确性和稳定性.首先,利用分而治之与URV分解确定块三对角矩阵的逆,并给出产生于计算中的逆的舍入误差.该算法其计算复杂性和误差均比块高斯-若尔当消元法的计算复杂性和误差好.其次,若原矩阵是I-块对角占优(Ⅱ-块对角占优)块三对角矩阵,则约化矩阵保持同样的性质.给出几类比较矩阵的一些性质并分析块三对角矩阵BLU分解的向后误差.再次,利用Amodio和Mazzia提出的一般矩阵LU分解的块表达式给出块三对角矩阵,块三对角块H-矩阵和复对称块三对角矩阵BLU分解的块表达式,研究其准确性和向后误差并体现其优越性.最后,为了避免分解因子不是三角形式的状况,研究块三对角的分块LU分解并分析其准确性和稳定性.
基于泰勒级数的一阶表达式和矩阵-向量方程,研究LU和QR分解因子的扰动理论.对于一般非奇异矩阵的LU分解,考虑Chang和Paige提出的问题并研究完全选主元LU分解的相关问题.另外,因严格t-对角占优可以刻画对角病态,所以考虑非对称严格t-对角占优线性方程组的LU和QR分解的向后误差.
基于Amodio和Mazzia提出的增长因子定义,研究块三对角块H-矩阵和复对称块三对角矩阵在BLU分解下的增长因子.同时,提出广义的Buckley矩阵的定义,并基于Wilkinson提出的增长因子的定义,讨论该矩阵及其一个相关的矩阵在高斯消去下的增长因子.另外,算法的更新和删减在很多方面均有涉及.分析极分解的更新和秩r更新的扰动理论.
研究鞍点问题的扰动理论.首先,探讨当扰动系统仍是鞍点系统时的扰动理论.由于系统解的准确性表明其依靠于鞍点矩阵的子矩阵,为了消除其影响而提出一些方案.为了提高解的准确性,提出一个尺度化.其次,当原系统是鞍点系统而扰动系统是广义鞍点系统时,考虑块LDLT分解的灵敏性并讨论此情形下的扰动理论.再次,由于鞍点问题与最小二乘问题的密切联系,给出一个解拥有不定对角加权矩阵的加权线性最小二乘问题的可证向后稳定算法.然而,当加权矩阵是广义鞍点矩阵时,类似的算法不一定是向后稳定的,所以最后讨论此时算法是可证向后稳定的情形.