|
Consider a square matrix c. If:
|
 |
[1] |
for some scalar
and some vector v
0, then we call
an eigenvalue and v an
eigenvector of c. This means
that an eigenvector of a matrix is any vector for which multiplication by
the matrix is not different from multiplication by a scalar (the
eigenvalue).
Suppose
is an eigenvalue and v is a corresponding eigenvector of
c. For any scalar a, the vector av is also
an eigenvector corresponding to eigenvalue
.
This follows because
|
c(av)
= a(cv) = a( v)
=
(av) |
[2] |
Accordingly, eigenvectors are uniquely determined only up
to scalar multiplication. If a set of eigenvectors are linearly
independent, we say they are distinct.
To determine eigenvalues and eigenvectors of a matrix, we
focus first on the eigenvalues. Rearranging [1], we
obtain
|
(c –
I)v
= 0 |
[3] |
where I is the identity matrix. This
equation will hold for some nonzero vector v if and only if
the matrix (c –
I)
is singular. Accordingly, we seek values
for which the matrix (c –
I)
has a determinant of 0.
Consider matrix
|
 |
[4] |
for which
|
 |
[5] |
This has determinant
|
 |
[6] |
which is a third-order polynomial. It has roots
= –1, 2, and 3. We find corresponding eigenvectors v by
substituting the eigenvalues into [3] and solving. For
example, with
= –1, [3] becomes:
|
 |
[7] |
By inspection, a solution is v = (1, 3, –3).
Obviously, any multiple of this is also a solution. We repeat the same
analysis for the other eigenvalues. Results are indicated in Exhibit 1.
|
|
 |
|
eigenvalue |
eigenvector |
|
–1 |
(1, 3, –3) |
|
2 |
(1, 0, 0) |
|
3 |
(3, 1, 3) |
|
|
Eigenvalues and eigenvectors of matrix [4]. |
The approach we employed in our example is useful for
deriving an important result. Consider an arbitrary n n
matrix c. To find its eigenvalues, we construct the
determinant of and set it equal to 0. This results in an
-order
polynomial equation. By the
fundamental theorem of algebra, it has n solutions. We conclude
that every matrix has n eigenvalues. Of course, some may be
complex. Others may be
repeated.
In practical applications, eigenvalues are not calculated
in this manner. Although setting the determinant of (c –
I)
equal to 0 and solving is theoretically useful, there are more efficient
algorithms, which are implemented in various software packages.
Eigenvalues have a number of convenient properties. A
matrix and its transpose both have the same eigenvalues. If
is an eigenvalue of a nonsingular matrix, then 1/
is an eigenvalue of its inverse. The product of the eigenvalues of a
matrix equals its determinant.
Consider an intuitive example. A sphere of unit radius is positioned at
the center of a three-dimensional coordinate system. It is rotating about
the
.
The matrix:
|
 |
[8] |
describes a one-eighth (45 )
rotation of the sphere. For example, multiplying c by the
vector (1,0,0) yields the vector (.7071, .7071, 0), which is rotated 45 .
This is depicted in Exhibit 2.
Intuitively, what might be an eigenvector of the matrix
c? Is there a point on the unit sphere that a 45
rotation transforms into a multiple of itself? Of course! Consider the
point at the north pole. It is the point (0, 0, 1), and it is transformed
into itself. We conclude that an eigenvector of c is the
vector (0, 0, 1). The corresponding eigenvalue is 1. Because it is a 3 3
matrix, c has two other eigenvalues, but they are both
complex numbers.
|
|