[Numpy-discussion] Eigenvectors in Matlab vs. Numpy
Saket
saketn@gmail....
Sun Jun 29 21:00:22 CDT 2008
Hmm... so the relationship Ax = Lx should hold for every eigenvalue
and corresponding eigenvector of A, right? But, consider the first
eigenvalue,eigenvector pair:
for i,eval in enumerate(d):
print abs(numpy.dot(A,v[i]) - numpy.dot(eval,v[i])).max()
return
Outputs: 1.928
I thought maybe the ith eigenvector corresponds to a different (not
the ith) eigenvalue, but there doesn't seem to be any eigenvalue which
corresponds to the ith eigenvector such that the relationship holds...
Thanks again.
Saket
atlas_threads_info:
NOT AVAILABLE
blas_opt_info:
libraries = ['f77blas', 'cblas', 'atlas']
library_dirs = ['C:\\local\\lib\\yop\\sse2']
define_macros = [('ATLAS_INFO', '"\\"?.?.?\\""')]
language = c
atlas_blas_threads_info:
NOT AVAILABLE
lapack_opt_info:
libraries = ['lapack', 'f77blas', 'cblas', 'atlas']
library_dirs = ['C:\\local\\lib\\yop\\sse2']
define_macros = [('ATLAS_INFO', '"\\"?.?.?\\""')]
language = f77
atlas_info:
libraries = ['lapack', 'f77blas', 'cblas', 'atlas']
library_dirs = ['C:\\local\\lib\\yop\\sse2']
language = f77
lapack_mkl_info:
NOT AVAILABLE
blas_mkl_info:
NOT AVAILABLE
atlas_blas_info:
libraries = ['f77blas', 'cblas', 'atlas']
library_dirs = ['C:\\local\\lib\\yop\\sse2']
language = c
mkl_info:
NOT AVAILABLE
On Sun, Jun 29, 2008 at 9:15 PM, Charles R Harris
<charlesr.harris@gmail.com> wrote:
>
>
> On Sun, Jun 29, 2008 at 6:47 PM, Saket <saketn@gmail.com> wrote:
>>
>> Hi,
>>
>> I'm having this weird problem when computing eigenvalues/vectors with
>> Numpy. I have the following symmetric matrix, B:
>>
>> -0.3462 0.6538 0.5385 -0.4615 0.6538 -0.3462 -0.3462
>> -0.3462
>> 0.6538 -0.3462 0.5385 -0.4615 0.6538 -0.3462 -0.3462
>> -0.3462
>> 0.5385 0.5385 -0.6154 0.3846 0.5385 -0.4615 -0.4615
>> -0.4615
>> -0.4615 -0.4615 0.3846 -0.6154 -0.4615 0.5385 0.5385
>> 0.5385
>> 0.6538 0.6538 0.5385 -0.4615 -0.3462 -0.3462 -0.3462
>> -0.3462
>> -0.3462 -0.3462 -0.4615 0.5385 -0.3462 -0.3462 0.6538
>> 0.6538
>> -0.3462 -0.3462 -0.4615 0.5385 -0.3462 0.6538 -0.3462
>> 0.6538
>> -0.3462 -0.3462 -0.4615 0.5385 -0.3462 0.6538 0.6538
>> -0.3462
>>
>> I compute the eigenvalues and eigenvectors of B using
>> numpy.linalg.eig(B). I get the following eigenvalues:
>>
>> [ 2.79128785e+00 -1.79128785e+00 1.64060486e-16 -3.07692308e-01
>> -1.00000000e+00 -1.00000000e+00 -1.00000000e+00 -1.00000000e+00]
>>
>> I do the same thing in Matlab and get the SAME eigenvalues. However,
>> my eigenVECTORS in Matlab versus numpy are different. It makes no
>> sense to me. In general, the following relationship should hold: Bx =
>> Lx, where B is my matrix, x is an eigenvector, and L is the
>> corresponding eigenvalue. For the eigenvectors that Matlab returns, I
>> have confirmed that the relationship does hold. But for the Numpy
>> eigenvectors, it doesn't!
>>
>> Any idea why this might be happening? I did some computations myself
>> and it looks like the Matlab output is correct. Just seems like the
>> eigenvectors that Numpy is returning are wrong...
>>
>> Thanks for any suggestions.
>
> Also note that the -1 eigenvalue has multiplicity 4. This means that any set
> of orthogonal vectors spanning the same eigenspace will do for eigenvectors,
> i.e., they aren't unique and roundoff error is likely to have a large effect
> on what you end up with.
>
> Chuck
>
>
>
> _______________________________________________
> Numpy-discussion mailing list
> Numpy-discussion@scipy.org
> http://projects.scipy.org/mailman/listinfo/numpy-discussion
>
>
More information about the Numpy-discussion
mailing list