The formulas for the perturbation analysis with modes of close eigenvalues are derived in this paper. Emphasis is made on the consistency of the straightforward perturbation process, given the complete terms of perturbations in the zeroth-order, which is a form of Rayleigh quotient, and in the higher-orders. By dividing the perturbation of eigenvector into two parts, the first-order perturbation with respect to the modes of close eigenvalues is moved into the zeroth-order perturbation. The normality condition is employed to compute the higher-order perturbations of eigenvector. The algorithm can be condensed to a single mode with a distinct eigenvalue, and this can accelerate the convergence of the perturbation analysis. The example confirms that the perturbation approximation obtained from the suggested procedure is in a good accuracy on the eigenvalues, eigenvectors, and normality.