Mathematica 9 is now available
Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
1997
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*November
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 1997

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Re: Numerical Left Eigenvectors

  • To: mathgroup at smc.vnet.net
  • Subject: [mg6406] Re: [mg6392] Re: [mg6369] Numerical Left Eigenvectors
  • From: "w.meeussen" <w.meeussen.vdmcc at vandemoortele.be>
  • Date: Sun, 16 Mar 1997 19:25:22 -0500 (EST)
  • Sender: owner-wri-mathgroup at wolfram.com

At 01:21 15-03-97 -0500, seanross at worldnet.att.net wrote:
>Gregory Dwyer wrote:
>> 
>> Greetings -
>> 
>> When I use the "eigenvector" function to numerically calculate
>> eigenvectors, Mathematica seems to assume that I always want right
>> eigenvectors.  Is there some way to numerically calculate left eigenvectors?
>> 
>> Thanks.
>> 
>> Greg Dwyer.
>> Entomology, UMASS Amherst
>> dwyer at ent.umass.edu
>
>
>What is a left eigenvector?  The only eigenvectors I know of are ones 
>that satisfy the following equation:
>
>Matrix . vector = constant vector.  Is this a right eigenvector?  Does 
>left eigenvector mean
>vector. Matrix = constant vector?
>
>If so, have you looked for a way to express Matrix . column vector in 
>terms of row vector . Matrix?  For a 2x2 case, it looks like 
>Matrix.columnvector = rowvector.TransposeMatrix
>
>I am not sure how commutable square matrices and vectors are supposed to 
>be.  Perhaps someone else out there with experience in something like 
>group theory or who has a copy of MathTensor etc. could enlighten you 
>further.
>
>
>

hi,

straight from The (electronic) Book under 'Eigensystem':

In[102]:=
it
Out[102]=
{{2,3,4,5},{1,2,7,8},{6,9,2,1}}
In[104]:=
mit=Join[it,{{3,3,9,9}}]
Out[104]=
{{2,3,4,5},{1,2,7,8},{6,9,2,1},{3,3,9,9}}
In[135]:=
{vals,vecs}=Eigensystem[N at mit];
vals
vecs//MatrixForm
Out[135]=
{19.1216,-2.27556+1.89411 I,-2.27556-1.89411 I,0.429553}
Out[136]//MatrixForm=
 *** stuff clipped ***
In[137]:=
mit.Transpose[vecs] == Transpose[vecs].DiagonalMatrix[vals]
Out[137]=
True
In[138]:=
(vecs.Transpose[mit] //Chop)==(DiagonalMatrix[vals].vecs//Chop)
Out[138]=
True

beware : without Chop no equality! different order of operations gives tiny
rounding errors. Except for pure Integers of course.

wouter.

Dr. Wouter L. J. MEEUSSEN
eu000949 at pophost.eunet.be
w.meeussen.vdmcc at vandemoortele.be



  • Prev by Date: Re: Numerical Left Eigenvectors
  • Next by Date: Re: Where is Terry Robb ?
  • Previous by thread: Re: Numerical Left Eigenvectors
  • Next by thread: Re: Numerical Left Eigenvectors