Take the 2-minute tour ×
MathOverflow is a question and answer site for professional mathematicians. It's 100% free, no registration required.

Suppose I have the symmetric tridiagonal matrix:

$ \begin{pmatrix} a & b_{1} & 0 & ... & 0 \\\ b_{1} & a & b_{2} & & ... \\\ 0 & b_{2} & a & ... & 0 \\\ ... & & ... & & b_{n-1} \\\ 0 & ... & 0 & b_{n-1} & a \end{pmatrix} $

All of the entries can be taken to be positive real numbers and all of the $a_{i}$ are equal. I know that when the $b_{i}$'s are equal (the matrix is uniform), there are closed-form expressions for the eigenvalues and eigenvectors in terms of cosine and sine functions. Additionally, I know of the recurrence relation:

$det(A_{n}) = a\cdot det(A_{n-1}) - b_{n-1}^{2}\cdot det(A_{n-2})$

Additionally, since my matrix is real-symmetric, I know that its eigenvalues are real.

Is there anything else I can determine about the eigenvalues? Furthermore, is there a closed-form expression for them?

share|improve this question
    
you might as well take $a=0$, since this is just an additive constant for each eigenvalue –  Carlo Beenakker May 22 '13 at 23:30
    
There are older questions on MO dealing with eigenvalues of symmetric tridiagonal; no closed form. You might want to ask on the scicomp stackexchange. –  Suvrit May 23 '13 at 0:02

3 Answers 3

up vote 4 down vote accepted

The type of matrix you have written down is called Jacobi matrix and people are still discovering new things about them basically their properties fill entire bookcases at a mathematics library. One of the reasons is the connection to orthogonal polynomials. Basically, if $\{p_n(x)\}_{n\geq 0}$ is a family of orthogonal polynomials, then they obey a recursion relation of the form $$ b_n p_{n+1}(x) + (a_n- x) p_n(x) + b_{n-1} p_{n-1}(x) = 0. $$ You should be able to recognize the form of your matrix from this.

As far as general properties of the eigenvalues, let me mention two:

  1. The eigenvalues are simple. In fact one has $\lambda_j - \lambda_{j-1} \geq e^{-c n}$, where $c$ is some constant that depends on the $b_j$.

  2. The eigenvalues of $A$ and $A_{n-1}$ interlace.

share|improve this answer
    
Do you have any references for the first property? Thanks for your help! –  FlamingWilderbeest May 23 '13 at 13:58
    
No. The proofs are simple enough that I don't know where they can be found. math.caltech.edu/Szego.html might contain these things. The main reason I'm hesitating posting a proof is that most proofs are heavy on notation, e.g. introduce orthogonal polynomials. –  Helge May 23 '13 at 17:11

Amongst the polynomials that can arise as characteristic polynomials of tridiagonal matrices with zero diagonal, one finds the Hermite polynomials. Schur showed that Hermite polynomials of even degree are irreducible and that their Galois groups are not solvable. Hence there can be no closed form expression for the zeros in terms of the $b_i$'s in general.

share|improve this answer
    
I wasn't aware of this. Thanks a lot! –  FlamingWilderbeest May 23 '13 at 1:28
1  
Another argument is: it's easy to reduce any symmetric matrix to tridiagonal with similarity transforms. So if this problem were easy to solve, all symmetric eigenproblems would be. –  Federico Poloni May 23 '13 at 6:40
    
Well, I didn't mean to imply that this would be easy to solve. It is too hopeful to expect a closed-form expression for the eigenvalues in terms of the bi's. However, knowing anything about what they look like would be tremendously useful to me –  FlamingWilderbeest May 23 '13 at 14:57
    
Closed form expression built from arithmetic operations and radicals. When all the $b_i$ are equal we also do not get a closed form expressions of this type or do we? –  Vít Tuček Jun 17 at 1:22

Mathematica gives to you the closed form that you want. All you have to do is use de recurrence package of the program

share|improve this answer
2  
Chris Godsil's answer seems to suggest otherwise –  Yemon Choi Jun 16 at 22:05

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.