|
| 1 | +import { vec3, vec4, mat4, quat } from 'https://esm.sh/gl-matrix'; |
| 2 | + |
| 3 | +export default class Controller{ |
| 4 | + constructor(){ |
| 5 | + } |
| 6 | + evaluate_step(state){ |
| 7 | + const q_des = [0, 0, 0, 1]; |
| 8 | + const stateDict = JSON.parse(state.get_state()); |
| 9 | + const params = JSON.parse(state.get_parameters()); |
| 10 | + |
| 11 | + const p = vec3.fromValues(...stateDict['position']); |
| 12 | + const q = quat.fromValues( |
| 13 | + stateDict['orientation'][1], |
| 14 | + stateDict['orientation'][2], |
| 15 | + stateDict['orientation'][3], |
| 16 | + stateDict['orientation'][0] |
| 17 | + ); |
| 18 | + const v = vec3.fromValues(...stateDict['linear_velocity']); |
| 19 | + const omega = vec3.fromValues(...stateDict['angular_velocity']); |
| 20 | + |
| 21 | + const m = params['dynamics']['mass']; |
| 22 | + const g = -params['dynamics']['gravity'][2]; |
| 23 | + const l = 0.028; |
| 24 | + const k_q = params['dynamics']['rotor_torque_constants'][0]; |
| 25 | + const [a, b, c] = params['dynamics']['rotor_thrust_coefficients'][0]; |
| 26 | + |
| 27 | + const k_p = 3.0; |
| 28 | + const k_d = 2.0; |
| 29 | + const k_R = 1.0; |
| 30 | + const k_omega = 0.1; |
| 31 | + |
| 32 | + const p_des = vec3.fromValues(0, 0, 0); |
| 33 | + const v_des = vec3.fromValues(0, 0, 0); |
| 34 | + |
| 35 | + const e_p = vec3.create(); |
| 36 | + vec3.subtract(e_p, p, p_des); |
| 37 | + const e_v = vec3.create(); |
| 38 | + vec3.subtract(e_v, v, v_des); |
| 39 | + |
| 40 | + const a_des = vec3.create(); |
| 41 | + vec3.scale(a_des, e_p, -k_p); |
| 42 | + const temp = vec3.create(); |
| 43 | + vec3.scale(temp, e_v, -k_d); |
| 44 | + vec3.add(a_des, a_des, temp); |
| 45 | + |
| 46 | + const T = m * (a_des[2] + g); |
| 47 | + |
| 48 | + const q_des_conj = quat.create(); |
| 49 | + quat.conjugate(q_des_conj, q_des); |
| 50 | + const q_err = quat.create(); |
| 51 | + quat.multiply(q_err, q_des_conj, q); |
| 52 | + |
| 53 | + if (q_err[3] < 0){ |
| 54 | + quat.scale(q_err, q_err, -1); |
| 55 | + } |
| 56 | + |
| 57 | + const e_R = vec3.fromValues(q_err[0], q_err[1], q_err[2]); |
| 58 | + vec3.scale(e_R, e_R, 2); |
| 59 | + |
| 60 | + const tau = vec3.create(); |
| 61 | + vec3.scale(tau, e_R, -k_R); |
| 62 | + const temp2 = vec3.create(); |
| 63 | + vec3.scale(temp2, omega, -k_omega); |
| 64 | + vec3.add(tau, tau, temp2); |
| 65 | + |
| 66 | + const A_inv = mat4.fromValues( |
| 67 | + 1/4, -1/(4*l), -1/(4*l), -1/(4*k_q), |
| 68 | + 1/4, -1/(4*l), 1/(4*l), 1/(4*k_q), |
| 69 | + 1/4, 1/(4*l), 1/(4*l), -1/(4*k_q), |
| 70 | + 1/4, 1/(4*l), -1/(4*l), 1/(4*k_q) |
| 71 | + ); |
| 72 | + const controlInputs = vec4.fromValues(T, tau[0], tau[1], tau[2]); |
| 73 | + const f = vec4.create(); |
| 74 | + vec4.transformMat4(f, controlInputs, A_inv); |
| 75 | + |
| 76 | + const f_clipped = [ |
| 77 | + Math.max(0, f[0]), |
| 78 | + Math.max(0, f[1]), |
| 79 | + Math.max(0, f[2]), |
| 80 | + Math.max(0, f[3]) |
| 81 | + ]; |
| 82 | + |
| 83 | + function solveRPM(fi) { |
| 84 | + const discriminant = b * b - 4 * a * (c - fi); |
| 85 | + if (discriminant < 0 || fi < c) { |
| 86 | + return 0.0; |
| 87 | + } |
| 88 | + const rpm = (-b + Math.sqrt(discriminant)) / (2 * a); |
| 89 | + return Math.max(0.0, Math.min(1.0, rpm)); |
| 90 | + } |
| 91 | + |
| 92 | + const action = f_clipped.map(solveRPM); |
| 93 | + |
| 94 | + return action; |
| 95 | + } |
| 96 | + reset() { |
| 97 | + } |
| 98 | +} |
| 99 | + |
0 commit comments