Eigenvalue Perturbation and Generalized Krylov Subspace Method

In this paper, we study the computational aspect of eigenvalue perturbation theory. In previous research, high order perturbation terms were often derived from Taylor series expansion. Computations based on such an approach can be both unstable and highly complicated. We present here with an approach based on the differential formulation of perturbation theory where the high order perturbation can be naturally obtained. The high order perturbation can be interpreted as a generalized Krylov subspace approximation and its convergence rate can be analyzed accordingly. This approach provides a simple and stable method to compute a few eigenvalues of a slightly modified system. 1 Introduction Perturbed eigenvalue problems are frequently encountered in physics and engineering applications. If a new system is constructed by slightly modifying from the original system, it is desirable to use the information from the original system to compute the eigensolutions of the new system. The traditio..