Similarity Transformation Doesn't seem to work

  • Context: Graduate 
  • Thread starter Thread starter I_am_learning
  • Start date Start date
  • Tags Tags
    Transformation Work
Click For Summary
SUMMARY

The discussion centers on the behavior of similarity transformations for square matrices, specifically addressing the case where matrix A has all eigenvalues equal to zero. The transformation A' = B-1AB, where B is a non-singular matrix, is expected to yield the same eigenvalues for A and A'. However, the user encountered discrepancies when A = [0 1 0; 0 0 0; 1 0 0], leading to eigenvalues that appear non-zero due to numerical errors. The conclusion emphasizes the importance of considering numerical stability when dealing with small eigenvalues.

PREREQUISITES
  • Understanding of eigenvalues and eigenvectors
  • Familiarity with similarity transformations in linear algebra
  • Knowledge of numerical methods and their implications on matrix computations
  • Experience with MATLAB or similar computational tools for matrix operations
NEXT STEPS
  • Investigate numerical stability in eigenvalue computations
  • Learn about the implications of matrix conditioning on eigenvalue accuracy
  • Explore MATLAB functions for eigenvalue analysis, such as eig and inv
  • Study the properties of non-singular matrices and their impact on similarity transformations
USEFUL FOR

Mathematicians, linear algebra students, and computational scientists interested in matrix theory and numerical analysis.

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
Inline13.gif
=0
Inline13.gif
= 0
Inline15.gif
= 0
Inline18.gif
= 0
Inline21.gif
= 0
Inline24.gif
=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.
 

Similar threads

  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 6 ·
Replies
6
Views
5K
  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 28 ·
Replies
28
Views
5K
  • · Replies 6 ·
Replies
6
Views
3K
  • · Replies 5 ·
Replies
5
Views
6K
  • · Replies 8 ·
Replies
8
Views
3K
  • · Replies 7 ·
Replies
7
Views
3K
  • · Replies 1 ·
Replies
1
Views
670
  • · Replies 6 ·
Replies
6
Views
4K