Hermitian半正定线性系统中定长迭代法收敛的必要条件总结.docx

Hermitian半正定线性系统中定长迭代法收敛的必要条件 Andreas Frommer, Daniel B. Szyld 关键词 线性系统 奇异线性系统 定长迭代法 半模 收敛性分析 摘要 在之前的文章中我们得到了定常迭代收敛的能量半模的充分条件，且该迭代是对于系数矩阵是Hermitian的且半正定的线性系统的。在这篇文章中，我们将说明在什么情况下这些条件依旧必要，以及在什么情况下不然。 介绍 我们考虑一个（奇异）线性系统 （1） 假设其中系数矩阵是Hermitian的且半正定的。 如果是大型稀疏矩阵，迭代思想是解决问题（1）的标准方法。在本文中，我们关注定常迭代法，例如，代数多重网格法，加性和乘性施瓦茨法。有时，通过把它们作为Kryov子空间法的预调节器，这些迭代法是加速的，就像共轭梯度一样。当我们不考虑后者的任何细节，只涉及其中通常假设预调节器是收敛的作为前提条件的情况，所以我们的工作在这种情况中也是相关的。 对于定常迭代法，我们通常考虑用非奇异矩阵分离，得到一个迭代矩阵,并且得到迭代式 . （2） 众所周知，当是非奇异的，对任意起始向量当且仅当时，迭代（2）收敛到一个（1）的唯一解，其中是的谱半径，这样的矩阵被称为收敛的；见例[3,19]。我们之前提到当时非奇异时，给定矩阵和一个收敛矩阵，就存在唯一相应的非奇异矩阵使得[10]。 当是Hermitian的且正定时，在下面的结果[20.Satz 1,p.156]中给出了迭代法（2）收敛的充分条件（也见[7,p.111],[8.p.21]）,现在通常称之为-正则分解定理。如果 是正定的，则的分解（Hermitian正定矩阵）被称为-正则的。 定理1.令是Hermitian正定的，且是一个-正则分解。则 ，其中。 经典的证明使用了Stein的理论（见例，[14, 7.1.8]）。下面我们给出一个本质上相同的证明，但是强调不仅要求，并且要有其中算子范数的给出是通过能量模，其中 ， （3） 表示欧式内积。这个结果很值得说明，因为它表明当能量范数是系统（1）的正则范数时，迭代（2）的误差是单调递减的。进一步，定理1中的“能量范数”也是正则收敛的，也就是说表明这个分解是-正则的。 定理2. 令是Hermitian正定的。则，是一个-正则分解当且仅当。 证明. 可得 则，对任意向量有当且仅当是正定的。 我们提到定理1的其他收敛情况也是可能的。其中一种典型结果就是给定一个-正则分解和Hermitian矩阵，表明是正定的，见例，[9],[14,E71.9]，且进一步，可以得到当不是Hermitian矩阵时的结果，见例[1],[7.p.111]，并被引用其中。 许多作者都给出了半正定矩阵普遍形式下收敛的充分条件；我们将在第二部分看到，其中我们将强调基于下面定义给出的-半范数的条件。在第三部分，我们将回答以下问题：在那种情况中这些充分条件也是必要的。 半定的情况和-半范数 当半定时，（3）式定义了一个半范数，而不是范数。在这部分我们将调查诱导算子-半范数在基于分解的迭代下收敛的充分条件的细节。除了命题7以外，这部分的结果并不新鲜；然而，算子半范数的系统使用在结果的公式中提供了一种同一方法，是我们所感兴趣的。当讨论根据文章下面所述的定理6所得的各种收敛结果时我们会论证这一点。 的零空间和表示其范围，我们假设。这意味着（1）的解集非空且它给出了（1）的一些解的一个仿射空间 。 接着[5]，我们考虑一般情况下在（1）的哪种迭代矩阵是这样的形式 （4） 其中是一个有可能奇异的矩阵。由矩阵和诱导得到迭代 （5） 我们要注意这里，作为与非奇异相反的情况，对于给定的和，（4）中的矩阵并不唯一；参见[2]. 对于的最低假设要求它到是单射，因为不这样的话，如果，在中且，迭代式（5）的结果全都会迭代到，也就是迭代并不能收敛到（1）的解。 （4）中迭代算子的一般形式尤其能说明由分解所得到的迭代，其中非奇异，在这种情况下可被当做。（4）在时的迭代也是一些奇异矩阵的广义逆的伪逆解；见[4,11,12]，其中对这种迭代有过研究。这种情况尤其会出现在构造边界条件是诺依曼或罗宾型的施瓦茨迭代的分析中；见，例[13,16,18]。 根据下面的定义，迭代式（5）的收敛等价于是半收敛的；见，例[3,4,17]。 定义3.一个矩阵是半收敛的，如果，是系数1的唯一特征值且是的半单特征值，即就是，它的几何多重性等于代数多重性。 上的-半范数可以导出上的算子-半范数通过下式 （6） 在正定的情况下，表示而且因此迭代（5）是收敛的。在半正定的情况下，用算子-半范数可以有类似的结果。 定理4.令是Hermitian的并且是半正定的。令是迭代式（5）的迭代算子，且令。则 （i）在中是