Rayleigh quotient iteration is an eigenvalue algorithm which extends the idea of the inverse iteration by using the Rayleigh quotient to obtain increasingly accurate eigenvalue estimates.
Contents
Rayleigh quotient iteration is an iterative method, that is, it must be repeated until it converges to an answer (this is true for all eigenvalue algorithms). Fortunately, very rapid convergence is guaranteed and no more than a few iterations are needed in practice. The Rayleigh quotient iteration algorithm converges cubically for Hermitian or symmetric matrices, given an initial vector that is sufficiently close to an eigenvector of the matrix that is being analyzed.
Algorithm
The algorithm is very similar to inverse iteration, but replaces the estimated eigenvalue at the end of each iteration with the Rayleigh quotient. Begin by choosing some value
Calculate the next approximation of the eigenvector
where
To compute more than one eigenvalue, the algorithm can be combined with a deflation technique.
Note that for very small problems it is beneficial to replace the matrix inverse with the adjugate, which will yield the same iteration because it's equal to the inverse up to an irrelevant scale (the inverse of the determinant, specifically). The adjugate is easier to compute explicitly than the inverse (though the inverse is easier to apply to a vector for problems that aren't small), and is more numerically sound because it remains well defined as the eigenvalue converges.
Example
Consider the matrix
for which the exact eigenvalues are
(where
The largest eigenvalue is
We begin with an initial eigenvalue guess of
Then, the first iteration yields
the second iteration,
and the third,
from which the cubic convergence is evident.
Octave Implementation
The following is a simple implementation of the algorithm in Octave.