Dickson polynomial: Difference between revisions

From formulasearchengine
Jump to navigation Jump to search
en>Xnn
removed merge tag
 
en>R. J. Mathar
References: one more reference
Line 1: Line 1:
Hi there, my name is Leah. Participating in region audio is anything I will hardly ever give up. My [http://www.alexa.com/search?q=position&r=topsites_index&p=bigtop position] is a travel agent. [http://www.Britannica.com/search?query=Michigan Michigan] is wherever we have been dwelling for years. If you want to find out additional check out out my web page: http://www.deboeck-incoming.com/aspnet_client/goedkoop/nike-air-max-one-143285250.htm<br><br>
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:


Here is my blog post - [http://www.deboeck-incoming.com/aspnet_client/goedkoop/nike-air-max-one-143285250.htm nike air max one]
:<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:

A(λ)x=0,

where x is a vector (the nonlinear "eigenvector") and A is a matrix-valued function of the number λ (the nonlinear "eigenvalue"). (More generally, A(λ) 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 Bv=λv, where B is a square matrix, corresponds to A(λ)=BλI, 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:

A(λ)x=(A2λ2+A1λ+A0)x=0,

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 y=λx. In terms of x and y, the quadratic eigenvalue problem becomes:

(A000I)(xy)=λ(A1A2I0)(xy),

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).