计算矩阵符号函数和不变子空间的Padé逼近方法

来源 :浙江大学 | 被引量 : 0次 | 上传用户:c_zhang08
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文主要研究矩阵特征值问题中的一个典型问题:不变子空间的计算。历史上出现过的许多求解特征值问题的经典算法在求解此问题时,都碰到了很大的困难,有时甚至不能采用。近年来不变子空间问题引起了许多学者的关注,并产生了许多优秀的算法,矩阵符号函数技术就是其中的一种。 在本文中,我们首先介绍了矩阵符号函数的定义及性质,回顾了计算矩阵符号函数的经典算法,然后分析了历史上的一些改进算法,在此基础上进一步研究矩阵符号函数数值计算方法,及其在计算指定特征值范围的不变子空间上的应用。首先,我们考虑传统的计算符号函数的迭代法的一种有理形式的改进算法。这一改进形式并不增加计算量,但可以很好地改进算法的稳定性和计算精度。同时,对于Hermite情形,这一有理方法适用于借助于矩阵QR分解来避免算法中的求逆运算。其次,作为本论文的最主要部分,我们首次提出了利用Padé逼近计算矩阵符号函数的数值计算方法。一阶Padé逼近被用作为基本的迭代函数。与传统的计算符号函数的迭代法相比,一阶Padé迭代方法具有明显的优点:在基本不增加计算量的同时,它有更快的收敛性。我们分析并证明了其三阶收敛速度。进而,我们也初步讨论了用计算量较少的高阶Padé逼近,来近似地计算不变子空间。本文最后给出了数值例子,从矩阵的性态等方面说明所提出了三种算法的数值特点。数值例子表明,有理形式的改进算法以及一阶Padé迭代算法,对于解决不变子空间问题是行之有效的。而高阶Padé逼近,对于指定区域外部特征值与指定区域没有很好分离度的困难情形,有着引人的特性。这一部分有待于进一步地深入研究。
其他文献
早在18世纪初,人们开始对弦振动方程问题及其解法进行探讨研究,这标志着偏微分方程这门学科开始诞生.但是,它并未引起很多学者的重视,直到19世纪它才迅速发展起来.至今,国内外已经
本文给出债务关系的模糊矩阵表示,通过模糊矩阵幂运算实现了间接债务关系的基本要素-债务链,债务圈,债务量的刻画。为三角债的分析研究和决策提供依据。 给出划转的数学公式,第一次把圈和链统一在一个模式之下,并通过矩阵的变换得以实现。这样,划转可以从任何债务链着手,从根本上简化了准备过程。进一步的讨论从理论上证明了它具有不改变单位债务量和可实现性优点且易于做成软件系统,实现电算化,成为解决三角债问题
本文旨在研究应用十分广泛的随机优化算法及其收敛性问题.首先研究了第二阶段问题为二次规划的补偿随机规划问题,提出了该问题的基于Montc CarIo或Quasi-MontcCarlo 模拟的近
分数阶微积分在数学和工程方面已经成为人们特别熟知的概念,其是关于任意阶积分和微分的理论,推广了整数阶积分和任意阶导数,是目前非线性泛函分析最活跃的领域之一.其中,分
非线性泛函分析是数学领域中一门重要的学科,广泛地应用在各领域出现的非线性问题,如物理学、工程学、经济数学等.对于非线性泛函分析的深入研究,有着深刻的理论和应用价值.
本文首先利用试验函数方法研究了无界区域上两类非局部抛物型方程的初(边)值问题解的不存在性;另外,我们利用能量方法,研究了一类四阶波动方程的整体解的存在性和不存在性,以及整