program4

program4 - Program 4 Numerical Linear Algebra 1 Fall 2010...

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
Program 4 Numerical Linear Algebra 1 Fall 2010 Due date: via email by 11:59PM on Friday, 11/5/10 General Task Implement the Implicit QR algorithm for symmetric tridiagonal matrices. Explore the con- vergence behavior for various sizes of n and eigevalue distributions of T . Section 8.3.5 discusses this algorithm and has an example. Pay attention to the behavior of the subdi- agonal elements in positions other than n,n - 1, i.e., watch for deflation of the problem. You do not have to exploit this in your algorithm but if you do discuss how it helps your computational complexity. There are forms of tridiagonal matrices that have known eigenvalues in addition to ex- amples you can find in various text books. For example, if T R n and the diagonals are constant, i.e., a typical row has only β,α,β as the three nonzero elements then λ k = α + 2 β cos ± n + 1 ² (verify this identity if you use it). You can also use Gershgorin’s three theorems do generate a matrix with various distributions on eigenvalues. (One is in the text on p. 320). See
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.
Ask a homework question - tutors are online