A special mathematical apparatus for monomial (1,0,-1)–matrices–(4×4) is applied to the description of kinematics and kinetics of rotation in space of an asymmetric rigid body about a fixed point. Using matrix representations for vector operations, Rodrigues vector formula is converted to concise, ordered record in quaternionic matrices. Using transposition operations, forward and indirect coordinate transformations’ matrices when rotating in space are represented by a commutative matrix compositions formed by Rodrigues-Hamilton parameters. The universal algorithm for coordinate transformation matrices formation is offered for an arbitrary sequence of independent rotations. This algorithm was tested on a set of three independent rotations’ systems and their connection with Gibbs vector components was established. Using quaternionic matrices’ properties there were derived invariants of angular velocity, kinetic energy, inertia matrix when rotating. The formulas for time differentiation of vector-function (hodograph), second rank tensor (inertia matrix) in the absolute and local reference systems are represented by equivalent record in quaternionic matrices.