Paparazzi UAS  v5.8.2_stable-0-g6260b7c
Paparazzi is a free software Unmanned Aircraft System.
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Modules Pages
ahrs_int_cmpl_quat_wrapper.c
Go to the documentation of this file.
1 /*
2  * Copyright (C) 2015 Felix Ruess <felix.ruess@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, see
18  * <http://www.gnu.org/licenses/>.
19  */
20 
28 #include "subsystems/ahrs.h"
29 #include "subsystems/abi.h"
30 #include "state.h"
31 
32 #ifndef AHRS_ICQ_OUTPUT_ENABLED
33 #define AHRS_ICQ_OUTPUT_ENABLED TRUE
34 #endif
35 PRINT_CONFIG_VAR(AHRS_ICQ_OUTPUT_ENABLED)
36 
37 
41 
42 static void set_body_state_from_quat(void);
43 
44 #if PERIODIC_TELEMETRY
46 #include "mcu_periph/sys_time.h"
47 #include "state.h"
48 
49 static void send_quat(struct transport_tx *trans, struct link_device *dev)
50 {
51  struct Int32Quat *quat = stateGetNedToBodyQuat_i();
52  pprz_msg_send_AHRS_QUAT_INT(trans, dev, AC_ID,
58  &(quat->qi),
59  &(quat->qx),
60  &(quat->qy),
61  &(quat->qz),
62  &ahrs_icq_id);
63 }
64 
65 static void send_euler(struct transport_tx *trans, struct link_device *dev)
66 {
67  struct Int32Eulers ltp_to_imu_euler;
68  int32_eulers_of_quat(&ltp_to_imu_euler, &ahrs_icq.ltp_to_imu_quat);
69  struct Int32Eulers *eulers = stateGetNedToBodyEulers_i();
70  pprz_msg_send_AHRS_EULER_INT(trans, dev, AC_ID,
71  &ltp_to_imu_euler.phi,
72  &ltp_to_imu_euler.theta,
73  &ltp_to_imu_euler.psi,
74  &(eulers->phi),
75  &(eulers->theta),
76  &(eulers->psi),
77  &ahrs_icq_id);
78 }
79 
80 static void send_bias(struct transport_tx *trans, struct link_device *dev)
81 {
82  pprz_msg_send_AHRS_GYRO_BIAS_INT(trans, dev, AC_ID,
85 }
86 
87 static void send_geo_mag(struct transport_tx *trans, struct link_device *dev)
88 {
89  struct FloatVect3 h_float;
90  h_float.x = MAG_FLOAT_OF_BFP(ahrs_icq.mag_h.x);
91  h_float.y = MAG_FLOAT_OF_BFP(ahrs_icq.mag_h.y);
92  h_float.z = MAG_FLOAT_OF_BFP(ahrs_icq.mag_h.z);
93  pprz_msg_send_GEO_MAG(trans, dev, AC_ID,
94  &h_float.x, &h_float.y, &h_float.z, &ahrs_icq_id);
95 }
96 
97 static void send_filter_status(struct transport_tx *trans, struct link_device *dev)
98 {
99  uint8_t mde = 3;
100  uint16_t val = 0;
101  if (!ahrs_icq.is_aligned) { mde = 2; }
103  /* set lost if no new gyro measurements for 50ms */
104  if (t_diff > 50000) { mde = 5; }
105  pprz_msg_send_STATE_FILTER_STATUS(trans, dev, AC_ID, &ahrs_icq_id, &mde, &val);
106 }
107 #endif
108 
109 
113 #ifndef AHRS_ICQ_IMU_ID
114 #define AHRS_ICQ_IMU_ID ABI_BROADCAST
115 #endif
116 PRINT_CONFIG_VAR(AHRS_ICQ_IMU_ID)
120 #ifndef AHRS_ICQ_MAG_ID
121 #define AHRS_ICQ_MAG_ID AHRS_ICQ_IMU_ID
122 #endif
123 PRINT_CONFIG_VAR(AHRS_ICQ_MAG_ID)
131 
132 
133 static void gyro_cb(uint8_t __attribute__((unused)) sender_id,
134  uint32_t stamp, struct Int32Rates *gyro)
135 {
136  ahrs_icq_last_stamp = stamp;
137 #if USE_AUTO_AHRS_FREQ || !defined(AHRS_PROPAGATE_FREQUENCY)
138  PRINT_CONFIG_MSG("Calculating dt for AHRS_ICQ propagation.")
139  /* timestamp in usec when last callback was received */
140  static uint32_t last_stamp = 0;
141 
142  if (last_stamp > 0 && ahrs_icq.is_aligned) {
143  float dt = (float)(stamp - last_stamp) * 1e-6;
144  ahrs_icq_propagate(gyro, dt);
146  }
147  last_stamp = stamp;
148 #else
149  PRINT_CONFIG_MSG("Using fixed AHRS_PROPAGATE_FREQUENCY for AHRS_ICQ propagation.")
150  PRINT_CONFIG_VAR(AHRS_PROPAGATE_FREQUENCY)
152  const float dt = 1. / (AHRS_PROPAGATE_FREQUENCY);
153  ahrs_icq_propagate(gyro, dt);
155  }
156 #endif
157 }
158 
159 static void accel_cb(uint8_t __attribute__((unused)) sender_id,
160  uint32_t __attribute__((unused)) stamp,
161  struct Int32Vect3 *accel)
162 {
163 #if USE_AUTO_AHRS_FREQ || !defined(AHRS_CORRECT_FREQUENCY)
164  PRINT_CONFIG_MSG("Calculating dt for AHRS int_cmpl_quat accel update.")
165  static uint32_t last_stamp = 0;
166  if (last_stamp > 0 && ahrs_icq.is_aligned) {
167  float dt = (float)(stamp - last_stamp) * 1e-6;
168  ahrs_icq_update_accel(accel, dt);
170  }
171  last_stamp = stamp;
172 #else
173  PRINT_CONFIG_MSG("Using fixed AHRS_CORRECT_FREQUENCY for AHRS int_cmpl_quat accel update.")
174  PRINT_CONFIG_VAR(AHRS_CORRECT_FREQUENCY)
175  if (ahrs_icq.is_aligned) {
176  const float dt = 1. / (AHRS_CORRECT_FREQUENCY);
177  ahrs_icq_update_accel(accel, dt);
179  }
180 #endif
181 }
182 
183 static void mag_cb(uint8_t __attribute__((unused)) sender_id,
184  uint32_t __attribute__((unused)) stamp,
185  struct Int32Vect3 *mag)
186 {
187 #if USE_AUTO_AHRS_FREQ || !defined(AHRS_MAG_CORRECT_FREQUENCY)
188  PRINT_CONFIG_MSG("Calculating dt for AHRS int_cmpl_quat mag update.")
189  static uint32_t last_stamp = 0;
190  if (last_stamp > 0 && ahrs_icq.is_aligned) {
191  float dt = (float)(stamp - last_stamp) * 1e-6;
192  ahrs_icq_update_mag(mag, dt);
194  }
195  last_stamp = stamp;
196 #else
197  PRINT_CONFIG_MSG("Using fixed AHRS_MAG_CORRECT_FREQUENCY for AHRS int_cmpl_quat mag update.")
198  PRINT_CONFIG_VAR(AHRS_MAG_CORRECT_FREQUENCY)
199  if (ahrs_icq.is_aligned) {
200  const float dt = 1. / (AHRS_MAG_CORRECT_FREQUENCY);
201  ahrs_icq_update_mag(mag, dt);
203  }
204 #endif
205 }
206 
207 static void aligner_cb(uint8_t __attribute__((unused)) sender_id,
208  uint32_t stamp __attribute__((unused)),
209  struct Int32Rates *lp_gyro, struct Int32Vect3 *lp_accel,
210  struct Int32Vect3 *lp_mag)
211 {
212  if (!ahrs_icq.is_aligned) {
213  if (ahrs_icq_align(lp_gyro, lp_accel, lp_mag)) {
215  }
216  }
217 }
218 
219 static void body_to_imu_cb(uint8_t sender_id __attribute__((unused)),
220  struct FloatQuat *q_b2i_f)
221 {
223 }
224 
225 static void geo_mag_cb(uint8_t sender_id __attribute__((unused)), struct FloatVect3 *h)
226 {
228  MAG_BFP_OF_REAL(h->z));
229 }
230 
231 static void gps_cb(uint8_t sender_id __attribute__((unused)),
232  uint32_t stamp __attribute__((unused)),
233  struct GpsState *gps_s)
234 {
235  ahrs_icq_update_gps(gps_s);
236 }
237 
238 static bool_t ahrs_icq_enable_output(bool_t enable)
239 {
240  ahrs_icq_output_enabled = enable;
242 }
243 
245 static void set_body_state_from_quat(void)
246 {
247  if (ahrs_icq_output_enabled) {
248  /* Compute LTP to BODY quaternion */
249  struct Int32Quat ltp_to_body_quat;
250  struct Int32Quat *body_to_imu_quat = orientationGetQuat_i(&ahrs_icq.body_to_imu);
251  int32_quat_comp_inv(&ltp_to_body_quat, &ahrs_icq.ltp_to_imu_quat, body_to_imu_quat);
252  /* Set state */
253  stateSetNedToBodyQuat_i(&ltp_to_body_quat);
254 
255  /* compute body rates */
256  struct Int32Rates body_rate;
257  struct Int32RMat *body_to_imu_rmat = orientationGetRMat_i(&ahrs_icq.body_to_imu);
258  int32_rmat_transp_ratemult(&body_rate, body_to_imu_rmat, &ahrs_icq.imu_rate);
259  /* Set state */
260  stateSetBodyRates_i(&body_rate);
261  }
262 }
263 
265 {
266  ahrs_icq_output_enabled = AHRS_ICQ_OUTPUT_ENABLED;
267  ahrs_icq_init();
269 
270  /*
271  * Subscribe to scaled IMU measurements and attach callbacks
272  */
273  AbiBindMsgIMU_GYRO_INT32(AHRS_ICQ_IMU_ID, &gyro_ev, gyro_cb);
274  AbiBindMsgIMU_ACCEL_INT32(AHRS_ICQ_IMU_ID, &accel_ev, accel_cb);
275  AbiBindMsgIMU_MAG_INT32(AHRS_ICQ_MAG_ID, &mag_ev, mag_cb);
276  AbiBindMsgIMU_LOWPASSED(ABI_BROADCAST, &aligner_ev, aligner_cb);
277  AbiBindMsgBODY_TO_IMU_QUAT(ABI_BROADCAST, &body_to_imu_ev, body_to_imu_cb);
278  AbiBindMsgGEO_MAG(ABI_BROADCAST, &geo_mag_ev, geo_mag_cb);
279  AbiBindMsgGPS(ABI_BROADCAST, &gps_ev, gps_cb);
280 
281 #if PERIODIC_TELEMETRY
282  register_periodic_telemetry(DefaultPeriodic, PPRZ_MSG_ID_AHRS_QUAT_INT, send_quat);
283  register_periodic_telemetry(DefaultPeriodic, PPRZ_MSG_ID_AHRS_EULER_INT, send_euler);
284  register_periodic_telemetry(DefaultPeriodic, PPRZ_MSG_ID_AHRS_GYRO_BIAS_INT, send_bias);
285  register_periodic_telemetry(DefaultPeriodic, PPRZ_MSG_ID_GEO_MAG, send_geo_mag);
286  register_periodic_telemetry(DefaultPeriodic, PPRZ_MSG_ID_STATE_FILTER_STATUS, send_filter_status);
287 #endif
288 }
int32_t psi
in rad with INT32_ANGLE_FRAC
static abi_event geo_mag_ev
Event structure to store callbacks in a linked list.
Definition: abi_common.h:65
unsigned short uint16_t
Definition: types.h:16
struct AhrsIntCmplQuat ahrs_icq
void ahrs_icq_update_mag(struct Int32Vect3 *mag, float dt)
angular rates
static abi_event body_to_imu_ev
static struct Int32RMat * orientationGetRMat_i(struct OrientationReps *orientation)
Get vehicle body attitude rotation matrix (int).
struct OrientationReps body_to_imu
Dispatcher to register actual AHRS implementations.
struct Int32Vect3 mag_h
Generic transmission transport header.
Definition: transport.h:89
#define AHRS_PROPAGATE_FREQUENCY
Definition: hf_float.c:51
Periodic telemetry system header (includes downlink utility and generated code).
int32_t theta
in rad with INT32_ANGLE_FRAC
static void body_to_imu_cb(uint8_t sender_id, struct FloatQuat *q_b2i_f)
void ahrs_icq_update_gps(struct GpsState *gps_s)
#define VECT3_ASSIGN(_a, _x, _y, _z)
Definition: pprz_algebra.h:124
Main include for ABI (AirBorneInterface).
static void mag_cb(uint8_t sender_id, uint32_t stamp, struct Int32Vect3 *mag)
#define MAG_FLOAT_OF_BFP(_ai)
static abi_event mag_ev
static void accel_cb(uint8_t sender_id, uint32_t stamp, struct Int32Vect3 *accel)
static void gyro_cb(uint8_t sender_id, uint32_t stamp, struct Int32Rates *gyro)
static abi_event accel_ev
static uint8_t ahrs_icq_id
static abi_event gyro_ev
static bool_t ahrs_icq_enable_output(bool_t enable)
static struct Int32Quat * stateGetNedToBodyQuat_i(void)
Get vehicle body attitude quaternion (int).
Definition: state.h:1084
static void aligner_cb(uint8_t sender_id, uint32_t stamp, struct Int32Rates *lp_gyro, struct Int32Vect3 *lp_accel, struct Int32Vect3 *lp_mag)
void ahrs_icq_update_accel(struct Int32Vect3 *accel, float dt)
void ahrs_icq_propagate(struct Int32Rates *gyro, float dt)
void ahrs_icq_set_body_to_imu_quat(struct FloatQuat *q_b2i)
static uint32_t get_sys_time_usec(void)
Get the time in microseconds since startup.
Definition: sys_time_arch.h:39
Roation quaternion.
Paparazzi specific wrapper to run floating point complementary filter.
static bool_t ahrs_icq_output_enabled
if TRUE with push the estimation results to the state interface
int32_t r
in rad/s with INT32_RATE_FRAC
void ahrs_icq_register(void)
static void gps_cb(uint8_t sender_id, uint32_t stamp, struct GpsState *gps_s)
struct Int32Rates gyro_bias
#define AHRS_ICQ_IMU_ID
ABI binding for IMU data.
void ahrs_register_impl(AhrsEnableOutput enable)
Register an AHRS implementation.
Definition: ahrs.c:57
struct Int32Quat ltp_to_imu_quat
#define MAG_BFP_OF_REAL(_af)
Architecture independent timing functions.
data structure for GPS information
Definition: gps.h:67
euler angles
uint16_t val[TCOUPLE_NB]
static abi_event gps_ev
unsigned long uint32_t
Definition: types.h:18
PRINT_CONFIG_MSG("USE_INS_NAV_INIT defaulting to TRUE")
bool_t ahrs_icq_align(struct Int32Rates *lp_gyro, struct Int32Vect3 *lp_accel, struct Int32Vect3 *lp_mag)
#define DefaultPeriodic
Set default periodic telemetry.
Definition: telemetry.h:66
void ahrs_icq_init(void)
void int32_quat_comp_inv(struct Int32Quat *a2b, struct Int32Quat *a2c, struct Int32Quat *b2c)
Composition (multiplication) of two quaternions.
static void stateSetBodyRates_i(struct Int32Rates *body_rate)
Set vehicle body angular rate (int).
Definition: state.h:1144
static uint32_t ahrs_icq_last_stamp
static void stateSetNedToBodyQuat_i(struct Int32Quat *ned_to_body_quat)
Set vehicle body attitude from quaternion (int).
Definition: state.h:1046
#define AHRS_ICQ_OUTPUT_ENABLED
static const struct usb_device_descriptor dev
Definition: usb_ser_hw.c:69
void int32_eulers_of_quat(struct Int32Eulers *e, struct Int32Quat *q)
static abi_event aligner_ev
int32_t phi
in rad with INT32_ANGLE_FRAC
unsigned char uint8_t
Definition: types.h:14
API to get/set the generic vehicle states.
static struct Int32Quat * orientationGetQuat_i(struct OrientationReps *orientation)
Get vehicle body attitude quaternion (int).
static void geo_mag_cb(uint8_t sender_id, struct FloatVect3 *h)
static void set_body_state_from_quat(void)
Rotate angles and rates from imu to body frame and set state.
struct Int32Rates imu_rate
int32_t p
in rad/s with INT32_RATE_FRAC
#define ABI_BROADCAST
Broadcast address.
Definition: abi_common.h:56
#define AHRS_COMP_ID_ICQ
Definition: ahrs.h:35
rotation matrix
static struct Int32Eulers * stateGetNedToBodyEulers_i(void)
Get vehicle body attitude euler angles (int).
Definition: state.h:1096
void int32_rmat_transp_ratemult(struct Int32Rates *rb, struct Int32RMat *m_b2a, struct Int32Rates *ra)
rotate anglular rates by transposed rotation matrix.
#define AHRS_ICQ_MAG_ID
ABI binding for magnetometer data.
int32_t q
in rad/s with INT32_RATE_FRAC
Rotation quaternion.
int8_t register_periodic_telemetry(struct periodic_telemetry *_pt, uint8_t _id, telemetry_cb _cb)
Register a telemetry callback function.
Definition: telemetry.c:46
enum AhrsICQStatus status
status of the AHRS, AHRS_ICQ_UNINIT or AHRS_ICQ_RUNNING