LinearAlgebra - Maple Programming Help

Online Help

All Products    Maple    MapleSim


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

LinearAlgebra

  

EigenConditionNumbers

  

compute condition numbers of eigenvalues or eigenvectors of a numeric eigenvalue conditioning problem

 

Calling Sequence

Parameters

Description

Examples

Calling Sequence

EigenConditionNumbers(A, C, bal, out, options, outopts)

Parameters

A

-

Matrix; numeric Matrix whose eigenvalues' or eigenvectors' condition numbers are required

C

-

(optional) Matrix; numeric Matrix for the generalized numeric eigenvalue conditioning problem

bal

-

(optional) identical('balance')=symbol; specifies whether to allow scaling or permutation, where the symbol is any of 'both', 'permute', 'scale', or 'none'

out

-

(optional) identical('output')={name,list}; specifies whether to return condition numbers of the eigenvalues and eigenvectors or the eigenvalues and eigenvectors

options

-

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

outopts

-

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

Description

• 

The EigenConditionNumbers(A) command solves the simple eigenvalue conditioning problem by returning an expression sequence.  This expression sequence consists of two components.  The first is a Vector V1 of the floating-point condition numbers of the eigenvalues of A, and the second is a Vector V2 of the floating-point condition numbers of the eigenvectors of A.  Optionally, the floating-point eigenvalues or the floating-point eigenvectors are returned.

  

In the simple eigenvalue problem, if A is real-symmetric or complex hermitian, the returned eigenvalues are real.

• 

The EigenConditionNumbers(..) function solves the simple eigenvalue conditioning problem A·x=λ·x and the generalized eigenvalue conditioning 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 M.E1..1,i=Vi.E1..1,i, and for the generalized problem M.E1..1,i=Vi.C.E1..1,i or βi·M·E1..−1,i=αi·C·E1..−1,i.

• 

The EigenConditionNumbers(A, C) command solves the generalized numeric eigenvalue conditioning problem.

  

In the generalized eigenvalue conditioning problem, if A and C are real-symmetric or complex hermitian and C also has the positive_definite attribute, the eigenvalues are real.

  

In the generalized eigenvalue conditioning problem, Vi=βiαi for each finite eigenvalue Vi.

• 

The out option specifies what the command should return.  It can be one of vectors, values, conditionvalues, conditionvectors or a list of these.

• 

The options option provides additional information (readonly, shape, storage, order, datatype, and attributes) to the Matrix or Vector constructor that builds the result. 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).

• 

For extensive details of the relationships between condition numbers and error bounds, see

  

E. Anderson et al. LAPACK users' guide. 3d ed. SIAM, 1999.

Examples

withLinearAlgebra:

AMatrix3,3,1,3,6,3,5,6,3,3,4,'datatype'='float'

A:=1.3.6.3.5.6.3.3.4.

(1)

EigenConditionNumbersA

0.4714045207910320.4362983775049870.745873405008495,6.1.5543122344752210-151.3322676295501910-15

(2)

v,eEigenConditionNumbersA,'output'='values','vectors'

v,e:=4.00000000000000+0.I2.00000000000000+0.I2.+0.I,0.577350269189626+0.I0.0866100166393946+0.I0.808190323654766+0.I0.577350269189626+0.I0.907714476800681+0.I0.254240841587667+0.I0.577350269189626+0.I0.410552230080643+0.I0.531215582621217+0.I

(3)

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

0.173220033278790+0.I1.81542895360136+0.I0.821104460161287+0.I=0.173220033278789+0.I1.81542895360136+0.I0.821104460161286+0.I

(4)

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

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

EigenConditionNumbersM,N,output='vectors','values','conditionvalues','conditionvectors'

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,2.3604228759176110.66875895618884.91085988501676,3.455943437063792.298794840362603.02673623742847

(5)

EigenConditionNumbersM,N,output='values'

4.55651173111532+0.I1.56423606281299+0.I0.223044899071558+0.I

(6)

EigenConditionNumbersM,N,output='α','β'

2.99144878241616+0.I14.9822766669746+0.I1.38335073553634+0.I,0.6565216900437839.578015123901266.20211778567742

(7)

See Also

LinearAlgebra

LinearAlgebra[CharacteristicMatrix]

LinearAlgebra[Eigenvalues]

LinearAlgebra[Eigenvectors]

LinearAlgebra[IdentityMatrix]

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