A common approach to understanding and controlling robotic legged locomotion is the construction and analysis of simplified mathematical models that capture essential features of locomotor behaviours. However, the representational power of such simple mathematical models is inevitably limited due to the non-linear and complex nature of biological locomotor systems. Attempting to identify and explicitly incorporate key non-linearities into the model is challenging, increases complexity, and decreases the analytic utility of the resulting models. In this paper, we adopt a data-driven approach, with the goal of furnishing an input–output representation of a locomotor system. Our method is based on approximating the hybrid dynamics of a legged locomotion model around its limit cycle as a Linear Time Periodic (LTP) system. Perturbing inputs to the locomotor system with small chirp signals yield the input–output data necessary for the application of LTP system identification techniques, allowing us to estimate harmonic transfer functions (HTFs) associated with the local LTP approximation to the system dynamics around the limit cycle. We compare actual system responses with responses predicted by the HTF, providing evidence that data-driven system identification methods can be used to construct models for locomotor behaviours.