Paparazzi UAS  v7.0_unstable
Paparazzi is a free software Unmanned Aircraft System.
stabilization_attitude_euler_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 
28 #include "generated/airframe.h"
29 
32 
33 #include "std.h"
34 #include "paparazzi.h"
36 #include "state.h"
37 
38 #ifndef USE_ATT_REF
39 #define USE_ATT_REF 1
40 #endif
41 
44 
47 
48 float stabilization_att_fb_cmd[COMMANDS_NB];
49 float stabilization_att_ff_cmd[COMMANDS_NB];
50 
51 #if PERIODIC_TELEMETRY
53 
54 static void send_att(struct transport_tx *trans, struct link_device *dev)
55 {
56  struct FloatRates *body_rate = stateGetBodyRates_f();
57  struct FloatEulers *att = stateGetNedToBodyEulers_f();
58  float foo = 0.0;
59  pprz_msg_send_STAB_ATTITUDE_FLOAT(trans, dev, AC_ID,
60  &(body_rate->p), &(body_rate->q), &(body_rate->r),
61  &(att->phi), &(att->theta), &(att->psi),
68  &stabilization_att_fb_cmd[COMMAND_ROLL],
69  &stabilization_att_fb_cmd[COMMAND_PITCH],
70  &stabilization_att_fb_cmd[COMMAND_YAW],
71  &stabilization_att_ff_cmd[COMMAND_ROLL],
72  &stabilization_att_ff_cmd[COMMAND_PITCH],
73  &stabilization_att_ff_cmd[COMMAND_YAW],
74  &stabilization_cmd[COMMAND_ROLL],
75  &stabilization_cmd[COMMAND_PITCH],
76  &stabilization_cmd[COMMAND_YAW],
77  &foo, &foo, &foo);
78 }
79 
80 static void send_att_ref(struct transport_tx *trans, struct link_device *dev)
81 {
82  pprz_msg_send_STAB_ATTITUDE_REF_FLOAT(trans, dev, AC_ID,
95 }
96 #endif
97 
99 {
100 
102 
104  STABILIZATION_ATTITUDE_PHI_PGAIN,
105  STABILIZATION_ATTITUDE_THETA_PGAIN,
106  STABILIZATION_ATTITUDE_PSI_PGAIN);
107 
109  STABILIZATION_ATTITUDE_PHI_DGAIN,
110  STABILIZATION_ATTITUDE_THETA_DGAIN,
111  STABILIZATION_ATTITUDE_PSI_DGAIN);
112 
114  STABILIZATION_ATTITUDE_PHI_IGAIN,
115  STABILIZATION_ATTITUDE_THETA_IGAIN,
116  STABILIZATION_ATTITUDE_PSI_IGAIN);
117 
119  STABILIZATION_ATTITUDE_PHI_DDGAIN,
120  STABILIZATION_ATTITUDE_THETA_DDGAIN,
121  STABILIZATION_ATTITUDE_PSI_DDGAIN);
122 
124 
125 #if PERIODIC_TELEMETRY
126  register_periodic_telemetry(DefaultPeriodic, PPRZ_MSG_ID_STAB_ATTITUDE_FLOAT, send_att);
127  register_periodic_telemetry(DefaultPeriodic, PPRZ_MSG_ID_STAB_ATTITUDE_REF_FLOAT, send_att_ref);
128 #endif
129 }
130 
131 void stabilization_attitude_read_rc(bool in_flight, bool in_carefree, bool coordinated_turn)
132 {
133  stabilization_attitude_read_rc_setpoint_eulers_f(&stab_att_sp_euler, in_flight, in_carefree, coordinated_turn);
134 }
135 
137 {
138 
139  /* reset psi setpoint to current psi angle */
141 
143 
145 }
146 
148 {
149  stab_att_sp_euler.phi = 0.0;
150  stab_att_sp_euler.theta = 0.0;
152 }
153 
155 {
157 }
158 
160 {
161  struct FloatQuat quat_f;
162  QUAT_FLOAT_OF_BFP(quat_f, *quat);
164 }
165 
167 {
168  struct FloatVect2 cmd_f;
169  cmd_f.x = ANGLE_FLOAT_OF_BFP(cmd->x);
170  cmd_f.y = ANGLE_FLOAT_OF_BFP(cmd->y);
171 
172  /* Rotate horizontal commands to body frame by psi */
173  float psi = stateGetNedToBodyEulers_f()->psi;
174  float s_psi = sinf(psi);
175  float c_psi = cosf(psi);
176  stab_att_sp_euler.phi = -s_psi * cmd_f.x + c_psi * cmd_f.y;
177  stab_att_sp_euler.theta = -c_psi * cmd_f.x - s_psi * cmd_f.y;
179 }
180 
182 {
184 }
185 
186 #define MAX_SUM_ERR 200
187 
188 void stabilization_attitude_run(bool in_flight)
189 {
190 
191 #if USE_ATT_REF
192  static const float dt = (1./PERIODIC_FREQUENCY);
194 #else
198 #endif
199 
200  /* Compute feedforward */
201  stabilization_att_ff_cmd[COMMAND_ROLL] =
203  stabilization_att_ff_cmd[COMMAND_PITCH] =
205  stabilization_att_ff_cmd[COMMAND_YAW] =
207 
208  /* Compute feedback */
209  /* attitude error */
210  struct FloatEulers *att_float = stateGetNedToBodyEulers_f();
211  struct FloatEulers att_err;
212  EULERS_DIFF(att_err, att_ref_euler_f.euler, *att_float);
213  FLOAT_ANGLE_NORMALIZE(att_err.psi);
214 
215  if (in_flight) {
216  /* update integrator */
219  } else {
221  }
222 
223  /* rate error */
224  struct FloatRates *rate_float = stateGetBodyRates_f();
225  struct FloatRates rate_err;
226  RATES_DIFF(rate_err, att_ref_euler_f.rate, *rate_float);
227 
228  /* PID */
229 
230  stabilization_att_fb_cmd[COMMAND_ROLL] =
231  stabilization_gains.p.x * att_err.phi +
232  stabilization_gains.d.x * rate_err.p +
234 
235  stabilization_att_fb_cmd[COMMAND_PITCH] =
236  stabilization_gains.p.y * att_err.theta +
237  stabilization_gains.d.y * rate_err.q +
239 
240  stabilization_att_fb_cmd[COMMAND_YAW] =
241  stabilization_gains.p.z * att_err.psi +
242  stabilization_gains.d.z * rate_err.r +
244 
245 
246  stabilization_cmd[COMMAND_ROLL] =
247  (stabilization_att_fb_cmd[COMMAND_ROLL] + stabilization_att_ff_cmd[COMMAND_ROLL]);
248  stabilization_cmd[COMMAND_PITCH] =
249  (stabilization_att_fb_cmd[COMMAND_PITCH] + stabilization_att_ff_cmd[COMMAND_PITCH]);
250  stabilization_cmd[COMMAND_YAW] =
251  (stabilization_att_fb_cmd[COMMAND_YAW] + stabilization_att_ff_cmd[COMMAND_YAW]);
252 
253  /* bound the result */
254  BoundAbs(stabilization_cmd[COMMAND_ROLL], MAX_PPRZ);
255  BoundAbs(stabilization_cmd[COMMAND_PITCH], MAX_PPRZ);
256  BoundAbs(stabilization_cmd[COMMAND_YAW], MAX_PPRZ);
257 }
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
#define FLOAT_ANGLE_NORMALIZE(_a)
#define FLOAT_EULERS_ZERO(_e)
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 EULERS_BOUND_CUBE(_v, _min, _max)
Definition: pprz_algebra.h:317
#define EULERS_COPY(_a, _b)
Definition: pprz_algebra.h:268
#define EULERS_DIFF(_c, _a, _b)
Definition: pprz_algebra.h:295
#define VECT3_ASSIGN(_a, _x, _y, _z)
Definition: pprz_algebra.h:125
#define EULERS_ADD(_a, _b)
Definition: pprz_algebra.h:281
#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 FloatEulers * stateGetNedToBodyEulers_f(void)
Get vehicle body attitude euler angles (float).
Definition: state.h:1143
static struct FloatRates * stateGetBodyRates_f(void)
Get vehicle body angular rate (float).
Definition: state.h:1200
uint16_t foo
Definition: main_demo5.c:58
#define MAX_PPRZ
Definition: paparazzi.h:8
Paparazzi floating 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 FloatEulers stab_att_sp_euler
with INT32_ANGLE_FRAC
void stabilization_attitude_init(void)
void stabilization_attitude_enter(void)
void stabilization_attitude_set_quat_setpoint_i(struct Int32Quat *quat)
static void send_att_ref(struct transport_tx *trans, struct link_device *dev)
struct AttRefEulerFloat att_ref_euler_f
struct FloatAttitudeGains stabilization_gains
void stabilization_attitude_set_failsafe_setpoint(void)
float stabilization_att_ff_cmd[COMMANDS_NB]
void stabilization_attitude_set_stab_sp(struct StabilizationSetpoint *sp)
void stabilization_attitude_set_rpy_setpoint_i(struct Int32Eulers *rpy)
float stabilization_att_fb_cmd[COMMANDS_NB]
void stabilization_attitude_read_rc(bool in_flight, bool in_carefree, bool coordinated_turn)
struct FloatEulers stabilization_att_sum_err
void stabilization_attitude_set_earth_cmd_i(struct Int32Vect2 *cmd, int32_t heading)
void stabilization_attitude_run(bool in_flight)
static void send_att(struct transport_tx *trans, struct link_device *dev)
float stabilization_attitude_get_heading_f(void)
void stabilization_attitude_read_rc_setpoint_eulers_f(struct FloatEulers *sp, bool in_flight, bool in_carefree, bool coordinated_turn)
Read an attitude setpoint from the RC.
void attitude_ref_euler_float_enter(struct AttRefEulerFloat *ref, float psi)
void attitude_ref_euler_float_update(struct AttRefEulerFloat *ref, struct FloatEulers *sp_eulers, float dt)
void attitude_ref_euler_float_init(struct AttRefEulerFloat *ref)
Attitude reference state/output (euler 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
float heading
Definition: wedgebug.c:258