Hereditary property

In the mathematical discipline of linear algebra, eigendecomposition or sometimes spectral decomposition is the factorization of a matrix into a canonical form, whereby the matrix is represented in terms of its eigenvalues and eigenvectors. Only diagonalizable matrices can be factorized in this way.

Fundamental theory of matrix eigenvectors and eigenvalues

Mining Engineer (Excluding Oil ) Truman from Alma, loves to spend time knotting, largest property developers in singapore developers in singapore and stamp collecting. Recently had a family visit to Urnes Stave Church.

A (non-zero) vector v of dimension N is an eigenvector of a square (N×N) matrix A if and only if it satisfies the linear equation

$\mathbf{A}\mathbf{v}=\lambda \mathbf{v}$

where λ is a scalar, termed the eigenvalue corresponding to v. That is, the eigenvectors are the vectors which the linear transformation A merely elongates or shrinks, and the amount that they elongate/shrink by is the eigenvalue. The above equation is called the eigenvalue equation or the eigenvalue problem.

This yields an equation for the eigenvalues

$p\left(\lambda \right):=\det (\mathbf{A}-\lambda \mathbf{I})=0.$

We call p(λ) the characteristic polynomial, and the equation, called the characteristic equation, is an Nth order polynomial equation in the unknown λ. This equation will have N_λ distinct solutions, where 1 ≤ N_λ ≤ N . The set of solutions, i.e. the eigenvalues, is sometimes called the spectrum of A.

We can factor p as

$p\left(\lambda \right)=(\lambda -{\lambda }_1{)}^{n_1}(\lambda -{\lambda }_2{)}^{n_2}\cdots (\lambda -{\lambda }_k{)}^{n_k}=0.$

The integer n_i is termed the algebraic multiplicity of eigenvalue λ_i. The algebraic multiplicities sum to N:

$\sum _{i=1}^{N_{\lambda }}{n}_i=N.$

For each eigenvalue, λ_i, we have a specific eigenvalue equation

$(\mathbf{A}-{\lambda }_i\mathbf{I})\mathbf{v}=0.$

There will be 1 ≤ m_i ≤ n_i linearly independent solutions to each eigenvalue equation. The m_i solutions are the eigenvectors associated with the eigenvalue λ_i. The integer m_i is termed the geometric multiplicity of λ_i. It is important to keep in mind that the algebraic multiplicity n_i and geometric multiplicity m_i may or may not be equal, but we always have m_i ≤ n_i. The simplest case is of course when m_i = n_i = 1. The total number of linearly independent eigenvectors, N_v, can be calculated by summing the geometric multiplicities

$\sum _{i=1}^{N_{\lambda }}{m}_i=N_{\mathbf{v}}.$

The eigenvectors can be indexed by eigenvalues, i.e. using a double index, with v_i,j being the j^th eigenvector for the i^th eigenvalue. The eigenvectors can also be indexed using the simpler notation of a single index v_k, with k = 1, 2, ..., N_v.

Eigendecomposition of a matrix

Let A be a square (N×N) matrix with N linearly independent eigenvectors, $q_i(i=1,\dots ,N).$ Then A can be factorized as

$\mathbf{A}=\mathbf{Q}\mathrm{\Lambda }{\mathbf{Q}}^{-1}$

where Q is the square (N×N) matrix whose i^th column is the eigenvector $q_i$ of A and Λ is the diagonal matrix whose diagonal elements are the corresponding eigenvalues, i.e., ${\mathrm{\Lambda }}_{ii}={\lambda }_i$. Note that only diagonalizable matrices can be factorized in this way. For example, the defective matrix $\left(\begin{array}{cc}1& 1\\ 0& 1\\ \end{array}\right)$ cannot be diagonalized.

The eigenvectors $q_i(i=1,\dots ,N)$ are usually normalized, but they need not be. A non-normalized set of eigenvectors, $v_i(i=1,\dots ,N),$ can also be used as the columns of Q. That can be understood by noting that the magnitude of the eigenvectors in Q gets canceled in the decomposition by the presence of Q⁻¹.

Example

Taking a 2 × 2 real matrix $\mathbf{A}=\left[\begin{array}{cc}1& 0\\ 1& 3\\ \end{array}\right]$ as an example to be decomposed into a diagonal matrix through multiplication of a non-singular matrix $\mathbf{B}=\left[\begin{array}{cc}a& b\\ c& d\\ \end{array}\right],[a,b,c,d]\in \mathbb{ℝ}$.

Then

${\left[\begin{array}{cc}a& b\\ c& d\\ \end{array}\right]}^{-1}\left[\begin{array}{cc}1& 0\\ 1& 3\\ \end{array}\right]\left[\begin{array}{cc}a& b\\ c& d\\ \end{array}\right]=\left[\begin{array}{cc}x& 0\\ 0& y\\ \end{array}\right]$, for some real diagonal matrix $\left[\begin{array}{cc}x& 0\\ 0& y\\ \end{array}\right]$.

Shifting $\mathbf{B}$ to the right hand side:

$\left[\begin{array}{cc}1& 0\\ 1& 3\\ \end{array}\right]\left[\begin{array}{cc}a& b\\ c& d\\ \end{array}\right]=\left[\begin{array}{cc}a& b\\ c& d\\ \end{array}\right]\left[\begin{array}{cc}x& 0\\ 0& y\\ \end{array}\right]$

The above equation can be decomposed into 2 simultaneous equations:

$\{\begin{array}{l}\left[\begin{array}{cc}1& 0\\ 1& 3\end{array}\right]\left[\begin{array}{c}a\\ c\end{array}\right]=\left[\begin{array}{c}ax\\ cx\end{array}\right]\\ \left[\begin{array}{cc}1& 0\\ 1& 3\end{array}\right]\left[\begin{array}{c}b\\ d\end{array}\right]=\left[\begin{array}{c}by\\ dy\end{array}\right]\end{array}$

Factoring out the eigenvalues $x$ and $y$:

$\{\begin{array}{l}\left[\begin{array}{cc}1& 0\\ 1& 3\end{array}\right]\left[\begin{array}{c}a\\ c\end{array}\right]=x\left[\begin{array}{c}a\\ c\end{array}\right]\\ \left[\begin{array}{cc}1& 0\\ 1& 3\end{array}\right]\left[\begin{array}{c}b\\ d\end{array}\right]=y\left[\begin{array}{c}b\\ d\end{array}\right]\end{array}$

Letting $\overrightarrow{a}=\left[\begin{array}{c}a\\ c\end{array}\right],\overrightarrow{b}=\left[\begin{array}{c}b\\ d\end{array}\right]$, this gives us two vector equations:

$\{\begin{array}{l}A\overrightarrow{a}=x\overrightarrow{a}\\ A\overrightarrow{b}=y\overrightarrow{b}\end{array}$

And can be represented by a single vector equation involving 2 solutions as eigenvalues:

$\mathbf{A}\mathbf{u}=\lambda \mathbf{u}$

where $\lambda$ represents the two eigenvalues $x$ and $y$, $\mathbf{u}$ represents the vectors $\overrightarrow{a}$ and $\overrightarrow{b}$.

Shifting $\lambda \mathbf{u}$ to the left hand side and factorizing $\mathbf{u}$ out

$(\mathbf{A}-\lambda \mathbf{I})\mathbf{u}=0$

Since $\mathbf{B}$ is non-singular, it is essential that $\mathbf{u}$ is non-zero. Therefore,

$(\mathbf{A}-\lambda \mathbf{I})=\mathbf{0}$

Considering the determinant of $(\mathbf{A}-\lambda \mathbf{I})$,

$\left[\begin{array}{cc}1-\lambda & 0\\ 1& 3-\lambda \end{array}\right]=0$

Thus

$(1-\lambda )(3-\lambda )=0$

Giving us the solutions of the eigenvalues for the matrix $\mathbf{A}$ as $\lambda =1$ or $\lambda =3$, and the resulting diagonal matrix from the eigendecomposition of $\mathbf{A}$ is thus $\left[\begin{array}{cc}1& 0\\ 0& 3\end{array}\right]$.

Putting the solutions back into the above simultaneous equations

$\{\begin{array}{l}\left[\begin{array}{cc}1& 0\\ 1& 3\end{array}\right]\left[\begin{array}{c}a\\ c\end{array}\right]=1\left[\begin{array}{c}a\\ c\end{array}\right]\\ \left[\begin{array}{cc}1& 0\\ 1& 3\end{array}\right]\left[\begin{array}{c}b\\ d\end{array}\right]=3\left[\begin{array}{c}b\\ d\end{array}\right]\end{array}$

Solving the equations, we have $a=-2c,a\in \mathbb{ℝ}$ and $b=0,d\in \mathbb{ℝ}$

Thus the matrix $\mathbf{B}$ required for the eigendecomposition of $\mathbf{A}$ is $\left[\begin{array}{cc}-2c& 0\\ c& d\end{array}\right],[c,d]\in \mathbb{ℝ}$. i.e. :

${\left[\begin{array}{cc}-2c& 0\\ c& d\\ \end{array}\right]}^{-1}\left[\begin{array}{cc}1& 0\\ 1& 3\\ \end{array}\right]\left[\begin{array}{cc}-2c& 0\\ c& d\\ \end{array}\right]=\left[\begin{array}{cc}1& 0\\ 0& 3\\ \end{array}\right],[c,d]\in \mathbb{ℝ}$

Matrix inverse via eigendecomposition

Mining Engineer (Excluding Oil ) Truman from Alma, loves to spend time knotting, largest property developers in singapore developers in singapore and stamp collecting. Recently had a family visit to Urnes Stave Church.

If matrix A can be eigendecomposed and if none of its eigenvalues are zero, then A is nonsingular and its inverse is given by

${\mathbf{A}}^{-1}=\mathbf{Q}{\mathrm{\Lambda }}^{-1}{\mathbf{Q}}^{-1}$

Furthermore, because Λ is a diagonal matrix, its inverse is easy to calculate:

${\left[{\mathrm{\Lambda }}^{-1}\right]}_{ii}=\frac{1}{{\lambda }_i}$

Practical implications

When eigendecomposition is used on a matrix of measured, real data, the inverse may be less valid when all eigenvalues are used unmodified in the form above. This is because as eigenvalues become relatively small, their contribution to the inversion is large. Those near zero or at the "noise" of the measurement system will have undue influence and could hamper solutions (detection) using the inverse.

Two mitigations have been proposed: 1) truncating small/zero eigenvalues, 2) extending the lowest reliable eigenvalue to those below it.

The first mitigation method is similar to a sparse sample of the original matrix, removing components that are not considered valuable. However, if the solution or detection process is near the noise level, truncating may remove components that influence the desired solution.

The second mitigation extends the eigenvalue so that lower values have much less influence over inversion, but do still contribute, such that solutions near the noise will still be found.

The reliable eigenvalue can be found by assuming that eigenvalues of extremely similar and low value are a good representation of measurement noise (which is assumed low for most systems).

If the eigenvalues are rank-sorted by value, then the reliable eigenvalue can be found by minimization of the Laplacian of the sorted eigenvalues:

$\min \left|{\mathrm{\nabla }}^2{\lambda }_s\right|$

where the eigenvalues are subscripted with an 's' to denote being sorted. The position of the minimization is the lowest reliable eigenvalue. In measurement systems, the square root of this reliable eigenvalue is the average noise over the components of the system.

Functional calculus

The eigendecomposition allows for much easier computation of power series of matrices. If f(x) is given by

$f(x)=a_0+a_{1}x+a_2{x}^2+\cdots$

then we know that

$f\left(\mathbf{A}\right)=\mathbf{Q}f\left(\mathrm{\Lambda }\right){\mathbf{Q}}^{-1}$

Because Λ is a diagonal matrix, functions of Λ are very easy to calculate:

${\left[f\left(\mathrm{\Lambda }\right)\right]}_{ii}=f\left({\lambda }_i\right)$

The off-diagonal elements of f(Λ) are zero; that is, f(Λ) is also a diagonal matrix. Therefore, calculating f(A) reduces to just calculating the function on each of the eigenvalues .

A similar technique works more generally with the holomorphic functional calculus, using

${\mathbf{A}}^{-1}=\mathbf{Q}{\mathrm{\Lambda }}^{-1}{\mathbf{Q}}^{-1}$

from above. Once again, we find that

${\left[f\left(\mathrm{\Lambda }\right)\right]}_{ii}=f\left({\lambda }_i\right)$

Examples

${\mathbf{A}}^2=(\mathbf{Q}\mathrm{\Lambda }{\mathbf{Q}}^{-1})(\mathbf{Q}\mathrm{\Lambda }{\mathbf{Q}}^{-1})=\mathbf{Q}\mathrm{\Lambda }({\mathbf{Q}}^{-1}\mathbf{Q})\mathrm{\Lambda }{\mathbf{Q}}^{-1}=\mathbf{Q}{\mathrm{\Lambda }}^2{\mathbf{Q}}^{-1}$

${\mathbf{A}}^n=\mathbf{Q}{\mathrm{\Lambda }}^n{\mathbf{Q}}^{-1}$

Decomposition for special matrices

Normal matrices

A complex normal matrix ($A^{*}A=A{A}^{*}$) has an orthogonal eigenvector basis, so a normal matrix can be decomposed as

$\mathbf{A}=\mathbf{U}\mathrm{\Lambda }{\mathbf{U}}^{*}$

where U is a unitary matrix. Further, if A is Hermitian ($A=A^{*}$), which implies that it is also complex normal, the diagonal matrix Λ has only real values, and if A is unitary, Λ takes all its values on the complex unit circle.

Real symmetric matrices

As a special case, for every N×N real symmetric matrix, the eigenvectors can be chosen such that they are real, orthogonal to each other and have norm one. Thus a real symmetric matrix A can be decomposed as

$\mathbf{A}=\mathbf{Q}\mathrm{\Lambda }{\mathbf{Q}}^T$

where Q is an orthogonal matrix, and Λ is a diagonal matrix whose entries are the eigenvalues of A.

Useful facts

Useful facts regarding eigenvalues

• The product of the eigenvalues is equal to the determinant of A

$\det \left(\mathbf{A}\right)=\prod _{i=1}^{N_{\lambda }}{\lambda }_i^{n_i}$

Note that each eigenvalue is raised to the power n_i, the algebraic multiplicity.

• The sum of the eigenvalues is equal to the trace of A

$\mathrm{tr}\left(\mathbf{A}\right)=\sum _{i=1}^{N_{\lambda }}{n}_i{\lambda }_i$

Note that each eigenvalue is multiplied by n_i, the algebraic multiplicity.

• If the eigenvalues of A are λ_i, and A is invertible, then the eigenvalues of A⁻¹ are simply λ_i⁻¹.
• If the eigenvalues of A are λ_i, then the eigenvalues of f(A) are simply f(λ_i), for any holomorphic function f.

Useful facts regarding eigenvectors

• If A is Hermitian and full-rank, the basis of eigenvectors may be chosen to be mutually orthogonal. The eigenvalues are real.
• The eigenvectors of A⁻¹ are the same as the eigenvectors of A.

Useful facts regarding eigendecomposition

• A can be eigendecomposed if and only if

$N_{\mathbf{v}}=N$

• If p(λ) has no repeated roots, i.e. N_λ = N, then A can be eigendecomposed.

• The statement "A can be eigendecomposed" does not imply that A has an inverse.

• The statement "A has an inverse" does not imply that A can be eigendecomposed.

Useful facts regarding matrix inverse

• $\mathbf{A}$ can be inverted if and only if

${\lambda }_i\ne 0\mathrm{\forall }i$

• If ${\lambda }_i\ne 0\mathrm{\forall }i$ and $N_{\mathbf{v}}=N$, the inverse is given by

${\mathbf{A}}^{-1}=\mathbf{Q}{\mathrm{\Lambda }}^{-1}{\mathbf{Q}}^{-1}$

Numerical computations

DTZ gives a comprehensive integrated property and services administration resolution for buyers, corporate house for sale In singapore owners, management firms and occupiers of property whatever their needs with the only goal of optimising and enhancing the investment worth of their property. We at the moment make use of a staff of more than 70 skilled staffs who are well-trained and dedicated to collectively achieving our purchasers' objectives.

Actual estate agency specialising in non-public condos and landed properties island vast. 10 Winstedt Highway, District 10, #01-thirteen, Singapore 227977. Property providers for enterprise relocation. Situated at 371 Beach Street, #19-10 KeyPoint, Singapore 199597. Property agents for homes, town houses, landed property, residences and condominium for sales and rentals of properties. Administration letting services for property homeowners. is there a single authority in singapore who regulates real property agents that i can file a complaint with for unethical behaviour? or is CASE is simply route? The 188 pages of Secrets and techniques of Singapore Property Gurus are full of professional knowledge and life altering wisdom. Asian industrial property market outlook Property Listing Supervisor Property Advertising Services

Should sellers go along with an agent who claims to specialize in your space? His experience might turn out to be useful, but he is probably additionally advertising a number of models within the neighbourhood – and so they're all your rivals. Within the worst-case state of affairs, your house may be used as a "showflat" as house owner YS Liang found. "Weekend after weekend, our agent would convey a stream of individuals to speed-go to our apartment, leaving within minutes. She did not even try to promote our condominium. It felt like we were just one of the many tour stops for her clients," he complains.

Step one in direction of conducting enterprise as an actual property company in Singapore is to include an organization, or if you happen to're going the partnership or sole-proprietorship route, register your Limited Legal responsibility Partnership or sole-proprietorship with the ACRA (Accounting and Company Regulatory Authority of Singapore) Whether or not you might be considering to promote, let, hire or buy a new industrial property, we're right here to assist. Search and browse our commercial property section. Possess not less than 3 years of working expertise below a Singapore licensed real-property agency; Sale, letting and property administration and taxation companies. three Shenton Means, #10-08 Shenton Home, Singapore 068805. Real property agents for purchasing, promoting, leasing, and renting property. Caveat Search

Firstly, the events might take into account to rescind the sale and buy agreement altogether. This avenue places the contracting events to a position as if the contract didn't happen. It's as if the contract was terminated from the start and events are put back into place that they were before the contract. Any items or monies handed are returned to the respective original house owners. As the worldwide real property market turns into extra refined and worldwide real property investments will increase, the ERA real estate network is well equipped to offer professional recommendation and guidance to our shoppers in making critical actual estate decisions. Relocationg, leasing and sales of properties for housing, food and beverage, retail and workplace wants.

Pasir Panjang, Singapore - $5,000-6,000 per 30 days By likelihood one among our buddies here in Singapore is an agent and we made contact for her to help us locate an residence, which she did. days from the date of execution if the doc is signed in Singapore; Be a Singapore Citizen or PR (Permanent Resident); The regulations also prohibit property agents from referring their shoppers to moneylenders, to discourage irresponsible shopping for. Brokers are additionally prohibited from holding or dealing with money on behalf of any party in relation to the sale or purchase of any property situated in Singapore, and the lease of HDB property. - Negotiate To Close A Sale together with sale and lease of HDB and private properties) Preparing your house for sale FEATURED COMMERCIAL AGENTS Property Guides

i) registered as a patent agent or its equal in any nation or territory, or by a patent workplace, specified within the Fourth Schedule; The business-specific tips for the true property agency and telecommunication sectors have been crafted to address considerations about scenarios that particularly apply to the two sectors, the PDPC stated. Mr Steven Tan, Managing Director of OrangeTee real property company, nonetheless, felt that it was a matter of "practising until it becomes part of our knowledge". "After a while, the agents ought to know the spirit behind the (Act)," he stated. Rising office sector leads real property market efficiency, while prime retail and enterprise park segments moderate and residential sector continues in decline Please choose an attendee for donation.

Numerical computation of eigenvalues

Suppose that we want to compute the eigenvalues of a given matrix. If the matrix is small, we can compute them symbolically using the characteristic polynomial. However, this is often impossible for larger matrices, in which case we must use a numerical method.

In practice, eigenvalues of large matrices are not computed using the characteristic polynomial. Computing the polynomial becomes expensive in itself, and exact (symbolic) roots of a high-degree polynomial can be difficult to compute and express: the Abel–Ruffini theorem implies that the roots of high-degree (5 or above) polynomials cannot in general be expressed simply using nth roots. Therefore, general algorithms to find eigenvectors and eigenvalues are iterative.

Iterative numerical algorithms for approximating roots of polynomials exist, such as Newton's method, but in general it is impractical to compute the characteristic polynomial and then apply these methods. One reason is that small round-off errors in the coefficients of the characteristic polynomial can lead to large errors in the eigenvalues and eigenvectors: the roots are an extremely ill-conditioned function of the coefficients.

A simple and accurate iterative method is the power method: a random vector $v$ is chosen and a sequence of unit vectors is computed as

$\frac{Av}{\Vert Av\Vert },\frac{A^{2}v}{\Vert A^{2}v\Vert },\frac{A^{3}v}{\Vert A^{3}v\Vert },\dots$

This sequence will almost always converge to an eigenvector corresponding to the eigenvalue of greatest magnitude, provided that v has a nonzero component of this eigenvector in the eigenvector basis (and also provided that there is only one eigenvalue of greatest magnitude). This simple algorithm is useful in some practical applications; for example, Google uses it to calculate the page rank of documents in their search engine. Also, the power method is the starting point for many more sophisticated algorithms. For instance, by keeping not just the last vector in the sequence, but instead looking at the span of all the vectors in the sequence, one can get a better (faster converging) approximation for the eigenvector, and this idea is the basis of Arnoldi iteration. Alternatively, the important QR algorithm is also based on a subtle transformation of a power method.

Numerical computation of eigenvectors

Once the eigenvalues are computed, the eigenvectors could be calculated by solving the equation

$(\mathbf{A}-{\lambda }_i\mathbf{I}){\mathbf{v}}_{i,j}=0$

using Gaussian elimination or any other method for solving matrix equations.

However, in practical large-scale eigenvalue methods, the eigenvectors are usually computed in other ways, as a byproduct of the eigenvalue computation. In power iteration, for example, the eigenvector is actually computed before the eigenvalue (which is typically computed by the Rayleigh quotient of the eigenvector). In the QR algorithm for a Hermitian matrix (or any normal matrix), the orthonormal eigenvectors are obtained as a product of the Q matrices from the steps in the algorithm. (For more general matrices, the QR algorithm yields the Schur decomposition first, from which the eigenvectors can be obtained by a backsubstitution procedure.) For Hermitian matrices, the Divide-and-conquer eigenvalue algorithm is more efficient than the QR algorithm if both eigenvectors and eigenvalues are desired.

Additional topics

Generalized eigenspaces

Recall that the geometric multiplicity of an eigenvalue can be described as the dimension of the associated eigenspace, the nullspace of λI − A. The algebraic multiplicity can also be thought of as a dimension: it is the dimension of the associated generalized eigenspace (1st sense), which is the nullspace of the matrix (λI − A)^k for any sufficiently large k. That is, it is the space of generalized eigenvectors (1st sense), where a generalized eigenvector is any vector which eventually becomes 0 if λI − A is applied to it enough times successively. Any eigenvector is a generalized eigenvector, and so each eigenspace is contained in the associated generalized eigenspace. This provides an easy proof that the geometric multiplicity is always less than or equal to the algebraic multiplicity.

This usage should not be confused with the generalized eigenvalue problem described below.

Conjugate eigenvector

A conjugate eigenvector or coneigenvector is a vector sent after transformation to a scalar multiple of its conjugate, where the scalar is called the conjugate eigenvalue or coneigenvalue of the linear transformation. The coneigenvectors and coneigenvalues represent essentially the same information and meaning as the regular eigenvectors and eigenvalues, but arise when an alternative coordinate system is used. The corresponding equation is

$Av=\lambda v^{*}.$

For example, in coherent electromagnetic scattering theory, the linear transformation A represents the action performed by the scattering object, and the eigenvectors represent polarization states of the electromagnetic wave. In optics, the coordinate system is defined from the wave's viewpoint, known as the Forward Scattering Alignment (FSA), and gives rise to a regular eigenvalue equation, whereas in radar, the coordinate system is defined from the radar's viewpoint, known as the Back Scattering Alignment (BSA), and gives rise to a coneigenvalue equation.

Generalized eigenvalue problem

A generalized eigenvalue problem (2nd sense) is the problem of finding a vector v that obeys

$A\mathbf{v}=\lambda B\mathbf{v}$

where A and B are matrices. If v obeys this equation, with some λ, then we call v the generalized eigenvector of A and B (in the 2nd sense), and λ is called the generalized eigenvalue of A and B (in the 2nd sense) which corresponds to the generalized eigenvector v. The possible values of λ must obey the following equation

$\det (A-\lambda B)=0.$

In the case we can find $n\in \mathbb{\mathbb{N}}$ linearly independent vectors $\{{\mathbf{v}}_{\mathbf{1}},\dots ,{\mathbf{v}}_{\mathbf{n}}\}$ so that for every $i\in \{1,\dots ,n\}$, $A{\mathbf{v}}_i={\lambda }_{i}B{\mathbf{v}}_i$, where ${\lambda }_i\in \mathbb{𝔽}$ then we define the matrices P and D such that

$P=\left(\begin{array}{ccc}|& & |\\ {\mathbf{v}}_{\mathbf{1}}& \cdots & {\mathbf{v}}_{\mathbf{n}}\\ |& & |\\ \end{array}\right)$ ≡ $\left(\begin{array}{ccc}({\mathbf{v}}_{\mathbf{1}}{)}_1& \cdots & ({\mathbf{v}}_{\mathbf{n}}{)}_1\\ ⋮& & ⋮\\ ({\mathbf{v}}_{\mathbf{1}}{)}_n& \cdots & ({\mathbf{v}}_{\mathbf{n}}{)}_n\\ \end{array}\right)$

$(D{)}_{ij}=\{\begin{array}{ll}{\lambda }_i,& \text{if }i=j\\ 0,& \text{else}\end{array}$

Then the following equality holds

$\mathbf{A}=\mathbf{B}\mathbf{P}\mathbf{D}{\mathbf{P}}^{\mathbf{-}\mathbf{1}}$

And the proof is

$\mathbf{A}\mathbf{P}=\mathbf{A}\left(\begin{array}{ccc}|& & |\\ {\mathbf{v}}_{\mathbf{1}}& \cdots & {\mathbf{v}}_{\mathbf{n}}\\ |& & |\\ \end{array}\right)=\left(\begin{array}{ccc}|& & |\\ A{\mathbf{v}}_{\mathbf{1}}& \cdots & A{\mathbf{v}}_{\mathbf{n}}\\ |& & |\\ \end{array}\right)=$

$\left(\begin{array}{ccc}|& & |\\ {\lambda }_{1}B{\mathbf{v}}_{\mathbf{1}}& \cdots & {\lambda }_{n}B{\mathbf{v}}_{\mathbf{n}}\\ |& & |\\ \end{array}\right)=\left(\begin{array}{ccc}|& & |\\ B{\mathbf{v}}_{\mathbf{1}}& \cdots & B{\mathbf{v}}_{\mathbf{n}}\\ |& & |\\ \end{array}\right)\mathbf{D}=\mathbf{B}\mathbf{P}\mathbf{D}$

And since P is invertible, we multiply the equation from the right by its inverse and QED.

The set of matrices of the form $A-\lambda B$, where $\lambda$ is a complex number, is called a pencil; the term matrix pencil can also refer to the pair (A,B) of matrices. If B is invertible, then the original problem can be written in the form

$B^{-1}Av=\lambda v$

which is a standard eigenvalue problem. However, in most situations it is preferable not to perform the inversion, but rather to solve the generalized eigenvalue problem as stated originally. This is especially important if A and B are Hermitian matrices, since in this case $B^{-1}A$ is not generally Hermitian and important properties of the solution are no longer apparent.

If A and B are Hermitian and B is a positive-definite matrix, the eigenvalues λ are real and eigenvectors v₁ and v₂ with distinct eigenvalues are B-orthogonal ($v_1^{*}B{v}_2=0$). Also, in this case it is guaranteed that there exists a basis of generalized eigenvectors (it is not a defective problem). This case is sometimes called a Hermitian definite pencil or definite pencil.

See also

• Matrix decomposition
• List of matrices
• Eigenvalue, eigenvector and eigenspace
• Spectral theorem
• Householder transformation
• Frobenius covariant
• Sylvester's formula

References

43 year old Petroleum Engineer Harry from Deep River, usually spends time with hobbies and interests like renting movies, property developers in singapore new condominium and vehicle racing. Constantly enjoys going to destinations like Camino Real de Tierra Adentro.

Bibliography

• 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.
  
  My blog: http://www.primaboinca.com/view_profile.php?userid=5889534
• 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.
  
  My blog: http://www.primaboinca.com/view_profile.php?userid=5889534
• 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.
  
  My blog: http://www.primaboinca.com/view_profile.php?userid=5889534
• 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.
  
  My blog: http://www.primaboinca.com/view_profile.php?userid=5889534
• 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.
  
  My blog: http://www.primaboinca.com/view_profile.php?userid=5889534

External links

• Interactive program & tutorial of Spectral Decomposition.