Eigenvectors and eigenvalues of nonsymmetric Tridiagonal matrix - MathOverflow most recent 30 from http://mathoverflow.net 2013-05-18T13:04:20Z http://mathoverflow.net/feeds/question/105225 http://www.creativecommons.org/licenses/by-nc/2.5/rdf http://mathoverflow.net/questions/105225/eigenvectors-and-eigenvalues-of-nonsymmetric-tridiagonal-matrix Eigenvectors and eigenvalues of nonsymmetric Tridiagonal matrix Cherep 2012-08-22T10:57:15Z 2012-08-26T17:15:43Z <p>Hi, the question is following: We have one matrix $$\begin{pmatrix} -\beta &amp; \Delta &amp; 0 &amp; 0 &amp;\cdots &amp; 0 &amp; 0 &amp; 0 \newline \beta &amp; -(\beta+\Delta) &amp; \Delta &amp; 0 &amp;\cdots &amp; 0 &amp; 0 &amp; 0 \newline 0 &amp; \beta &amp; -(\beta + \Delta) &amp; \Delta &amp;\cdots &amp; 0 &amp; 0 &amp; 0 \newline \vdots &amp; \vdots &amp; \vdots &amp; \vdots &amp; \ddots &amp; \vdots &amp; \vdots&amp; \vdots \newline 0 &amp; 0 &amp; 0 &amp; 0 &amp;\cdots &amp; \beta &amp; -(\beta + \Delta) &amp; \Delta \newline 0 &amp; 0 &amp; 0 &amp; 0 &amp; \cdots &amp; 0 &amp; \beta &amp; -(\beta + \Delta) \end{pmatrix}.$$</p> <p>Is it possible to find analytically all eigenvalues and eigenvectors of this matrix? What I found <a href="http://mathoverflow.net/questions/91161/eigenvectors-and-eigenvalues-of-tridiagonal-matrix" rel="nofollow">here</a> is quite similar, but not exactly the same. The element in the last row and last column $-(\beta +\Delta)$ can be replaced by $-\Delta$ if it simplifies the solution. $\beta>0$, $\Delta>0$.</p> <p>Thanks</p> http://mathoverflow.net/questions/105225/eigenvectors-and-eigenvalues-of-nonsymmetric-tridiagonal-matrix/105256#105256 Answer by Robert Israel for Eigenvectors and eigenvalues of nonsymmetric Tridiagonal matrix Robert Israel 2012-08-22T17:59:57Z 2012-08-26T17:15:43Z <p>As Denis remarked, the characteristic polynomial $P_n(\lambda)$ of the $n \times n$ matrix does satisfy a three-term linear recurrence $$P_{n+2}(\lambda) = (\lambda + \beta + \Delta) P_{n+1}(\lambda) - \beta \Delta P_n(\lambda)$$ with initial conditions $P_1(\lambda) = \lambda+\beta$, $P_2(\lambda) = \lambda^2+(2 \beta+\Delta)\lambda +\beta^2$, and then $$P_n(\lambda) = \frac{(P_2 - r_2 P_1)}{r_1 (r_1 - r_2)} r_1^n + \frac{(P_2 - r_1 P_1)}{r_2 (r_2 - r_1)} r_2^n$$ where $r_1$ and $r_2$ are the roots of $r^2 - (\lambda + \beta + \Delta) r + \beta \Delta$. However, I don't see how this implies a "closed form" for the eigenvalues. </p> <p>EDIT: By scaling, we may as well assume $\Delta=1$. I don't know if this leads to a "closed form", but the eigenvalue $\lambda_0$ that is analytic in a neighbourhood of $\beta=0$ has some interesting regularities in its Maclaurin series:</p> <p>$$\matrix{n=2 &amp; -\beta^{2}+2 \beta^{3}-5 \beta^{4}+14 \beta^{5}-42 \beta^{6}\cr&amp;+132 \beta^{7}-429 \beta^{8}+1430 \beta^{9}-4862 \beta^{10}\cr n=3 &amp; -\beta^{3}+2 \beta^{4}-\beta^{5}-6 \beta^{6}+20 \beta^{7}\cr&amp;-22 \beta^{8}-49 \beta^{9}+260 \beta^{10}-441 \beta^{11}\cr&amp;-320 \beta^{12}+3652 \beta^{13}\cr n=4 &amp; -\beta^{4}+2 \beta^{5}-\beta^{6}-8 \beta^{8}+26 \beta^{9}\cr&amp;-28 \beta^{10}+10 \beta^{11}-100 \beta^{12}+442 \beta^{13}\cr&amp;-729 \beta^{14}+532 \beta^{15}-1641 \beta^{16}\cr n=5 &amp; -\beta^{5}+2 \beta^{6}-\beta^{7}-10 \beta^{10}+32 \beta^{11}\cr&amp;-34 \beta^{12}+12 \beta^{13}-155 \beta^{15}+672 \beta^{16}\cr&amp;-1089 \beta^{17}+782 \beta^{18}-210 \beta^{19}\cr n=6 &amp; -\beta^{6}+2 \beta^{7}-\beta^{8}-12 \beta^{12}+38 \beta^{13}\cr&amp;-40 \beta^{14}+14 \beta^{15}-222 \beta^{18}+950 \beta^{19}\cr&amp;-1521 \beta^{20}+1080 \beta^{21}-287 \beta^{22}\cr n=7 &amp; -\beta^{7}+2 \beta^{8}-\beta^{9}-14 \beta^{14}+44 \beta^{15}\cr&amp;-46 \beta^{16}+16 \beta^{17}-301 \beta^{21}+1276 \beta^{22}\cr&amp;-2025 \beta^{23}+1426 \beta^{24}-376 \beta^{25}\cr n=8 &amp; -\beta^{8}+2 \beta^{9}-\beta^{10}-16 \beta^{16}+50 \beta^{17}\cr&amp;-52 \beta^{18}+18 \beta^{19}-392 \beta^{24}+1650 \beta^{25}\cr&amp;-2601 \beta^{26}+1820 \beta^{27}-477 \beta^{28}\cr n=9 &amp; -\beta^{9}+2 \beta^{10}-\beta^{11}-18 \beta^{18}+56 \beta^{19}\cr&amp;-58 \beta^{20}+20 \beta^{21}-495 \beta^{27}+2072 \beta^{28}\cr&amp;-3249 \beta^{29}+2262 \beta^{30}-590 \beta^{31}\cr n=10 &amp; -\beta^{10}+2 \beta^{11}-\beta^{12}-20 \beta^{20}+62 \beta^{21}\cr&amp;-64 \beta^{22}+22 \beta^{23}-610 \beta^{30}+2542 \beta^{31}\cr&amp;-3969 \beta^{32}+2752 \beta^{33}-715 \beta^{34}\cr}$$</p> <p>All coefficients are integers, and it starts with $-\beta^n + 2 \beta^{n+1} - \beta^{n+2} - 2 n \beta^{2n} + (6n+2) \beta^{2n+1} - (6n+4) \beta^{2n+2} + (2n+2) \beta^{2n+3} - (6n+1)n \beta^{3n} \ldots$</p> http://mathoverflow.net/questions/105225/eigenvectors-and-eigenvalues-of-nonsymmetric-tridiagonal-matrix/105398#105398 Answer by Felix Goldberg for Eigenvectors and eigenvalues of nonsymmetric Tridiagonal matrix Felix Goldberg 2012-08-24T17:02:19Z 2012-08-24T17:02:19Z <p>There are such formulas, but beware: they are rather intimidating.</p> <p>Look at this paper:</p> <p>"EXPLICIT EIGENVALUES AND INVERSES OF TRIDIAGONAL TOEPLITZ MATRICES WITH FOUR PERTURBED CORNERS" by WEN-CHYUAN YUEH and SUI SUN CHENG2, ANZIAM J. 49(2008), 361–387.</p> http://mathoverflow.net/questions/105225/eigenvectors-and-eigenvalues-of-nonsymmetric-tridiagonal-matrix/105428#105428 Answer by Tristram Shandy for Eigenvectors and eigenvalues of nonsymmetric Tridiagonal matrix Tristram Shandy 2012-08-24T22:15:54Z 2012-08-24T23:48:41Z <p>Various commenters have pointed out that one can easily get asymptotic estimates for the eigenvalues as $n \rightarrow \infty$ (specifically, with $\beta, \Delta > 0$ they converge to $-2 (\beta + \Delta) \cos(\pi k/2 n)^2$ for $k = 1$ to $n$). OTOH, you specify in the comments that you require the eigenvalues in exact form. This is not possible. Even for $n = 5$, $\beta = 2$, and $\Delta = 1$, the characteristic polynomial of the matrix above is $$x^5 + 14 x^4 + 70 x^3 + 150 x^2 + 129 x + 32$$ whose splitting field has Galois group $S_5$, and hence there cannot be any exact formula (in radicals) for the roots by a theorem of Abel. This suggests that it is highly unlikely that there exists any "exact" formula for the roots, although "exact formula" is a somewhat nebulous notion, since of course one can <i> define </i> a function $f_{n,k}(\beta,\Delta)$ to be the $k$th smallest root of the corresponding characteristic polynomial.</p>