Similarity Transformation Doesn't seem to work

  • Level: Graduate 
  • Thread starter Thread starter I_am_learning
  • Start date Start date
  • Tags Tags
    Transformation Work
Join the discussion
Registration is free. Start your own thread to ask a follow-up.
3 replies · 1K views · last reply
I_am_learning
Messages
681
Reaction score
16
if A is a square matrix, and A' = B-1AB is its similarity transform (with a non-singular similarity transformation matrix B), then the eigenvalues of A and A' are supposed to be the same. I can verify this for all most all cases of A. But, it doesn't seem to work, when the eigen values of A are all zero. For example, it doesn't work when A = [0 1 0;0 0 0; 1 0 0].
Yet I can't see why, it shouldn't work.
from http://mathworld.wolfram.com/SimilarityTransformation.html
the eigen values of A' is given as a solution of
Expression B^-1 A B - λI, where B^-1 is the inverse of B, followed by A, B and λI, shown as a matrix similarity form minus lambda times the identity.
=0
Inline formula: B^-1 A B - λI, where B^-1 is the inverse of B, A and B are matrices, and λ is an eigenvalue.
= 0
Equation B^-1 A B - B^-1 λ I B = 0 written in typeset form, with B^-1 on the left and B on the right in each term.
= 0
Inline expression B^-1 (A - λI) B, with B^-1 as a superscript-labelled inverse matrix, A minus lambda times the identity, and B on the right.
= 0
Equation: |B^{-1}| |A - λI| |B|, where B^{-1} is the inverse of B, the determinant bars denote determinants, and λ is the eigenvalue variable.
= 0
Plain text expression: |A - λI|, the determinant of A minus lambda times the identity matrix, written with vertical bars.
=0
which is the solution for eigen values of A.
 
Last edited by a moderator:
Physics news on Phys.org
for what matrix B did you find that the eigenvalues of B-1AB to be nonzero?
 
BruceW said:
for what matrix B did you find that the eigenvalues of B-1AB to be nonzero?
I created a random matrix T,
T = rand(3)
T =
0.4509 0.7447 0.1835
0.5470 0.1890 0.3685
0.2963 0.6868 0.6256
and, the eigen values I got are:
eig(inv(T)*A*T)
ans =
1.0e-05 *
0.4418 + 0.0000i
-0.2209 + 0.3826i
-0.2209 - 0.3826i
Oh! now that I post it, I realize that, it may have to do with numerical errors, the values are very small anyway.
Thanks
 
yeah, since the eigenvalues are of order 10^-5, it's likely they are nonzero just due to numerical error.