In mathematics, the quadratic eigenvalue problem (QEP), is to find scalar eigenvalues λ {\displaystyle \lambda } , left eigenvectors y {\displaystyle y} and right eigenvectors x {\displaystyle x} such that
Q ( λ ) x = 0 and y ∗ Q ( λ ) = 0 , {\displaystyle Q(\lambda )x=0~{\text{ and }}~y^{\ast }Q(\lambda )=0,}
where Q ( λ ) = λ 2 M + λ C + K {\displaystyle Q(\lambda )=\lambda ^{2}M+\lambda C+K} , with matrix coefficients M , C , K ∈ C n × n {\displaystyle M,\,C,K\in \mathbb {C} ^{n\times n}} and we require that M ≠ 0 {\displaystyle M\,\neq 0} , (so that we have a nonzero leading coefficient). There are 2 n {\displaystyle 2n} eigenvalues that may be infinite or finite, and possibly zero. This is a special case of a nonlinear eigenproblem. Q ( λ ) {\displaystyle Q(\lambda )} is also known as a quadratic polynomial matrix.
Spectral theory A QEP is said to be regular if det ( Q ( λ ) ) ≢ 0 {\displaystyle {\text{det}}(Q(\lambda ))\not \equiv 0} identically. The coefficient of the λ 2 n {\displaystyle \lambda ^{2n}} term in det ( Q ( λ ) ) {\displaystyle {\text{det}}(Q(\lambda ))} is det ( M ) {\displaystyle {\text{det}}(M)} , implying that the QEP is regular if M {\displaystyle M} is nonsingular. Eigenvalues at infinity and eigenvalues at 0 may be exchanged by considering the reversed polynomial, λ 2 Q ( λ − 1 ) = λ 2 K + λ C + M {\displaystyle \lambda ^{2}Q(\lambda ^{-1})=\lambda ^{2}K+\lambda C+M} . As there are 2 n {\displaystyle 2n} eigenvectors in a n {\displaystyle n} dimensional space, the eigenvectors cannot be orthogonal. It is possible to have the same eigenvector attached to different eigenvalues.
Applications
Systems of differential equations Quadratic eigenvalue problems arise naturally in the solution of systems of second order linear differential equations without forcing:
M q ″ ( t ) + C q ′ ( t ) + K q ( t ) = 0 {\displaystyle Mq''(t)+Cq'(t)+Kq(t)=0}
Where q ( t ) ∈ R n {\displaystyle q(t)\in \mathbb {R} ^{n}} , and M , C , K ∈ R n × n {\displaystyle M,C,K\in \mathbb {R} ^{n\times n}} . If all quadratic eigenvalues of Q ( λ ) = λ 2 M + λ C + K {\displaystyle Q(\lambda )=\lambda ^{2}M+\lambda C+K} are distinct, then the solution can be written in terms of the quadratic eigenvalues and right quadratic eigenvectors as
q ( t ) = ∑ j = 1 2 n α j x j e λ j t = X e Λ t α {\displaystyle q(t)=\sum _{j=1}^{2n}\alpha _{j}x_{j}e^{\lambda _{j}t}=Xe^{\Lambda t}\alpha }
… excerpt ends here. Continue reading the full article.
