论文部分内容阅读
大型辛矩阵特征值的计算在子结构链的动(静)力分析中十分重要,在离散时间最优控制大系统分析及金融数学中都有重要应用。在保证Hamilton结构始终不变的原则下来求解Hamilton矩阵的特征值是保证计算结果正确的最有效方式。常用的求解特征值的数值解法,只考虑数值精度,不考虑保证结构不变。 保守体系有辛规律,Hamilton体系是保守体系,因此有必要探讨应用辛算法来计算各类Ham ilton系统中的矩阵特征值。求解方式是多种多样的。Van Loan提出的平方约化法,保持了Hamilton结构,克服了常用的QR算法不能恰好保证每个半平面上都能求得n个特征值的缺陷。Benner提出了一个求解大型矩阵特征值问题的Lanczos方法,钟万勰建立了求解哈密尔顿矩阵特征值的共轭辛子空间逆迭代法和大型辛矩阵特征问题的逆迭代法。Bunse提出了求解实系数的Riccati方程的辛QR算法。本文的创新工作正是基于上述工作展开的。 文中第1-2章是学术背景研究。第1章系统地介绍了Hamilton体系,并研究了Ham ilton矩阵与辛矩阵的特性。第2章则介绍了矩阵特征值的常用数值计算方法。文中第3-4章主要是作者取得的一系列的创新成果,包括:从理论上建立了辛SL算法,分析了有效性和收敛性,以及如何用辛SL算法求解辛矩阵特征值。数值计算结果令人满意。 文中的创新成果主要有: ①就三种特殊形式的辛矩阵,建立了相应的求解特征值的算法。 ②针对上述三种矩阵,还建立了矩阵特征值的SL求解算法,并证明了算法的收敛性。 ③就反Hamilton矩阵,证明了也能应用SL算法有效求解特征值。 本文的第5章,是作者按照“上海交通大学数学系硕士研究生毕业要求”的条例完成的,是在阅读、理解大量科技文献之后经思考提炼而撰写的综合报告。主要综述了辛算法在弹性力学、波动方程以及DNA弹性杆力学分析中的应用及隐式辛算法的稳定性分析。将不同背景下的实际问题转化为Hamilton系统,选取合适的辛差分格式求解。数值算例表明:Hamilton系统的辛算法数值解是十分可靠的,数值结果是收敛的。