求解某些特别稀疏线性系统的数值解法.pdfVIP

中文摘要 摘要 大型稀疏线性系统来源于很多应用领域，譬如流体动力学，结构分析，电磁场 计算等等．将描述自然现象的偏微分方程离散后，通常就会得到一个稀疏的线性 系统．这样一来，实时高效的求解大型的稀疏线性系统对整个应用问题的解决有 着至关重要的作用．因此，近年来无论国内还是国外，大规模稀疏线性系统的求解 算法的研究已成为大规模科学与工程计算的一个重要研究领域．进一步，由于许 多实际问题产生的大规模稀疏线性系统，其系数矩阵往往都是具有某种特殊形式 或者某种特殊结构，因此本文主要研究的是一些特殊形式的稀疏线性系统快速、 有效的数值求解方法．全文共分为五章． 第一章介绍了大规模稀疏线性系统问题的来源、历史、发展现状以及本文所 涉及的几种特殊稀疏线性系统． 在第二章，我们给出了求解大规模稀疏位移线性系统的一种增广的重新启动 GMRES方法：每次重新启动时，将母系统所得到的多个误差向量添加到求解的 Krylov子空间中去，在新的增广的空间中求解母系统，而子系统的解通过强行使 其残量和母系统的残量平行得到．这样不仅能使我们在同一个空间中求解子母系 统，还能加速求解位移线性系统重新启动GMRES方法的收敛速度，数值试验也 表明这种方法的高效性． 第三章针对块三对角系统，给出了一种切频率过滤预条件子的变形．新的预 条件子是建立在块三对角矩阵的一种组合分解基础上，并满足特定的过滤性质得 到的，新的预条件子有着天然的并行性．我们简单分析了新的预条件子的一些性 质，在实际运用中，我们将所得的新的预条件子与传统的ILU(O)按照某种乘法 的形式结合起来使用．数值试验详细比较了这种新的预条件子与传统的预条件子 的数值效果，给出了这种预条件子的优势和缺陷． 第四章我们给出了对于求解Sylvester方程的一种预条件的梯度迭代方法，预 条件通过合理的选择两个辅助矩阵实现．这一想法可以看做为一般化线性系统的 分裂迭代到Sylvester方程中来．我们在数值试验中比较了这种迭代格式和原始 的迭代法，结果表明预条件的梯度迭代法在求解Sylvester方程时收敛得要更快， 另外我们也通过试验数值上分析了步长参数对于算法收敛的影响． 中文摘要 ii 在第五章，我们提出并且分析了对于一般鞍点问题的一种预条件子，这种预 Bai 条件子是建立矩阵分裂和最近提出的一种双参数的分裂迭代技术[Z．Z and G．H．Golub，IMA 了预条件后矩阵谱的性质，并且通过数值试验验证了我们的理论和这种预条件子 的效率． 关键词：Krylov子空间方法；位移线性系统；块三对角矩阵；Sylvester方程； 鞍点问题． 英文摘要 111 Abstract scale linear arisein asfluid Large sparse systemsmanyarea，such dynamics， structural calculationof SO on． analysis，numericalelectromagneticfields，and When thePDEs whichdescribethe obtain discreting phenomena，wegenerally linear these isofvital sparse