Current research on discrete and rhythmic movements differs in both experimental procedures and theory, despite the ubiquitous overlap between discrete and rhythmic components in everyday behaviors. Models of rhythmic movements usually use oscillatory systems mimicking central pattern generators (CPGs). In contrast, models of discrete movements often employ optimization principles, thereby reflecting the higher-level cortical resources involved in the generation of such movements. This letter proposes a unified model for the generation of both rhythmic and discrete movements. We show that a physiologically motivated model of a CPG can not only generate simple rhythmic movements with only a small set of parameters, but can also produce discrete movements if the CPG is fed with an exponentially decaying phasic input. We further show that a particular coupling between two of these units can reproduce main findings on in-phase and antiphase stability. Finally, we propose an integrated model of combined rhythmic and discrete movements for the two hands. These movement classes are sequentially addressed in this letter with increasing model complexity. The model variations are discussed in relation to the degree of recruitment of the higher-level cortical resources, necessary for such movements.

The question of how best to model rhythmic movements at self-selected amplitude-frequency combinations, and their variability, is a long-standing issue. This study presents a systematic analysis of a coupled oscillator system that has successfully accounted for the experimental result that humans' preferred oscillation frequencies closely correspond to the linear resonance frequencies of the biomechanical limb systems, a phenomenon known as resonance tuning or frequency scaling. The dynamics of the coupled oscillator model is explored by numerical integration in different areas of its parameter space, where a period doubling route to chaotic dynamics is discovered. It is shown that even in the regions of the parameter space with chaotic solutions, the model still effectively scales to the biomechanical oscillator's natural frequency. Hence, there is a solution providing for frequency scaling in the presence of chaotic variability. The implications of these results for interpreting variability as fundamentally stochastic or chaotic are discussed.