Thanks. (and also to Dan A.) I guess I owe a little bit of explanation as to why I even cared. I was looking for the simplest case of a 4-d rotation that would not "trivially" be reducible to a lower dimensional case. (Interestingly, it is exactly this difference that seems to throw Macsyma for a loop on this). Now I see (thanks to some offline comments from Dan) that any rotation in R^n is representable as some number of orthogonal 2-d rotations (leaving out many details). One thing I found interesting to try to understand is the meaning of the elements of a skew-symmetric matrix in relation to the rotation matrix obtained by exponentiating it. In 3-d, they represent a vector giving the axis of rotation, and the amount of rotation (the length of that vector). But in higher dimensions, there is no axis. It appears the (i,j) element represents the amount of rotation in the (i,j) coordinate plane. As there are (n choose 2) (i,j) planes, this is an easy way to see that the dimension of SO(n) is (n choose 2). The total rotation doesn't necessarily take the coordinate planes into themselves, due to interaction between rotations in coordinate planes with a common coordinate. You could consider a skew-symmetric matrix as the sum of a number of skew-symmetric matrices with only one pair of non-zero elements each. These each represent a simple rotation in one coordinate plane. These component rotations are combined when exponentiating to get the corresponding rotation matrix, using (also from Dan): exp(a+b+...) = lim[n->oo] (exp(a)^(1/n) exp(b)^(1/n) ...) ^n, or (I suppose) using the messy BCH stuff. Either way it is a commutative combination of the component rotations that, e.g, in 3-d is equivalent to vector addition of axial vectors of the component rotations. So, with some change of basis, an arbitrary rotation in R^n should be representable by a skew-symmetric matrix with only max(floor(n/2)) non-zero components, arranged so that the rotations are in independent coordinate planes, i.e. they have no shared coordinate. Shel At 01:05 PM 7/14/2003 -0700, Eugene Salamin wrote:
Consider, more generally,
[ 0 a -b c] [-a 0 d e] = M. [ b -d 0 f] [-c -e -f 0]
Maple calculates the characteristic polynomial to be
x^4 + (aa+bb+cc+dd+ee+ff) x^2 + (af+cd+be)^2 = 0.
The eigenvalues are +iA, -iA, +iB, -iB, with A and B messy real-valued expressions (assuming a, ..., f are real). R = exp(M) is a 4-dimensional rotation. There is a pair of orthogonal planes P1 and P2 (determined by the eigenvectors of M) such that R is the commutative product of a rotation by angle A in P1 with P2 fixed, and a rotation by angle B in P2 with P1 fixed.
--- Shel Kaphan <sjk@kaphan.org> wrote:
Is anyone's symbolic math package able to take the matrix exponential of
[ 0 a 0 -d ] [ -a 0 b 0 ] [ 0 -b 0 c ] [ d 0 -c 0 ]
?
I set Macsyma to work on it, and it seems only to be trying to melt my laptop and burn the house down after long enough to suspect it isn't going to finish.
If a result of that is available, I'd also like to know the resulting eigenvalues and eigenvectors.
Thanks in advance,
Shel
_______________________________________________ math-fun mailing list math-fun@mailman.xmission.com http://mailman.xmission.com/cgi-bin/mailman/listinfo/math-fun
__________________________________ Do you Yahoo!? SBC Yahoo! DSL - Now only $29.95 per month! http://sbc.yahoo.com
_______________________________________________ math-fun mailing list math-fun@mailman.xmission.com http://mailman.xmission.com/cgi-bin/mailman/listinfo/math-fun