线性响应特征值问题的数值方法研究

来源 :华中师范大学 | 被引量 : 0次 | 上传用户:cyalil
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文将主要讨论如下形式的线性响应特征值问题:Hz=[(0)MK(0)][yx]=λ[yx]=λz这里的K和M是n×n实对称矩阵,并且其中一个是正定的.这种形式的特征值问题广泛应用于时间相关的密度函数理论的线性响应扰动分析.因此,称之为线性响应特征值问题,或随机相位近似(RPA)特征值问题.本文主要讨论线性响应特征值问题的数值方法以及相应的收敛性分析。全文共分为五章。  第一章主要介绍了特征值问题的背景和求解特征值问题的投影类方法,以及线性响应特征值问题的背景和发展现状。  第二章讨论了线性响应特征值的问题的一些基本理论.在这些定理中,有一些是线性响应特征值问题的基本性质,而有一些则是最近的一些新的理论结果.另外,在这一章的末尾,我们还简单介绍了Bai和Li(SIAMJ.MatrixAnal.Appl.,toappear)提出的求解线性响应特征值问题的局部最优的块预条件的4-维共轭梯度法(LOBP4DCG)。  我们主要工作集中在第三章和第四章.在第三章中,我们分析了线性响应特征值问题的两类Lanczos类型方法的收敛性。第一种是求解对称特征值问题的经典Lanczos方法的自然推广,而第二种是则由Tisper(JETPLetters,70(1999),PP.751-755)提出的.通过收敛性分析,我们得到了这两种方法的特征值和特征向量的收敛性定理.这些收敛性定理暗示着第一种Lanczos方法要比Tisper的Lanczos方法快的多.数值实验也证实了我们的结论。  在第四章中,我们提出了求解线性响应特征值问题的块Chebyshev-Davidson方法,并用以求其最小正的特征值以及所对应的特征向量.Chebyshev-Davidson算法的关键是需要一个有效的特征值上界.为此,我们提出了一个基于合理假设条件下的上界估计算法,同时我们还给出了Chebyshev-Davidson的自适应策略.当这个假设条件不成立的情况时,自适应的策略可以自适应的调整上界,从而保证了Chebyshev-Davidson方法仍然是有效的.另外,我们还估计了Chebyshev-Davidson算法的收敛率.最后的数值实验了也说明了Chebyshev-Davidson方法是有效的.在最后一章中,我们总结了本文的主要工作,以及对进一步研究线性响应特征值问题的展望。
其他文献
声波和电磁波传输问题在逆传输理论中扮演着重要的角色。这是由于材料的折射率不能够由远场数据唯一确定,即使多频数据能够被得到。最近,一种新的质量方法出现了,该方法是利用传
3一李代数在数学和数学物理的相关领域有着广泛的应用,特别是度量3一李代数在弦理论中的应用,吸引了越来越多的数学和物理工作者的关注,本文主要研究3一李双代数的结构分类.基于
The competitive removal of copper and cadmium from aqueous solutions using scoria has been investigated. Scoria was characterized by various methods, such as XR
本文主要针对一类线性Sobolev方程和一类非线性Sobolev方程分别构造了不同的数值计算格式.第一章简单介绍了数值求解Sobolev方程的发展现状。  第二章研究了一类线性Sobole
概率极限理论是概率统计学科的一个主要分支,也是概率论的其它分支和数理统计的重要基础.人们在实践中认识到事件发生的频率具有渐进稳定性,即随着试验次数的增加,事件发生的频
仿切触度量几何作为仿复几何的奇数维对偶,无论在理论数学还是数学物理方面都具有重要的研究价值.自从上世纪七十年代这类流形被提出以来,一直被众多几何学家与物理学家所关注.
设G是一个n阶有限群,令Ψ(G)=∑x∈Go(x),其中o(x)表示元素x的阶数.在文献[3]中,H.Amiri和I.M.Isaacs等人证明了在所有的n阶群中,元素阶数之和最大的群一定是循环群Zn.如果对于任
近几年来,对超有限Ⅱ1型因子R中算子的研究非常广泛.本文研究了超有限Ⅱ1型因子中一类算子uf(v),其中f是单位圆周上S1上有界的勒贝格可测函数,并且u,v是R的两个生成元满足u*u=v*v
本文通过对荣华二采区10
设G是一个图,A是一个阿贝尔群,对G通过连续收缩非平凡的A-连通的子图,直到没有非平凡的A-连通的子图剩余为止,得到的图记为G*,我们就说G能A-可收缩到G*.K4通过增加一个顶点v且点v