二阶双正交过程和二次特征值问题求解

来源 :复旦大学 | 被引量 : 0次 | 上传用户:ShakeBabyYTT
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
在该文中,我们提出了二阶双正交过程的算法,通过这一过程可以得到二阶Krylov子空间的一组双正交基.有别于标准的Krylov子空间,二阶Krylov子空间的构造需要基于一对矩阵和初试向量.给定一组双正交基之后,我们采用斜投影技术来求解二次特征值问题.由于在二阶双正交过程中只涉及到矩阵向量乘法运算,同时通过揭示其中的矩阵关系,我们对存储量进行了优化;又因为在算法中隐式的构造了线性化问题的Krylov子空间,使得算法不仅在存储空间上而且在运算时间上达到较优的效果,所以这个方法适用于大规模和稀疏问题的求解.我们给出的数值例子也证实,通过这一算法得到的Ritz值和对应的Titz向量可以快速和准确的收敛到真实解,并且通常是那些具有较大模的特征值.
其他文献
在粗糙集理论的代数观点下,约简关注的只有数据库中相容示例的部分;而在粗糙集理论的信息论观点下,约简开始关注库中不相容的部分,试图保证数据库中所有数据的不同分辨关系不
大气非线性重力波演变和相互作用过程对重大灾害性天气的形成具有重要的触发作用,重力波的演变囿于Schrodinger动力学.研究非线性重力波的演化方程、解的存在性和数值格式不
该文工作属于模糊矩阵理论和应用研究,具体工作针对一般(非相似)模糊矩阵传递闭包的计算、简化与应用.该文首先给出了网络最大路的概念,并把其归结为求模糊矩阵的传递闭包,这
本文运用Galerkin方法证明了一类Kirchhoff型非线性阻尼耦合梁方程组{ utt+u(4)-M(‖u(1)‖2+‖v(1)‖2)u(2)+f(u)+g(ut)=h1(x,t)vtt+v(4)-M(‖u(1)‖2+‖v(1)‖2)u(2)+f(v)+
近年来,微分方程边值问题已经成为非线性微分方程理论研究中一个比较活跃且具有现实意义的问题,在许多领域中都提出了由微分方程边值问题来描述的数学模型.四阶微分方程起源于