大型线性方程组求解技术及在计算电磁学中应用研究.pdf

摘要 摘要 计算数学与科学工程计算涉及到航空航天、现代生物与医学、石油勘探、环 境科学、隐身器件设计等国民经济与国防建设的各个方面，其中往往需要求解一 个或一系列大型线性系统．而且，随着问题规模的增大，相应线性系统的未知数个 数大大增加，动辄上百万、千万，甚至上亿．这些大型线性方程组的求解是求解整 个问题的基础和关键所在，其计算量也占整个计算过程非常大的比重，有的甚至达 到80％以上．超大规模线性系统求解能力的缺乏成为解决某些实际问题的瓶颈．大 型线性方程组求解研究是现代科学计算的重要课题和焦点之一，高效、快速的求解 方法研究既有理论意义又有实际价值．本文旨在深入研究求解大型线性系统的高性 能算法，并特别针对电磁计算中产生的大型线性系统，构造有效的算法． 构造对称正定矩阵的不完全分解预条件子．不完全分解预条件迭代法的通常做 法是，先对系数矩阵进行重排，再不完全分解，然后再进行预条件迭代求解．本文将 近似最小度排序嵌入不完全分解过程，提出了一类新型的高效预条件技术．给出了 基于IKJ高斯消元的不完全Cholesky分解，并讨论了一些实现细节，然后分析了精确 计算和近似计算分解过程中节点度的方法．结合CG方法，提出的预处理技术能很大 地加快迭代法收敛速度．最后用丰富的数值实验表明与最小度排序耦合的不完全分 解预条件技术的效率，要高于按照一般做法得到的预条件子的效率． (AHSS)方法．这类方法包含两步迭代，本文从理论上研究了其收敛性以及最优参数 的选择问题．在LHSS迭代和AHSS迭代的每一步，都需要求解两个线性方程组．如果 用直接法求解这两个方程组，其每个方程组的计算量与求解原方程组相当，那么这 里的两步迭代就失去了意义．为了加快两步迭代的收敛速度，在求解这两个方程组 时，不准确求出其解，而是用迭代法求其近似解，得到了所谓的ILHSS与IAHSS方法． 给出的数值例子验证了本章提出方法的有效性． 基于选主元策略研究对称不定矩阵的较为稳定的预条件技术．该类预条件子修 预条件子一般也具有不定性，期望与原方程组的系数矩阵更相近．通过数值例子表 明，作用于SQMR方法时，该类预条件子显示出了很好的效果． 针对非对称矩阵，研究基于Givens变换的不完全正交预处理技术．直接将得到 摘要 的不完全QR分解作为原方程组的预条件子，并通过数值实验验证了提出的预条件 方法的高效性．然后着重研究了排序对Givens旋转的QR分解的影响，给出了一种减 少QR分解中所作Givens旋转次数的排序算法． 针对开域电磁计算问题中出现的大型线性方程组，研究适用有效的预条件子． 先讨论用FEM方法求解散射问题时遇到的线性方程组的求解，提出并实现了MIC分 解预条件子，然后研究FEM／MoM混合方法求解辐射和散射问题中产生的复杂线性 方程组的求解． 关键词：迭代法，预条件技术，Krylov子空间方法，电磁计算，不完全分解，排序技术， 选主元，非Hermitian正定矩阵 Ⅱ ABSTRACT ABSTRACT mathematicsandscientific aleinvolvedin Computational computing manyappli- cations，such觞acromeehanics，modembiomedicine，oilreservoir,stealth technology． ale tOsolveoneoraseriesof oflineal required Usually,We systems equations．Wi血 theincrementofthesizeofthe sizeofthe linea