Steklov 特征值问题基于移位反迭代的多水平校 正方案-计算数学专业论文.docxVIP

I万方数据 I 万方数据 摘 要 Steklov 特征值问题有很强的物理背景及广泛的应用, 尤其在流体力学方面, 因此 越来越多的学者开始关注 Steklov 特征值问题的有限元方法. Lin 和 Xie 介绍了一种基 于一步校正的多水平方法. 移位反迭代是解矩阵特征值的基本方法, Yang 和 Bi 建立了 基于移位反迭代的二网格方法. 组合 Lin 和 Xie 的校正技巧和移位反迭代, 本文提出了 Steklov 特征值问题基于移位反迭代的多水平校正方案, 并给出了严格的理论分析. 这个 多水平方案包括了在越来越细的有限元空间上用移位反迭代来求解边值问题，即用预 先得到的特征值作为位移, 特征函数作为迭代初值, 以及在粗有限元空间上求解 Steklov 特征值问题. 可以证明, 通过每一个一步校正之后, 都能够提高预先给定近似特征对的精 度. 最后, 用 Matlab 在单位方形区域和 L-型区域上给出数值实验来验证我们所提出的方 案是有效的. 此外, 随着迭代次数的增加，近似特征值越来越接近准确特征值，这导致 我们求解的方程越来越接近奇异. 虽然数值求解上没有困难, 但是在一步校正中可以固 定位移, 从而得到一个新的方案, 并且给出了相应的数值实验. 理论分析和数值实验都表 明本文提出的方案不仅适用于单特征值, 而且也适用于重特征值. 关键词: Steklov 特征值问题; 多水平方案; 移位反迭代; 一步校正. II II 万方数据 Abstract Steklov eigenvalue problems have important physical background and wide applica- tions especially in fluid, thus finite element methods for Steklov eigenvalue problems have attracted the attention of more and more scholars. A multilevel method based on the cor- rection step was proposed by Lin and Xie. The shifted inverse iteration method is a basic approach for solving matrix eigenvalue problems. Now, two-grid methods based on the shifted inverse iteration have been established by Yang and Bi. Combining the correction technique proposed by Lin and Xie and the shifted inverse iteration, a multilevel correc- tion scheme for the Steklov eigenvalue problem is proposed in this paper and the rigorous theoretical analysis is given. This type of multilevel method includes solving boundary value problems by the shifted inverse iteration in finer and finer finite element spaces , i.e., using the previously obtained eigenvalue as the shift and eigenfunction as the itera- tion initial value, and Steklov eigenvalue problems in coarse finite element spaces. And it is proved that the accuracy of the resulting eigenpair approximation can be improved after each correction step. Finally, numerical experiments with Matlab on the unit square and L-shape domain are reported to validate the efficiency of the proposed method. In addition, with the iteratio