LinearAlgebra - Maple Programming Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Mathematics : Linear Algebra : LinearAlgebra Package : Eigenvalues : LinearAlgebra/Eigenvectors

LinearAlgebra

  

Eigenvectors

  

compute the eigenvectors of a Matrix

 

Calling Sequence

Parameters

Description

Examples

Calling Sequence

Eigenvectors(A, C, imp, out, options, outopts)

Parameters

A

-

Matrix; Matrix whose Eigenvectors are required

C

-

(optional) Matrix; Matrix for the generalized eigenvector problem

imp

-

(optional) BooleanOpt(implicit); specify that roots of the characteristic polynomial contain implicit RootOfs or radicals

out

-

(optional) equation of the form output = obj where obj is one of 'values', 'vectors', or 'list', or a list containing one or more of these names; selects format of the result object

options

-

(optional); constructor options for the result object(s)

outopts

-

(optional) equation(s) of the form outputoptions[o]=list where o is one of 'values', 'vectors', or 'list'; constructor options for the specified result object

Description

• 

The Eigenvectors(A) command solves the simple eigenvector problem by returning an expression sequence whose first member is a Vector V of the eigenvalues of A, and whose second member is a Matrix E whose columns are the eigenvectors of A.  The ith column of E is an eigenvector associated with the ith eigenvalue of the returned Vector V.

  

In the simple floating-point eigenvector problem, if A has either the symmetric or the hermitian indexing function then the returned eigenvalue object has float[8] or sfloat datatype. Otherwise the returned eigenvalue object has complex[8] or complexsfloat datatype.

• 

The Eigenvectors(..) function solves the simple eigenvector problem A·x=λ·x and the generalized eigenvalue problem A·x=λ·C·x.  The eigenvalues V[i] and the nonzero eigenvectors E[1 .. -1, i] are such that for the simple eigenvalue problem A.E1..1,i=Vi.E1..1,i, and for the generalized problem A.E1..1,i=Vi.C.E1..1,i.

• 

The Eigenvectors(A, C) command solves the generalized eigenvector problem.

  

In the generalized floating-point eigenvector problem, if A and C have either symmetric or hermitian indexing functions and C also has the positive_definite attribute then the returned eigenvalue object has float[8] or sfloat datatype. Otherwise the returned eigenvalue object has complex[8] or complexsfloat datatype.

• 

With an eigenvalue of multiplicity k>1, there may be fewer than k linearly independent eigenvectors. In this case, the matrix is called defective.  By design, the returned matrix always has full column dimension.  Therefore, in the defective case, some of the columns that are returned are zero.  Thus, they are not eigenvectors.  With the option, output=list, only eigenvectors are returned.  For more information, see LinearAlgebra[JordanForm] and LinearAlgebra[SchurForm].

• 

The format in which the Eigenvectors of A are returned is determined by parameter out.  By default, an expression sequence is returned as described above. In the case that the output is specified as list, a list of lists is returned.  The first element of each sublist is an eigenvalue, the second element is its multiplicity, and the third element is a set of linearly independent eigenvectors corresponding to that eigenvalue.

• 

The options option provides additional information (readonly, shape, storage, order, datatype, and attributes) to the Matrix or Vector constructor that builds the result(s). These options may also be provided in the form outputoptions[o]=[...], where [...] represents a Maple list.  If a constructor option is provided in both the calling sequence directly and in an outputoptions[o] option, the latter takes precedence (regardless of the order).

• 

This function is part of the LinearAlgebra package, and so it can be used in the form Eigenvectors(..) only after executing the command with(LinearAlgebra). However, it can always be accessed through the long form of the command by using LinearAlgebra[Eigenvectors](..).

Examples

withLinearAlgebra:

Find the eigenvalues and eigenvectors of the Matrix A.

A1,3,6|3,5,6|3,3,4

A:=133353664

(1)

v,eEigenvectorsA

v,e:=422,12111201110

(2)

Verify A·x=λx for the second eigenvalue and second eigenvector.

A.e1..1,2=v2.e1..1,2

202=202

(3)

Find the eigenvectors of the Matrix B using the option output=list.

B1,2,3|2,4,6|5,10,15

B:=12524103615

(4)

EigenvectorsB,output='list'

0,2,501,210,20,1,13231

(5)

Find the generalized eigenvectors of M and N.

MMatrix6.,8.,5.,8.,8.,9.,5.,9.,6.,datatype=float:

NMatrix6.,3.,3.,3.,8.,8.,3.,8.,9.,datatype=float:

v,eEigenvectorsM,N,output='vectors','values'

v,e:=0.139566030823037+0.I0.940687706820450+0.I1.+0.I1.+0.I1.+0.I0.0726922335169216+0.I0.911459942195583+0.I0.399204400339096+0.I0.953350738117875+0.I,4.55651173111532+0.I1.56423606281299+0.I0.223044899071558+0.I

(6)

A simple floating-point eigenvector problem.

N1.0,4.0I|4.0I,1.0

N:=1.04.0I4.0I1.0

(7)

evalsN,evecsNEigenvectorsN

evalsN,evecsN:=5.00000000000000+0.I3.+0.I,0.0.707106781186547I0.707106781186547+0.I0.707106781186548+0.I0.0.707106781186547I

(8)

By default, the resulting eigenvalues are given as complex numbers.

VectorOptionsevalsN,datatype,MatrixOptionsevecsN,datatype

complex8,complex8

(9)

However, if N has the hermitian indexing function, then the resulting eigenvalues are given as real floating-point numbers.

evalsN,evecsNEigenvectorsMatrixN,shape=hermitian

evalsN,evecsN:=3.5.,0.0.707106781186547I0.0.707106781186547I0.7071067811865470.I0.707106781186547+0.I

(10)

VectorOptionsevalsN,datatype,MatrixOptionsevecsN,datatype

float8,complex8

(11)

See Also

LinearAlgebra[CharacteristicMatrix]

LinearAlgebra[Eigenvalues]

LinearAlgebra[IdentityMatrix]

LinearAlgebra[JordanForm]

LinearAlgebra[SchurForm]

Matrix

Vector

 


Download Help Document

Was this information helpful?



Please add your Comment (Optional)
E-mail Address (Optional)
What is ? This question helps us to combat spam