In this paper, we propose a variation-based method to linearize the nonlinear dynamics of robotic systems, whose configuration spaces contain the manifolds S and SO(3), along dynamically-feasible reference trajectories. The proposed variation-based linearization results in an implicitly time-varying linear system, representing the error dynamics, that is globally valid. We illustrate this metho...