一种新的适用于JacobianFree+NewtonGMRESm方法的自适应预处理算法.pdfVIP

一种新的适用于Jaeobian-Free Newton-GMRES(m)方法的自适应预处理算法 一种新的适用于Jacobian-Free 自适应预处理算法 沈沉陈颖梅生伟 (清华大学电机系，邮编100084) 摘 要 本文介绍了一种新的适用于Jaeobian-free 测试结果表明，采用这种新的预处理算法后，Newton—GMRES(m)方法的收敛性得到了极大的提高。将新的方 应用于求解电力系统的潮流方程。通过与常规Newton法对比，可以明显地看到新算法在收敛特性方面的优越性。 关键词：Ncwton-GMRES(m)：预处理；潮流 1引言 Jacobian-free 微分方程和常微分方程的数值计算方法。这类方法综合了Newton算法和GMRES算法的特点， 具备超线性的收敛特性【ll丌。由于采用了有限差分技术近似Jaeobian矩阵与向量的乘积，在计算 是生成开销非常大的情况下具有重要的实用价值。特别是在分布式计算和并行计算中，Jacobin．free 的方法更是得到了大量应用。 为了加快JFNG(m)方法的收敛速度，方法之一是降低GMRES方法的叠代次数。为了达到这 一目的，可以采用预处理技术。目前，各种预处理技术得到了广泛的研究和应用【2-181。很多预处 进行ILU分解I闭，或者是基于待求非线性方程组本身的一些具体特征∞】获得预处理矩阵。另外 的一些工作则集中在利用Amoldi过程中获得的信息构建自适应预处理矩阵，这其中包括增维 Krylov子空间的一类方法t9-141和自适应更新预处理矩阵的一类方法…61。以上的这些方法都具有一 个共同点，即都需要在预处理过程中预估雅科比矩阵的特征值或特征向量，从而增加了算法的复 杂性和计算开销。 ． 本文针对JFNG(m)方法提出了一种新的自适应预处理算法。这个算法的特点是预处理矩阵的 形成不需要预估雅科比矩阵或是其特征值，而是利用Alnoldi过程中生成的各Krylov子空间信息 不断修正预处理矩阵。这些投影间的内在联系保证了预处理矩阵始终是雅科比矩阵的一个很好的 二 近似。 具体的做法是，参照文献【4】提出的利用低阶Broyden修正方法在牛顿叠代过程的每一步更新 预处理矩阵的思路，在。眦S(m)叠代过程的每一步更新相应的预处理矩阵。采用的也是Broyden 修正方法，只是修正依据来源于Alnoldi过程中生成的Krylov子空间信息。 为了叙述方便，我们将在第二二节简要介绍标准的Newton．GMRES(m)右预处理方法，秩1修 正方法(Rank-one update)，以及文献141的主要工作。在本文的第三部分详细介绍本文提出的 第十届全国电工数学学术年会论文集 GMRES(m)叠代过程中的自适应预处理算法。第四部分中我们将把提出的方法应用予电力系统 潮流问题的求解并给出测试结果。第五部分为结论。 2基本算法 利用JFNG(m)方法求解形如f(x)=o，x∈R“f∈R”的非线性方程的一般算法如下： (1)令七=-1，选择叠代初始点‰。重复以下过程直至l厂GtⅪ2errtol|Ⅳ (2)k=k+1 (3)计算残差R^=f(x女) (4)进入GMREs(m)叠代求解修正方程：f’GI逾^=一凡=一厂(x七) a)ro=一尺。，z=1，P=11,oll：，V，=roerrtolG=@oll：0 b) 如果PerrtolG并且，m，计算，=，+1 v川=f’X々， ～f(x。+nⅣ，)一厂G。) ∥ 正交化巧+l=Ivl，V2，…，V『+IJ J 求解Hessenberg阵H=(hg 生成并应用GiVe瞄旋转p=0(g)，+．0 c)令‘．』