可对称化不定非线性问题的不精确Newton法

本文讨论了处理可对称化不定问题的不精确Newton方法,并针对问题的特殊结构提出了不精确Newton-PSMINRES算法。理论分析与数值试验表明,Newton-PSMINRES算法优于其它处理可对称化不定问题的不精确Newton-Krylov算法。

对称不定问题的不精确Newton法

1.引 言 非线性方程组F（x）=0的数值求解,经典的算法是Newton迭代;xk+1=xk+sk,k=0,1,2,…,（1.1）其中的sk满足F’（xk）sk=-F（xk）;k=0,1,2,….（1.2）这里x0为迭代的初始点,{xk}称为Newton迭代序列.当变量个数比较多时,每一步Newton迭代中计算Jacobi矩阵F’（xk）和求解线性方程组（1.2）的代价非常高;特别当xk远离方程组的解x＊时,高精度地求解线性方程组（1.2）