stmbl/shared/comps/stp.c

86 lines
1.7 KiB
C

#include "hal.h"
#include "math.h"
#include "defines.h"
#include "angle.h"
HAL_COMP(stp);
HAL_PIN(target);
HAL_PIN(pos);
HAL_PIN(pos_in);
HAL_PIN(pos_out);
HAL_PIN(vel_in);
HAL_PIN(vel_out);
HAL_PIN(acc_out);
HAL_PIN(max_vel);
HAL_PIN(max_acc);
HAL_PIN(dtg);
struct stp_ctx_t {
float p0;
float p;
float v0;
float pold;
float vold;
};
static void rt_func(float period, volatile void *ctx_ptr, volatile hal_pin_inst_t *pin_ptr) {
struct stp_ctx_t *ctx = (struct stp_ctx_t *)ctx_ptr;
struct stp_pin_ctx_t *pins = (struct stp_pin_ctx_t *)pin_ptr;
//float p = PIN(pos_in);
//float ctx->v0 = PIN(vel_in);
//ctx->p0 += minus(p, ctx->pold);
//ctx->pold = p;
float p1 = PIN(target);
float vmax = PIN(max_vel);
float amax = MAX(PIN(max_acc), 0.1);
float vel = SIGN(p1 - ctx->p0) * amax * sqrtf(2.0 * ABS(p1 - ctx->p0) / amax);
vel = LIMIT(vel, vmax);
vel = CLAMP(vel, ctx->v0 - amax * period, ctx->v0 + amax * period);
if(ABS(vel) > ABS(p1 - ctx->p0) / period) {
vel = (p1 - ctx->p0) / period;
}
float acc = (vel - ctx->vold) / period;
ctx->vold = vel;
ctx->v0 = vel;
ctx->p0 += ctx->v0 * period;
ctx->p += ctx->p0 - ctx->pold;
ctx->pold = ctx->p0;
//p += vel * period;
PIN(pos_out) = mod(ctx->p);
PIN(vel_out) = vel;
PIN(acc_out) = acc;
PIN(dtg) = p1 - ctx->p0;
}
const hal_comp_t stp_comp_struct = {
.name = "stp",
.nrt = 0,
.rt = rt_func,
.frt = 0,
.nrt_init = 0,
.hw_init = 0,
.rt_start = 0,
.frt_start = 0,
.rt_stop = 0,
.frt_stop = 0,
.ctx_size = sizeof(struct stp_ctx_t),
.pin_count = sizeof(struct stp_pin_ctx_t) / sizeof(struct hal_pin_inst_t),
};