Dickson polynomial: Difference between revisions
en>Xnn removed merge tag |
en>R. J. Mathar →References: one more reference |
||
Line 1: | Line 1: | ||
A '''nonlinear eigenproblem''' is a generalization of an ordinary [[Eigenvalue, eigenvector and eigenspace|eigenproblem]] to equations that depend [[nonlinearly]] on the eigenvalue. Specifically, it refers to equations of the form: | |||
:<math>A(\lambda) \mathbf{x} = 0 , \,</math> | |||
where '''x''' is a [[vector (mathematics)|vector]] (the nonlinear "eigenvector") and ''A'' is a [[matrix (mathematics)|matrix]]-valued [[function (mathematics)|function]] of the number <math>\lambda</math> (the nonlinear "eigenvalue"). (More generally, <math>A(\lambda)</math> could be a [[linear map]], but most commonly it is a finite-dimensional, usually square, matrix.) ''A'' is usually required to be a [[holomorphic]] function of <math>\lambda</math> (in some [[domain (mathematics)|domain]]). | |||
For example, an ordinary linear eigenproblem <math>B\mathbf{v} = \lambda \mathbf{v}</math>, where ''B'' is a square matrix, corresponds to <math>A(\lambda) = B - \lambda I</math>, where ''I'' is the [[identity matrix]]. | |||
One common case is where ''A'' is a [[polynomial matrix]], which is called a '''polynomial eigenvalue problem'''. In particular, the specific case where the polynomial has [[degree of a polynomial|degree]] two is called a [[quadratic eigenvalue problem]], and can be written in the form: | |||
:<math>A(\lambda) \mathbf{x} = ( A_2 \lambda^2 + A_1 \lambda + A_0) \mathbf{x} = 0 , \,</math> | |||
in terms of the constant square matrices ''A''<sub>0,1,2</sub>. This can be converted into an ordinary linear generalized eigenproblem of twice the size by defining a new vector <math>\mathbf{y} = \lambda \mathbf{x}</math>. In terms of '''x''' and '''y''', the quadratic eigenvalue problem becomes: | |||
:<math>\begin{pmatrix} -A_0 & 0 \\ 0 & I \end{pmatrix} \begin{pmatrix} \mathbf{x} \\ \mathbf{y} \end{pmatrix} = \lambda | |||
\begin{pmatrix} A_1 & A_2 \\ I & 0 \end{pmatrix} \begin{pmatrix} \mathbf{x} \\ \mathbf{y} \end{pmatrix} | |||
, </math> | |||
where ''I'' is the identity matrix. More generally, if ''A'' is a matrix polynomial of degree ''d'', then one can convert the nonlinear eigenproblem into a linear (generalized) eigenproblem of ''d'' times the size. | |||
Besides converting them to ordinary eigenproblems, which only works if ''A'' is polynomial, there are other methods of solving nonlinear eigenproblems based on the [[Jacobi-Davidson algorithm]] or based on [[Newton's method]] (related to [[inverse iteration]]). | |||
==References== | |||
* [[Françoise Tisseur]] and Karl Meerbergen, "The quadratic eigenvalue problem," ''SIAM Review'' '''43''' (2), 235-286 (2001). | |||
* Gene H. Golub and Henk A. van der Vorst, "Eigenvalue computation in the 20th century," ''Journal of Computational and Applied Mathematics'' '''123''', 35-65 (2000). | |||
* Philippe Guillaume, "Nonlinear eigenproblems," ''SIAM J. Matrix. Anal. Appl.'' '''20''' (3), 575-595 (1999). | |||
* Axel Ruhe, "Algorithms for the nonlinear eigenvalue problem," ''SIAM Journal on Numerical Analysis'' '''10''' (4), 674-689 (1973). | |||
[[Category:Linear algebra]] |
Revision as of 18:29, 13 December 2013
A nonlinear eigenproblem is a generalization of an ordinary eigenproblem to equations that depend nonlinearly on the eigenvalue. Specifically, it refers to equations of the form:
where x is a vector (the nonlinear "eigenvector") and A is a matrix-valued function of the number (the nonlinear "eigenvalue"). (More generally, could be a linear map, but most commonly it is a finite-dimensional, usually square, matrix.) A is usually required to be a holomorphic function of (in some domain).
For example, an ordinary linear eigenproblem , where B is a square matrix, corresponds to , where I is the identity matrix.
One common case is where A is a polynomial matrix, which is called a polynomial eigenvalue problem. In particular, the specific case where the polynomial has degree two is called a quadratic eigenvalue problem, and can be written in the form:
in terms of the constant square matrices A0,1,2. This can be converted into an ordinary linear generalized eigenproblem of twice the size by defining a new vector . In terms of x and y, the quadratic eigenvalue problem becomes:
where I is the identity matrix. More generally, if A is a matrix polynomial of degree d, then one can convert the nonlinear eigenproblem into a linear (generalized) eigenproblem of d times the size.
Besides converting them to ordinary eigenproblems, which only works if A is polynomial, there are other methods of solving nonlinear eigenproblems based on the Jacobi-Davidson algorithm or based on Newton's method (related to inverse iteration).
References
- Françoise Tisseur and Karl Meerbergen, "The quadratic eigenvalue problem," SIAM Review 43 (2), 235-286 (2001).
- Gene H. Golub and Henk A. van der Vorst, "Eigenvalue computation in the 20th century," Journal of Computational and Applied Mathematics 123, 35-65 (2000).
- Philippe Guillaume, "Nonlinear eigenproblems," SIAM J. Matrix. Anal. Appl. 20 (3), 575-595 (1999).
- Axel Ruhe, "Algorithms for the nonlinear eigenvalue problem," SIAM Journal on Numerical Analysis 10 (4), 674-689 (1973).