外文翻译word版(共7页).doc
《外文翻译word版(共7页).doc》由会员分享,可在线阅读,更多相关《外文翻译word版(共7页).doc(7页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、精选优质文档-倾情为你奉上RICHARD W.COTTLE斯坦福大学对于矩阵schur补的表现性质 (基于1974年11月12日的会议)摘要在本文中,笔者关注的是Schur补可以在数值线性代数中使用的一些方法。 首先概述的是变量消除和块枢转:商属性是为了下一个用于测试的矩阵的前主子式是否是非零的特定标志。另一个有用的应用是在计算实对称矩阵的惯性:惯性有助于检查这种矩阵正(半)定性。 这可以在以数学规划问题被利用检查非凸二次函数的准凸(拟凸)的非负定限。1.-介绍在1917年,I.schur18 发表了带有Schur引理的论文,这篇论文后来被称为Schur的行列式公式,它讲的是,如果M是一个方阵
2、,并且A是非奇异的M的主子阵,然后,则 (1) 。 最近,E.Haynsworth和其他人已经引入术语Schur补作为一个矩阵名称,比如,它出现在(1)中。参见12。更确切的说, 当,时,就把称作在M中A的schur补。习惯的写成 (2) 之所以这种表示法吸引人,是因为(1),则(3) 。在某些不参考这些性质是公知的或在数学科学的其他分支被使用的情况下,有一些论文中已经出现了对Schur补的两个主要理论特性的开发。本说明的目的是回顾schur补的产生以及如何它们的属性可以在数值线性代数被利用的方式。这个问题的更广泛的理解可以在作者的论文中找到8。 2.-一般定理大家在解线性方程组的过程中,消去
3、一个或多个变量都遇到了schur补。比如,思考Mz=0这个线性方程组, 因此(4) (5) 如果 ,则消除变量x可得出下式 (6) 另一个常见的例子是在块旋转中,有这样一个系统(在均一的情况下)(7) (8) 其中,则这样的操作被称为一个块主体枢21,17,23,或回转10。 在上面的两个例子中,M必须是方阵。当M是方阵且它的主子阵A和D是可逆矩阵时,在这种情况下,和的逆在公式中扮演了一个突出的角色。Schur补被Carlson,,Hayns-worth,,和Markham在摩尔 - 彭罗斯逆中推广,W. N. Anderson, Jr开展了另一种类型的发展,使用在短路操作代替Schur补方面
4、,如想深入了解电路网络中的短路操作理论可参考3,4,11。3.-商属性已知M,A和E为方阵,且 则矩阵Schur补的商属性可描述为:(10) Crabtree 和 Haynsworth 9提出了商公式(10),接着Ostrowski16对它提出了质疑,但是它早已在数理统计(在计算多元条件正态分布的协方差矩阵)和数学规划关键代数中所熟知。见1。商属性有一个好的应用是:判断一个矩阵的主子式是否为奇异的。如果被测试的矩阵是M,且它的主子式的矩阵顺序为,则的起始输入为的起始输入,特别的有. 如果通过了测试,即它是非零的或者是所被要求的,那么的起始输入是定义好的,它的起始输入可被作为或者。因此,可得到一
5、个信息,关于从到,如果通过了测试,然后可通过计算得到 但是通过商属性可得: 所以可获得关于det 的一个信息,这个过程可一直重复知道M的所有主子式被测试完或者找到一个未能拥有这种性质的主子式。4-惯性公式矩阵Schur补的第二个性质是在惯性的概念上是有用的。当前的讨论是对于实对称矩阵,比如,对于矩阵M来说,它的惯性是一个三维的 其中,分别对应矩阵M的正,负,和零特征值的数量,可包含多重的。因此,对于一个阶为n的矩阵来说 这样,可得出关于秩和符号之间的数量关系 显然,从矩阵M的秩,符号,和阶可算出矩阵M的惯性。 给出分块矩阵 ,其中A是非退化的矩阵,易得出(11) 这种关系被叫做惯性公式,由Ha
6、yns-wortht 12所证明。这个公式被Haynsworth 和 Ostrowski计算分块矩阵的惯性,对于这个分块矩阵的一个特殊形式是 ,其中A是任意且对称的,B是阶为k的非退化的矩阵,Haynsworth 和Ostrowski在前引中给出了结果的备用证明, 由于Carlson和 Schneider 6中所描述的那样(12) 在 8 中给出了 (12) 的归纳证明,先使用了是公式的成立,然后使矩阵Schur补M中的一个2 X 2的矩阵 得出了一个和M有相同性质结构的矩阵。5.-对于计算惯性的一个算法对于惯性公式(11)和特殊情况(12), k=1在 计算实对称矩阵M的惯性中扮演了一个重要
7、角色。注意到M的对角线上,如果M的对角线上有一个非零的输入,如,则通过(11)可得 (13) 其中 如果M是非零的且在对角线上存在零元素,则存在非零输入,如,并且 ,则有 (14) 如果M是一个阶为n的零矩阵,则 我们注意到当M是对称矩阵时,也是对称矩阵。 矩阵M必定有下列性质之一:(P1) ;(P2) , (P3) 。我们采用了下列公约,让成为空矩阵并设置让表示要计算的的Schur补,让表示枢纽(块),则 则一般地算法的步骤仅仅对于,它的阶从nk. k从0开始。步骤 0 :到目前为止,我们有(15) 如果有性质(),则到步骤。步骤 1 :我们可以假设,让 然后 如果,则停止,否则返回到步骤0
8、,并用k+1代替k。步骤 2 :我们可以假设,让 , 则 如果,则停止,否则返回到步骤0,并用k+2代替k。步骤 3 :在这种情况下,因为它的阶是nk,我们有 ,停止,此时M的惯性已被找到。这个算法是基于确定的算术上的。然而当在数字电路计算机中使用有限的准确性时,它缺少数字稳定性。这个缺陷由Bunch 和 Parlett 5所克服,对于获得一个分解,它给出了数字性地稳定对角线枢纽步骤。其中L为单位下三角,D是1 X 1和2 X 2矩阵的直接总和,且。或许惯性概念最熟练的应用被发现是在测试实对称矩阵是否为正定(半)矩阵方面。它在本质上是确定的的问题而不是一个寻找特征值本身,以便对发现他们是否是负
9、的。测试非凸二次函数对于准凸(或伪凸)在上的非负卦限是在二次规划中。这可以通过进行两个试验来完成。第一是验证 是一个非正矩阵,第二是查明。如果当的(15)中的第二成分超过1以上,上述用于计算M的惯性算法可以总是被终止。参考文献1 ANDERSON T. W., An Introduction to Multivariate Statistical Analysis.Wiley, New York (1958).2 ANDERSON W. N. Jr., Shorted operators. SIAM J. Appl. Math. 20 (1971),520-525.3 ANDERSON W.
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 外文 翻译 word
限制150内