Divergent components of motion
Abstract
Some nonlinear control systems admit an exponential dichotomy (Coppel, 1966), that is to say, their dynamics can be decomposed into (exponentially) stable and unstable components. Walking robots fall into this category, and we call their unstable components divergent components of motion (DCM). The concept of DCM has been fruitfully applied to the linear inverted pendulum (LIP) for both walking pattern generation and balance feedback control. But DCMs can be found for other models as well! In this talk, we will discuss DCMs for the variableheight inverted pendulum (VHIP), an extension of the LIP where the controller can add height variations. Ideally, we would like our robot to behave as a LIP (nominal height) unless some perturbation occurs and the robot resorts to the heightvariation strategy, if it has to. Deciding when to use or not this strategy may seem "smart" or predictive, but we will see that it can be implemented straightforwardly as linear feedback over a 4D DCM.
Videos
Chatty Presentation
TL;DW
Content
Slides 
References
Feedback control of a 4D DCM for the variableheight inverted pendulum  
Walking trajectory generation with height variations 
Discussion
Thanks to all those who have contributed to the conversation so far. Feel free to leave a reply using the form below, or subscribe to the Discussion's atom feed to stay tuned.

Attendee #1
Posted on
What matrices did you use to generate the figure on slide 4?

Stéphane
Posted on
This figure corresponds to:
 \(A = \begin{bmatrix} +2 & 1 \\ 0 & +1 \end{bmatrix}\) for \(\mathrm{eig}(A) = \{2, 1\}\)
 \(A = \begin{bmatrix} 2 & 1 \\ 0 & +1 \end{bmatrix}\) for \(\mathrm{eig}(A) = \{2, 1\}\)
 \(A = \begin{bmatrix} 2 & 1 \\ 0 & 1 \end{bmatrix}\) for \(\mathrm{eig}(A) = \{2, 1\}\)


Attendee #2
Posted on
Why did you seem to doubt that \(\omega\) is a DCM, isn't it clearly divergent?

Stéphane
Posted on
Yes, the point I had doubts on is about the "of motion" part. Previously, when the DCM was directly computed by linear combination of the CoM position and velocity, it was clear that "it diverges" and "it is a component of motion" imply that it is a DCM. But here, \(\omega\) appears as a technical choice we make in order to diagonalize the statetransition matrix after changing variable.


Attendee #3
Posted on
How do you choose the remaining proportional gain on slide 12?

Stéphane
Posted on
On the real robot, it will depend on your control cycle and in particular on the bandwidth of the force control loop (admittance control on our robots, see slide 26). DCM and force control gains are coupled when the two are run at roughly the same frequency, as is the case here, and we are not modeling this coupling. For practical advice, check out this note on tuning stabilizer gains.


Attendee #4
Posted on
How did you select the poles in the final leastsquares formulation?

Stéphane
Posted on
In general we could have four gains on the diagonal of the closedloop statetransition matrix, but in practice we often use the same gains for different directions (for instance, the same gain for both sagittal and lateral DCM feedback in the LIP). I followed this practice, using a single normalized gain \(k > 1\) and scaling it on each axis by a factor consistent with the equations of motion.


Attendee #5
Posted on
Is there unicity of the DCMs or exponential dichotomy?

Stéphane
Posted on
No! For instance, in a previous work we had used a different DCM for the VHIP whose formula included the ZMP as well. Multiplying a DCM by a nonzero scalar also yields a DCM, there may be "classes" of equivalent DCMs for some equivalence relation, but I wonder what it could be...
