A guidance mode based on Incremental Nonlinear Dynamic Inversion. More...

#include "generated/airframe.h"
#include "firmwares/rotorcraft/guidance/guidance_indi.h"
#include "modules/radio_control/radio_control.h"
#include "firmwares/rotorcraft/stabilization.h"
#include "firmwares/rotorcraft/guidance/guidance_h.h"
#include "firmwares/rotorcraft/guidance/guidance_v.h"
#include "firmwares/rotorcraft/autopilot_rc_helpers.h"
#include "mcu_periph/sys_time.h"
#include "state.h"
#include "autopilot.h"
#include "filters/low_pass_filter.h"
#include "modules/core/abi.h"
#include "modules/datalink/telemetry.h"

Include dependency graph for guidance_indi.c:

Go to the source code of this file.

Macros
#define	GUIDANCE_INDI_ACCEL_SP_ID ABI_BROADCAST

#define	GUIDANCE_INDI_FILTER_CUTOFF 3.0

Functions
static void	accel_sp_cb (uint8_t sender_id, uint8_t flag, struct FloatVect3 *accel_sp)
	ABI callback that obtains the acceleration setpoint from telemetry flag: 0 -> 2D, 1 -> 3D. More...

static void	guidance_indi_propagate_filters (struct FloatEulers *eulers)
	Low pass the accelerometer measurements to remove noise from vibrations. More...

static void	guidance_indi_calcG (struct FloatMat33 *Gmat)

static void	guidance_indi_calcG_yxz (struct FloatMat33 Gmat, struct FloatEulers euler_yxz)

static void	send_indi_guidance (struct transport_tx trans, struct link_device dev)

void	guidance_indi_init (void)
	Init function. More...

void	guidance_indi_enter (void)
	Call upon entering indi guidance. More...

struct StabilizationSetpoint	guidance_indi_run (struct FloatVect3 *accel_sp, float heading_sp)

struct StabilizationSetpoint	guidance_indi_run_mode (bool in_flight UNUSED, struct HorizontalGuidance gh, struct VerticalGuidance gv, enum GuidanceIndi_HMode h_mode, enum GuidanceIndi_VMode v_mode)

Variables
float	guidance_indi_pos_gain = 0.5

float	guidance_indi_speed_gain = 1.8

abi_event	accel_sp_ev

struct FloatVect3	indi_accel_sp = {0.0f, 0.0f, 0.0f}

bool	indi_accel_sp_set_2d = false

bool	indi_accel_sp_set_3d = false

struct FloatVect3	speed_sp = {0.0f, 0.0f, 0.0f}

struct FloatVect3	sp_accel = {0.0f, 0.0f, 0.0f}

float	thrust_dyn = 0.f

float	thrust_act = 0.f

Butterworth2LowPass	filt_accel_ned [3]

Butterworth2LowPass	roll_filt

Butterworth2LowPass	pitch_filt

Butterworth2LowPass	thrust_filt

struct FloatMat33	Ga

struct FloatMat33	Ga_inv

struct FloatVect3	control_increment

float	filter_cutoff = GUIDANCE_INDI_FILTER_CUTOFF

float	guidance_indi_max_bank = GUIDANCE_H_MAX_BANK

float	time_of_accel_sp_2d = 0.0

float	time_of_accel_sp_3d = 0.0

struct FloatEulers	guidance_euler_cmd

struct ThrustSetpoint	thrust_sp

float	thrust_in

Detailed Description

A guidance mode based on Incremental Nonlinear Dynamic Inversion.

Based on the papers: Cascaded Incremental Nonlinear Dynamic Inversion Control for MAV Disturbance Rejection https://www.researchgate.net/publication/312907985_Cascaded_Incremental_Nonlinear_Dynamic_Inversion_Control_for_MAV_Disturbance_Rejection

Gust Disturbance Alleviation with Incremental Nonlinear Dynamic Inversion https://www.researchgate.net/publication/309212603_Gust_Disturbance_Alleviation_with_Incremental_Nonlinear_Dynamic_Inversion

Definition in file guidance_indi.c.

Macro Definition Documentation

◆ GUIDANCE_INDI_ACCEL_SP_ID

#define GUIDANCE_INDI_ACCEL_SP_ID ABI_BROADCAST

Definition at line 64 of file guidance_indi.c.

◆ GUIDANCE_INDI_FILTER_CUTOFF

#define GUIDANCE_INDI_FILTER_CUTOFF 3.0

Definition at line 101 of file guidance_indi.c.

Function Documentation

◆ accel_sp_cb()

static void accel_sp_cb	(	uint8_t	sender_id,
		uint8_t	flag,
		struct FloatVect3 *	accel_sp
	)

static

ABI callback that obtains the acceleration setpoint from telemetry flag: 0 -> 2D, 1 -> 3D.

Definition at line 460 of file guidance_indi.c.

References get_sys_time_float(), indi_accel_sp, indi_accel_sp_set_2d, indi_accel_sp_set_3d, time_of_accel_sp_2d, time_of_accel_sp_3d, FloatVect3::x, FloatVect3::y, and FloatVect3::z.

Referenced by guidance_indi_init().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ guidance_indi_calcG()

UNUSED void guidance_indi_calcG ( struct FloatMat33 * Gmat )

static

Parameters

Gmat	array to write the matrix to [3x3]

Calculate the matrix of partial derivatives of the roll, pitch and thrust. w.r.t. the NED accelerations for ZYX eulers ddx = G*[dtheta,dphi,dT]

Definition at line 431 of file guidance_indi.c.

References FloatEulers::phi, FloatEulers::psi, RMAT_ELMT, stateGetNedToBodyEulers_f(), mesonh.mesonh_atmosphere::T, and FloatEulers::theta.

Here is the call graph for this function:

◆ guidance_indi_calcG_yxz()

void guidance_indi_calcG_yxz	(	struct FloatMat33 *	Gmat,
		struct FloatEulers *	euler_yxz
	)

static

Parameters

Gmat	array to write the matrix to [3x3]

Calculate the matrix of partial derivatives of the pitch, roll and thrust. w.r.t. the NED accelerations for YXZ eulers ddx = G*[dtheta,dphi,dT]

Definition at line 403 of file guidance_indi.c.

References FloatEulers::phi, RMAT_ELMT, mesonh.mesonh_atmosphere::T, and FloatEulers::theta.

◆ guidance_indi_enter()

void guidance_indi_enter ( void )

Call upon entering indi guidance.

Definition at line 168 of file guidance_indi.c.

References Stabilization::cmd, filt_accel_ned, filter_cutoff, RotorcraftNavigation::heading, init_butterworth_2_low_pass(), nav, pitch_filt, FloatEulers::psi, roll_filt, stabilization, stateGetNedToBodyEulers_f(), thrust_act, thrust_dyn, thrust_filt, and thrust_in.

Here is the call graph for this function:

◆ guidance_indi_init()

void guidance_indi_init ( void )

Init function.

Definition at line 153 of file guidance_indi.c.

References accel_sp_cb(), accel_sp_ev, DefaultPeriodic, FLOAT_EULERS_ZERO, guidance_euler_cmd, GUIDANCE_INDI_ACCEL_SP_ID, register_periodic_telemetry(), send_indi_guidance(), thrust_sp, and THRUST_SP_SET_ZERO.

Here is the call graph for this function:

◆ guidance_indi_propagate_filters()

void guidance_indi_propagate_filters ( struct FloatEulers * eulers )

static

Low pass the accelerometer measurements to remove noise from vibrations.

The roll and pitch also need to be filtered to synchronize them with the acceleration

Definition at line 385 of file guidance_indi.c.

References filt_accel_ned, FloatEulers::phi, pitch_filt, roll_filt, stateGetAccelNed_f(), FloatEulers::theta, update_butterworth_2_low_pass(), NedCoor_f::x, NedCoor_f::y, and NedCoor_f::z.

Here is the call graph for this function:

◆ guidance_indi_run()

struct StabilizationSetpoint guidance_indi_run	(	struct FloatVect3 *	accel_sp,
		float	heading_sp
	)

Parameters

accel_sp	accel setpoint in NED frame [m/s^2]
heading_sp	the desired heading [rad]

Returns: stabilization setpoint structure

main indi guidance function

Definition at line 168 of file guidance_indi.c.

◆ guidance_indi_run_mode()

struct StabilizationSetpoint guidance_indi_run_mode	(	bool in_flight	UNUSED,
		struct HorizontalGuidance *	gh,
		struct VerticalGuidance *	gv,
		enum GuidanceIndi_HMode	h_mode,
		enum GuidanceIndi_VMode	v_mode
	)