Loading [MathJax]/jax/output/HTML-CSS/jax.js

 

 

 

Eigenvalues and Lanczos' method

Basic features with a real symmetric matrix (and normally huge n>106 and sparse) ˆA of dimension n×n:

  • Lanczos' algorithm generates a sequence of real tridiagonal matrices Tk of dimension k×k with kn, with the property that the extremal eigenvalues of Tk are progressively better estimates of ˆA' extremal eigenvalues.* The method converges to the extremal eigenvalues.
  • The similarity transformation is
ˆT=ˆQTˆAˆQ, with the first vector ˆQˆe1=ˆq1.

We are going to solve iteratively ˆT=ˆQTˆAˆQ, with the first vector ˆQˆe1=ˆq1. We can write out the matrix ˆQ in terms of its column vectors ˆQ=[ˆq1ˆq2ˆqn].