@@ -50,9 +50,7 @@ def triangle(iterations, dt):
5050 reversal_vals = np .array (reversal_vals )[idx ]
5151 reversal_ts = t [reversal_idxs ]
5252
53- x = np .interp (t , reversal_ts , reversal_vals )
54-
55- return x
53+ return np .interp (t , reversal_ts , reversal_vals )
5654
5755
5856def effective_wheel_radius (v = 20 ):
@@ -164,19 +162,13 @@ def run(duration=4, dt=0.01):
164162 t = np .arange (0 , duration , dt )
165163 iterations = len (t )
166164
167- # hills
168- disturbances = np .array (np .zeros ([3 , iterations + 1 ]))
165+ disturbances = np .array (np .zeros ([3 , iterations + 1 ])) # disturbance from hills
169166 h = hills (iterations + 1 , dt , factor = 0.5 * duration / 2 )
170167 disturbances [2 , :] = h [0 :disturbances .shape [1 ]]
171168
172- # controls
173169 controls = np .array ([0 ])
174-
175- # initial condition
176- state_vals = np .array ([[0 ], [0 ], [0 ]])
177-
178- # desired vel
179- v_r = desired_velocity (iterations , factor = 0.5 * iterations * dt / 2 )
170+ state_vals = np .array ([[0 ], [0 ], [0 ]]) # initial condition
171+ v_r = desired_velocity (iterations , factor = 0.5 * iterations * dt / 2 ) # desired, reference velocity
180172
181173 for i in range (1 , iterations + 1 ):
182174 new_state , u = step_forward (state_vals , disturbances [:, 0 :i ], v_r [0 :i ], dt )
0 commit comments