Maths - Skew Symmetric Matrix

A matrix is skew symmetric if its elements meet the following rule:

a_ij = - a_ji

The leading diagonal terms must be zero since in this case a= -a which is only true when a=0.

We can express this as:

[A]^t = -[A]

Representing Vector cross Multiplication

We want to find a 3x3 matrix which is equivalent to vector cross multiplication (as described here). Note: cross multiplication only applies to 3 dimensional vectors and represents a vector which is perpendicular to both the vectors being multiplied.

In other words, if,

C = A × B

where:

• A,B and C are 3 dimensional vectors
• × is the cross product

Then we want to find:

C = [~A] B

where:

[~A] = 3x3 Skew Symmetric Matrix which we want to find

The components of the cross product are (as described here):

Cx = Ay * Bz - By * Az
Cy = Az * Bx - Bz * Ax
Cz = Ax * By - Bx * Ay

So [~A] is equivalent to A x

[~A] =

0 -Az Ay Az 0 -Ax -Ay Ax 0

Properties of Skew Symmetric Matrix

Jacobis theorem

There is no inverse of skew symmetric matrix in the form used to represent cross multiplication (or any odd dimension skew symmetric matrix), if there were then we would be able to get an inverse for the vector cross product but this is not possible.

To demonstrate this we can try to calculate its inverse as follows:

We can use the expression for inverse given here and apply it to skew symmetric matrix to get:

[~A] =1/det[m] *

Ax*Ax Ay*Ax Az*Ax Ax*Ay Ay*Ay Az*Ay Ax*Az Ay*Az Az*Az

where det = Ay*Az*Ax -Az*Ax*Ay = 0

So all terms tend towards 1/0 = infinity

'tilde' matrix

If we have a vector cross multiplication, for instance:

= x

We can replace this expession by an equivilant matrix transform:

= [~r]

Where [~r] is a square matrix derived as follows, if we expand out the terms for vector cross multiplation we get:

v_x = r_y * w_z - w_y * r_z
v_y = r_z * w_x - w_z * r_x
v_z = r_x * w_y - w_x * r_y

we can rearange these terms to be in matrix form as follows:

v_x = 0 * w_x - r_z* w_y + r_y * w_z
v_y = r_z * w_x + 0 * w_y - r_x * w_z
v_z =-r_y * w_x + r_x * w_y + 0 * w_z

which gives:

=

0 - rz ry rz 0 - rx -ry rx 0

so

[~r] =

0 - rz ry rz 0 - rx -ry rx 0

If instead of we want to turn into the matrix, then we can use the following identiy for vectors (reversing the order of operands changes the sign) :

x = - ( x )

Therefore if = x

then = -[~w], where:

-[~w] =

0 wz -wy -wz 0 wx wy -wx 0

Vector Cross Product

Linear Functions

Standards