不可约M-矩阵的一种判别法.docVIP

不可约M-矩阵的一种判别法 许晓玲 (闽江学院数学系，福建，福州350108) 关晋瑞 (厦门大学数学科学学院，福建，厦门361005) 摘要：本文中我们提出了一个判定不可约M-矩阵的实用算法，给出了相应的理论分析，并用数值算例展示了该算法的有效性和优越性。 关键词：M-矩阵；判别法；不可约 中图分类号：O151.21 1 引言 M-矩阵是一类很重要的特殊矩阵，自1937年由Ostrowski提出之后，由于它的重要性和优美的性质，得到了深入的研究和广泛的应用。从那时起，新的性质和等价条件不断被发现，1977年Plemmons在[11]中总结的非奇异M-矩阵的等价条件已有40个，在后来的专著[2]中又扩充到多达50个。另一方面，M-矩阵的应用十分广泛，数学上应用在矩阵理论，微分方程数值解，Markov链，线性互补问题，线性方程组迭代法等问题中，其他学科如物理，生物，经济中也有着广泛的应用。这方面的详细内容可参考[2][7][8][15]。 下面我们给出M-矩阵的定义及一些基本性质，主要来自[2] [15]。 记。 对任意的，我们总可以将表示为，其中为常数，是一个非负矩阵。若，我们称是一个M-矩阵。特别的，当时称是一个非奇异M-矩阵，当时称是一个奇异M-矩阵。 关于非奇异M-矩阵我们有下面几个常见的等价条件。 定理1.1 设，则下列各条件等价： (1) 是一个非奇异的M-矩阵； (2) 可逆，且； (3) 存在，使得； (4) 的任意特征值都有正实部； (5) 的对角元素都为正的，且存在正的对角矩阵，使得是严格对角占优矩阵。 在实际应用中我们最关心的问题是，一个给定的矩阵是否为M-矩阵，以便作进一步的研究。因此，判别一个给定的矩阵是否为M-矩阵是一个关键。不过首先要说明的是，用M-矩阵的定义或者等价条件去验证并不是一个好的途径，因为要验证这些条件并非易事。从定理1.1我们可以看到，若用等价条件(1)去验证需要计算，对于一般的矩阵这并不容易求出。用(2)验证需要计算，运算量大且容易受误差影响。用(3)验证需要找到合适的正向量，然后并没有合适的途径去寻找。用(4)验证需要计算所有的特征值，更是不可行。 因此，寻求M-矩阵更简洁实用的判定条件是一个值得研究的课题。 近几十年来国内外很多数学工作者都对这个问题作了深入的研究，得到了许多优美而实用的判别条件，文献[3][4][5][6][9][10][12]是其中一些比较好的结果。然而在实际应用中发现，这些判别方法尽管各具特色却都有一些不足之处。总的来说存在这样一些问题：一是判别条件难以验证，这在应用中很不方便；二是判别范围过窄，这些判别条件大都是一些充分条件，只能判别出M-矩阵的一个子集，对于一般的M-矩阵却束手无策；三是对于奇异的M-矩阵则不能或不易判定；四是这些判别方法不易在计算机上实现。 针对上述判别法存在的这些问题，我们在本文中提出一种新的迭代判别法，该判别法的特点是：首先，判定范围广泛，非奇异或奇异的M-矩阵都可判别出来，只要要求判定矩阵不可约即可；其次，计算量小，每步迭代只需要对判定矩阵的某个列乘以一个非零常数即可；此外，易于在计算机上实现，在应用中很方便。 下面我们在第2节中给出迭代判别法的理论分析及算法，在第3节中给出一些数值例子，最后一节给出一个简单的总结。 2 主要结果 以下我们记，对于，记，，。在推导算法前我们需要用到下面的一些结论，其中引理2.1用M-矩阵的定义很容易验证，引理2.2来自[15，p377]，2.3来自[15，p376]，引理2.4来自[14，p38]。 引理2.1 设，则是M-矩阵当且仅当是M-矩阵，其中为正对角矩阵。 引理2.2 设是一个不可约M-矩阵，则，对任意。 引理2.3 是不可约M-矩阵，当且仅当存在向量，使得。 将在第行展开化简可得，这说明是广义(非严格)对角占优矩阵。关于广义对角占优矩阵的知识可以参考[15，Chapter 3]。上述引理也可以表述为：是不可约M-矩阵，当且仅当是广义(非严格)对角占优矩阵。 引理2.4 设矩阵，向量，若存在，使得，，则有 ，。 下面定理2.1中的矩阵我们称之为不可约对角占劣矩阵，该定理结论说明，不可约对角占劣矩阵是广义严格对角占劣矩阵。 定理2.1 设不可约，对角线元素非零，，对任意，且至少有一个严格不等式成立，则存在向量，，使得 ，对任意。 证：记。显然，若，令 。 由于不可约，非空。若，令 。 由于不可约，非空。继续上述过程，由于有限，因此存在使得。 对所有的，选取。然后对所有的用乘以的第列，得到新的矩阵。 注意此时在行仍然是对角占劣的，且对于任意的， 。 即在行也是对角占劣的。 对所有的，选取。 然后对所有的用乘以的第列，得到新的矩阵。类似上述讨论，得在行都是对角占劣的。依次下去，