Exercise 10.
QR-Algorithm
A much better method for approximating the eigenvalues of an invertible matrix A depends on the factor-
izatio (using the Gram-Schmidt algorithm) of A in the form
A = QR
where Q is orthogonal and R is invertible and upper triangular
The QR-algorithm
uses this repeatedly to create a sequence of matrices A? = A, A?, A?, ..., as follows:
1. Define A? = A and factor it as A? = Q?R?.
2. Define A? = R?Q? and factor it as A? = Q?R?.
3. Define A? = R?Q? and factor it as A? = Q?R?.
...
In general, A? is factored as A? = Q?R? and we define A??? = R?Q?. Then A??? is similar to A? (in fact,
A??? = R?Q? = (Q??¹A?Q?)Q?, and hence each A? has the same eigenvalues as A. If the eigenvalues of A are
real and have distinct absolute values, the remarkable thing is that the sequence of matrices A?, A?, A?, ...
converges to an upper triangular matrix with these eigenvalues on the main diagonal. [See below for the
case of complex eigenvalues.]
If A = $\begin{bmatrix} 1 & 1\\ 2 & 0 \end{bmatrix}$ use the QR-algorithm to approximate the eigenvalues.