1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135

50

In general, there is no guarantee that the discrete system Jacobian will be invertible for each step.

However, with well-chosen choices of RVs and LPPs, singularities generally occur only for

trajectories far enough from useful limit cycles to be of little interest.

Typical singularities for the

human model, occur when it is lying on its side and when the current stance leg does not contact

the ground at all during the current step.


The general form of the Jacobian for a 2D system is:

IMAGE Imgs/thesis.final.w683.gif

J=

(3.10)

Using finite differences, a minimum of N+1 samples are required to construct the linear model of

h(K) for a control system of dimension N: one for the nominal operating point and one for each

control dimension, each of which yields a column in the Jacobian.

In this case, each "sample"

consists of a simulation of one step with a different value of K.

samples may be required.

In practice, a greater number of

For the two dimensions of our bipedal control system, we use four sample simulations, two for

each dimension.

An additional simulation computes the final motion for the step after the final

PCG scalings have been chosen using the model.

choose to work with two simplified forms:

Rather than using the complete Jacobian, we

  1. Assume independent control dimensions and observe only the primary RV.

This

corresponds to the assumption of a diagonal Jacobian (i.e. ðqi / ðkj= 0, i[!]j).

refer to this form as superposition(SP) sampling.


We will

  1. Assume near-independent control dimensions but allow the operating point to move as

each final perturbation scaling is determined. This corresponds to a form of triangular

[CONVERTED BY MYRMIDON]