Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Generalized eigenvectors/eigenvalues

  1. Dec 12, 2011 #1
    Mathematica has this command "Eigensystem[{m,a}]", which (to quote their documentation) "gives the generalized eigenvalues and eigenvectors of m with respect to a." I have never encountered this concept before, ever - that there can be eigenvectors of matrices with respect to other matrices. All I have ever come across is that [itex]\lambda[/itex] is a generalized eigenvalue of [itex]A[/itex] with generalized eigenvector [itex]\vec x[/itex] if there exists some [itex]p \in \mathbb N[/itex] such that [itex](A-\lambda I)^p\vec x = 0[/itex].

    Can someone please explain what it *means* to be a "generalized eigenvalue or eigenvector" of m with respect to a? Maybe it is related to the concept I mentioned above, but if so, I don't see it.
  2. jcsd
  3. Dec 13, 2011 #2
    Maybe it is answered in their documentation??
  4. Dec 13, 2011 #3
    The only thing they provide is some dodgy example: If you let

    A = \begin{pmatrix} 1. & 2. \\ 3. & 4. \end{pmatrix}, \quad B = \begin{pmatrix} 1. & 4. \\ 9. & 16. \end{pmatrix},

    (and they *HAVE* to be decimal entries; if you just put 1 (instead of 1.), you get an error, which just adds to the mystery), then "Eigensystem[{A,B}]" returns the following two complex eigenvalues and two vectors

    &\{0.25 + 0.193649 I,
    0.25 - 0.193649 I\}, \\
    &[-0.848472 + 0.0858378 I, 0.498043 + 0.157099 I],\\
    &[-0.848472 - 0.0858378 I, 0.498043 - 0.157099 I]

    Last edited by a moderator: Dec 13, 2011
  5. Dec 13, 2011 #4


    User Avatar
    Homework Helper

    If I recall correctly a generalized eigensystem is one for which

    [tex]Ax = \lambda Bx[/tex]

    for given matrices A and B.

    EDIT: See here.
  6. Dec 13, 2011 #5

    D H

    User Avatar
    Staff Emeritus
    Science Advisor

  7. Dec 13, 2011 #6


    User Avatar
    Science Advisor
    Homework Helper

    A practical application of all this is the vibration of a flexible structure, where A and B represent the stiffness and mass properties. If the matrices are large and sparse (and often also Hermitian and positive semi-definite), solving the generalized eigenproblem is much more efficient than solving the equivalent problem [itex]B^{-1}Ax = \lambda x[/itex], (assuming B is invertible) because [itex]B^{-1}A[/itex] looks like an arbitrary full non-symmetric matrix with no obvious "special properties" to leading to a more efficient solution.

    FWIW there are solution procedures that represent [itex]\lambda[/itex] as a ratio of two numbers, with conventions to represent the "indeterminate" or "infinte" eigenvalues and their corresponding vectors. The vectors are well defined and meaningful as the basis vectors of subspaces, even if the corresponding eigenvalues are not so well defined.
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook