[Date Index]
[Thread Index]
[Author Index]
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**
| |