Jon Freeman / Mbed 2 deprecated Loco_10_17_SD

Electric Locomotive control system. Touch screen driver control, includes regenerative braking, drives 4 brushless motors, displays speed MPH, system volts and power

Dependencies:   BSP_DISCO_F746NG FastPWM LCD_DISCO_F746NG SD_DISCO_F746NG TS_DISCO_F746NG mbed

costab.cpp@0:23cc72b18e74, 2017-11-12 (annotated)

Committer:
JonFreeman
Date:
Sun Nov 12 06:26:29 2017 +0000
Revision:
0:23cc72b18e74
Child:
1:8ef34deb5177
Electric loco controller up to Nov 2017 includes SD odometer

Who changed what in which revision?

UserRevisionLine numberNew contents of line
JonFreeman 0:23cc72b18e74 1 #include "mbed.h"
JonFreeman 0:23cc72b18e74 2 static const double HALF_PI = 2.0 * atan(1.0);
JonFreeman 0:23cc72b18e74 3 static const double TWO_PI = 8.0 * atan(1.0);
JonFreeman 0:23cc72b18e74 4
JonFreeman 0:23cc72b18e74 5 double jcos (double angle) {
JonFreeman 0:23cc72b18e74 6 //At costabgen with 180 points
JonFreeman 0:23cc72b18e74 7 static const double costab[] = {
JonFreeman 0:23cc72b18e74 8 +1.000000, +0.999962, +0.999848, +0.999657, +0.999391, +0.999048, +0.998630, +0.998135,
JonFreeman 0:23cc72b18e74 9 +0.997564, +0.996917, +0.996195, +0.995396, +0.994522, +0.993572, +0.992546, +0.991445,
JonFreeman 0:23cc72b18e74 10 +0.990268, +0.989016, +0.987688, +0.986286, +0.984808, +0.983255, +0.981627, +0.979925,
JonFreeman 0:23cc72b18e74 11 +0.978148, +0.976296, +0.974370, +0.972370, +0.970296, +0.968148, +0.965926, +0.963630,
JonFreeman 0:23cc72b18e74 12 +0.961262, +0.958820, +0.956305, +0.953717, +0.951057, +0.948324, +0.945519, +0.942641,
JonFreeman 0:23cc72b18e74 13 +0.939693, +0.936672, +0.933580, +0.930418, +0.927184, +0.923880, +0.920505, +0.917060,
JonFreeman 0:23cc72b18e74 14 +0.913545, +0.909961, +0.906308, +0.902585, +0.898794, +0.894934, +0.891007, +0.887011,
JonFreeman 0:23cc72b18e74 15 +0.882948, +0.878817, +0.874620, +0.870356, +0.866025, +0.861629, +0.857167, +0.852640,
JonFreeman 0:23cc72b18e74 16 +0.848048, +0.843391, +0.838671, +0.833886, +0.829038, +0.824126, +0.819152, +0.814116,
JonFreeman 0:23cc72b18e74 17 +0.809017, +0.803857, +0.798636, +0.793353, +0.788011, +0.782608, +0.777146, +0.771625,
JonFreeman 0:23cc72b18e74 18 +0.766044, +0.760406, +0.754710, +0.748956, +0.743145, +0.737277, +0.731354, +0.725374,
JonFreeman 0:23cc72b18e74 19 +0.719340, +0.713250, +0.707107, +0.700909, +0.694658, +0.688355, +0.681998, +0.675590,
JonFreeman 0:23cc72b18e74 20 +0.669131, +0.662620, +0.656059, +0.649448, +0.642788, +0.636078, +0.629320, +0.622515,
JonFreeman 0:23cc72b18e74 21 +0.615661, +0.608761, +0.601815, +0.594823, +0.587785, +0.580703, +0.573576, +0.566406,
JonFreeman 0:23cc72b18e74 22 +0.559193, +0.551937, +0.544639, +0.537300, +0.529919, +0.522499, +0.515038, +0.507538,
JonFreeman 0:23cc72b18e74 23 +0.500000, +0.492424, +0.484810, +0.477159, +0.469472, +0.461749, +0.453990, +0.446198,
JonFreeman 0:23cc72b18e74 24 +0.438371, +0.430511, +0.422618, +0.414693, +0.406737, +0.398749, +0.390731, +0.382683,
JonFreeman 0:23cc72b18e74 25 +0.374607, +0.366501, +0.358368, +0.350207, +0.342020, +0.333807, +0.325568, +0.317305,
JonFreeman 0:23cc72b18e74 26 +0.309017, +0.300706, +0.292372, +0.284015, +0.275637, +0.267238, +0.258819, +0.250380,
JonFreeman 0:23cc72b18e74 27 +0.241922, +0.233445, +0.224951, +0.216440, +0.207912, +0.199368, +0.190809, +0.182236,
JonFreeman 0:23cc72b18e74 28 +0.173648, +0.165048, +0.156434, +0.147809, +0.139173, +0.130526, +0.121869, +0.113203,
JonFreeman 0:23cc72b18e74 29 +0.104528, +0.095846, +0.087156, +0.078459, +0.069756, +0.061049, +0.052336, +0.043619,
JonFreeman 0:23cc72b18e74 30 +0.034899, +0.026177, +0.017452, +0.008727, +0.000000, 0.0 } ;
JonFreeman 0:23cc72b18e74 31 //End of costab
JonFreeman 0:23cc72b18e74 32 const int costab_points_i = (sizeof(costab) / sizeof(double)) - 1;
JonFreeman 0:23cc72b18e74 33 const double costab_points_d = (double) costab_points_i;
JonFreeman 0:23cc72b18e74 34 bool isneg = false;
JonFreeman 0:23cc72b18e74 35 int offset;
JonFreeman 0:23cc72b18e74 36 while (angle > TWO_PI) angle -= TWO_PI;
JonFreeman 0:23cc72b18e74 37 while (angle < 0.0) angle += TWO_PI;
JonFreeman 0:23cc72b18e74 38 // now got 0.0 <= angle <= 2PI
JonFreeman 0:23cc72b18e74 39 offset = (int)(angle * costab_points_d / HALF_PI); // ang2 is 0.0 <= ang2 <= 2.0 * PI
JonFreeman 0:23cc72b18e74 40 if (offset > 3 * costab_points_i) { // quadrant = 3;
JonFreeman 0:23cc72b18e74 41 offset = (4 * costab_points_i) - offset - 1;
JonFreeman 0:23cc72b18e74 42 }
JonFreeman 0:23cc72b18e74 43 else if (offset > 2 * costab_points_i) { // quadrant = 2;
JonFreeman 0:23cc72b18e74 44 offset -= 2 * costab_points_i;
JonFreeman 0:23cc72b18e74 45 isneg = true;
JonFreeman 0:23cc72b18e74 46 }
JonFreeman 0:23cc72b18e74 47 else if (offset > 1 * costab_points_i) { // quadrant = 1;
JonFreeman 0:23cc72b18e74 48 offset = (2 * costab_points_i) - offset - 1;
JonFreeman 0:23cc72b18e74 49 isneg = true;
JonFreeman 0:23cc72b18e74 50 } // else quadrant = 0;
JonFreeman 0:23cc72b18e74 51 if (offset < 0) offset = 0;
JonFreeman 0:23cc72b18e74 52 if (offset > costab_points_i) offset = costab_points_i;
JonFreeman 0:23cc72b18e74 53 if (isneg) return - costab[offset];
JonFreeman 0:23cc72b18e74 54 else return costab[offset];
JonFreeman 0:23cc72b18e74 55 }
JonFreeman 0:23cc72b18e74 56
JonFreeman 0:23cc72b18e74 57 double jsin (double angle) {
JonFreeman 0:23cc72b18e74 58 return -jcos (angle + HALF_PI);
JonFreeman 0:23cc72b18e74 59 }
JonFreeman 0:23cc72b18e74 60
JonFreeman 0:23cc72b18e74 61