Jacobi Eigenvalue Algorithm - Description

Description

Let S be a symmetric matrix, and G = G(i,j,θ) be a Givens rotation matrix. Then:

is symmetric and similar to S.

Furthermore, S′ has entries:

\begin{align} S'_{ii} &= c^2\, S_{ii} - 2\, s c \,S_{ij} + s^2\, S_{jj} \\ S'_{jj} &= s^2 \,S_{ii} + 2 s c\, S_{ij} + c^2 \, S_{jj} \\ S'_{ij} &= S'_{ji} = (c^2 - s^2 ) \, S_{ij} + s c \, (S_{ii} - S_{jj} ) \\ S'_{ik} &= S'_{ki} = c \, S_{ik} - s \, S_{jk} & k \ne i,j \\ S'_{jk} &= S'_{kj} = s \, S_{ik} + c \, S_{jk} & k \ne i,j \\ S'_{kl} &= S_{kl} &k,l \ne i,j
\end{align}

where s = sin(θ) and c = cos(θ).

Since G is orthogonal, S and S′ have the same Frobenius norm ||·||F (the square-root sum of squares of all components), however we can choose θ such that S′ij = 0, in which case S′ has a larger sum of squares on the diagonal:

Set this equal to 0, and rearrange:

if

In order to optimize this effect, Sij should be the largest off-diagonal component, called the pivot.

The Jacobi eigenvalue method repeatedly performs rotations until the matrix becomes almost diagonal. Then the elements in the diagonal are approximations of the (real) eigenvalues of S.

Read more about this topic:  Jacobi Eigenvalue Algorithm

Famous quotes containing the word description:

    “The next Augustan age will dawn on the other side of the Atlantic. There will, perhaps, be a Thucydides at Boston, a Xenophon at New York, and, in time, a Virgil at Mexico, and a Newton at Peru. At last, some curious traveller from Lima will visit England and give a description of the ruins of St. Paul’s, like the editions of Balbec and Palmyra.”
    —Horace Walpole (1717–1797)

    “It is possible—indeed possible even according to the old conception of logic—to give in advance a description of all ‘true’ logical propositions. Hence there can never be surprises in logic.”
    —Ludwig Wittgenstein (1889–1951)

    “It [Egypt] has more wonders in it than any other country in the world and provides more works that defy description than any other place.”
    —Herodotus (c. 484–424 B.C.)