We introduce the use of harmonic analysis to decompose the state space of symmetric robotic systems into orthogonal isotypic subspaces. These are lower-dimensional spaces that capture distinct, symmetric, and synergistic motions. For linear dynamics, we characterize how this decomposition leads to a subdivision of the dynamics into independent linear systems on each subspace, a property we term dynamics harmonic analysis (DHA). To exploit this property, we use Koopman operator theory to propose an equivariant deep-learning architecture that leverages the properties of DHA to learn a global linear model of the system dynamics. Our architecture, validated on synthetic systems and the dynamics of locomotion of a quadrupedal robot, exhibits enhanced generalization, sample efficiency, and interpretability, with fewer trainable parameters and computational costs.
Next to each animation we plot the kinetic energy of each of the configurations on an isotypic subspace (i.e., the kinetic energy of each background robot). This allows us to quantify how relevant each subspace is for the generation of the motion of interest. For instance, for this trotting gait, we can clearly identify that the dynamics evolve primarily in the second and third isotypic subspaces. That is, in a 6-dimensional subspace, instead of the 12-dimensional joint-space.