LinearAlgebra
Eigenvectors
compute the eigenvectors of a Matrix
Calling Sequence
Parameters
Description
Examples
Eigenvectors(A, C, imp, out, options, outopts)
A
-
Matrix; Matrix whose Eigenvectors are required
C
(optional) Matrix; Matrix for the generalized eigenvector problem
imp
(optional) equation of the form implicit=true or false; 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
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 float8 or sfloat datatype. Otherwise the returned eigenvalue object has complex8 or complex⁡sfloat 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 float8 or sfloat datatype. Otherwise the returned eigenvalue object has complex8 or complex⁡sfloat 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](..).
with⁡LinearAlgebra:
Find the eigenvalues and eigenvectors of the Matrix A.
A≔−1,−3,−6|3,5,6|−3,−3,−4
A≔−13−3−35−3−66−4
v,e≔Eigenvectors⁡A
v,e≔22−4,1−1121012011
Verify A·x=λ⁢x for the second eigenvalue and second eigenvector.
A·e1..−1,2=v2·e1..−1,2
−202=−202
Find the eigenvectors of the Matrix B using the option output=list.
B≔1,2,3|2,4,6|5,10,15
B≔12524103615
Eigenvectors⁡B,output=list
20,1,13231,0,2,−501,−210
Find the generalized eigenvectors of M and N.
M≔Matrix⁡6.,8.,5.,8.,8.,9.,5.,9.,6.,datatype=float:
N≔Matrix⁡6.,3.,3.,3.,8.,8.,3.,8.,9.,datatype=float:
v,e≔Eigenvectors⁡M,N,output=vectors,values
v,e≔−0.139566030823037+0.⁢I−0.940687706820450+0.⁢I−1.+0.⁢I1.+0.⁢I−1.+0.⁢I0.0726922335169212+0.⁢I−0.911459942195583+0.⁢I0.399204400339097+0.⁢I0.953350738117875+0.⁢I,−4.55651173111531+0.⁢I1.56423606281299+0.⁢I0.223044899071558+0.⁢I
A simple floating-point eigenvector problem.
N≔1.0,4.0⁢I|−4.0⁢I,1.0
N≔1.0−4.0⁢I4.0⁢I1.0
evalsN,evecsN≔Eigenvectors⁡N
evalsN,evecsN≔5.00000000000000+0.⁢I−3.+0.⁢I,−0.−0.707106781186547⁢I0.707106781186547+0.⁢I0.707106781186548+0.⁢I0.−0.707106781186547⁢I
By default, the resulting eigenvalues are given as complex numbers.
VectorOptions⁡evalsN,datatype,MatrixOptions⁡evecsN,datatype
complex8,complex8
However, if N has the hermitian indexing function, then the resulting eigenvalues are given as real floating-point numbers.
evalsN,evecsN≔Eigenvectors⁡Matrix⁡N,shape=hermitian
evalsN,evecsN≔−3.5.,0.−0.707106781186547⁢I0.−0.707106781186547⁢I−0.707106781186547−0.⁢I0.707106781186547+0.⁢I
float8,complex8
See Also
LinearAlgebra[CharacteristicMatrix]
LinearAlgebra[Eigenvalues]
LinearAlgebra[IdentityMatrix]
LinearAlgebra[JordanForm]
LinearAlgebra[SchurForm]
Matrix
Vector
Download Help Document