2011-03-19 07:20:11 -03:00
/// -*- tab-width: 4; Mode: C++; c-basic-offset: 4; indent-tabs-mode: nil -*-
2011-12-03 21:54:38 -04:00
2012-10-01 02:02:49 -03:00
static void
2011-11-07 18:24:32 -04:00
get_stabilize_roll(int32_t target_angle)
2010-12-19 12:40:33 -04:00
{
2012-08-21 23:19:50 -03:00
// angle error
2013-03-28 23:14:31 -03:00
target_angle = wrap_180_cd(target_angle - ahrs.roll_sensor);
2011-01-25 01:53:36 -04:00
2012-08-21 23:19:50 -03:00
// limit the error we're feeding to the PID
target_angle = constrain(target_angle, -4500, 4500);
2011-12-11 03:25:52 -04:00
2012-11-26 22:02:41 -04:00
// convert to desired Rate:
2012-08-21 23:19:50 -03:00
int32_t target_rate = g.pi_stabilize_roll.get_p(target_angle);
2011-12-11 03:25:52 -04:00
2012-08-21 23:19:50 -03:00
int16_t i_stab;
if(labs(ahrs.roll_sensor) < 500) {
target_angle = constrain(target_angle, -500, 500);
i_stab = g.pi_stabilize_roll.get_i(target_angle, G_Dt);
}else{
i_stab = g.pi_stabilize_roll.get_integrator();
}
2012-06-05 20:41:31 -03:00
2012-10-01 02:02:49 -03:00
// set targets for rate controller
2012-10-14 05:47:46 -03:00
set_roll_rate_target(target_rate+i_stab, EARTH_FRAME);
2011-04-25 02:12:59 -03:00
}
2012-10-01 02:02:49 -03:00
static void
2011-11-07 18:24:32 -04:00
get_stabilize_pitch(int32_t target_angle)
2011-01-25 01:53:36 -04:00
{
2012-08-21 23:19:50 -03:00
// angle error
2013-03-28 23:14:31 -03:00
target_angle = wrap_180_cd(target_angle - ahrs.pitch_sensor);
2011-02-17 05:36:33 -04:00
2012-08-21 23:19:50 -03:00
// limit the error we're feeding to the PID
target_angle = constrain(target_angle, -4500, 4500);
2011-12-11 03:25:52 -04:00
2012-08-21 23:19:50 -03:00
// convert to desired Rate:
int32_t target_rate = g.pi_stabilize_pitch.get_p(target_angle);
2011-12-03 21:54:38 -04:00
2012-08-21 23:19:50 -03:00
int16_t i_stab;
if(labs(ahrs.pitch_sensor) < 500) {
target_angle = constrain(target_angle, -500, 500);
i_stab = g.pi_stabilize_pitch.get_i(target_angle, G_Dt);
}else{
i_stab = g.pi_stabilize_pitch.get_integrator();
}
2012-10-01 02:02:49 -03:00
// set targets for rate controller
2012-10-14 05:47:46 -03:00
set_pitch_rate_target(target_rate + i_stab, EARTH_FRAME);
2012-01-29 02:00:05 -04:00
}
2010-12-19 12:40:33 -04:00
2012-10-01 02:02:49 -03:00
static void
2012-01-29 02:00:05 -04:00
get_stabilize_yaw(int32_t target_angle)
{
2012-08-21 23:19:50 -03:00
int32_t target_rate,i_term;
int32_t angle_error;
2012-10-01 02:02:49 -03:00
int32_t output = 0;
2012-03-25 03:55:49 -03:00
2012-08-21 23:19:50 -03:00
// angle error
2013-03-28 23:14:31 -03:00
angle_error = wrap_180_cd(target_angle - ahrs.yaw_sensor);
2011-07-10 21:47:08 -03:00
2012-08-21 23:19:50 -03:00
// limit the error we're feeding to the PID
2012-11-26 22:02:41 -04:00
2012-08-21 23:19:50 -03:00
angle_error = constrain(angle_error, -4500, 4500);
2011-12-29 14:23:18 -04:00
2012-08-21 23:19:50 -03:00
// convert angle error to desired Rate:
target_rate = g.pi_stabilize_yaw.get_p(angle_error);
i_term = g.pi_stabilize_yaw.get_i(angle_error, G_Dt);
2011-12-29 14:23:18 -04:00
2012-08-21 23:19:50 -03:00
// do not use rate controllers for helicotpers with external gyros
2012-03-25 03:55:49 -03:00
#if FRAME_CONFIG == HELI_FRAME
2012-10-06 01:46:19 -03:00
if(motors.ext_gyro_enabled) {
g.rc_4.servo_out = constrain((target_rate + i_term), -4500, 4500);
2012-08-21 23:19:50 -03:00
}
2011-12-08 08:30:47 -04:00
#endif
2012-03-25 03:55:49 -03:00
2012-03-27 01:33:21 -03:00
#if LOGGING_ENABLED == ENABLED
2012-08-21 23:19:50 -03:00
// log output if PID logging is on and we are tuning the yaw
2013-01-08 01:45:12 -04:00
if( g.log_bitmask & MASK_LOG_PID && g.radio_tuning == CH6_YAW_KP ) {
pid_log_counter++;
if( pid_log_counter >= 10 ) { // (update rate / desired output rate) = (100hz / 10hz) = 10
pid_log_counter = 0;
2012-08-21 23:19:50 -03:00
Log_Write_PID(CH6_YAW_KP, angle_error, target_rate, i_term, 0, output, tuning_value);
}
}
2012-03-27 01:33:21 -03:00
#endif
2012-03-25 03:55:49 -03:00
2012-10-01 02:02:49 -03:00
// set targets for rate controller
2012-10-14 05:47:46 -03:00
set_yaw_rate_target(target_rate+i_term, EARTH_FRAME);
}
2012-10-01 02:02:49 -03:00
static void
2012-02-09 01:09:24 -04:00
get_acro_roll(int32_t target_rate)
{
2012-08-21 23:19:50 -03:00
target_rate = target_rate * g.acro_p;
2012-10-01 02:02:49 -03:00
// set targets for rate controller
2012-10-14 05:47:46 -03:00
set_roll_rate_target(target_rate, BODY_FRAME);
2012-02-09 01:09:24 -04:00
}
2012-10-01 02:02:49 -03:00
static void
2012-02-09 01:09:24 -04:00
get_acro_pitch(int32_t target_rate)
{
2012-08-21 23:19:50 -03:00
target_rate = target_rate * g.acro_p;
2012-10-01 02:02:49 -03:00
// set targets for rate controller
2012-10-14 05:47:46 -03:00
set_pitch_rate_target(target_rate, BODY_FRAME);
2012-02-09 01:09:24 -04:00
}
2012-10-01 02:02:49 -03:00
static void
2012-02-11 02:22:19 -04:00
get_acro_yaw(int32_t target_rate)
2012-07-14 16:26:13 -03:00
{
2012-10-14 05:47:46 -03:00
target_rate = target_rate * g.acro_p;
2012-10-01 02:02:49 -03:00
// set targets for rate controller
2012-10-14 05:47:46 -03:00
set_yaw_rate_target(target_rate, BODY_FRAME);
2012-07-14 16:26:13 -03:00
}
2012-10-18 11:24:34 -03:00
// Roll with rate input and stabilized in the earth frame
static void
get_roll_rate_stabilized_ef(int32_t stick_angle)
{
int32_t angle_error = 0;
// convert the input to the desired roll rate
2012-10-23 09:30:50 -03:00
int32_t target_rate = stick_angle * g.acro_p - (roll_axis * g.acro_balance_roll)/100;
2012-10-18 11:24:34 -03:00
// convert the input to the desired roll rate
roll_axis += target_rate * G_Dt;
2013-03-28 23:14:31 -03:00
roll_axis = wrap_180_cd(roll_axis);
2012-10-18 11:24:34 -03:00
// ensure that we don't reach gimbal lock
2013-03-31 23:41:52 -03:00
if (labs(roll_axis) > 4500 && g.acro_trainer_enabled) {
2012-10-18 11:24:34 -03:00
roll_axis = constrain(roll_axis, -4500, 4500);
2013-03-28 23:14:31 -03:00
angle_error = wrap_180_cd(roll_axis - ahrs.roll_sensor);
2012-10-18 11:24:34 -03:00
} else {
// angle error with maximum of +- max_angle_overshoot
2013-03-28 23:14:31 -03:00
angle_error = wrap_180_cd(roll_axis - ahrs.roll_sensor);
2012-10-18 11:24:34 -03:00
angle_error = constrain(angle_error, -MAX_ROLL_OVERSHOOT, MAX_ROLL_OVERSHOOT);
}
2013-03-16 05:14:21 -03:00
if (motors.armed() == false || ((g.rc_3.control_in == 0) && !ap.failsafe_radio)) {
2012-10-18 11:24:34 -03:00
angle_error = 0;
}
// update roll_axis to be within max_angle_overshoot of our current heading
2013-03-28 23:14:31 -03:00
roll_axis = wrap_180_cd(angle_error + ahrs.roll_sensor);
2012-10-18 11:24:34 -03:00
// set earth frame targets for rate controller
// set earth frame targets for rate controller
set_roll_rate_target(g.pi_stabilize_roll.get_p(angle_error) + target_rate, EARTH_FRAME);
}
// Pitch with rate input and stabilized in the earth frame
static void
get_pitch_rate_stabilized_ef(int32_t stick_angle)
{
int32_t angle_error = 0;
// convert the input to the desired pitch rate
2012-10-23 09:30:50 -03:00
int32_t target_rate = stick_angle * g.acro_p - (pitch_axis * g.acro_balance_pitch)/100;
2012-10-18 11:24:34 -03:00
// convert the input to the desired pitch rate
pitch_axis += target_rate * G_Dt;
2013-03-28 23:14:31 -03:00
pitch_axis = wrap_180_cd(pitch_axis);
2012-10-18 11:24:34 -03:00
// ensure that we don't reach gimbal lock
2012-12-10 09:27:46 -04:00
if (labs(pitch_axis) > 4500) {
2012-10-18 11:24:34 -03:00
pitch_axis = constrain(pitch_axis, -4500, 4500);
2013-03-28 23:14:31 -03:00
angle_error = wrap_180_cd(pitch_axis - ahrs.pitch_sensor);
2012-10-18 11:24:34 -03:00
} else {
// angle error with maximum of +- max_angle_overshoot
2013-03-28 23:14:31 -03:00
angle_error = wrap_180_cd(pitch_axis - ahrs.pitch_sensor);
2012-10-18 11:24:34 -03:00
angle_error = constrain(angle_error, -MAX_PITCH_OVERSHOOT, MAX_PITCH_OVERSHOOT);
}
2013-03-16 05:14:21 -03:00
if (motors.armed() == false || ((g.rc_3.control_in == 0) && !ap.failsafe_radio)) {
2012-10-18 11:24:34 -03:00
angle_error = 0;
}
// update pitch_axis to be within max_angle_overshoot of our current heading
2013-03-28 23:14:31 -03:00
pitch_axis = wrap_180_cd(angle_error + ahrs.pitch_sensor);
2012-10-18 11:24:34 -03:00
// set earth frame targets for rate controller
set_pitch_rate_target(g.pi_stabilize_pitch.get_p(angle_error) + target_rate, EARTH_FRAME);
}
2012-10-14 09:30:16 -03:00
// Yaw with rate input and stabilized in the earth frame
static void
get_yaw_rate_stabilized_ef(int32_t stick_angle)
{
int32_t angle_error = 0;
// convert the input to the desired yaw rate
int32_t target_rate = stick_angle * g.acro_p;
2012-10-18 11:24:34 -03:00
// convert the input to the desired yaw rate
2012-10-14 09:30:16 -03:00
nav_yaw += target_rate * G_Dt;
2013-03-28 23:14:31 -03:00
nav_yaw = wrap_360_cd(nav_yaw);
2012-10-14 09:30:16 -03:00
2012-10-16 22:18:24 -03:00
// calculate difference between desired heading and current heading
2013-03-28 23:14:31 -03:00
angle_error = wrap_180_cd(nav_yaw - ahrs.yaw_sensor);
2012-10-14 09:30:16 -03:00
2012-10-16 22:18:24 -03:00
// limit the maximum overshoot
2012-10-18 11:24:34 -03:00
angle_error = constrain(angle_error, -MAX_YAW_OVERSHOOT, MAX_YAW_OVERSHOOT);
2012-10-14 09:30:16 -03:00
2013-03-16 05:14:21 -03:00
if (motors.armed() == false || ((g.rc_3.control_in == 0) && !ap.failsafe_radio)) {
2012-10-14 09:30:16 -03:00
angle_error = 0;
}
2012-10-18 11:24:34 -03:00
// update nav_yaw to be within max_angle_overshoot of our current heading
2013-03-28 23:14:31 -03:00
nav_yaw = wrap_360_cd(angle_error + ahrs.yaw_sensor);
2012-10-14 09:30:16 -03:00
// set earth frame targets for rate controller
set_yaw_rate_target(g.pi_stabilize_yaw.get_p(angle_error)+target_rate, EARTH_FRAME);
}
2012-10-14 05:47:46 -03:00
// set_roll_rate_target - to be called by upper controllers to set roll rate targets in the earth frame
void set_roll_rate_target( int32_t desired_rate, uint8_t earth_or_body_frame ) {
rate_targets_frame = earth_or_body_frame;
if( earth_or_body_frame == BODY_FRAME ) {
roll_rate_target_bf = desired_rate;
}else{
roll_rate_target_ef = desired_rate;
}
}
// set_pitch_rate_target - to be called by upper controllers to set pitch rate targets in the earth frame
void set_pitch_rate_target( int32_t desired_rate, uint8_t earth_or_body_frame ) {
rate_targets_frame = earth_or_body_frame;
if( earth_or_body_frame == BODY_FRAME ) {
pitch_rate_target_bf = desired_rate;
}else{
pitch_rate_target_ef = desired_rate;
}
}
// set_yaw_rate_target - to be called by upper controllers to set yaw rate targets in the earth frame
void set_yaw_rate_target( int32_t desired_rate, uint8_t earth_or_body_frame ) {
rate_targets_frame = earth_or_body_frame;
if( earth_or_body_frame == BODY_FRAME ) {
yaw_rate_target_bf = desired_rate;
}else{
yaw_rate_target_ef = desired_rate;
}
}
// update_rate_contoller_targets - converts earth frame rates to body frame rates for rate controllers
void
update_rate_contoller_targets()
{
if( rate_targets_frame == EARTH_FRAME ) {
// convert earth frame rates to body frame rates
2012-10-28 16:12:04 -03:00
roll_rate_target_bf = roll_rate_target_ef - sin_pitch * yaw_rate_target_ef;
pitch_rate_target_bf = cos_roll_x * pitch_rate_target_ef + sin_roll * cos_pitch_x * yaw_rate_target_ef;
yaw_rate_target_bf = cos_pitch_x * cos_roll_x * yaw_rate_target_ef - sin_roll * pitch_rate_target_ef;
2012-10-14 05:47:46 -03:00
}
}
2012-10-14 09:30:16 -03:00
2012-10-01 02:02:49 -03:00
// run roll, pitch and yaw rate controllers and send output to motors
// targets for these controllers comes from stabilize controllers
void
run_rate_controllers()
{
2012-10-06 01:46:19 -03:00
#if FRAME_CONFIG == HELI_FRAME // helicopters only use rate controllers for yaw and only when not using an external gyro
if(!motors.ext_gyro_enabled) {
2012-11-26 22:02:41 -04:00
g.rc_1.servo_out = get_heli_rate_roll(roll_rate_target_bf);
g.rc_2.servo_out = get_heli_rate_pitch(pitch_rate_target_bf);
g.rc_4.servo_out = get_heli_rate_yaw(yaw_rate_target_bf);
2012-10-06 01:46:19 -03:00
}
#else
2012-10-01 02:02:49 -03:00
// call rate controllers
2012-10-14 05:47:46 -03:00
g.rc_1.servo_out = get_rate_roll(roll_rate_target_bf);
g.rc_2.servo_out = get_rate_pitch(pitch_rate_target_bf);
g.rc_4.servo_out = get_rate_yaw(yaw_rate_target_bf);
2012-10-06 01:46:19 -03:00
#endif
2012-11-23 02:57:49 -04:00
2012-12-03 10:05:14 -04:00
// run throttle controller if accel based throttle controller is enabled and active (active means it has been given a target)
if( g.throttle_accel_enabled && throttle_accel_controller_active ) {
2012-11-24 00:41:17 -04:00
set_throttle_out(get_throttle_accel(throttle_accel_target_ef), true);
2012-11-23 02:57:49 -04:00
}
2012-10-01 02:02:49 -03:00
}
2012-11-26 22:02:41 -04:00
#if FRAME_CONFIG == HELI_FRAME
2012-12-08 15:49:38 -04:00
// init_rate_controllers - set-up filters for rate controller inputs
void init_rate_controllers()
{
// initalise low pass filters on rate controller inputs
// 1st parameter is time_step, 2nd parameter is time_constant
2013-01-23 01:50:44 -04:00
rate_roll_filter.set_cutoff_frequency(0.01f, 2.0f);
rate_pitch_filter.set_cutoff_frequency(0.01f, 2.0f);
// rate_yaw_filter.set_cutoff_frequency(0.01f, 2.0f);
2012-12-08 15:49:38 -04:00
// other option for initialisation is rate_roll_filter.set_cutoff_frequency(<time_step>,<cutoff_freq>);
}
2012-11-26 22:02:41 -04:00
static int16_t
get_heli_rate_roll(int32_t target_rate)
{
2012-12-04 15:42:04 -04:00
int32_t p,i,d,ff; // used to capture pid values for logging
2012-11-26 22:02:41 -04:00
int32_t current_rate; // this iteration's rate
int32_t rate_error; // simply target_rate - current_rate
int32_t output; // output from pid controller
// get current rate
current_rate = (omega.x * DEGX100);
// filter input
2012-11-26 22:50:39 -04:00
current_rate = rate_roll_filter.apply(current_rate);
2012-11-26 22:02:41 -04:00
// call pid controller
2012-12-04 15:42:04 -04:00
rate_error = target_rate - current_rate;
2012-11-26 22:02:41 -04:00
p = g.pid_rate_roll.get_p(rate_error);
if (motors.flybar_mode == 1) { // Mechanical Flybars get regular integral for rate auto trim
if (target_rate > -50 && target_rate < 50){ // Frozen at high rates
i = g.pid_rate_roll.get_i(rate_error, G_Dt);
} else {
i = g.pid_rate_roll.get_integrator();
}
} else {
i = g.pid_rate_roll.get_leaky_i(rate_error, G_Dt, RATE_INTEGRATOR_LEAK_RATE); // Flybarless Helis get huge I-terms. I-term controls much of the rate
}
d = g.pid_rate_roll.get_d(rate_error, G_Dt);
2012-12-04 15:42:04 -04:00
ff = g.heli_roll_ff * target_rate;
2012-11-26 22:02:41 -04:00
2012-12-04 15:42:04 -04:00
output = p + i + d + ff;
2012-11-26 22:02:41 -04:00
// constrain output
output = constrain(output, -4500, 4500);
#if LOGGING_ENABLED == ENABLED
// log output if PID logging is on and we are tuning the rate P, I or D gains
if( g.log_bitmask & MASK_LOG_PID && (g.radio_tuning == CH6_RATE_KP || g.radio_tuning == CH6_RATE_KI || g.radio_tuning == CH6_RATE_KD) ) {
2013-01-08 01:45:12 -04:00
pid_log_counter++;
if( pid_log_counter >= 10 ) { // (update rate / desired output rate) = (100hz / 10hz) = 10
pid_log_counter = 0;
2012-11-26 22:02:41 -04:00
Log_Write_PID(CH6_RATE_KP, rate_error, p, i, d, output, tuning_value);
}
}
#endif
// output control
return output;
}
static int16_t
get_heli_rate_pitch(int32_t target_rate)
{
2012-12-04 15:42:04 -04:00
int32_t p,i,d,ff; // used to capture pid values for logging
2012-11-26 22:02:41 -04:00
int32_t current_rate; // this iteration's rate
int32_t rate_error; // simply target_rate - current_rate
int32_t output; // output from pid controller
// get current rate
current_rate = (omega.y * DEGX100);
// filter input
2012-11-26 22:50:39 -04:00
current_rate = rate_pitch_filter.apply(current_rate);
2012-11-26 22:02:41 -04:00
// call pid controller
2012-12-04 15:42:04 -04:00
rate_error = target_rate - current_rate;
2012-11-26 22:02:41 -04:00
p = g.pid_rate_pitch.get_p(rate_error); // Helicopters get huge feed-forward
if (motors.flybar_mode == 1) { // Mechanical Flybars get regular integral for rate auto trim
if (target_rate > -50 && target_rate < 50){ // Frozen at high rates
i = g.pid_rate_pitch.get_i(rate_error, G_Dt);
} else {
i = g.pid_rate_pitch.get_integrator();
}
} else {
i = g.pid_rate_pitch.get_leaky_i(rate_error, G_Dt, RATE_INTEGRATOR_LEAK_RATE); // Flybarless Helis get huge I-terms. I-term controls much of the rate
}
d = g.pid_rate_pitch.get_d(rate_error, G_Dt);
2012-12-04 15:42:04 -04:00
ff = g.heli_pitch_ff*target_rate;
2012-11-26 22:02:41 -04:00
2012-12-04 15:42:04 -04:00
output = p + i + d + ff;
2012-11-26 22:02:41 -04:00
// constrain output
output = constrain(output, -4500, 4500);
#if LOGGING_ENABLED == ENABLED
// log output if PID logging is on and we are tuning the rate P, I or D gains
if( g.log_bitmask & MASK_LOG_PID && (g.radio_tuning == CH6_RATE_KP || g.radio_tuning == CH6_RATE_KI || g.radio_tuning == CH6_RATE_KD) ) {
2013-01-08 01:45:12 -04:00
if( pid_log_counter == 0 ) { // relies on get_heli_rate_roll to update pid_log_counter
2012-11-26 22:02:41 -04:00
Log_Write_PID(CH6_RATE_KP+100, rate_error, p, i, 0, output, tuning_value);
}
}
#endif
// output control
return output;
}
static int16_t
get_heli_rate_yaw(int32_t target_rate)
{
2012-12-04 15:42:04 -04:00
int32_t p,i,d,ff; // used to capture pid values for logging
2012-11-26 22:02:41 -04:00
int32_t current_rate; // this iteration's rate
int32_t rate_error;
int32_t output;
// get current rate
current_rate = (omega.z * DEGX100);
// filter input
// current_rate = rate_yaw_filter.apply(current_rate);
// rate control
rate_error = target_rate - current_rate;
// separately calculate p, i, d values for logging
p = g.pid_rate_yaw.get_p(rate_error);
i = g.pid_rate_yaw.get_i(rate_error, G_Dt);
d = g.pid_rate_yaw.get_d(rate_error, G_Dt);
2012-12-04 15:42:04 -04:00
ff = g.heli_yaw_ff*target_rate;
2012-11-26 22:02:41 -04:00
2012-12-04 15:42:04 -04:00
output = p + i + d + ff;
2012-11-26 22:02:41 -04:00
output = constrain(output, -4500, 4500);
#if LOGGING_ENABLED == ENABLED
// log output if PID loggins is on and we are tuning the yaw
2013-01-08 01:45:12 -04:00
if( g.log_bitmask & MASK_LOG_PID && (g.radio_tuning == CH6_YAW_RATE_KP || g.radio_tuning == CH6_YAW_RATE_KD) ) {
pid_log_counter++;
if( pid_log_counter >= 10 ) { // (update rate / desired output rate) = (100hz / 10hz) = 10
pid_log_counter = 0;
2012-11-26 22:02:41 -04:00
Log_Write_PID(CH6_YAW_RATE_KP, rate_error, p, i, d, output, tuning_value);
}
}
#endif
// output control
return output;
}
#endif // HELI_FRAME
#if FRAME_CONFIG != HELI_FRAME
2012-07-14 16:26:13 -03:00
static int16_t
2012-08-21 23:19:50 -03:00
get_rate_roll(int32_t target_rate)
2012-02-11 02:22:19 -04:00
{
2012-10-28 16:12:04 -03:00
int32_t p,i,d; // used to capture pid values for logging
int32_t current_rate; // this iteration's rate
int32_t rate_error; // simply target_rate - current_rate
int32_t output; // output from pid controller
2012-07-14 16:26:13 -03:00
2012-08-21 23:19:50 -03:00
// get current rate
current_rate = (omega.x * DEGX100);
2012-07-14 16:26:13 -03:00
2012-08-21 23:19:50 -03:00
// call pid controller
2012-10-28 16:12:04 -03:00
rate_error = target_rate - current_rate;
p = g.pid_rate_roll.get_p(rate_error);
2012-10-10 03:54:46 -03:00
// freeze I term if we've breached roll-pitch limits
if( motors.reached_limit(AP_MOTOR_ROLLPITCH_LIMIT) ) {
2012-10-28 16:12:04 -03:00
i = g.pid_rate_roll.get_integrator();
2012-10-10 03:54:46 -03:00
}else{
2012-10-28 16:12:04 -03:00
i = g.pid_rate_roll.get_i(rate_error, G_Dt);
2012-10-10 03:54:46 -03:00
}
2012-10-28 16:12:04 -03:00
d = g.pid_rate_roll.get_d(rate_error, G_Dt);
output = p + i + d;
2012-07-14 16:26:13 -03:00
2012-08-21 23:19:50 -03:00
// constrain output
output = constrain(output, -5000, 5000);
2012-07-14 01:24:22 -03:00
#if LOGGING_ENABLED == ENABLED
2012-08-21 23:19:50 -03:00
// log output if PID logging is on and we are tuning the rate P, I or D gains
if( g.log_bitmask & MASK_LOG_PID && (g.radio_tuning == CH6_RATE_KP || g.radio_tuning == CH6_RATE_KI || g.radio_tuning == CH6_RATE_KD) ) {
2013-01-08 01:45:12 -04:00
pid_log_counter++; // Note: get_rate_pitch pid logging relies on this function to update pid_log_counter so if you change the line above you must change the equivalent line in get_rate_pitch
if( pid_log_counter >= 10 ) { // (update rate / desired output rate) = (100hz / 10hz) = 10
pid_log_counter = 0;
Log_Write_PID(CH6_RATE_KP, rate_error, p, i, d, output, tuning_value);
2012-08-21 23:19:50 -03:00
}
}
2012-07-14 01:24:22 -03:00
#endif
2012-08-21 23:19:50 -03:00
// output control
return output;
2012-02-11 02:22:19 -04:00
}
2012-02-28 08:26:37 -04:00
static int16_t
2012-08-21 23:19:50 -03:00
get_rate_pitch(int32_t target_rate)
2011-01-25 01:53:36 -04:00
{
2012-08-21 23:19:50 -03:00
int32_t p,i,d; // used to capture pid values for logging
int32_t current_rate; // this iteration's rate
int32_t rate_error; // simply target_rate - current_rate
int32_t output; // output from pid controller
2012-04-07 22:19:20 -03:00
2012-08-21 23:19:50 -03:00
// get current rate
current_rate = (omega.y * DEGX100);
2012-04-07 22:19:20 -03:00
2012-08-21 23:19:50 -03:00
// call pid controller
rate_error = target_rate - current_rate;
2012-10-28 16:12:04 -03:00
p = g.pid_rate_pitch.get_p(rate_error);
2012-10-10 03:54:46 -03:00
// freeze I term if we've breached roll-pitch limits
if( motors.reached_limit(AP_MOTOR_ROLLPITCH_LIMIT) ) {
2012-10-28 16:12:04 -03:00
i = g.pid_rate_pitch.get_integrator();
2012-10-10 03:54:46 -03:00
}else{
2012-10-28 16:12:04 -03:00
i = g.pid_rate_pitch.get_i(rate_error, G_Dt);
2012-10-10 03:54:46 -03:00
}
2012-10-28 16:12:04 -03:00
d = g.pid_rate_pitch.get_d(rate_error, G_Dt);
output = p + i + d;
2012-08-21 23:19:50 -03:00
// constrain output
output = constrain(output, -5000, 5000);
2012-04-07 22:19:20 -03:00
#if LOGGING_ENABLED == ENABLED
2012-08-21 23:19:50 -03:00
// log output if PID logging is on and we are tuning the rate P, I or D gains
if( g.log_bitmask & MASK_LOG_PID && (g.radio_tuning == CH6_RATE_KP || g.radio_tuning == CH6_RATE_KI || g.radio_tuning == CH6_RATE_KD) ) {
2013-01-08 01:45:12 -04:00
if( pid_log_counter == 0 ) { // relies on get_rate_roll having updated pid_log_counter
Log_Write_PID(CH6_RATE_KP+100, rate_error, p, i, d, output, tuning_value);
2012-08-21 23:19:50 -03:00
}
}
2012-04-07 22:19:20 -03:00
#endif
2012-08-21 23:19:50 -03:00
// output control
return output;
2012-01-29 02:00:05 -04:00
}
2012-02-28 08:26:37 -04:00
static int16_t
2012-01-29 02:00:05 -04:00
get_rate_yaw(int32_t target_rate)
{
2012-08-21 23:19:50 -03:00
int32_t p,i,d; // used to capture pid values for logging
int32_t rate_error;
int32_t output;
2012-03-25 03:55:49 -03:00
2012-08-21 23:19:50 -03:00
// rate control
rate_error = target_rate - (omega.z * DEGX100);
2012-03-25 03:55:49 -03:00
2012-08-21 23:19:50 -03:00
// separately calculate p, i, d values for logging
p = g.pid_rate_yaw.get_p(rate_error);
2012-10-10 03:54:46 -03:00
// freeze I term if we've breached yaw limits
if( motors.reached_limit(AP_MOTOR_YAW_LIMIT) ) {
i = g.pid_rate_yaw.get_integrator();
}else{
i = g.pid_rate_yaw.get_i(rate_error, G_Dt);
}
2012-08-21 23:19:50 -03:00
d = g.pid_rate_yaw.get_d(rate_error, G_Dt);
2012-03-25 03:55:49 -03:00
2012-08-21 23:19:50 -03:00
output = p+i+d;
output = constrain(output, -4500, 4500);
2012-03-25 03:55:49 -03:00
2012-03-27 01:33:21 -03:00
#if LOGGING_ENABLED == ENABLED
2012-08-21 23:19:50 -03:00
// log output if PID loggins is on and we are tuning the yaw
2013-01-08 01:45:12 -04:00
if( g.log_bitmask & MASK_LOG_PID && g.radio_tuning == CH6_YAW_RATE_KP ) {
pid_log_counter++;
if( pid_log_counter >= 10 ) { // (update rate / desired output rate) = (100hz / 10hz) = 10
pid_log_counter = 0;
2012-08-21 23:19:50 -03:00
Log_Write_PID(CH6_YAW_RATE_KP, rate_error, p, i, d, output, tuning_value);
}
}
2012-03-27 01:33:21 -03:00
#endif
2012-03-25 03:55:49 -03:00
2012-11-26 22:02:41 -04:00
#if FRAME_CONFIG == TRI_FRAME
2012-08-21 23:19:50 -03:00
// constrain output
return output;
2012-11-26 22:02:41 -04:00
#else // !TRI_FRAME
2012-08-21 23:19:50 -03:00
// output control:
int16_t yaw_limit = 2200 + abs(g.rc_4.control_in);
2012-10-01 02:02:49 -03:00
2012-08-21 23:19:50 -03:00
// smoother Yaw control:
return constrain(output, -yaw_limit, yaw_limit);
2012-11-26 22:02:41 -04:00
#endif // TRI_FRAME
2010-12-19 12:40:33 -04:00
}
2012-11-26 22:02:41 -04:00
#endif // !HELI_FRAME
2010-12-19 12:40:33 -04:00
2012-01-25 08:55:14 -04:00
// calculate modified roll/pitch depending upon optical flow calculated position
2012-01-09 00:53:54 -04:00
static int32_t
2012-08-17 23:48:12 -03:00
get_of_roll(int32_t input_roll)
2012-01-09 00:53:54 -04:00
{
2012-12-10 11:09:45 -04:00
#if OPTFLOW == ENABLED
2012-08-21 23:19:50 -03:00
static float tot_x_cm = 0; // total distance from target
2012-02-28 08:26:37 -04:00
static uint32_t last_of_roll_update = 0;
2012-08-21 23:19:50 -03:00
int32_t new_roll = 0;
int32_t p,i,d;
// check if new optflow data available
if( optflow.last_update != last_of_roll_update) {
last_of_roll_update = optflow.last_update;
// add new distance moved
tot_x_cm += optflow.x_cm;
// only stop roll if caller isn't modifying roll
if( input_roll == 0 && current_loc.alt < 1500) {
p = g.pid_optflow_roll.get_p(-tot_x_cm);
2013-01-23 01:50:44 -04:00
i = g.pid_optflow_roll.get_i(-tot_x_cm,1.0f); // we could use the last update time to calculate the time change
d = g.pid_optflow_roll.get_d(-tot_x_cm,1.0f);
2012-08-21 23:19:50 -03:00
new_roll = p+i+d;
}else{
g.pid_optflow_roll.reset_I();
tot_x_cm = 0;
p = 0; // for logging
i = 0;
d = 0;
}
// limit amount of change and maximum angle
of_roll = constrain(new_roll, (of_roll-20), (of_roll+20));
#if LOGGING_ENABLED == ENABLED
// log output if PID logging is on and we are tuning the rate P, I or D gains
if( g.log_bitmask & MASK_LOG_PID && (g.radio_tuning == CH6_OPTFLOW_KP || g.radio_tuning == CH6_OPTFLOW_KI || g.radio_tuning == CH6_OPTFLOW_KD) ) {
2013-01-08 01:45:12 -04:00
pid_log_counter++; // Note: get_of_pitch pid logging relies on this function updating pid_log_counter so if you change the line above you must change the equivalent line in get_of_pitch
if( pid_log_counter >= 5 ) { // (update rate / desired output rate) = (100hz / 10hz) = 10
pid_log_counter = 0;
2012-08-21 23:19:50 -03:00
Log_Write_PID(CH6_OPTFLOW_KP, tot_x_cm, p, i, d, of_roll, tuning_value);
}
}
#endif // LOGGING_ENABLED == ENABLED
}
// limit max angle
2012-01-25 08:55:14 -04:00
of_roll = constrain(of_roll, -1000, 1000);
2012-04-21 08:17:12 -03:00
2012-08-17 23:48:12 -03:00
return input_roll+of_roll;
2012-01-09 00:53:54 -04:00
#else
2012-08-17 23:48:12 -03:00
return input_roll;
2012-11-23 02:57:49 -04:00
#endif
2012-01-09 00:53:54 -04:00
}
static int32_t
2012-08-17 23:48:12 -03:00
get_of_pitch(int32_t input_pitch)
2012-01-09 00:53:54 -04:00
{
2012-12-10 11:09:45 -04:00
#if OPTFLOW == ENABLED
2012-01-25 08:55:14 -04:00
static float tot_y_cm = 0; // total distance from target
2012-02-28 08:26:37 -04:00
static uint32_t last_of_pitch_update = 0;
2012-08-21 23:19:50 -03:00
int32_t new_pitch = 0;
int32_t p,i,d;
// check if new optflow data available
if( optflow.last_update != last_of_pitch_update ) {
last_of_pitch_update = optflow.last_update;
// add new distance moved
tot_y_cm += optflow.y_cm;
// only stop roll if caller isn't modifying pitch
if( input_pitch == 0 && current_loc.alt < 1500 ) {
p = g.pid_optflow_pitch.get_p(tot_y_cm);
2013-01-23 01:50:44 -04:00
i = g.pid_optflow_pitch.get_i(tot_y_cm, 1.0f); // we could use the last update time to calculate the time change
d = g.pid_optflow_pitch.get_d(tot_y_cm, 1.0f);
2012-08-21 23:19:50 -03:00
new_pitch = p + i + d;
}else{
tot_y_cm = 0;
g.pid_optflow_pitch.reset_I();
p = 0; // for logging
i = 0;
d = 0;
}
// limit amount of change
of_pitch = constrain(new_pitch, (of_pitch-20), (of_pitch+20));
#if LOGGING_ENABLED == ENABLED
// log output if PID logging is on and we are tuning the rate P, I or D gains
if( g.log_bitmask & MASK_LOG_PID && (g.radio_tuning == CH6_OPTFLOW_KP || g.radio_tuning == CH6_OPTFLOW_KI || g.radio_tuning == CH6_OPTFLOW_KD) ) {
2013-01-08 01:45:12 -04:00
if( pid_log_counter == 0 ) { // relies on get_of_roll having updated the pid_log_counter
2012-08-21 23:19:50 -03:00
Log_Write_PID(CH6_OPTFLOW_KP+100, tot_y_cm, p, i, d, of_pitch, tuning_value);
}
}
#endif // LOGGING_ENABLED == ENABLED
}
// limit max angle
of_pitch = constrain(of_pitch, -1000, 1000);
2012-04-21 08:17:12 -03:00
2012-08-17 23:48:12 -03:00
return input_pitch+of_pitch;
2012-01-09 00:53:54 -04:00
#else
2012-08-17 23:48:12 -03:00
return input_pitch;
2012-11-23 02:57:49 -04:00
#endif
}
2012-12-08 01:23:32 -04:00
/*************************************************************
* yaw controllers
*************************************************************/
2013-03-22 05:38:07 -03:00
// get_look_at_yaw - updates bearing to look at center of circle or do a panorama
// should be called at 100hz
static void get_circle_yaw()
{
static uint8_t look_at_yaw_counter = 0; // used to reduce update rate to 10hz
// if circle radius is zero do panorama
if( g.circle_radius == 0 ) {
// slew yaw towards circle angle
nav_yaw = get_yaw_slew(nav_yaw, ToDeg(circle_angle)*100, AUTO_YAW_SLEW_RATE);
}else{
look_at_yaw_counter++;
if( look_at_yaw_counter >= 10 ) {
look_at_yaw_counter = 0;
yaw_look_at_WP_bearing = pv_get_bearing_cd(inertial_nav.get_position(), yaw_look_at_WP);
}
// slew yaw
nav_yaw = get_yaw_slew(nav_yaw, yaw_look_at_WP_bearing, AUTO_YAW_SLEW_RATE);
}
// call stabilize yaw controller
get_stabilize_yaw(nav_yaw);
}
// get_look_at_yaw - updates bearing to location held in look_at_yaw_WP and calls stabilize yaw controller
// should be called at 100hz
static void get_look_at_yaw()
{
static uint8_t look_at_yaw_counter = 0; // used to reduce update rate to 10hz
look_at_yaw_counter++;
if( look_at_yaw_counter >= 10 ) {
look_at_yaw_counter = 0;
yaw_look_at_WP_bearing = pv_get_bearing_cd(inertial_nav.get_position(), yaw_look_at_WP);
}
// slew yaw and call stabilize controller
nav_yaw = get_yaw_slew(nav_yaw, yaw_look_at_WP_bearing, AUTO_YAW_SLEW_RATE);
get_stabilize_yaw(nav_yaw);
}
2012-12-08 01:23:32 -04:00
static void get_look_ahead_yaw(int16_t pilot_yaw)
{
// Commanded Yaw to automatically look ahead.
if (g_gps->fix && g_gps->ground_course > YAW_LOOK_AHEAD_MIN_SPEED) {
nav_yaw = get_yaw_slew(nav_yaw, g_gps->ground_course, AUTO_YAW_SLEW_RATE);
2013-03-28 23:14:31 -03:00
get_stabilize_yaw(wrap_360_cd(nav_yaw + pilot_yaw)); // Allow pilot to "skid" around corners up to 45 degrees
2012-12-08 01:23:32 -04:00
}else{
nav_yaw += pilot_yaw * g.acro_p * G_Dt;
2013-03-28 23:14:31 -03:00
nav_yaw = wrap_360_cd(nav_yaw);
2012-12-08 01:23:32 -04:00
get_stabilize_yaw(nav_yaw);
}
}
2012-11-23 02:57:49 -04:00
/*************************************************************
* throttle control
****************************************************************/
// update_throttle_cruise - update throttle cruise if necessary
static void update_throttle_cruise(int16_t throttle)
{
// ensure throttle_avg has been initialised
if( throttle_avg == 0 ) {
throttle_avg = g.throttle_cruise;
}
2012-11-24 09:50:09 -04:00
// calc average throttle if we are in a level hover
if (throttle > g.throttle_min && abs(climb_rate) < 60 && labs(ahrs.roll_sensor) < 500 && labs(ahrs.pitch_sensor) < 500) {
2013-01-10 14:42:24 -04:00
throttle_avg = throttle_avg * 0.99f + (float)throttle * 0.01f;
2012-11-23 02:57:49 -04:00
g.throttle_cruise = throttle_avg;
}
}
#if FRAME_CONFIG == HELI_FRAME
2012-11-24 00:41:17 -04:00
// get_angle_boost - returns a throttle including compensation for roll/pitch angle
2012-11-23 02:57:49 -04:00
// throttle value should be 0 ~ 1000
2012-11-24 00:41:17 -04:00
// for traditional helicopters
static int16_t get_angle_boost(int16_t throttle)
2012-11-23 02:57:49 -04:00
{
float angle_boost_factor = cos_pitch_x * cos_roll_x;
2013-01-23 01:50:44 -04:00
angle_boost_factor = 1.0f - constrain(angle_boost_factor, .5f, 1.0f);
2012-11-23 02:57:49 -04:00
int16_t throttle_above_mid = max(throttle - motors.throttle_mid,0);
2012-11-24 00:41:17 -04:00
// to allow logging of angle boost
angle_boost = throttle_above_mid*angle_boost_factor;
return throttle + angle_boost;
2012-11-23 02:57:49 -04:00
}
2012-11-24 00:41:17 -04:00
#else // all multicopters
// get_angle_boost - returns a throttle including compensation for roll/pitch angle
// throttle value should be 0 ~ 1000
static int16_t get_angle_boost(int16_t throttle)
{
float temp = cos_pitch_x * cos_roll_x;
int16_t throttle_out;
2013-01-23 01:50:44 -04:00
temp = constrain(temp, 0.5f, 1.0f);
2013-03-01 21:28:44 -04:00
// reduce throttle if we go inverted
temp = constrain(9000-max(labs(ahrs.roll_sensor),labs(ahrs.pitch_sensor)), 0, 3000) / (3000 * temp);
// apply scale and constrain throttle
2012-11-24 00:41:17 -04:00
throttle_out = constrain((float)(throttle-g.throttle_min) * temp + g.throttle_min, g.throttle_min, 1000);
// to allow logging of angle boost
angle_boost = throttle_out - throttle;
return throttle_out;
}
#endif // FRAME_CONFIG == HELI_FRAME
2012-11-23 02:57:49 -04:00
2012-12-03 10:05:14 -04:00
// set_throttle_out - to be called by upper throttle controllers when they wish to provide throttle output directly to motors
2012-11-23 02:57:49 -04:00
// provide 0 to cut motors
2012-11-24 00:41:17 -04:00
void set_throttle_out( int16_t throttle_out, bool apply_angle_boost )
2012-11-23 02:57:49 -04:00
{
2012-11-24 00:41:17 -04:00
if( apply_angle_boost ) {
g.rc_3.servo_out = get_angle_boost(throttle_out);
}else{
g.rc_3.servo_out = throttle_out;
// clear angle_boost for logging purposes
angle_boost = 0;
}
2013-02-27 03:57:04 -04:00
// update compass with throttle value
compass.set_throttle((float)g.rc_3.servo_out/1000.0f);
2012-11-23 02:57:49 -04:00
}
// set_throttle_accel_target - to be called by upper throttle controllers to set desired vertical acceleration in earth frame
void set_throttle_accel_target( int16_t desired_acceleration )
{
2012-12-03 10:05:14 -04:00
if( g.throttle_accel_enabled ) {
throttle_accel_target_ef = desired_acceleration;
throttle_accel_controller_active = true;
}else{
// To-Do log dataflash or tlog error
cliSerial->print_P(PSTR("Err: target sent to inactive acc thr controller!\n"));
}
}
// disable_throttle_accel - disables the accel based throttle controller
// it will be re-enasbled on the next set_throttle_accel_target
// required when we wish to set motors to zero when pilot inputs zero throttle
void throttle_accel_deactivate()
{
throttle_accel_controller_active = false;
2012-11-23 02:57:49 -04:00
}
// get_throttle_accel - accelerometer based throttle controller
// returns an actual throttle output (0 ~ 1000) to be sent to the motors
static int16_t
get_throttle_accel(int16_t z_target_accel)
{
2012-12-17 01:36:05 -04:00
static float z_accel_error = 0; // The acceleration error in cm.
2013-01-08 03:41:07 -04:00
static uint32_t last_call_ms = 0; // the last time this controller was called
2012-12-17 01:36:05 -04:00
int32_t p,i,d; // used to capture pid values for logging
int16_t output;
float z_accel_meas;
2013-01-08 03:41:07 -04:00
uint32_t now = millis();
2012-11-23 02:57:49 -04:00
// Calculate Earth Frame Z acceleration
2013-01-01 22:52:53 -04:00
z_accel_meas = -(ahrs.get_accel_ef().z + GRAVITY_MSS) * 100;
2012-11-23 02:57:49 -04:00
2013-01-08 03:41:07 -04:00
// reset target altitude if this controller has just been engaged
if( now - last_call_ms > 100 ) {
// Reset Filter
z_accel_error = 0;
} else {
// calculate accel error and Filter with fc = 2 Hz
2013-01-10 14:42:24 -04:00
z_accel_error = z_accel_error + 0.11164f * (constrain(z_target_accel - z_accel_meas, -32000, 32000) - z_accel_error);
2013-01-08 03:41:07 -04:00
}
last_call_ms = now;
2012-11-23 02:57:49 -04:00
// separately calculate p, i, d values for logging
p = g.pid_throttle_accel.get_p(z_accel_error);
// freeze I term if we've breached throttle limits
if( motors.reached_limit(AP_MOTOR_THROTTLE_LIMIT) ) {
i = g.pid_throttle_accel.get_integrator();
}else{
2013-01-23 01:50:44 -04:00
i = g.pid_throttle_accel.get_i(z_accel_error, .01f);
2012-11-23 02:57:49 -04:00
}
2013-01-23 01:50:44 -04:00
d = g.pid_throttle_accel.get_d(z_accel_error, .01f);
2012-11-23 02:57:49 -04:00
//
// limit the rate
output = constrain(p+i+d+g.throttle_cruise, g.throttle_min, g.throttle_max);
#if LOGGING_ENABLED == ENABLED
// log output if PID loggins is on and we are tuning the yaw
2013-01-08 01:45:12 -04:00
if( g.log_bitmask & MASK_LOG_PID && (g.radio_tuning == CH6_THR_ACCEL_KP || g.radio_tuning == CH6_THR_ACCEL_KI || g.radio_tuning == CH6_THR_ACCEL_KD) ) {
pid_log_counter++;
if( pid_log_counter >= 10 ) { // (update rate / desired output rate) = (50hz / 10hz) = 5hz
pid_log_counter = 0;
Log_Write_PID(CH6_THR_ACCEL_KP, z_accel_error, p, i, d, output, tuning_value);
2012-11-23 02:57:49 -04:00
}
}
#endif
return output;
}
2013-01-30 11:25:41 -04:00
// get_pilot_desired_throttle - transform pilot's throttle input to make cruise throttle mid stick
// used only for manual throttle modes
// returns throttle output 0 to 1000
#define THROTTLE_IN_MIDDLE 500 // the throttle mid point
static int16_t get_pilot_desired_throttle(int16_t throttle_control)
{
int16_t throttle_out;
// exit immediately in the simple cases
if( throttle_control == 0 || g.throttle_mid == 500) {
return throttle_control;
}
// ensure reasonable throttle values
throttle_control = constrain(throttle_control,0,1000);
g.throttle_mid = constrain(g.throttle_mid,300,700);
// check throttle is above, below or in the deadband
if (throttle_control < THROTTLE_IN_MIDDLE) {
// below the deadband
2013-01-31 10:21:18 -04:00
throttle_out = g.throttle_min + ((float)(throttle_control-g.throttle_min))*((float)(g.throttle_mid - g.throttle_min))/((float)(500-g.throttle_min));
2013-01-30 11:25:41 -04:00
}else if(throttle_control > THROTTLE_IN_MIDDLE) {
// above the deadband
throttle_out = g.throttle_mid + ((float)(throttle_control-500))*(float)(1000-g.throttle_mid)/500.0f;
}else{
// must be in the deadband
throttle_out = g.throttle_mid;
}
return throttle_out;
}
2012-11-23 02:57:49 -04:00
// get_pilot_desired_climb_rate - transform pilot's throttle input to
// climb rate in cm/s. we use radio_in instead of control_in to get the full range
// without any deadzone at the bottom
#define THROTTLE_IN_DEADBAND 100 // the throttle input channel's deadband in PWM
#define THROTTLE_IN_DEADBAND_TOP (THROTTLE_IN_MIDDLE+THROTTLE_IN_DEADBAND) // top of the deadband
#define THROTTLE_IN_DEADBAND_BOTTOM (THROTTLE_IN_MIDDLE-THROTTLE_IN_DEADBAND) // bottom of the deadband
static int16_t get_pilot_desired_climb_rate(int16_t throttle_control)
{
int16_t desired_rate = 0;
// throttle failsafe check
2013-03-16 05:14:21 -03:00
if( ap.failsafe_radio ) {
2012-11-23 02:57:49 -04:00
return 0;
}
// ensure a reasonable throttle value
throttle_control = constrain(throttle_control,0,1000);
// check throttle is above, below or in the deadband
if (throttle_control < THROTTLE_IN_DEADBAND_BOTTOM) {
// below the deadband
2012-12-21 23:52:49 -04:00
desired_rate = (int32_t)g.pilot_velocity_z_max * (throttle_control-THROTTLE_IN_DEADBAND_BOTTOM) / (THROTTLE_IN_MIDDLE - THROTTLE_IN_DEADBAND);
2012-11-23 02:57:49 -04:00
}else if (throttle_control > THROTTLE_IN_DEADBAND_TOP) {
// above the deadband
2012-12-21 23:52:49 -04:00
desired_rate = (int32_t)g.pilot_velocity_z_max * (throttle_control-THROTTLE_IN_DEADBAND_TOP) / (THROTTLE_IN_MIDDLE - THROTTLE_IN_DEADBAND);
2012-11-23 02:57:49 -04:00
}else{
// must be in the deadband
desired_rate = 0;
}
2012-11-24 03:45:28 -04:00
// desired climb rate for logging
desired_climb_rate = desired_rate;
2012-11-23 02:57:49 -04:00
return desired_rate;
}
// get_pilot_desired_acceleration - transform pilot's throttle input to a desired acceleration
// default upper and lower bounds are 500 cm/s/s (roughly 1/2 a G)
// returns acceleration in cm/s/s
static int16_t get_pilot_desired_acceleration(int16_t throttle_control)
{
int32_t desired_accel = 0;
// throttle failsafe check
2013-03-16 05:14:21 -03:00
if( ap.failsafe_radio ) {
2012-11-23 02:57:49 -04:00
return 0;
}
// ensure a reasonable throttle value
throttle_control = constrain(throttle_control,0,1000);
// check throttle is above, below or in the deadband
if (throttle_control < THROTTLE_IN_DEADBAND_BOTTOM) {
// below the deadband
2012-11-25 05:07:33 -04:00
desired_accel = (int32_t)ACCELERATION_MAX_Z * (throttle_control-THROTTLE_IN_DEADBAND_BOTTOM) / (THROTTLE_IN_MIDDLE - THROTTLE_IN_DEADBAND);
2012-11-23 02:57:49 -04:00
}else if(throttle_control > THROTTLE_IN_DEADBAND_TOP) {
// above the deadband
2012-11-25 05:07:33 -04:00
desired_accel = (int32_t)ACCELERATION_MAX_Z * (throttle_control-THROTTLE_IN_DEADBAND_TOP) / (THROTTLE_IN_MIDDLE - THROTTLE_IN_DEADBAND);
2012-11-23 02:57:49 -04:00
}else{
// must be in the deadband
desired_accel = 0;
}
return desired_accel;
}
// get_pilot_desired_direct_alt - transform pilot's throttle input to a desired altitude
// return altitude in cm between 0 to 10m
static int32_t get_pilot_desired_direct_alt(int16_t throttle_control)
{
int32_t desired_alt = 0;
2013-03-16 05:14:21 -03:00
// radio failsafe check
if( ap.failsafe_radio ) {
2012-11-23 02:57:49 -04:00
return 0;
}
// ensure a reasonable throttle value
throttle_control = constrain(throttle_control,0,1000);
desired_alt = throttle_control;
return desired_alt;
}
2013-03-17 01:24:49 -03:00
// get_initial_alt_hold - get new target altitude based on current altitude and climb rate
static int32_t
get_initial_alt_hold( int32_t alt_cm, int16_t climb_rate_cms)
{
int32_t target_alt;
int32_t linear_distance; // half the distace we swap between linear and sqrt and the distace we offset sqrt.
int32_t linear_velocity; // the velocity we swap between linear and sqrt.
linear_velocity = ALT_HOLD_ACCEL_MAX/g.pi_alt_hold.kP();
if (abs(climb_rate_cms) < linear_velocity) {
target_alt = alt_cm + climb_rate_cms/g.pi_alt_hold.kP();
} else {
linear_distance = ALT_HOLD_ACCEL_MAX/(2*g.pi_alt_hold.kP()*g.pi_alt_hold.kP());
if (climb_rate_cms > 0){
target_alt = alt_cm + linear_distance + (int32_t)climb_rate_cms*(int32_t)climb_rate_cms/(2*ALT_HOLD_ACCEL_MAX);
} else {
target_alt = alt_cm - ( linear_distance + (int32_t)climb_rate_cms*(int32_t)climb_rate_cms/(2*ALT_HOLD_ACCEL_MAX) );
}
}
return constrain(target_alt, alt_cm - ALT_HOLD_INIT_MAX_OVERSHOOT, alt_cm + ALT_HOLD_INIT_MAX_OVERSHOOT);
}
2012-11-23 02:57:49 -04:00
// get_throttle_rate - calculates desired accel required to achieve desired z_target_speed
// sets accel based throttle controller target
static void
2013-04-02 09:10:34 -03:00
get_throttle_rate(float z_target_speed)
2012-11-23 02:57:49 -04:00
{
2013-01-08 03:41:07 -04:00
static uint32_t last_call_ms = 0;
2012-12-17 01:36:05 -04:00
static float z_rate_error = 0; // The velocity error in cm.
2013-04-01 04:10:44 -03:00
static float z_target_speed_last = 0; // The requested speed from the previous iteration
2012-11-23 02:57:49 -04:00
int32_t p,i,d; // used to capture pid values for logging
2013-04-02 09:10:34 -03:00
int32_t output; // the target acceleration if the accel based throttle is enabled, otherwise the output to be sent to the motors
2013-01-08 03:41:07 -04:00
uint32_t now = millis();
2012-11-23 02:57:49 -04:00
2013-01-08 03:41:07 -04:00
// reset target altitude if this controller has just been engaged
if( now - last_call_ms > 100 ) {
// Reset Filter
z_rate_error = 0;
2013-04-01 04:10:44 -03:00
output = 0;
2013-01-08 03:41:07 -04:00
} else {
// calculate rate error and filter with cut off frequency of 2 Hz
2013-01-10 14:42:24 -04:00
z_rate_error = z_rate_error + 0.20085f * ((z_target_speed - climb_rate) - z_rate_error);
2013-04-01 04:10:44 -03:00
// feed forward acceleration based on change in desired speed.
2013-04-02 09:10:34 -03:00
output = (z_target_speed - z_target_speed_last) * 50.0f; // To-Do: replace 50 with dt
2013-01-08 03:41:07 -04:00
}
last_call_ms = now;
2012-11-23 02:57:49 -04:00
2013-04-01 04:10:44 -03:00
// store target speed for next iteration
z_target_speed_last = z_target_speed;
2012-11-23 02:57:49 -04:00
// separately calculate p, i, d values for logging
p = g.pid_throttle.get_p(z_rate_error);
// freeze I term if we've breached throttle limits
if(motors.reached_limit(AP_MOTOR_THROTTLE_LIMIT)) {
i = g.pid_throttle.get_integrator();
}else{
i = g.pid_throttle.get_i(z_rate_error, .02);
}
d = g.pid_throttle.get_d(z_rate_error, .02);
2013-04-02 09:10:34 -03:00
// consolidate and constrain target acceleration
2013-04-01 04:10:44 -03:00
output += p+i+d;
2013-04-02 09:10:34 -03:00
output = constrain_int32(output, -32000, 32000);
2012-11-23 02:57:49 -04:00
#if LOGGING_ENABLED == ENABLED
// log output if PID loggins is on and we are tuning the yaw
2012-12-15 00:06:24 -04:00
if( g.log_bitmask & MASK_LOG_PID && (g.radio_tuning == CH6_THROTTLE_KP || g.radio_tuning == CH6_THROTTLE_KI || g.radio_tuning == CH6_THROTTLE_KD) ) {
2013-01-08 01:45:12 -04:00
pid_log_counter++;
if( pid_log_counter >= 10 ) { // (update rate / desired output rate) = (50hz / 10hz) = 5hz
pid_log_counter = 0;
2012-12-03 10:05:14 -04:00
Log_Write_PID(CH6_THROTTLE_KP, z_rate_error, p, i, d, output, tuning_value);
2012-11-23 02:57:49 -04:00
}
}
#endif
2012-12-03 10:05:14 -04:00
// send output to accelerometer based throttle controller if enabled otherwise send directly to motors
if( g.throttle_accel_enabled ) {
// set target for accel based throttle controller
set_throttle_accel_target(output);
}else{
set_throttle_out(g.throttle_cruise+output, true);
}
2012-11-24 09:50:09 -04:00
2013-04-01 03:45:23 -03:00
// limit loiter & waypoint navigation from causing too much lean
// To-Do: ensure that this limit is cleared when this throttle controller is not running so that loiter is not left constrained for Position mode
wp_nav.set_angle_limit(4500 - constrain((z_rate_error - 100) * 10, 0, 3500));
2012-11-24 09:50:09 -04:00
// update throttle cruise
2012-12-05 21:27:38 -04:00
// TO-DO: this may not be correct because g.rc_3.servo_out has not been updated for this iteration
2012-11-24 09:50:09 -04:00
if( z_target_speed == 0 ) {
update_throttle_cruise(g.rc_3.servo_out);
}
2012-11-23 02:57:49 -04:00
}
// get_throttle_althold - hold at the desired altitude in cm
// updates accel based throttle controller targets
// Note: max_climb_rate is an optional parameter to allow reuse of this function by landing controller
static void
2012-12-21 23:52:49 -04:00
get_throttle_althold(int32_t target_alt, int16_t min_climb_rate, int16_t max_climb_rate)
2012-11-23 02:57:49 -04:00
{
2013-01-08 01:45:12 -04:00
int32_t alt_error;
2013-03-17 01:24:49 -03:00
float desired_rate;
int32_t linear_distance; // half the distace we swap between linear and sqrt and the distace we offset sqrt.
2012-12-17 01:36:05 -04:00
// calculate altitude error
2013-01-08 01:45:12 -04:00
alt_error = target_alt - current_loc.alt;
2012-12-17 01:36:05 -04:00
2012-12-26 08:18:16 -04:00
// check kP to avoid division by zero
if( g.pi_alt_hold.kP() != 0 ) {
2013-03-17 01:24:49 -03:00
linear_distance = ALT_HOLD_ACCEL_MAX/(2*g.pi_alt_hold.kP()*g.pi_alt_hold.kP());
2013-01-08 01:45:12 -04:00
if( alt_error > 2*linear_distance ) {
2013-03-17 01:24:49 -03:00
desired_rate = safe_sqrt(2*ALT_HOLD_ACCEL_MAX*(alt_error-linear_distance));
2013-01-08 01:45:12 -04:00
}else if( alt_error < -2*linear_distance ) {
2013-03-17 01:24:49 -03:00
desired_rate = -safe_sqrt(2*ALT_HOLD_ACCEL_MAX*(-alt_error-linear_distance));
2012-12-26 08:18:16 -04:00
}else{
2013-01-08 01:45:12 -04:00
desired_rate = g.pi_alt_hold.get_p(alt_error);
2012-12-26 08:18:16 -04:00
}
2012-12-22 22:17:26 -04:00
}else{
2012-12-26 08:18:16 -04:00
desired_rate = 0;
2012-12-17 01:36:05 -04:00
}
2012-11-23 02:57:49 -04:00
2012-12-21 23:52:49 -04:00
desired_rate = constrain(desired_rate, min_climb_rate, max_climb_rate);
2012-11-23 02:57:49 -04:00
// call rate based throttle controller which will update accel based throttle controller targets
get_throttle_rate(desired_rate);
2013-01-08 01:45:12 -04:00
// update altitude error reported to GCS
altitude_error = alt_error;
2012-11-23 02:57:49 -04:00
// TO-DO: enabled PID logging for this controller
}
2013-01-08 03:41:07 -04:00
// get_throttle_althold_with_slew - altitude controller with slew to avoid step changes in altitude target
// calls normal althold controller which updates accel based throttle controller targets
static void
2013-03-16 23:14:35 -03:00
get_throttle_althold_with_slew(int32_t target_alt, int16_t min_climb_rate, int16_t max_climb_rate)
2013-01-08 03:41:07 -04:00
{
// limit target altitude change
2013-01-23 01:50:44 -04:00
controller_desired_alt += constrain(target_alt-controller_desired_alt, min_climb_rate*0.02f, max_climb_rate*0.02f);
2013-01-08 03:41:07 -04:00
// do not let target altitude get too far from current altitude
controller_desired_alt = constrain(controller_desired_alt,current_loc.alt-750,current_loc.alt+750);
get_throttle_althold(controller_desired_alt, min_climb_rate-250, max_climb_rate+250); // 250 is added to give head room to alt hold controller
}
2012-12-21 23:52:49 -04:00
// get_throttle_rate_stabilized - rate controller with additional 'stabilizer'
// 'stabilizer' ensure desired rate is being met
// calls normal throttle rate controller which updates accel based throttle controller targets
static void
get_throttle_rate_stabilized(int16_t target_rate)
{
2013-01-10 14:42:24 -04:00
controller_desired_alt += target_rate * 0.02f;
2012-12-21 23:52:49 -04:00
// do not let target altitude get too far from current altitude
2013-01-08 03:41:07 -04:00
controller_desired_alt = constrain(controller_desired_alt,current_loc.alt-750,current_loc.alt+750);
2012-12-21 23:52:49 -04:00
2013-04-08 23:58:01 -03:00
// update target altitude for reporting purposes
set_target_alt_for_reporting(controller_desired_alt);
2012-12-21 23:52:49 -04:00
2013-01-08 03:41:07 -04:00
get_throttle_althold(controller_desired_alt, -g.pilot_velocity_z_max-250, g.pilot_velocity_z_max+250); // 250 is added to give head room to alt hold controller
2012-12-21 23:52:49 -04:00
}
2012-11-23 02:57:49 -04:00
// get_throttle_land - high level landing logic
// sends the desired acceleration in the accel based throttle controller
// called at 50hz
static void
get_throttle_land()
{
2012-12-29 00:51:14 -04:00
// if we are above 10m and the sonar does not sense anything perform regular alt hold descent
2013-01-08 03:41:07 -04:00
if (current_loc.alt >= LAND_START_ALT && !(g.sonar_enabled && sonar_alt_health >= SONAR_ALT_HEALTH_MAX)) {
2013-04-18 02:52:21 -03:00
get_throttle_althold_with_slew(LAND_START_ALT, -wp_nav.get_descent_velocity(), -abs(g.land_speed));
2012-11-23 02:57:49 -04:00
}else{
2012-11-24 09:50:09 -04:00
get_throttle_rate_stabilized(-abs(g.land_speed));
2012-11-29 08:08:19 -04:00
// detect whether we have landed by watching for minimum throttle and now movement
if (abs(climb_rate) < 20 && (g.rc_3.servo_out <= get_angle_boost(g.throttle_min) || g.pid_throttle_accel.get_integrator() <= -150)) {
2012-11-24 09:50:09 -04:00
if( land_detector < LAND_DETECTOR_TRIGGER ) {
land_detector++;
}else{
set_land_complete(true);
2013-03-16 05:14:21 -03:00
if( g.rc_3.control_in == 0 || ap.failsafe_radio ) {
2012-11-24 09:50:09 -04:00
init_disarm_motors();
}
}
}else{
// we've sensed movement up or down so decrease land_detector
if (land_detector > 0 ) {
land_detector--;
}
}
2012-11-23 02:57:49 -04:00
}
}
2012-12-29 00:51:14 -04:00
// get_throttle_surface_tracking - hold copter at the desired distance above the ground
// updates accel based throttle controller targets
static void
get_throttle_surface_tracking(int16_t target_rate)
{
static float target_sonar_alt = 0; // The desired altitude in cm above the ground
static uint32_t last_call_ms = 0;
2013-01-08 03:41:07 -04:00
float distance_error;
float sonar_induced_slew_rate;
2012-12-29 00:51:14 -04:00
uint32_t now = millis();
// reset target altitude if this controller has just been engaged
2013-01-08 03:41:07 -04:00
if( now - last_call_ms > 200 ) {
target_sonar_alt = sonar_alt + controller_desired_alt - current_loc.alt;
2012-12-29 00:51:14 -04:00
}
2013-01-08 03:41:07 -04:00
last_call_ms = now;
2012-12-29 00:51:14 -04:00
2013-01-10 14:42:24 -04:00
target_sonar_alt += target_rate * 0.02f;
2012-12-29 00:51:14 -04:00
2013-01-08 03:41:07 -04:00
distance_error = (target_sonar_alt-sonar_alt);
sonar_induced_slew_rate = constrain(fabs(THR_SURFACE_TRACKING_P * distance_error),0,THR_SURFACE_TRACKING_VELZ_MAX);
2012-12-29 00:51:14 -04:00
// do not let target altitude get too far from current altitude above ground
// Note: the 750cm limit is perhaps too wide but is consistent with the regular althold limits and helps ensure a smooth transition
target_sonar_alt = constrain(target_sonar_alt,sonar_alt-750,sonar_alt+750);
2013-01-08 03:41:07 -04:00
controller_desired_alt = current_loc.alt+(target_sonar_alt-sonar_alt);
2013-04-08 23:58:01 -03:00
// update target altitude for reporting purposes
set_target_alt_for_reporting(controller_desired_alt);
2012-12-29 00:51:14 -04:00
2013-01-08 03:41:07 -04:00
get_throttle_althold_with_slew(controller_desired_alt, target_rate-sonar_induced_slew_rate, target_rate+sonar_induced_slew_rate); // VELZ_MAX limits how quickly we react
2012-12-29 00:51:14 -04:00
}
2012-11-23 02:57:49 -04:00
/*
* reset all I integrators
*/
static void reset_I_all(void)
{
reset_rate_I();
reset_stability_I();
reset_throttle_I();
reset_optflow_I();
// This is the only place we reset Yaw
g.pi_stabilize_yaw.reset_I();
}
static void reset_rate_I()
{
g.pid_rate_roll.reset_I();
g.pid_rate_pitch.reset_I();
g.pid_rate_yaw.reset_I();
}
static void reset_optflow_I(void)
{
g.pid_optflow_roll.reset_I();
g.pid_optflow_pitch.reset_I();
of_roll = 0;
of_pitch = 0;
}
static void reset_throttle_I(void)
{
// For Altitude Hold
g.pi_alt_hold.reset_I();
g.pid_throttle.reset_I();
2012-12-17 01:36:05 -04:00
g.pid_throttle_accel.reset_I();
2012-11-23 02:57:49 -04:00
}
2013-01-31 03:30:03 -04:00
static void set_accel_throttle_I_from_pilot_throttle(int16_t pilot_throttle)
2013-01-11 23:20:37 -04:00
{
// shift difference between pilot's throttle and hover throttle into accelerometer I
2013-01-31 03:30:03 -04:00
g.pid_throttle_accel.set_integrator(pilot_throttle-g.throttle_cruise);
2013-01-11 23:20:37 -04:00
}
2012-11-23 02:57:49 -04:00
static void reset_stability_I(void)
{
// Used to balance a quad
// This only needs to be reset during Auto-leveling in flight
g.pi_stabilize_roll.reset_I();
g.pi_stabilize_pitch.reset_I();
2012-01-09 00:53:54 -04:00
}