用矩阵分解和Krylov子空间方法求解线性矩阵方程

来源 :中国科学院研究生院 中国科学院数学与系统科学研究院 中国科学院大学 | 被引量 : 0次 | 上传用户:hzwn001
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
矩阵方程的求解问题在控制论、结构工程、振动理论、系统参数识别和经济领域等方面有广泛的应用背景.近几十年来,一直受到国内外学者的广泛关注,并取得了一系列重要的成果.这一问题的研究是数值代数研究领域的一个十分活跃的分支.对于大型的矩阵方程(如Sylvester线性矩阵方程,Lyapunov线性矩阵方程,和非线性代数Riccati方程等)的研究,在国际国内已有一定的基础。例如,Erik Elmroth(瑞典),Isak Gustavson(瑞典),Bo Kagstrom(瑞典)和Fred Gustavson(美国)利用矩阵分解和递归方法给出了Sylvester型线性矩阵方程求解的数值方法,并设计出了高效的计算软件包(见SIAM Riew,Vol.26,No.1,2004),Yousef Saad(美国)对大型稀疏特征值问题进行了深入的研究,并大量地应用于实际之中。P.Lancaster(加拿大)和G.H.Golub(美国)等也分别研究了Lyapunov方程和Sylvester方程的相容性条件,最小二乘解和求解的算法。八十年代以来,国内的张磊,孙继广,胡锡炎和戴华等相继考虑了方程AX=B在约束条件下的解和最小二乘解及矩阵的最佳逼近问题,K.G.Woodgate(美国)和J.C.Allwright(美国)在对称半正定矩阵类上导出了方程AX=B的最小二乘解存在的条件,并给出了数值方法。  作者利用较多的数值代数知识,如矩阵的分解、广义逆和迭代方法等解决了几大类约束矩阵方程的解的相容性条件,通解表达式,最小二乘解和解的最佳逼近问题,特别地,将解代数方程系统的Krylov子空间方法,如CG,CGNE和CGNR方法应用到求解线性矩阵方程当中,其中包括AX=B,AXB=C,AX+XB=C(Sylvester方程),AX+ XA*=C(Lyapunov方程),AXB+CXD=E,(AX, XB)=(C,D)(左右逆特征值问题)等.为了适应大规模科学计算的需要,将矩阵分解方法和Krylov子空间方法结合起来,在存储量和计算时间上都获得了令人满意的结果。  报告的主要内容分为五章.第一章列出了通篇所用的基本符号,矩阵的范数和子空间理论,并列举了线性矩阵方程的基本类型及某些背景范例。  第二章介绍了矩阵分解理论,并利用矩阵分解理论:(1)解决了矩阵方程AX=B的两类反问题;(2)得到了矩阵方程AXB=C的极小范数Hermitian解和skew-Hermitian解;(3)得到了矩阵方程AXAT+BYBT=C的几类约束条件下的解。  第三章将求解线性代数方程组的共轭梯度法(CG方法)推广到求解对称正定型(SPD)线性矩阵方程AXB=C和AX+XB=C(Sylvester方程)中去。  第四章将求解线性代数方程组的共轭梯度法方程法(CGNE方法)推广到求解非对称正定型(SPD)线性矩阵方程AXB=C和(AX,XB)=(C,D)(左右逆特征值问题)中去,用迭代方法得到了这两类矩阵方程Hermitian极小范数数值解。  最后一章将矩阵分解的方法和Krylov子空间方法结合起来,给出了用递归法求大型Sylvester矩阵方程的解和大型矩阵方程AX=B的Hermitian解的思路。
其他文献
二硫化钼(MoS_2)纳米材料作为析氢反应(HER)催化剂Pt的替代材料得到了研究人员的广泛关注。其中,边缘结构对MoS_2纳米带催化析氢活性的影响甚为关键。通过施加应力和引入缺陷等方式能够显著地改变MoS_2纳米带体系的几何构型、电子性质和催化析氢活性。MoS_2纳米带的催化析氢活性及其调控的理论研究,对进一步提高MoS_2材料催化析氢活性、设计新型高效的电化学析氢反应催化剂有着重要的意义。因此
光学全息存储技术的不断发展对存储材料提出了严格要求,评价一种材料的全息性能需要综合考虑动态范围、光敏度、散射噪声、制备与处理过程等多方面因素。与无机晶体相比,有机聚
作为我国重要保护树种的杜仲树,既是重要的药材、又是提制优质橡胶的原材料,还是优美的园林绿化树种,倍受人们喜爱,近年来种植面积越来越大。但是,在杜仲栽培过程中,由于种
里德堡原子的实验研究取得突破性进展,激发了人们对扩展J-C模型研究的兴趣,其主要方法是通过腔耗散或加Kerr介质来修正J-C模型。把J-C模型同时扩展到腔耗散和加Kerr介质的情况,
取九个玉米自交系:合二、331、406,MO17、洛85、苏37、273、B87和292组成完全双列杂交(不包括反交)用Griffing分析法(方法2,模型Ⅰ)对株高、单株叶数及叶面积、粒重和籽粒密
茶树果实虽含有较高的脂肪,可供工业用,而且经过脱涩处理后还可供食用,但其经济价值远不如茶叶与水果,未被人们所重视。尤其对其形态与分类的研究更是冷门。笔者在多年对茶
近年来,随着我国基础工程建设的不断推进和铁路数据的通信网陆建设的深化,与铁路工程有关的通信业务,越来越多的寄托于网络数据通网上。本文主要对铁路数据通信网络和IP地址
在量子物理学创立于发展过程中,单光子源担任十分重要的角色,总是引起研究学者的关注。一方面,单光子源作为一个非经典光源可以被用来演示量子力学的基本原理,展现量子的奇异性。
化石燃料的大量开采和使用,导致能源短缺,环境污染以及全球气候变暖等一系列问题。而太阳能资源具有取之不尽用之不竭,清洁无污染,可靠安全等优点,使得太阳能转化技术成为当前研究
荧光探针能够与目标分析物(离子、小分子等)发生作用,通过荧光光谱的变化对目标分析物进行检测分析的一种化学物质,由荧光团、连接臂和识别基团三部分组成。其中,有机荧光团一般具有很大的共轭体系、刚性平面及能够被进一步修饰等特点。由于氧杂蒽染料具有独特的结构及良好的光学性质,因此在荧光探针的设计中常被用作荧光团。本论文以新型氧杂蒽染料为荧光团,设计出不同的荧光探针,用于检测识别氟离子和金属离子(铜离子和汞