forked from kektobiologist/motion-simulation
-
Notifications
You must be signed in to change notification settings - Fork 0
/
tracking.cpp
63 lines (57 loc) · 1.89 KB
/
tracking.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
#include "tracking.hpp"
#include <math.h>
#include <QDebug>
#include <assert.h>
#define PREDICTION_PACKET_DELAY 4
double sgn(double x) {
return (x > 0) - (0 > x);
}
MiscData Tracker::genControls(Pose s, int &vl, int &vr, int prevVl, int prevVr, double t) {
Q_UNUSED(prevVl);
Q_UNUSED(prevVr);
Pose ref(traj->x(t)*fieldXConvert, traj->y(t)*fieldXConvert, traj->theta(t));
//qDebug() << ref.x() << endl;
double ur1 = traj->v(t);
double ur2 = traj->thetad(t);
// err coordinates are in cm!
Error err(ref, s);
double zeta = 1, omegan = 10, g = .1;
double k1 = 2*zeta*omegan;
double k3 = k1;
double k2;
if (fabs(ur1) < 1) {
k2 = g*fabs(ur1);
} else {
k2 = (omegan*omegan - ur2*ur2)/fabs(ur1);
}
// NOTE: hardcoding k3 = 0!
// k3 = 0;
// v = K.e
// err.e3 = 0;
double v1 = -k1*err.e1;
double v2 = -sgn(ur1)*k2*err.e2 -k3*err.e3;
double v = ur1*cos(err.e3) - v1;
double w = ur2-v2;
v/=Constants::ticksToCmS;
vl = (v - Constants::d*w/2);
vr = (v + Constants::d*w/2);
// vl, vr transform now
double vsat_ticks = Constants::vsat/Constants::ticksToCmS;
if (vl > vsat_ticks)
vl = vsat_ticks;
else if (vl < -vsat_ticks)
vl = -vsat_ticks;
if (vr > vsat_ticks)
vr = vsat_ticks;
else if (vr < -vsat_ticks)
vr = -vsat_ticks;
// just for testing, calculate reference vl and vr also
double vl_ref = (ur1/Constants::ticksToCmS - Constants::d*ur2/2);
double vr_ref = (ur1/Constants::ticksToCmS + Constants::d*ur2/2);
return MiscData(ur1, ur2, v1, v2, t, v, w, vl, vr, vl_ref, vr_ref);
}
Pose Tracker::getNewStartPose(double t){
double timeLMs = 16.;
//qDebug() << traj->y(t)*fieldXConvert << "Dqa "<< endl ;
return Pose(traj->x(t + 4*timeLMs*0.001)*fieldXConvert, traj->y(t + 4*timeLMs*0.001)*fieldXConvert, traj->theta(t + 3*timeLMs*0.001));
}