linalg(deprecated)/eigenvectors - Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : linalg(deprecated)/eigenvectors

linalg(deprecated)

  

eigenvectors

  

find the eigenvectors of a matrix

 

Calling Sequence

Parameters

Description

Examples

Calling Sequence

eigenvectors(A)

eigenvectors(A, 'radical')

eigenvectors(A, 'implicit')

Parameters

A

-

square matrix

Description

• 

Important: The linalg package has been deprecated. Use the superseding command LinearAlgebra[Eigenvectors], instead.

  

- For information on migrating linalg code to the new packages, see examples/LinearAlgebraMigration.

• 

The procedure eigenvectors computes the eigenvalues and eigenvectors of A.  That is, for each eigenvalue lambda of A it solves the linear system IlambdaAX=0 for X.

• 

The result returned is a sequence of lists of the form &lsqb;ei,mi,{v&lsqb;1,i&rsqb;,tripledotplaceholderv&lsqb;ni,i&rsqb;}&rsqb;, where the ei are the eigenvalues, mi their algebraic multiplicities, {v&lsqb;1,i&rsqb;,...,v&lsqb;ni,i&rsqb;} is a set of basis vectors for the eigenspace corresponding to ei, and 1<=ni<=mi is the dimension of the eigenspace. The first example below shows how to extract the eigenvalues and their eigenvectors from the result.

• 

Numeric Case: If the matrix A contains any floating-point (decimal) numbers, the eigenvectors are computed numerically.  A standard numerical algorithm is used.  All floating point arithmetic is done at Digits digits of precision.  Note that the matrix entries on input must all be all of type numeric or complex(numeric).

• 

Symbolic Case: Otherwise the eigenvalues and eigenvectors are computed symbolically (exactly). First the characteristic polynomial is computed and solved for its roots (the eigenvalues) lambda[i] symbolically.  Then for each eigenvalue, a basis for its eigenspace is computed by computing the null space of the characteristic matrix AlambdaiI.

• 

Although eigenvalues and eigenvectors can in principle be computed for a matrix over any field F, one must be able to compute effectively with the roots of a univariate polynomial over F. For this reason, this routine is limited at present to the field of algebraic numbers or algebraic functions.  Therefore, the input A must be a matrix of rationals, or polynomials, or a matrix of algebraic numbers or algebraic functions represented using Maple's RootOf or radical notation.

• 

If the optional second argument is 'radical' is specified, Maple will try to express the eigenvalues and eigenvectors in terms of exact radicals. Note this may fail if any factor of the characteristic polynomial has degree 5 or higher. If the optional second argument is 'implicit' then Maple will express the eigenvalues in terms of an algebraic extension of any nonlinear factors of the characteristic polynomial using Maple's RootOf notation. This is always possible. If no optional second argument is specified, then Maple will use radicals as the default.  The examples below will make these differences clear.

• 

The command with(linalg,eigenvectors) allows the use of the abbreviated form of this command.

Examples

Important: The linalg package has been deprecated. Use the superseding command LinearAlgebra[Eigenvectors], instead.

withlinalg&colon;

Amatrix3&comma;3&comma;1&comma;3&comma;3&comma;3&comma;5&comma;3&comma;6&comma;6&comma;4

A133353664

(1)

eeigenvaluesA

e4&comma;2&comma;2

(2)

veigenvectorsA

v4&comma;1&comma;112&comma;2&comma;2&comma;101&comma;110

(3)

v11

4

(4)

v12

1

(5)

v13

112

(6)

v21

2

(7)

v22

2

(8)

Here is a numerical example (A is a symmetric Toeplitz matrix)

Atoeplitz1.0&comma;2.0&comma;3.0

A1.02.03.02.01.02.03.02.01.0

(9)

eigenvaluesA

2.00000000000000&comma;0.701562118716425&comma;5.70156211871643

(10)

eigenvectorsA

2.000000002&comma;1&comma;0.70710678111.4710-100.7071067809&comma;5.701562115&comma;1&comma;0.60591280040.51549913480.6059128001&comma;0.7015621183&comma;1&comma;0.36451293300.85689009960.3645129334

(11)

Here is the same example done with exact arithmetic

Atoeplitz1&comma;2&comma;3

A123212321

(12)

eigenvaluesA

2&comma;521241&comma;52&plus;1241

(13)

eigenvectorsA

2&comma;1&comma;101&comma;52&plus;1241&comma;1&comma;134&plus;14411&comma;521241&comma;1&comma;13414411

(14)

By default the eigenvalues and eigenvectors are computed using radicals. To compute implicitly in terms of RootOf's, use the 'implicit' option

eigenvectorsA&comma;&apos;implicit&apos;

RootOf_Z25_Z4&comma;1&comma;112RootOf_Z25_Z421&comma;2&comma;1&comma;101

(15)

Here is an example of a symbolic matrix

Ttoeplitza&comma;b&comma;c

Tabcbabcba

(16)

eigenvaluesT

ac&comma;a&plus;12c&plus;128b2&plus;c2&comma;a&plus;12c128b2&plus;c2

(17)

eigenvectorsT&comma;&apos;implicit&apos;

ac&comma;1&comma;101&comma;RootOf_Z2&plus;2ac_Z&plus;a2&plus;ca2b2&comma;1&comma;1RootOf_Z2&plus;2ac_Z&plus;a2&plus;ca2b2&plus;a&plus;cb1

(18)

The sqrt(2) in the following example makes the computation more difficult.

Ttoeplitz2&comma;a&comma;3

T2a3a2a3a2

(19)

eigenvectorsT

32&plus;2&plus;128a2&plus;9&comma;1&comma;r&comma;32&plus;2128a2&plus;9&comma;1&comma;r&comma;23&comma;1&comma;r

(20)

This final example shows a case where the algebraic multiplicity of the eigenvalue u is 4 but the dimension of the eigenspace is only 3.

Bmatrix4&comma;4&comma;0&comma;1&comma;0&comma;0&comma;u2&comma;2u&comma;0&comma;0&comma;us&comma;s&comma;u&comma;0&comma;ut&comma;t&comma;0&comma;u

B0100u22u00ussu0utt0u

(21)

eigenvectorsB

u&comma;4&comma;0001&comma;0010&comma;1u100

(22)

See Also

linalg(deprecated)[charmat]

linalg(deprecated)[eigenvalues]

linalg(deprecated)[nullspace]

linalg(deprecated)[vector]

LinearAlgebra

RootOf