复Hermite矩阵特征值问题的实IsotropicLanczos算法-福州大学学报.PDF

第４１卷 第３期 福州大学学报（自然科学版） Ｖｏｌ．４１Ｎｏ．３ ２０１３年６月 ＪｏｕｒｎａｌｏｆＦｕｚｈｏｕＵｎｉｖｅｒｓｉｔｙ（ＮａｔｕｒａｌＳｃｉｅｎｃｅＥｄｉｔｉｏｎ） Ｊｕｎ．２０１３ ＤＯＩ：１０．７６３１／ｉｓｓｎ．１０００－２２４３．２０１３．０３．０２６１ 文章编号：１０００－２２４３（２０１３）０３－０２６１－０５ 复Ｈｅｒｍｉｔｅ矩阵特征值问题的实ＩｓｏｔｒｏｐｉｃＬａｎｃｚｏｓ算法 １ ２ 赖降周 ，卢琳璋 （１．福州大学数学与计算机科学学院，福州 福建　３５０１１６；２．厦门大学数学科学学院，厦门福建　３６１００５） 摘要：研究了复Ｈｅｒｍｉｔｅ矩阵经Ｗｉｌｋｉｎｓｏｎ实数转化后矩阵的性质，利用其对称性和反 Ｈａｍｉｌｔｏｎｉａｎ结构，给出 了特征值问题的隐式重启ＩｓｏｔｒｏｐｉｃＬａｎｃｚｏｓ保结构算法．数值试验表明，这种方法求解出的特征对残量很小， 具有较高的精度． 关键词：Ｈｅｒｍｉｔｅ矩阵；特征值问题；Ｈａｍｉｌｔｏｎｉａｎ结构；ＩｓｏｔｒｏｐｉｃＬａｎｃｚｏｓ算法 中图分类号：Ｏ２４１　 文献标识码：Ａ ＡｒｅａｌｉｓｏｔｒｏｐｉｃＬａｎｃｚｏｓａｌｇｏｒｉｔｈｍｆｏｒｅｉｇｅｎｐｒｏｂｌｅｍｏｆｃｏｍｐｌｅｘＨｅｒｍｉｔｅｍａｔｒｉｘ １ ２ ＬＡＩＪｉａｎｇ－ｚｈｏｕ，ＬＵＬｉｎ－ｚｈａｎｇ （１．ＣｏｌｌｅｇｅｏｆＭａｔｈｅｍａｔｉｃａｎｄＣｏｍｐｕｔｅｒＳｃｉｅｎｃｅ，ＦｕｚｈｏｕＵｎｉｖｅｒｓｉｔｙ，Ｆｕｚｈｏｕ，Ｆｕｊｉａｎ３５０１１６Ｃｈｉｎａ； ２．ＣｏｌｌｅｇｅｏｆＭａｔｈｅｍａｔｉｃＳｃｉｅｎｃｅ，ＸｉａｍｅｎＵｎｉｖｅｒｓｉｔｙ，Ｘｉａｍｅｎ，Ｆｕｊｉａｎ，３６１００５，Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ：ＴｈｉｓａｒｔｉｃｌｅｄｅａｌｓｗｉｔｈｔｈｅｅｉｇｅｎｐｒｏｂｌｅｍｏｆｃｏｍｐｌｅｘＨｅｒｍｉｔｅｍａｔｒｉｘ．Ｆｉｒｓｔｌｙｔｈｅｐｒｏｐｏｓｉｔｉｏｎ ｏｆｔｈｅｍａｔｒｉｘｖｉａＷｉｌｋｉｎｓｏｎ’ｓｔｒａｎｓｆｏｒｍｗａｓｒｅｓｅａｒｃｈｅｄ，ｔｈｅｎａｎｉｍｐｌｉｃｉｔｒｅｓｔａｒｔＩｓｏｔｒｏｐｉｃＬａｎｃｚｏｓｍｅｔｈ ｏｄｗａｓａｌｓｏｐｒｅｓｅｎｔｅｄａｉｄｉｔｓｓｙｍｍｅｔｒｙａｎｄｓｋｅｗ－Ｈａｍｉｌｔｏｎｉａｎｐｒｏｐｏｓｉｔｉｏｎ．Ｆｉｎａｌｌｙｎｕｍｅｒｉｃａｌｒｅｓｕｌｔｓ ｓｈｏｗｅｄｔｈａｔｔｈｉｓａｌｇｏｒｉｔｈｍｃａｎｅｆｆｉｃｉｅｎｔｌｙｓｏｌｖｅｔｈｅｅｉｇｅｎｖａｌｕｅｏｆｃｏｍｐｌｅｘＨｅｒｍｉｔｅｍａｔｒｉｘｗｉｔｈｈｉｇｈｐｒｅ ｃｉｓｉｏｎ． Ｋｅｙｗｏｒｄｓ：ｅｉｇｅｎｖａｌｕｅｐｒｏｂｌｅｍ；Ｈａｍｉｌｔｏｎｉａｎｓｔｒｕｃｔｕｒｅ；ＩｓｏｔｒｏｐｉｃＬａｎｃｚｏｓｍｅｔｈｏｄ ０　引言 对称矩阵和Ｈｅｒｍｉｔｅ矩阵有许多共同的地方，以至于大家似乎通常只习惯于讨论对称矩阵而有忽视 ［１］ Ｈｅｒｍｉｔｅ矩阵的一种倾向 ．事实上，从计算的角度来说，这两类矩阵最不同的地方是：Ｈｅｒｍｉｔｅ矩阵通常 要涉及到复运算，而实对称阵不要．具体到特征值问题来说，复Ｈｅｒｍｉｔｅ矩阵的特征值都是实的，但特征 向量却肯定是复的．由于复运算的工作量通常是实运算的４倍，因此，如何减少计算Ｈｅｒｍｉｔｅ矩阵特征对 的复运算量、非常值得研究． ［２］ Ｗｉｌｋｉｎｓｏｎ 考虑了这样的问题，他把一个Ｎ阶复Ｈｅｒｍｉｔｅ矩阵特征对计算问题转化为一个２Ｎ阶实对 称阵特征对的计算问题，开辟了一条通过变换方式将它的特征值问题转化为实对称阵的特征值问题来处 理之路．具体来说，假定 Ｈ ＝Ａ＋ｉＢ