Eigensystem
 Usage
Notes
 Further Examples
These are the eigenvalues of an exact  matrix m1.
In[1]:=
|
Out[1]=
|
These are the eigenvectors of m1.
In[2]:=
|
Out[2]=
|
This lists the eigenvalues and the eigenvectors for m1. The eigenvalues are listed first.
In[3]:=
|
Out[3]=
|
This verifies the result.
In[4]:=
|
Out[4]=
|
Here are two inexact  matrices.
In[5]:=
|
This gives the generalized eigenvalues and eigenvectors of m2 with respect to m3.
In[6]:=
|
Out[6]=
|
This checks the result.
In[7]:=
|
Out[7]=
|
Here are two more matrices.
In[8]:=
|
One of the generalized eigenvalues is  .
In[9]:=
|
Out[9]=
|
This checks the values of the two finite eigenvalues.
In[10]:=
|
Out[10]=
|
The reciprocals of the eigenvalues of  are equal to the eigenvalues of  .
In[11]:=
|
Out[11]=
|
In[12]:=
|
Using Eigensystem to diagonalize a matrix We can use Eigensystem to diagonalize a matrix that has distinct eigenvalues.
In[13]:=
|
In[14]:=
|
The eigenvalues give the diagonalized form.
In[15]:=
|
Out[15]//MatrixForm=
|
The columns of the change of basis matrix ch are the eigenvectors.
In[16]:=
|
We check the diagonalization.
In[17]:=
|
Out[17]//MatrixForm=
|
In[18]:=
|
Using Eigenvalues to find the roots of a polynomial Solve uses the Jenkins-Traub algorithm to find approximate numeric roots to a polynomial in one variable. Another method uses the eigenvalues of the companion matrix for that polynomial. This gives the companion matrix of a polynomial.
In[19]:=
|
Here is a polynomial and its companion matrix.
In[20]:=
|
In[21]:=
|
Out[21]//MatrixForm=
|
This sorts the eigenvalues of the polynomial's companion matrix.
In[22]:=
|
Out[22]=
|
Here is a comparison with the Mathematica default technique of finding the roots of a polynomial.
In[23]:=
|
Out[23]=
|
In[24]:=
|
When a matrix is defective The  matrix m1 is defective since the number of linearly independent eigenvectors is less than  .
In[25]:=
|
This lists the eigenvalues and the eigenvectors for m1. The eigenvalues are listed first.
In[26]:=
|
Out[26]=
|
Although the vectors in vectors are the eigenvectors of the matrix m1, they provide only one linearly independent vector,  . The other computed vectors are linearly dependent on  . It follows that m1 is not diagonalizable.
In[27]:=
|
Out[27]=
|
Note that the computations are inexact, and the eigenvectors are also inexact. Here it is evident that the terms of order  are roundoff errors. However, it is not always correct to use Chop since Chop may remove terms which are not roundoff errors.
In[28]:=
|
Out[28]=
|
In[29]:=
|
| 
