...
where L1-L12 are the link lengths, phuman(tk) is the ankle position of human gait motion capture data at time frame k, pjansen(tk,L) is the kinematic end point of Theo Jansen mechanism at time frame k, Φ Φ is the rotation angle of the ground and τ is the phase lag.
The optimization results are shown below.