Eigenvalues and eigenspaces of orthogonal (or rotation) matrices

Anonymous

Given an orthogonal transformation of finite order, e.g.

Matrix([[0,0,0,-1],[1,0,0,-1],[0,1,0,-1],[0,0,1,-1]])

Its eigenvalues are going to be of the form

exp(I*pi/5),exp(2*I*pi/5),...,exp(2*I*pi*m),...

corresponding to a splitting of the matrix into rotation (and reflection) matrices. I'd like to extract these fractions m (mod ZZ) and study the corresponding (real rotation and reflection) eigenspaces.

My impression is that Sage isn't suitable for doing this directly, but that I should use e.g. the Maxima or Mathematica interface? Any suggestions for the most suitable method?

Comments

Well... here is a start:

sage: M=matrix(SR,[[0,0,0,-1],[1,0,0,-1],[0,1,0,-1],[0,0,1,-1]])
sage: list(map(lambda u:u.maxima_methods().polarform(), M.eigenvalues()))
[1/4*sqrt((sqrt(5) + 1)^2 - 2*sqrt(5) + 10)*e^(-I*pi - I*arctan(-2*sqrt(-1/2*sqrt(5) + 5/2)/(sqrt(5) + 1))),
 1/4*sqrt((sqrt(5) + 1)^2 - 2*sqrt(5) + 10)*e^(I*pi + I*arctan(-2*sqrt(-1/2*sqrt(5) + 5/2)/(sqrt(5) + 1))),
 1/4*sqrt((sqrt(5) - 1)^2 + 2*sqrt(5) + 10)*e^(-I*arctan(2*sqrt(1/2*sqrt(5) + 5/2)/(sqrt(5) - 1))),
 1/4*sqrt((sqrt(5) - 1)^2 + 2*sqrt(5) + 10)*e^(I*arctan(2*sqrt(1/2*sqrt(5) + 5/2)/(sqrt(5) - 1)))]

Emmanuel Charpentier ( 5 years ago )

add a comment

Comments

Much, much better that my attempt ...

Emmanuel Charpentier ( 5 years ago )

Thanks, this is excellent! (I need the actual fractions (1/5,2/5,3/5,4/5 in the example) which I can now deduce inductively from your code using the multiplicities, and then I suppose I can obtain the eigenspaces by first comparing those eigenvalues with E(5)^i in UniversalCyclotomicField, and then taking the real and complex parts of the matching eigenvectors gives the rotation planes. In other words, this is very doable now but if you know some smart tricks for this again then please do let me know!)

Bob67846 ( 5 years ago )

add a comment

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Eigenvalues and eigenspaces of orthogonal (or rotation) matrices

Comments

1 Answer

Comments

Your Answer

Question Tools

Stats

Related questions

Eigenvalues and eigenspaces of orthogonal (or rotation) matrices savecancel

Comments

1 Answer

Comments

Your Answer

Question Tools

Stats

Related questions

Eigenvalues and eigenspaces of orthogonal (or rotation) matrices