一类Krylov子空间方法的加速算法研究

一类Krylov子空间方法的加速算法研究

论文摘要

在某些特殊情况下,Krylov子空间方法的残量下降曲线会出现不稳定的现象,如停滞现象,不规则振荡现象等,从而严重影响到Krylov子空间方法的收敛速度。本文以带多右端项的大规模稀疏线性系统为研究对象,针对上述不稳定问题提出两类用于求解对称正定线性系统的新型算法:带自适应重新启动过程的免中断块共轭梯度(Adaptive Restart procedures for the Breakdown-Free Block Conjugate Gradient,简称AR-BFBCG)算法,基于投影的收缩块共轭梯度(Projected variant of the Deflated Block Conjugate Gradient,简称PD-BCG)算法。详细的研究内容与主要成果如下:提出了一种自适应重启的AR-BFBCG方法。基于Powell在1977年针对共轭梯度算法提出的重新启动过程(后命名为:Powell’s restart),以及Dai等人于2004年在Powell’s restart基础上的工作,我们将这种重新启动过程推广并应用到求解带多右端项线性方程组的块共轭梯度(Block Conjugate Gradient,简称BCG)算法中。同时,结合Ji等人在2017年提出的免中断块共轭梯度(Breakdown-Free Block Conjugate Gradient,简称BFBCG)算法,提出AR-BFBCG方法用于求解带多右端项的对称正定线性系统。该方法继承了Powell’s restart的优点,可通过重启帮助消除BCG方法中出现的残差收敛行为不规则的现象,进而提高BCG的收敛速度。同时保留了BFBCG方法中免秩亏损的优点,可避免带重启的BCG方法在执行中可能会出现由秩亏损而导致的算法中断问题。数值试验表明,这类方法可有效“打破”残差收敛行为中的不规则现象,从而在一定程度上减少计算量,提高BFBCG算法的收敛速度。此外,在求解病态系统,或带秩亏损的多右端项线性系统时,该方法的数值优越性更加明显。设计了一种基于投影的PD-BCG方法。Chen在2011年提出收缩块共轭梯度(Deflated Block Conjugate Gradient,简称D-BCG)算法,该方法可通过收缩掉极端特征值(比如小特征值)来减小预处理矩阵的条件数,从而提高BCG算法收敛速度。但在有限精度算法(finite arithmetic)框架下,残差空间与收缩空间之间的正交性会随着算法的执行而逐渐丢失,而这种正交性的丢失会明显地影响算法的稳定性,导致残量下降曲线出现长时间停滞,甚至因正交性丢失过于严重而根本无法收敛。重正交(reorthogonalization)方法可用于修复这种丢失的正交性,但是其代价相当大。基于此,我们提出了PD-BCG算法来“延缓”这种正交性的丢失,并且理论表明这种“延缓”不需要任何额外代价,同时理论还表明PD-BCG与DBCG在数学上是等价的。最后数值实验表明,PD-BCG算法能够将这种正交性保持在一定程度内,进而有效解决了收敛停滞甚至不收敛的问题,验证了其有效性与数值稳定性。

论文目录

  • 摘要
  • ABSTRACT
  • 第一章 绪论
  •   1.1 研究工作背景和意义
  •   1.2 块Krylov子空间方法的发展概述
  •     1.2.1 块Krylov子空间算法
  •     1.2.2 块共轭梯度算法
  •   1.3 本文主要内容和创新点
  •   1.4 本文结构安排
  •   1.5 符号说明
  • 第二章 求解带多右端线性系统的AR-BFBCG算法
  •   2.1 引言
  •   2.2 AR-BFBCG算法
  •     2.2.1 BFBCG算法
  •     2.2.2 动态重启策略
  •     2.2.3 AR-BFBCG算法的构造
  •     2.2.4 AR-BFBCG算法的性质
  •   2.3 数值实验
  •   2.4 本章小结
  • 第三章 求解带多右端线性系统的PD-BCG算法
  •   3.1 引言
  •   3.2 PD-BCG算法
  •     3.2.1 D-BCG算法
  •     3.2.2 PD-BCG算法的构造
  •     3.2.3 PD-BCG算法的性质
  •   3.3 自适应重正交策略
  •   3.4 数值实验
  •   3.5 本章小结
  • 第四章 总结与展望
  •   4.1 总结
  •   4.2 后续工作展望
  • 致谢
  • 参考文献
  • 攻硕期间取得的研究成果
  • 文章来源

    类型: 硕士论文

    作者: 向燕菲

    导师: 荆燕飞

    关键词: 免中断块共轭梯度方法,收缩块共轭梯度方法,投影,收缩特征值,自适应重启

    来源: 电子科技大学

    年度: 2019

    分类: 基础科学

    专业: 数学,数学

    单位: 电子科技大学

    分类号: O241.6

    总页数: 53

    文件大小: 17089K

    下载量: 43

    相关论文文献

    • [1].一种全局收敛的杂交谱共轭梯度方法[J]. 广东技术师范学院学报 2015(05)
    • [2].非线性共轭梯度反演法在铁路特长隧道勘查中的应用[J]. 铁道科学与工程学报 2020(05)
    • [3].一种求解矩阵填充问题的交替共轭梯度最小化法[J]. 太原师范学院学报(自然科学版) 2020(03)
    • [4].一种最大共轭梯度连续泛函的网络峰值预测[J]. 信阳师范学院学报(自然科学版) 2015(02)
    • [5].大地电磁场非线性共轭梯度三维反演[J]. 矿床地质 2012(S1)
    • [6].基于共轭梯度分解算法的电网谐波估计[J]. 计算机技术与发展 2013(02)
    • [7].2.5维海洋非线性共轭梯度反演应用效果研究[J]. 海洋技术学报 2018(02)
    • [8].一个新的解非线性对称方程组的非单调共轭梯度方法[J]. 广西科学 2009(02)
    • [9].大地电磁非线性共轭梯度反演试验——以内蒙古扎鲁特盆地数据为例[J]. 地球物理学进展 2014(05)
    • [10].改进的约束共轭梯度闪光照相图像重建算法[J]. 强激光与粒子束 2011(08)
    • [11].大型线性方程组的迭代求解(英文)[J]. 物理化学学报 2009(10)
    • [12].一类具有充分下降性的共轭梯度方法[J]. 广东技术师范学院学报 2010(06)
    • [13].三维地质模型中地震波共轭梯度非线性走时反演[J]. 地震地磁观测与研究 2018(04)
    • [14].基于共轭梯度的极速学习机[J]. 计算机应用 2015(10)
    • [15].一种新的非线性共轭梯度方法及其收敛性(英文)[J]. 数学杂志 2013(06)
    • [16].大地电磁非线性共轭梯度一维反演[J]. 石油物探 2014(06)
    • [17].一个三项LS共轭梯度方法[J]. 广西科学 2013(04)
    • [18].基于共轭梯度对数分解的大数据分类模型[J]. 科技通报 2015(12)
    • [19].修正LS共轭梯度方法及其收敛性[J]. 西南师范大学学报(自然科学版) 2016(07)
    • [20].基于混合随机共轭梯度的最小二乘逆时偏移[J]. 石油地球物理勘探 2018(06)
    • [21].概率功能度量求解的共轭梯度步长调节法[J]. 计算力学学报 2018(06)
    • [22].基于共轭梯度迭代算法受控AR模型的参数辨识[J]. 吉首大学学报(自然科学版) 2016(06)
    • [23].用稳定双共轭梯度方法数值求解球坐标系下的Poisson方程[J]. 中国科学技术大学学报 2017(08)
    • [24].基于二维数据的加权非线性共轭梯度三维反演[J]. 工程地球物理学报 2016(05)
    • [25].求解稀疏线性方程组的预处理共轭梯度并行算法[J]. 纺织高校基础科学学报 2014(04)
    • [26].一个求解非线性不等式约束优化问题的带有共轭梯度参数的广义梯度投影算法(英文)[J]. 应用数学 2011(02)
    • [27].一类修正的共轭梯度方法及其全局收敛性[J]. 云南大学学报(自然科学版) 2008(S1)
    • [28].基于加权奇异值分解截断共轭梯度的电容层析图像重建[J]. 光学精密工程 2010(03)
    • [29].修正的Dai-Liao三项共轭梯度方法[J]. 四川理工学院学报(自然科学版) 2019(05)
    • [30].三项修正LS共轭梯度方法及其收敛性研究[J]. 应用数学学报 2017(06)

    标签:;  ;  ;  ;  ;  

    一类Krylov子空间方法的加速算法研究
    下载Doc文档

    猜你喜欢