几类结构矩阵谱问题.docVIP

厦门大学 硕士学位论文 几类结构矩阵的谱问题 姓名：曾祝明 申请学位级别：硕士 专业：计算数学 指导教师：卢琳璋 几类结构矩阵的谱问题  Ｕｌ  摘  要  众所周知，在工程计算和实际应用中有许多问题最终都归结为矩阵计算问题，而 且不同的应用会导出—些具有特殊结构的矩降ｉ博．最常见的—些结构矩阵有Ｔｏｅｐ－－ ｌｉｔｚ矩阵ｌａ／一ｊ】，Ｈａｎｋｅｌ矩阵【毗＋Ｊ】，Ｔｏｅｐｌｉｔｚ－ｐｌｕｓ－Ｈａｎｋｅｌ矩阵，Ｃａｕｃｈｙ矩阵【：盈１ ｉ】 等等．处理与这些结构矩阵有关的矩阵计算问题（例如计算特征值、求解线性方程组 等），若矩阵的阶数较小时，通常的经典算法是可行的（例如ＬＵ分解算法、ＱＲ算法 等）．然而，在许多实际应用当中，矩阵的阶数仃很大（ｎ—１０６—１０９）或某个线性方 程组需要多次计算直到得到—个满意的结果（例如迭代法时），此时这些经典的算法 由于代价太大而失去了实际意义． 因此，针对这些结构矩阵的特点而设计—些能利用它们的结构的，数值稳定的快 速算法，具有非常重要的意义．正因为结构矩阵在实际应用中所具有的重要意义，国 内外众多的学者将目光投入到这—领域．结构矩阵的快速算法中最著名的莫过于快 速傅里叶变换（即ＦＦＴ），有许多胬彭享法均是由快速傅里叶变换导出的．因此，著 名数学家Ｃｈａｒｌｅｓ Ｖａｎ Ｌｏａｎ曾这样评价快速傅里叶变换算法： “从计算的角度看， 快速傅里叶变换是本世纪最杰出的成就之一，毫不夸张地说，快速傅里叶变换改变 了科学与工程计算的面貌，如果没有它，生活将会是另—种景象”． 本论文主要研究了实Ｈａｎｋｅｌ．ｃｉｒｃｕｌａｎｔ和Ｈａｎｋｅｌ．ｓｋｅｗ．ｃｉｒｃｕｌａｎｔ矩阵的奇异值分 解，给出了对称Ｔｏｅｐｌｉｔｚ－ｐｌｕｓ－Ｈａｎｋｅｌ矩阵特征值的快速算法和这个计算矩阵特征值 算法的数值实验．理论和数值实验显示，这个侥速算法是行之有效的． 第—章，我们简单介绍了研究结构矩阵决速算法的现实意义，研究概况以及常 用的研究方法，同时也给出了与本论文有关的几．类结构矩阵的定义及其基本性质． 第二章，我们给出了ｎ阶对称Ｔｏｅｐｌｉｔｚ－ｐｌｕｓ－Ｈａｎｋｅｌ矩阵与—个佗维向量乘积的 快速算法；并利用他阶矩阵的对称陛，对其实施Ｌａｎｃｚｏｓ三对角化和ＱＲ对角化， 计算出矩阵的所有特征值．该算法的计算复杂度为Ｄ（ｎ２ ｌｏｇ ｎ）． 第三章，我们研究了实Ｈａｎｋｅｌ－ｃｉｒｃｕｌａｎｔ、Ｈａｎｋｅｌ－ｓｋｅｗ．ｃｉｒｃｕｌａｎｔ矩阵与Ｈａｎ－ ｋｅｌ矩阵的关系，并给出了它们的奇异值分解，为我们研究实Ｈａｎｋｅｌ－ｃｉｒｃｕｌａｎｔ矩阵、 Ｈａｎｋｅｌ－ｓｋｅｗ－ｃｉｒｃｕｌａｎｔ矩阵的奇异值分解提供了理论基础． 关键词：结构矩阵；快速算法；Ｔｏｅｐｌｉｔｚ矩阵；Ｈａｎｋｃｌ矩阵；对称Ｔｏｅｐｌｉｔｚ．ｐｌｕｓ－Ｈａｎｋｅｌ 矩阵；Ｈａｎｋｅｌ．ｃｉｒｃｕｌａｎｔ矩阵；Ｈａｎｋｅｌ－ｓｋｅｗ－ｃｉｒｃｕｌａｎｔ矩阵。 Ａｂｓｔｒａｃｔ  Ｍａｎｙ ｉｍｐｏｒｔａｎｔ ｐｒｏｂｌｅｍｓ ｉｎ ｅｎｇｉｎｅｅｒｉｎｇ ａｎｄ ａｐｐｌｉｅｄ  ｓｃｉｅｎｃｅｓ  ｃａｎ  ｂｅ ｆｉ－ ｎａｌｌｙ ｒｅｄｕｃｅｄ ｔｏ ｍａｔｒｉｘ ｃｏｍｐｕｔａｔｉｏｎ ｐｒｏｂｌｅｍｓ．Ｍｏｒｅｏｖｅｒ，ｖａｒｉｏｕｓ ａｐｐｌｉｃａｔｉｏｎｓ ｏｆｔｅｎ ｉｎｔｒｏｄｕｃｅ ｓｏｍｅ ｓｐｅｃｉａｌ ｓｔｒｕｃｔｕｒｅｓ ｉｎｔｏ ｔｈｅ ｃｏｒｒｅｓｐｏｎｄｉｎｇ ｍａｔｒｉｃｅｓ．Ａｍｏｎｇ ｃｌａｓｓｉｃａｌ ｅｘａｍｐｌｅｓ ａｒｅ Ｔｏｅｐｌｉｔｚ ｍａｔｒｉｃｅｓ【ａｉ－ｉ］，Ｈａｎｋｅｌ ｍａｔｒｉｃｅｓ【ａｉ州】，Ｔｏｅｐｌｉｔｚ－ ｐｌｕｓ－Ｈａｎｋｅｌ ｍａｔｒｉｃｅｓ，Ｃａｕｅｈｙ ｍａｔｒｉｃｅｓ【去】ａｎｄ ｏｔｈｅｒｓ．Ｗｈｅｎ ｗｅ ｄｅａｌ ｗｉｔｈ ｔｈｅ ｒｅｌａｔｅｄ ｍａｔｒｉｘ ｐｒｏｂｌｅｍ（ｓｕｃｈ ａｓ ｃｏｍｐｕｔｉｎｇ ｅｉｇｅｎｖａｌｕｅ，ｓｏｌｖｉｎｇ ｌｉｎｅａｒ ｓｙｓｔｅｍｓ ａｎｄ ８０ ｏｎ）ｗｉｔｈ ｓｐｅｃｉａｌ ｓｔｒｕｃｔｕｒｅ，ａｓ ｗｅ ｋｎｏｗ，ｔｈｅ ｓｔａｎｄａｒｄ ｌｉｎｅａｒ ａｌｇｅｂｒａ（ｓｕｃｈ ａｓ Ｇａｕｓｓｉａｎ ｅｌｉｍｉｎａｔｉｏｎ ａｌｇｏｒｉｔｈｍ，Ｑｒｔ ａｌｇｏｒｉｔｈｍ，ｅｔｃ）ａｒｅ，ｏｆ ｃｏｕｒｓｅ，ｒｅａｄｉｌｙ ａｖａｉｌ． ａｂｌｅ ｆｏｒ ｓｍａｌｌ－ｓｉｚｅ ｐｒｏｂｌｅｍｓ．Ｈｏｗｅｖｅｒ，ｉｎ ｍａｎｙ ｐｒａｃｔｉｃａｌ ａｐｐｌｉｃａｔｉｏｎｓ，ｔｈｅ ｏｒｄｅｒ ｏｆ ｍａｔｒｉｃｅｓ  ａｒｅ