Paparazzi UAS  v7.0_unstable
Paparazzi is a free software Unmanned Aircraft System.
stabilization_attitude_quat_float.c
Go to the documentation of this file.
1 /*
2  * Copyright (C) 2008-2009 Antoine Drouin <poinix@gmail.com>
3  *
4  * This file is part of paparazzi.
5  *
6  * paparazzi is free software; you can redistribute it and/or modify
7  * it under the terms of the GNU General Public License as published by
8  * the Free Software Foundation; either version 2, or (at your option)
9  * any later version.
10  *
11  * paparazzi is distributed in the hope that it will be useful,
12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14  * GNU General Public License for more details.
15  *
16  * You should have received a copy of the GNU General Public License
17  * along with paparazzi; see the file COPYING. If not, write to
18  * the Free Software Foundation, 59 Temple Place - Suite 330,
19  * Boston, MA 02111-1307, USA.
20  */
21 
26 #include "generated/airframe.h"
27 
31 
32 #include "std.h"
33 #include "paparazzi.h"
35 #include "math/pprz_algebra_int.h"
36 #include "state.h"
37 
39 
42 
44 
46 
48 struct FloatRates body_rate_d;
49 
50 float stabilization_att_fb_cmd[COMMANDS_NB];
51 float stabilization_att_ff_cmd[COMMANDS_NB];
52 
54 
55 static const float phi_pgain[] = STABILIZATION_ATTITUDE_PHI_PGAIN;
56 static const float theta_pgain[] = STABILIZATION_ATTITUDE_THETA_PGAIN;
57 static const float psi_pgain[] = STABILIZATION_ATTITUDE_PSI_PGAIN;
58 
59 static const float phi_dgain[] = STABILIZATION_ATTITUDE_PHI_DGAIN;
60 static const float theta_dgain[] = STABILIZATION_ATTITUDE_THETA_DGAIN;
61 static const float psi_dgain[] = STABILIZATION_ATTITUDE_PSI_DGAIN;
62 
63 static const float phi_igain[] = STABILIZATION_ATTITUDE_PHI_IGAIN;
64 static const float theta_igain[] = STABILIZATION_ATTITUDE_THETA_IGAIN;
65 static const float psi_igain[] = STABILIZATION_ATTITUDE_PSI_IGAIN;
66 
67 static const float phi_ddgain[] = STABILIZATION_ATTITUDE_PHI_DDGAIN;
68 static const float theta_ddgain[] = STABILIZATION_ATTITUDE_THETA_DDGAIN;
69 static const float psi_ddgain[] = STABILIZATION_ATTITUDE_PSI_DDGAIN;
70 
71 static const float phi_dgain_d[] = STABILIZATION_ATTITUDE_PHI_DGAIN_D;
72 static const float theta_dgain_d[] = STABILIZATION_ATTITUDE_THETA_DGAIN_D;
73 static const float psi_dgain_d[] = STABILIZATION_ATTITUDE_PSI_DGAIN_D;
74 
75 
76 #if defined COMMAND_ROLL_SURFACE && defined COMMAND_PITCH_SURFACE && defined COMMAND_YAW_SURFACE
77 #define HAS_SURFACE_COMMANDS 1
78 #endif
79 
80 #ifdef HAS_SURFACE_COMMANDS
81 static const float phi_pgain_surface[] = STABILIZATION_ATTITUDE_PHI_PGAIN_SURFACE;
82 static const float theta_pgain_surface[] = STABILIZATION_ATTITUDE_THETA_PGAIN_SURFACE;
83 static const float psi_pgain_surface[] = STABILIZATION_ATTITUDE_PSI_PGAIN_SURFACE;
84 
85 static const float phi_dgain_surface[] = STABILIZATION_ATTITUDE_PHI_DGAIN_SURFACE;
86 static const float theta_dgain_surface[] = STABILIZATION_ATTITUDE_THETA_DGAIN_SURFACE;
87 static const float psi_dgain_surface[] = STABILIZATION_ATTITUDE_PSI_DGAIN_SURFACE;
88 
89 static const float phi_igain_surface[] = STABILIZATION_ATTITUDE_PHI_IGAIN_SURFACE;
90 static const float theta_igain_surface[] = STABILIZATION_ATTITUDE_THETA_IGAIN_SURFACE;
91 static const float psi_igain_surface[] = STABILIZATION_ATTITUDE_PSI_IGAIN_SURFACE;
92 
93 static const float phi_ddgain_surface[] = STABILIZATION_ATTITUDE_PHI_DDGAIN_SURFACE;
94 static const float theta_ddgain_surface[] = STABILIZATION_ATTITUDE_THETA_DDGAIN_SURFACE;
95 static const float psi_ddgain_surface[] = STABILIZATION_ATTITUDE_PSI_DDGAIN_SURFACE;
96 #endif
97 
98 #define IERROR_SCALE 1024
99 
100 #if PERIODIC_TELEMETRY
102 
103 static void send_att(struct transport_tx *trans, struct link_device *dev)
104 {
105  struct FloatRates *body_rate = stateGetBodyRates_f();
106  struct FloatEulers *att = stateGetNedToBodyEulers_f();
107  pprz_msg_send_STAB_ATTITUDE_FLOAT(trans, dev, AC_ID,
108  &(body_rate->p), &(body_rate->q), &(body_rate->r),
109  &(att->phi), &(att->theta), &(att->psi),
116  &stabilization_att_fb_cmd[COMMAND_ROLL],
117  &stabilization_att_fb_cmd[COMMAND_PITCH],
118  &stabilization_att_fb_cmd[COMMAND_YAW],
119  &stabilization_att_ff_cmd[COMMAND_ROLL],
120  &stabilization_att_ff_cmd[COMMAND_PITCH],
121  &stabilization_att_ff_cmd[COMMAND_YAW],
122  &stabilization_cmd[COMMAND_ROLL],
123  &stabilization_cmd[COMMAND_PITCH],
124  &stabilization_cmd[COMMAND_YAW],
126 }
127 
128 static void send_att_ref(struct transport_tx *trans, struct link_device *dev)
129 {
130  pprz_msg_send_STAB_ATTITUDE_REF_FLOAT(trans, dev, AC_ID,
143 }
144 
145 static void send_ahrs_ref_quat(struct transport_tx *trans, struct link_device *dev)
146 {
147  struct Int32Quat *quat = stateGetNedToBodyQuat_i();
148  struct Int32Quat refquat;
150  pprz_msg_send_AHRS_REF_QUAT(trans, dev, AC_ID,
151  &refquat.qi,
152  &refquat.qx,
153  &refquat.qy,
154  &refquat.qz,
155  &(quat->qi),
156  &(quat->qx),
157  &(quat->qy),
158  &(quat->qz));
159 }
160 #endif
161 
163 {
164  /* setpoints */
167  /* reference */
170 
171  for (int i = 0; i < STABILIZATION_ATTITUDE_GAIN_NB; i++) {
177 #ifdef HAS_SURFACE_COMMANDS
178  VECT3_ASSIGN(stabilization_gains[i].surface_p, phi_pgain_surface[i], theta_pgain_surface[i], psi_pgain_surface[i]);
179  VECT3_ASSIGN(stabilization_gains[i].surface_d, phi_dgain_surface[i], theta_dgain_surface[i], psi_dgain_surface[i]);
180  VECT3_ASSIGN(stabilization_gains[i].surface_i, phi_igain_surface[i], theta_igain_surface[i], psi_igain_surface[i]);
181  VECT3_ASSIGN(stabilization_gains[i].surface_dd, phi_ddgain_surface[i], theta_ddgain_surface[i], psi_ddgain_surface[i]);
182 #endif
183  }
184 
188 
189 #if PERIODIC_TELEMETRY
190  register_periodic_telemetry(DefaultPeriodic, PPRZ_MSG_ID_STAB_ATTITUDE_FLOAT, send_att);
191  register_periodic_telemetry(DefaultPeriodic, PPRZ_MSG_ID_STAB_ATTITUDE_REF_FLOAT, send_att_ref);
193 #endif
194 }
195 
197 {
199  // This could be bad -- Just say no.
200  return;
201  }
202  gain_idx = idx;
204 }
205 
207 {
208 
209  /* reset psi setpoint to current psi angle */
211 
212  struct FloatQuat *state_quat = stateGetNedToBodyQuat_f();
213 
215 
217 }
218 
220 {
221  /* set failsafe to zero roll/pitch and current heading */
222  float heading2 = stabilization_attitude_get_heading_f() / 2;
223  stab_att_sp_quat.qi = cosf(heading2);
224  stab_att_sp_quat.qx = 0.0;
225  stab_att_sp_quat.qy = 0.0;
226  stab_att_sp_quat.qz = sinf(heading2);
227 }
228 
230 {
231  // copy euler setpoint for debugging
233 
235 }
236 
238 {
241 }
242 
244 {
245  struct FloatVect2 cmd_f;
246  cmd_f.x = ANGLE_FLOAT_OF_BFP(cmd->x);
247  cmd_f.y = ANGLE_FLOAT_OF_BFP(cmd->y);
248  float heading_f;
249  heading_f = ANGLE_FLOAT_OF_BFP(heading);
250 
251  quat_from_earth_cmd_f(&stab_att_sp_quat, &cmd_f, heading_f);
252 }
253 
255 {
258 }
259 
260 #ifndef GAIN_PRESCALER_FF
261 #define GAIN_PRESCALER_FF 1
262 #endif
263 static void attitude_run_ff(float ff_commands[], struct FloatAttitudeGains *gains, struct FloatRates *ref_accel)
264 {
265  /* Compute feedforward based on reference acceleration */
266 
267  ff_commands[COMMAND_ROLL] = GAIN_PRESCALER_FF * gains->dd.x * ref_accel->p;
268  ff_commands[COMMAND_PITCH] = GAIN_PRESCALER_FF * gains->dd.y * ref_accel->q;
269  ff_commands[COMMAND_YAW] = GAIN_PRESCALER_FF * gains->dd.z * ref_accel->r;
270 #ifdef HAS_SURFACE_COMMANDS
271  ff_commands[COMMAND_ROLL_SURFACE] = GAIN_PRESCALER_FF * gains->surface_dd.x * ref_accel->p;
272  ff_commands[COMMAND_PITCH_SURFACE] = GAIN_PRESCALER_FF * gains->surface_dd.y * ref_accel->q;
273  ff_commands[COMMAND_YAW_SURFACE] = GAIN_PRESCALER_FF * gains->surface_dd.z * ref_accel->r;
274 #endif
275 }
276 
277 #ifndef GAIN_PRESCALER_P
278 #define GAIN_PRESCALER_P 1
279 #endif
280 #ifndef GAIN_PRESCALER_D
281 #define GAIN_PRESCALER_D 1
282 #endif
283 #ifndef GAIN_PRESCALER_I
284 #define GAIN_PRESCALER_I 1
285 #endif
286 static void attitude_run_fb(float fb_commands[], struct FloatAttitudeGains *gains, struct FloatQuat *att_err,
287  struct FloatRates *rate_err, struct FloatRates *rate_err_d, struct FloatQuat *sum_err)
288 {
289  /* PID feedback */
290  fb_commands[COMMAND_ROLL] =
291  GAIN_PRESCALER_P * gains->p.x * att_err->qx +
292  GAIN_PRESCALER_D * gains->d.x * rate_err->p +
293  GAIN_PRESCALER_D * gains->rates_d.x * rate_err_d->p +
294  GAIN_PRESCALER_I * gains->i.x * sum_err->qx;
295 
296  fb_commands[COMMAND_PITCH] =
297  GAIN_PRESCALER_P * gains->p.y * att_err->qy +
298  GAIN_PRESCALER_D * gains->d.y * rate_err->q +
299  GAIN_PRESCALER_D * gains->rates_d.y * rate_err_d->q +
300  GAIN_PRESCALER_I * gains->i.y * sum_err->qy;
301 
302  fb_commands[COMMAND_YAW] =
303  GAIN_PRESCALER_P * gains->p.z * att_err->qz +
304  GAIN_PRESCALER_D * gains->d.z * rate_err->r +
305  GAIN_PRESCALER_D * gains->rates_d.z * rate_err_d->r +
306  GAIN_PRESCALER_I * gains->i.z * sum_err->qz;
307 
308 #ifdef HAS_SURFACE_COMMANDS
309  fb_commands[COMMAND_ROLL_SURFACE] =
310  GAIN_PRESCALER_P * gains->surface_p.x * att_err->qx +
311  GAIN_PRESCALER_D * gains->surface_d.x * rate_err->p +
312  GAIN_PRESCALER_I * gains->surface_i.x * sum_err->qx;
313 
314  fb_commands[COMMAND_PITCH_SURFACE] =
315  GAIN_PRESCALER_P * gains->surface_p.y * att_err->qy +
316  GAIN_PRESCALER_D * gains->surface_d.y * rate_err->q +
317  GAIN_PRESCALER_I * gains->surface_i.y * sum_err->qy;
318 
319  fb_commands[COMMAND_YAW_SURFACE] =
320  GAIN_PRESCALER_P * gains->surface_p.z * att_err->qz +
321  GAIN_PRESCALER_D * gains->surface_d.z * rate_err->r +
322  GAIN_PRESCALER_I * gains->surface_i.z * sum_err->qz;
323 #endif
324 }
325 
326 void stabilization_attitude_run(bool enable_integrator)
327 {
328 
329  /*
330  * Update reference
331  */
332  static const float dt = (1./PERIODIC_FREQUENCY);
334 
335  /*
336  * Compute errors for feedback
337  */
338 
339  /* attitude error */
340  struct FloatQuat att_err;
341  struct FloatQuat *att_quat = stateGetNedToBodyQuat_f();
342  float_quat_inv_comp(&att_err, att_quat, &att_ref_quat_f.quat);
343  /* wrap it in the shortest direction */
344  float_quat_wrap_shortest(&att_err);
345 
346  /* rate error */
347  struct FloatRates rate_err;
348  struct FloatRates *body_rate = stateGetBodyRates_f();
349  RATES_DIFF(rate_err, att_ref_quat_f.rate, *body_rate);
350  /* rate_d error */
351  RATES_DIFF(body_rate_d, *body_rate, last_body_rate);
352  RATES_COPY(last_body_rate, *body_rate);
353 
354 #define INTEGRATOR_BOUND 1.0
355  /* integrated error */
356  if (enable_integrator) {
363  } else {
364  /* reset accumulator */
366  }
367 
369 
372 
373  stabilization_cmd[COMMAND_ROLL] = stabilization_att_fb_cmd[COMMAND_ROLL] + stabilization_att_ff_cmd[COMMAND_ROLL];
374  stabilization_cmd[COMMAND_PITCH] = stabilization_att_fb_cmd[COMMAND_PITCH] + stabilization_att_ff_cmd[COMMAND_PITCH];
375  stabilization_cmd[COMMAND_YAW] = stabilization_att_fb_cmd[COMMAND_YAW] + stabilization_att_ff_cmd[COMMAND_YAW];
376 
377 #ifdef HAS_SURFACE_COMMANDS
378  stabilization_cmd[COMMAND_ROLL_SURFACE] = stabilization_att_fb_cmd[COMMAND_ROLL_SURFACE] +
379  stabilization_att_ff_cmd[COMMAND_ROLL_SURFACE];
380  stabilization_cmd[COMMAND_PITCH_SURFACE] = stabilization_att_fb_cmd[COMMAND_PITCH_SURFACE] +
381  stabilization_att_ff_cmd[COMMAND_PITCH_SURFACE];
382  stabilization_cmd[COMMAND_YAW_SURFACE] = stabilization_att_fb_cmd[COMMAND_YAW_SURFACE] +
383  stabilization_att_ff_cmd[COMMAND_YAW_SURFACE];
384 #endif
385 
386  /* bound the result */
387  BoundAbs(stabilization_cmd[COMMAND_ROLL], MAX_PPRZ);
388  BoundAbs(stabilization_cmd[COMMAND_PITCH], MAX_PPRZ);
389  BoundAbs(stabilization_cmd[COMMAND_YAW], MAX_PPRZ);
390 }
391 
392 void stabilization_attitude_read_rc(bool in_flight, bool in_carefree, bool coordinated_turn)
393 {
394 #if USE_EARTH_BOUND_RC_SETPOINT
395  stabilization_attitude_read_rc_setpoint_quat_earth_bound_f(&stab_att_sp_quat, in_flight, in_carefree, coordinated_turn);
396 #else
397  stabilization_attitude_read_rc_setpoint_quat_f(&stab_att_sp_quat, in_flight, in_carefree, coordinated_turn);
398 #endif
399  //float_quat_wrap_shortest(&stab_att_sp_quat);
400 }
float q
in rad/s
float phi
in radians
float p
in rad/s
float r
in rad/s
float theta
in radians
float psi
in radians
static void float_quat_identity(struct FloatQuat *q)
initialises a quaternion to identity
void float_quat_inv_comp(struct FloatQuat *b2c, struct FloatQuat *a2b, struct FloatQuat *a2c)
Composition (multiplication) of two quaternions.
#define FLOAT_EULERS_ZERO(_e)
static void float_quat_wrap_shortest(struct FloatQuat *q)
void float_quat_of_eulers(struct FloatQuat *q, struct FloatEulers *e)
quat of euler roation 'ZYX'
void float_eulers_of_quat(struct FloatEulers *e, struct FloatQuat *q)
euler rotation 'ZYX'
#define FLOAT_RATES_ZERO(_r)
euler angles
Roation quaternion.
angular rates
#define RATES_COPY(_a, _b)
Definition: pprz_algebra.h:337
#define QUAT_BFP_OF_REAL(_qi, _qf)
Definition: pprz_algebra.h:752
#define VECT3_ASSIGN(_a, _x, _y, _z)
Definition: pprz_algebra.h:125
#define RATES_DIFF(_c, _a, _b)
Definition: pprz_algebra.h:372
#define QUAT_FLOAT_OF_BFP(_qf, _qi)
Definition: pprz_algebra.h:745
#define EULERS_FLOAT_OF_BFP(_ef, _ei)
Definition: pprz_algebra.h:709
#define ANGLE_FLOAT_OF_BFP(_ai)
euler angles
Rotation quaternion.
static struct Int32Quat * stateGetNedToBodyQuat_i(void)
Get vehicle body attitude quaternion (int).
Definition: state.h:1113
static struct FloatEulers * stateGetNedToBodyEulers_f(void)
Get vehicle body attitude euler angles (float).
Definition: state.h:1143
static struct FloatQuat * stateGetNedToBodyQuat_f(void)
Get vehicle body attitude quaternion (float).
Definition: state.h:1131
static struct FloatRates * stateGetBodyRates_f(void)
Get vehicle body angular rate (float).
Definition: state.h:1200
static float p[2][2]
static uint32_t idx
float gains[MAX_SAMPLES_LEARNING]
#define MAX_PPRZ
Definition: paparazzi.h:8
Paparazzi floating point algebra.
Paparazzi fixed point algebra.
General attitude stabilization interface for rotorcrafts.
struct FloatEulers stab_sp_to_eulers_f(struct StabilizationSetpoint *sp)
int32_t stabilization_cmd[COMMANDS_NB]
Stabilization commands.
Definition: stabilization.c:34
struct FloatQuat stab_sp_to_quat_f(struct StabilizationSetpoint *sp)
struct FloatEulers stab_att_sp_euler
with INT32_ANGLE_FRAC
void stabilization_attitude_init(void)
stabilization_attitude_init
static void send_ahrs_ref_quat(struct transport_tx *trans, struct link_device *dev)
struct FloatAttitudeGains stabilization_gains[STABILIZATION_ATTITUDE_GAIN_NB]
static void attitude_run_ff(float ff_commands[], struct FloatAttitudeGains *gains, struct FloatRates *ref_accel)
static const float theta_ddgain[]
void stabilization_attitude_run(bool enable_integrator)
struct FloatRates last_body_rate
struct AttRefQuatFloat att_ref_quat_f
void stabilization_attitude_enter(void)
void stabilization_attitude_set_quat_setpoint_i(struct Int32Quat *quat)
static const float theta_igain[]
static const float phi_dgain_d[]
static void send_att_ref(struct transport_tx *trans, struct link_device *dev)
static const float theta_pgain[]
void stabilization_attitude_set_failsafe_setpoint(void)
float stabilization_att_ff_cmd[COMMANDS_NB]
static const float psi_igain[]
static const float theta_dgain_d[]
static const float phi_ddgain[]
static const float phi_igain[]
static const float psi_pgain[]
struct FloatQuat stabilization_att_sum_err_quat
static int gain_idx
struct FloatQuat stab_att_sp_quat
with INT32_QUAT_FRAC
static const float phi_pgain[]
void stabilization_attitude_set_stab_sp(struct StabilizationSetpoint *sp)
static const float theta_dgain[]
static const float psi_ddgain[]
static void attitude_run_fb(float fb_commands[], struct FloatAttitudeGains *gains, struct FloatQuat *att_err, struct FloatRates *rate_err, struct FloatRates *rate_err_d, struct FloatQuat *sum_err)
static const float phi_dgain[]
void stabilization_attitude_set_rpy_setpoint_i(struct Int32Eulers *rpy)
struct FloatRates body_rate_d
static const float psi_dgain[]
float stabilization_att_fb_cmd[COMMANDS_NB]
void stabilization_attitude_read_rc(bool in_flight, bool in_carefree, bool coordinated_turn)
void stabilization_attitude_gain_schedule(uint8_t idx)
void stabilization_attitude_set_earth_cmd_i(struct Int32Vect2 *cmd, int32_t heading)
static const float psi_dgain_d[]
#define INTEGRATOR_BOUND
static void send_att(struct transport_tx *trans, struct link_device *dev)
#define STABILIZATION_ATTITUDE_GAIN_IDX_DEFAULT
#define STABILIZATION_ATTITUDE_GAIN_NB
void quat_from_earth_cmd_f(struct FloatQuat *quat, struct FloatVect2 *cmd, float heading)
Quaternion transformation functions.
void stabilization_attitude_read_rc_setpoint_quat_earth_bound_f(struct FloatQuat *q_sp, bool in_flight, bool in_carefree, bool coordinated_turn)
void stabilization_attitude_read_rc_setpoint_quat_f(struct FloatQuat *q_sp, bool in_flight, bool in_carefree, bool coordinated_turn)
Read attitude setpoint from RC as quaternion Interprets the stick positions as axes.
float stabilization_attitude_get_heading_f(void)
Read an attitude setpoint from the RC.
void attitude_ref_quat_float_update(struct AttRefQuatFloat *ref, struct FloatQuat *sp_quat, float dt)
void attitude_ref_quat_float_enter(struct AttRefQuatFloat *ref, struct FloatQuat *state_quat)
void attitude_ref_quat_float_init(struct AttRefQuatFloat *ref)
void attitude_ref_quat_float_schedule(struct AttRefQuatFloat *ref, uint8_t idx)
Attitude reference models and state/output (float)
API to get/set the generic vehicle states.
static const struct usb_device_descriptor dev
Definition: usb_ser_hw.c:74
Stabilization setpoint.
Definition: stabilization.h:42
int8_t register_periodic_telemetry(struct periodic_telemetry *_pt, uint8_t _id, telemetry_cb _cb)
Register a telemetry callback function.
Definition: telemetry.c:51
Periodic telemetry system header (includes downlink utility and generated code).
#define DefaultPeriodic
Set default periodic telemetry.
Definition: telemetry.h:66
int int32_t
Typedef defining 32 bit int type.
Definition: vl53l1_types.h:83
unsigned char uint8_t
Typedef defining 8 bit unsigned char type.
Definition: vl53l1_types.h:98
float heading
Definition: wedgebug.c:258