大型稀疏鞍点问题的快速迭代算法研究

来源 :华南师范大学 | 被引量 : 0次 | 上传用户:xjdszcjl
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
鞍点问题在工程和科学计算上有着极其广泛的应用,如计算流体力学,偏微分方程的混合有限元离散,最优控制等各方面的许多问题都归结为大型稀疏鞍点问题的求解.随着问题规模的增大,用直接法求解这类问题的计算费用非常昂贵,迭代地求解该类问题成为唯一可行的选择.因此本论文的选题具有十分重要的现实意义. 分裂方法和Krylov子空间方法是人们所熟知的两类求解鞍点问题的迭代方法.分裂迭代方法算法结构简单,因而易于实施.然而,为了获得较快的收敛速度,这类方法一般都涉及到最优参数的选择.到目前为止,确定最优参数依然是一个非常困难的问题.Krylov子空间算法具有复杂的算法结构,且需要较大的储存空间.从而,这类方法用于求解小规模系统时较有效,但求解大型的鞍点问题时倾向于收敛很慢或失去效力.因而需要寻求好的预条件子,以改善此类算法的性能,拓宽其有效性的范围.所以,从事大型稀疏鞍点问题的快速迭代算法研究也非常必要. 本文首先提出了非对称块超松弛型(UBOR-型)分裂迭代方法,该方法带有三个松弛参数T,w1,w2,是文[54]中SOR-型方法的推广,算法结构上又完全不同于文[5]所述的GSOR方法.UBOR-型方法收敛的充分必要条件在实际应用中具有很强的可操作性.当参数满足关系式T=w1+w2-w1w2时,导出了最优参数的求解公式,并获得了相应的最优谱半径.理论分析表明,UBOR-型方法有比SOR-型方法更快的渐进收敛速度,且至少可以获得与GSOR方法相同的收敛速度.数值实验支持理论结果,且进一步揭示:依据迭代次数,这种新的方法能比GSOR方法有效得多.基于UBOR-型方法中系数矩阵的分裂,构造了块超松弛型预条件子,以加速Krylov子空间方法的收敛.由于系数矩阵是非对称的,这里选用一种非对称的Krylov子空间方法,即广义最小残量方法(GMRES).文中给出了预处理后系数矩阵的特征值分布.利用最小多项式的理论,考虑了块超松弛型预条件子中参数值的几种有效选择.数值实验进一步表明,通过适当选取参数,相应的预处理GMRES方法较未预处理的方法或UBOR-型方法具有快得多的收敛速度. 其次,基于系数矩阵的Hermite和斜-Hermite分裂及交替方向迭代的思想,本文提出了一种新的交替方向迭代方法(New-ADI),并给出了该算法收敛的充分必要条件及几个充分条件,这些条件在应用中易于检验.通过适当选取参数,新方法比最优情形的Uzawa算法及文[30]最近提出的的交替方向迭代法有更快的收敛速度和更少的计算耗费.数值实例进一步表明了新方法的有效性.而且,利用New-ADI算法中的矩阵分裂,构建了新的交替预条件子以提高GMRES方法的收敛速度.这种新的预条件子带有两个参数.继分析预处理后系数矩阵的特征值分布之后,考虑了参数的选择并获得了相应情形的拟最优参数.带有这些拟最优参数的预条件子大大改善了GMRES方法的收敛性能.关于线性化的Navier-Stokes方程的数值实验展示了新预条件子的有效性,对于大粘性系数的情形,这种有效性更明显. 最后,利用组合策略,构造了一种新的组合预条件子以改善鞍点系统的谱分布.通过定义一种新的内积,将不定的预处理后系数矩阵转换成关于该非标准内积为对称正定的矩阵,以应用非标准内积下的共轭梯度方法(CG)求解相应系统.继研究CG方法的可应用性条件后,基于预处理后系数矩阵的谱条件数分析,给出了较佳的组合参数值.带有如此的参数值,组合预条件子导致相应的CG方法获得了非常快的收敛速度.关于Stokes测试问题的数值实验验证了理论结果.
其他文献
多面体最佳逼近问题是逼近论中一个重要的分支,它的应用十分广泛,特别在约束条件下的统计推断领域,各类保序回归问题(包括广义保序回归和多维保序回归等)都可以化为多面体逼近问
请下载后查看,本文暂不支持在线获取查看简介。 Please download to view, this article does not support online access to view profile.
期刊
图论是一门古老而又有趣的学科。它主要研究用某种方式联系起来的若干事物之间的二元或者多元的关系,其中包括拓扑图论、代数图论、化学图论、算法图论、网络图论、模糊图论等
课堂提问是一种最直接的师生双边交流活动,是课堂教学的重要手段之一。它是教师课堂教学艺术的组成部分,也是课堂教学重要的一环,直接关系到课堂教学的成败。如何提高历史课
学位
863计划自动化领域在七五期间确实取得了很大成绩,我代表国家科委向所有参加863自动化领域的同志和各个单位、部门表示衷心的感谢。自动化领域从系统的观点出发,用集成的方法
本文六章主要分为三大部分。 第一部分是基础部分。为了研究四元数随机矩阵的多元分布函数,就必须给出四元数矩阵变换的Jacobi行列式。而本文的研究涉及到奇异的四元数矩阵
排队论(随机服务系统)是研究系统由于随机因素的干扰出现排队现象的规律性的一门学科,它适用于一切服务系统,包括通信系统、交通与运输系统、生成服务系统、存储与装卸系统、管理
作为一个大学生,我看到太多由于课后太过放松而产生材料力学挂科的现象.对此,我将对材料力学课后如何学习以提高学习效率提出我自己的观点.
倒向随机微分方程是近几年来迅速发展并引起人们极大兴趣的一门学科。它的理论不仅被广泛地认为是研究金融数学,例如期权和衍生证券定价问题的主要工具,而且也是研究随机控制、