Version 1: coming together. With PD controller but without inverse kinematics
Dependencies: Encoder MODSERIAL mbed
Fork of motoraansturing_met_EMG by
Revision 1:1c22ce9f370b, committed 2013-11-01
- Comitter:
- Tess
- Date:
- Fri Nov 01 10:46:19 2013 +0000
- Parent:
- 0:73dcc8fb9900
- Child:
- 2:55bff07c1058
- Commit message:
- Version 1: coming together. Without inverse kinematics but with switch
Changed in this revision
main.cpp | Show annotated file Show diff for this revision Revisions of this file |
--- a/main.cpp Fri Nov 01 08:24:02 2013 +0000 +++ b/main.cpp Fri Nov 01 10:46:19 2013 +0000 @@ -3,44 +3,44 @@ #include "MODSERIAL.h" //high pass filter constantes 15Hz cutoff 4e orde - #define NUM0 0.2754 // constante - #define NUM1 -1.1017 // z^-1 - #define NUM2 1.6525 // z^-2etc. - #define NUM3 -1.1017 - #define NUM4 0.2754 +#define NUM0 0.2754 // constante +#define NUM1 -1.1017 // z^-1 +#define NUM2 1.6525 // z^-2etc. +#define NUM3 -1.1017 +#define NUM4 0.2754 - #define DEN0 1 // constante - #define DEN1 -1.5704 - #define DEN2 1.2756 - #define DEN3 -0.4844 - #define DEN4 0.0762 +#define DEN0 1 // constante +#define DEN1 -1.5704 +#define DEN2 1.2756 +#define DEN3 -0.4844 +#define DEN4 0.0762 //lowpass filter constantes 40 Hz 4e orde - #define NUM0_2 0.4328 // constante - #define NUM1_2 1.7314 // z^-1 - #define NUM2_2 2.5971 // z^-2etc. - #define NUM3_2 1.7314 - #define NUM4_2 0.4328 +#define NUM0_2 0.4328 // constante +#define NUM1_2 1.7314 // z^-1 +#define NUM2_2 2.5971 // z^-2etc. +#define NUM3_2 1.7314 +#define NUM4_2 0.4328 - #define DEN0_2 1 // constante - #define DEN1_2 2.3695 - #define DEN2_2 2.3140 - #define DEN3_2 1.0547 - #define DEN4_2 0.1874 +#define DEN0_2 1 // constante +#define DEN1_2 2.3695 +#define DEN2_2 2.3140 +#define DEN3_2 1.0547 +#define DEN4_2 0.1874 //lowpass filter constantes 3Hz 4e orde - #define NUM0_3 0.0000624 // constante - #define NUM1_3 0.0002495 // z^-1 - #define NUM2_3 0.0003743 // z^-2etc. - #define NUM3_3 0.0002495 - #define NUM4_3 0.0000624 +#define NUM0_3 0.0000624 // constante +#define NUM1_3 0.0002495 // z^-1 +#define NUM2_3 0.0003743 // z^-2etc. +#define NUM3_3 0.0002495 +#define NUM4_3 0.0000624 - #define DEN0_3 1 // constante - #define DEN1_3 -3.5078 - #define DEN2_3 4.6409 - #define DEN3_3 -2.7427 - #define DEN4_3 0.6105 +#define DEN0_3 1 // constante +#define DEN1_3 -3.5078 +#define DEN2_3 4.6409 +#define DEN3_3 -2.7427 +#define DEN4_3 0.6105 /******************************************************************************* * * @@ -48,191 +48,249 @@ * * ********************************************************************************/ -// dit is voor schakelaar die aan en uit wil gaan -DigitalIn toggle(PTD7); - -void toggle_on() -{ -} - -void toggle_off() -{ - // do nothing -} - /** keep_in_range -> float in, and keep_in_range if less than min, or larger than max **/ void keep_in_range(float * in, float min, float max); -/** variable to show when a new loop can be started*/ -/** volatile means that it can be changed in an */ -/** interrupt routine, and that that change is vis-*/ -/** ible in the main loop. */ +/** variable to show when a new loop can be started */ +/** volatile means that it can be changed in an interrupt routine, and that that change is visible in the main loop. */ volatile bool looptimerflag; /** function called by Ticker "looptimer" */ -/** variable 'looptimerflag' is set to 'true' */ -/** each time the looptimer expires. */ +/** variable 'looptimerflag' is set to 'true' each time the looptimer expires.*/ void setlooptimerflag(void) { looptimerflag = true; } - //emg variabelen - float emg_value_biceps, emg_value_triceps, emg_value_flexoren, emg_value_extensoren, dy; - AnalogIn emg_biceps(PTB0); //Analog input - AnalogIn emg_triceps(PTB1); - AnalogIn emg_flexoren(PTB2); - AnalogIn emg_extensoren(PTB3); +//emg variabelen -> analoge input +float emg_value_biceps, emg_value_triceps, emg_value_flexoren, emg_value_extensoren, dy; +AnalogIn emg_biceps(PTB0); +AnalogIn emg_triceps(PTB1); +AnalogIn emg_flexoren(PTB2); +AnalogIn emg_extensoren(PTB3); /* DIT IS DE FILTER FUNCTIE! aanroepen door "filter(signaal nummer)" filter(1): biceps meten filter(2): triceps meten filter(3): flexoren meten - filter(3): extensoren meten + filter(3): extensoren meten */ -float filter(int sig_number){ +float filter(int sig_number) +{ float sig_out; - - // eerst variabelen definieren - + + // eerst variabelen definiëren + //biceps - //filter 1 + //filter 1 float in0_biceps =0; static float in1_biceps =0, in2_biceps = 0, in3_biceps = 0, in4_biceps = 0; static float out0_biceps = 0, out1_biceps = 0 , out2_biceps = 0, out3_biceps = 0, out4_biceps = 0; - //filter 2 + //filter 2 float in0_2_biceps =0; static float in1_2_biceps =0, in2_2_biceps = 0, in3_2_biceps = 0, in4_2_biceps = 0; static float out0_2_biceps = 0, out1_2_biceps = 0 , out2_2_biceps = 0, out3_2_biceps = 0, out4_2_biceps = 0; - //filter 3 + //filter 3 float in0_3_biceps =0; static float in1_3_biceps =0, in2_3_biceps = 0, in3_3_biceps = 0, in4_3_biceps = 0; static float out0_3_biceps = 0, out1_3_biceps = 0 , out2_3_biceps = 0, out3_3_biceps = 0, out4_3_biceps = 0; - + //triceps - //filter 1 + //filter 1 float in0_triceps =0; static float in1_triceps =0, in2_triceps = 0, in3_triceps = 0, in4_triceps = 0; static float out0_triceps = 0, out1_triceps = 0 , out2_triceps = 0, out3_triceps = 0, out4_triceps = 0; - //filter 2 + //filter 2 float in0_2_triceps =0; static float in1_2_triceps =0, in2_2_triceps = 0, in3_2_triceps = 0, in4_2_triceps = 0; static float out0_2_triceps = 0, out1_2_triceps = 0 , out2_2_triceps = 0, out3_2_triceps = 0, out4_2_triceps = 0; - //filter 3 + //filter 3 float in0_3_triceps =0; static float in1_3_triceps =0, in2_3_triceps = 0, in3_3_triceps = 0, in4_3_triceps = 0; static float out0_3_triceps = 0, out1_3_triceps = 0 , out2_3_triceps = 0, out3_3_triceps = 0, out4_3_triceps = 0; - + //flexoren - //filter 1 + //filter 1 float in0_flexoren =0; static float in1_flexoren =0, in2_flexoren = 0, in3_flexoren = 0, in4_flexoren = 0; static float out0_flexoren = 0, out1_flexoren = 0 , out2_flexoren = 0, out3_flexoren = 0, out4_flexoren = 0; - //filter 2 + //filter 2 float in0_2_flexoren =0; static float in1_2_flexoren =0, in2_2_flexoren = 0, in3_2_flexoren = 0, in4_2_flexoren = 0; static float out0_2_flexoren = 0, out1_2_flexoren = 0 , out2_2_flexoren = 0, out3_2_flexoren = 0, out4_2_flexoren = 0; - //filter 3 + //filter 3 float in0_3_flexoren =0; static float in1_3_flexoren =0, in2_3_flexoren = 0, in3_3_flexoren = 0, in4_3_flexoren = 0; static float out0_3_flexoren = 0, out1_3_flexoren = 0 , out2_3_flexoren = 0, out3_3_flexoren = 0, out4_3_flexoren = 0; - + //extensoren - //filter 1 + //filter 1 float in0_extensoren =0; static float in1_extensoren =0, in2_extensoren = 0, in3_extensoren = 0, in4_extensoren = 0; static float out0_extensoren = 0, out1_extensoren = 0 , out2_extensoren = 0, out3_extensoren = 0, out4_extensoren = 0; - //filter 2 + //filter 2 float in0_2_extensoren =0; static float in1_2_extensoren =0, in2_2_extensoren = 0, in3_2_extensoren = 0, in4_2_extensoren = 0; static float out0_2_extensoren = 0, out1_2_extensoren = 0 , out2_2_extensoren = 0, out3_2_extensoren = 0, out4_2_extensoren = 0; - //filter 3 + //filter 3 float in0_3_extensoren =0; static float in1_3_extensoren =0, in2_3_extensoren = 0, in3_3_extensoren = 0, in4_3_extensoren = 0; static float out0_3_extensoren = 0, out1_3_extensoren = 0 , out2_3_extensoren = 0, out3_3_extensoren = 0, out4_3_extensoren = 0; - - - switch(sig_number){ - case 1: + + + switch(sig_number) { + case 1: // signaal filteren op 15 Hz HIGHPASS - in4_biceps = in3_biceps; in3_biceps = in2_biceps; in2_biceps = in1_biceps; in1_biceps = in0_biceps; + in4_biceps = in3_biceps; + in3_biceps = in2_biceps; + in2_biceps = in1_biceps; + in1_biceps = in0_biceps; in0_biceps = emg_biceps.read(); - out4_biceps = out3_biceps; out3_biceps = out2_biceps; out2_biceps = out1_biceps; out1_biceps = out0_biceps; - out0_biceps = (NUM0*in0_biceps + NUM1*in1_biceps + NUM2*in2_biceps + NUM3*in3_biceps + NUM4*in4_biceps - DEN1*out1_biceps - DEN2*out2_biceps - DEN3*out3_biceps - DEN4*out4_biceps ) / DEN0; - + out4_biceps = out3_biceps; + out3_biceps = out2_biceps; + out2_biceps = out1_biceps; + out1_biceps = out0_biceps; + out0_biceps = (NUM0*in0_biceps + NUM1*in1_biceps + NUM2*in2_biceps + NUM3*in3_biceps + NUM4*in4_biceps - DEN1*out1_biceps - DEN2*out2_biceps - DEN3*out3_biceps - DEN4*out4_biceps ) / DEN0; + //signaal filteren op 40 HZ LOWPASS - in4_2_biceps = in3_2_biceps; in3_2_biceps = in2_2_biceps; in2_2_biceps = in1_2_biceps; in1_2_biceps = in0_2_biceps; + in4_2_biceps = in3_2_biceps; + in3_2_biceps = in2_2_biceps; + in2_2_biceps = in1_2_biceps; + in1_2_biceps = in0_2_biceps; in0_2_biceps = out0_biceps; - out4_2_biceps = out3_2_biceps; out3_2_biceps = out2_2_biceps; out2_2_biceps = out1_2_biceps; out1_2_biceps = out0_2_biceps; + out4_2_biceps = out3_2_biceps; + out3_2_biceps = out2_2_biceps; + out2_2_biceps = out1_2_biceps; + out1_2_biceps = out0_2_biceps; out0_2_biceps = (NUM0_2*in0_2_biceps + NUM1_2*in1_2_biceps + NUM2_2*in2_2_biceps + NUM3_2*in3_2_biceps + NUM4_2*in4_2_biceps - DEN1_2*out1_2_biceps - DEN2_2*out2_2_biceps - DEN3_2*out3_2_biceps - DEN4_2*out4_2_biceps ) / DEN0_2; - + //signaal filteren op 5Hz LOWPASS - in4_3_biceps = in3_3_biceps; in3_3_biceps = in2_3_biceps; in2_3_biceps = in1_3_biceps; in1_3_biceps = in0_3_biceps; + in4_3_biceps = in3_3_biceps; + in3_3_biceps = in2_3_biceps; + in2_3_biceps = in1_3_biceps; + in1_3_biceps = in0_3_biceps; in0_3_biceps = abs(out0_2_biceps); - out4_3_biceps = out3_3_biceps; out3_3_biceps = out2_3_biceps; out2_3_biceps = out1_3_biceps; out1_3_biceps = out0_3_biceps; - out0_3_biceps = (NUM0_3*in0_3_biceps + NUM1_3*in1_3_biceps + NUM2_3*in2_3_biceps + NUM3_3*in3_3_biceps + NUM4_3*in4_3_biceps - DEN1_3*out1_3_biceps - DEN2_3*out2_3_biceps - DEN3_3*out3_3_biceps - DEN4_3*out4_3_biceps ) / DEN0_3; + out4_3_biceps = out3_3_biceps; + out3_3_biceps = out2_3_biceps; + out2_3_biceps = out1_3_biceps; + out1_3_biceps = out0_3_biceps; + out0_3_biceps = (NUM0_3*in0_3_biceps + NUM1_3*in1_3_biceps + NUM2_3*in2_3_biceps + NUM3_3*in3_3_biceps + NUM4_3*in4_3_biceps - DEN1_3*out1_3_biceps - DEN2_3*out2_3_biceps - DEN3_3*out3_3_biceps - DEN4_3*out4_3_biceps ) / DEN0_3; sig_out = out0_3_biceps; break; case 2: // signaal filteren op 15 Hz HIGHPASS - in4_triceps = in3_triceps; in3_triceps = in2_triceps; in2_triceps = in1_triceps; in1_triceps = in0_triceps; + in4_triceps = in3_triceps; + in3_triceps = in2_triceps; + in2_triceps = in1_triceps; + in1_triceps = in0_triceps; in0_triceps = emg_triceps.read(); - out4_triceps = out3_triceps; out3_triceps = out2_triceps; out2_triceps = out1_triceps; out1_triceps = out0_triceps; - out0_triceps = (NUM0*in0_triceps + NUM1*in1_triceps + NUM2*in2_triceps + NUM3*in3_triceps + NUM4*in4_triceps - DEN1*out1_triceps - DEN2*out2_triceps - DEN3*out3_triceps - DEN4*out4_triceps ) / DEN0; - + out4_triceps = out3_triceps; + out3_triceps = out2_triceps; + out2_triceps = out1_triceps; + out1_triceps = out0_triceps; + out0_triceps = (NUM0*in0_triceps + NUM1*in1_triceps + NUM2*in2_triceps + NUM3*in3_triceps + NUM4*in4_triceps - DEN1*out1_triceps - DEN2*out2_triceps - DEN3*out3_triceps - DEN4*out4_triceps ) / DEN0; + //signaal filteren op 40 HZ LOWPASS - in4_2_triceps = in3_2_triceps; in3_2_triceps = in2_2_triceps; in2_2_triceps = in1_2_triceps; in1_2_triceps = in0_2_triceps; + in4_2_triceps = in3_2_triceps; + in3_2_triceps = in2_2_triceps; + in2_2_triceps = in1_2_triceps; + in1_2_triceps = in0_2_triceps; in0_2_triceps = out0_triceps; - out4_2_triceps = out3_2_triceps; out3_2_triceps = out2_2_triceps; out2_2_triceps = out1_2_triceps; out1_2_triceps = out0_2_triceps; + out4_2_triceps = out3_2_triceps; + out3_2_triceps = out2_2_triceps; + out2_2_triceps = out1_2_triceps; + out1_2_triceps = out0_2_triceps; out0_2_triceps = (NUM0_2*in0_2_triceps + NUM1_2*in1_2_triceps + NUM2_2*in2_2_triceps + NUM3_2*in3_2_triceps + NUM4_2*in4_2_triceps - DEN1_2*out1_2_triceps - DEN2_2*out2_2_triceps - DEN3_2*out3_2_triceps - DEN4_2*out4_2_triceps ) / DEN0_2; - + //signaal filteren op 5Hz LOWPASS - in4_3_triceps = in3_3_triceps; in3_3_triceps = in2_3_triceps; in2_3_triceps = in1_3_triceps; in1_3_triceps = in0_3_triceps; + in4_3_triceps = in3_3_triceps; + in3_3_triceps = in2_3_triceps; + in2_3_triceps = in1_3_triceps; + in1_3_triceps = in0_3_triceps; in0_3_triceps = abs(out0_2_triceps); - out4_3_triceps = out3_3_triceps; out3_3_triceps = out2_3_triceps; out2_3_triceps = out1_3_triceps; out1_3_triceps = out0_3_triceps; - out0_3_triceps = (NUM0_3*in0_3_triceps + NUM1_3*in1_3_triceps + NUM2_3*in2_3_triceps + NUM3_3*in3_3_triceps + NUM4_3*in4_3_triceps - DEN1_3*out1_3_triceps - DEN2_3*out2_3_triceps - DEN3_3*out3_3_triceps - DEN4_3*out4_3_triceps ) / DEN0_3; + out4_3_triceps = out3_3_triceps; + out3_3_triceps = out2_3_triceps; + out2_3_triceps = out1_3_triceps; + out1_3_triceps = out0_3_triceps; + out0_3_triceps = (NUM0_3*in0_3_triceps + NUM1_3*in1_3_triceps + NUM2_3*in2_3_triceps + NUM3_3*in3_3_triceps + NUM4_3*in4_3_triceps - DEN1_3*out1_3_triceps - DEN2_3*out2_3_triceps - DEN3_3*out3_3_triceps - DEN4_3*out4_3_triceps ) / DEN0_3; sig_out = out0_3_triceps; break; case 3: // signaal filteren op 15 Hz HIGHPASS - in4_flexoren = in3_flexoren; in3_flexoren = in2_flexoren; in2_flexoren = in1_flexoren; in1_flexoren = in0_flexoren; + in4_flexoren = in3_flexoren; + in3_flexoren = in2_flexoren; + in2_flexoren = in1_flexoren; + in1_flexoren = in0_flexoren; in0_flexoren = emg_flexoren.read(); - out4_flexoren = out3_flexoren; out3_flexoren = out2_flexoren; out2_flexoren = out1_flexoren; out1_flexoren = out0_flexoren; - out0_flexoren = (NUM0*in0_flexoren + NUM1*in1_flexoren + NUM2*in2_flexoren + NUM3*in3_flexoren + NUM4*in4_flexoren - DEN1*out1_flexoren - DEN2*out2_flexoren - DEN3*out3_flexoren - DEN4*out4_flexoren ) / DEN0; - + out4_flexoren = out3_flexoren; + out3_flexoren = out2_flexoren; + out2_flexoren = out1_flexoren; + out1_flexoren = out0_flexoren; + out0_flexoren = (NUM0*in0_flexoren + NUM1*in1_flexoren + NUM2*in2_flexoren + NUM3*in3_flexoren + NUM4*in4_flexoren - DEN1*out1_flexoren - DEN2*out2_flexoren - DEN3*out3_flexoren - DEN4*out4_flexoren ) / DEN0; + //signaal filteren op 40 HZ LOWPASS - in4_2_flexoren = in3_2_flexoren; in3_2_flexoren = in2_2_flexoren; in2_2_flexoren = in1_2_flexoren; in1_2_flexoren = in0_2_flexoren; + in4_2_flexoren = in3_2_flexoren; + in3_2_flexoren = in2_2_flexoren; + in2_2_flexoren = in1_2_flexoren; + in1_2_flexoren = in0_2_flexoren; in0_2_flexoren = out0_flexoren; - out4_2_flexoren = out3_2_flexoren; out3_2_flexoren = out2_2_flexoren; out2_2_flexoren = out1_2_flexoren; out1_2_flexoren = out0_2_flexoren; + out4_2_flexoren = out3_2_flexoren; + out3_2_flexoren = out2_2_flexoren; + out2_2_flexoren = out1_2_flexoren; + out1_2_flexoren = out0_2_flexoren; out0_2_flexoren = (NUM0_2*in0_2_flexoren + NUM1_2*in1_2_flexoren + NUM2_2*in2_2_flexoren + NUM3_2*in3_2_flexoren + NUM4_2*in4_2_flexoren - DEN1_2*out1_2_flexoren - DEN2_2*out2_2_flexoren - DEN3_2*out3_2_flexoren - DEN4_2*out4_2_flexoren ) / DEN0_2; - + //signaal filteren op 5Hz LOWPASS - in4_3_flexoren = in3_3_flexoren; in3_3_flexoren = in2_3_flexoren; in2_3_flexoren = in1_3_flexoren; in1_3_flexoren = in0_3_flexoren; + in4_3_flexoren = in3_3_flexoren; + in3_3_flexoren = in2_3_flexoren; + in2_3_flexoren = in1_3_flexoren; + in1_3_flexoren = in0_3_flexoren; in0_3_flexoren = abs(out0_2_flexoren); - out4_3_flexoren = out3_3_flexoren; out3_3_flexoren = out2_3_flexoren; out2_3_flexoren = out1_3_flexoren; out1_3_flexoren = out0_3_flexoren; - out0_3_flexoren = (NUM0_3*in0_3_flexoren + NUM1_3*in1_3_flexoren + NUM2_3*in2_3_flexoren + NUM3_3*in3_3_flexoren + NUM4_3*in4_3_flexoren - DEN1_3*out1_3_flexoren - DEN2_3*out2_3_flexoren - DEN3_3*out3_3_flexoren - DEN4_3*out4_3_flexoren ) / DEN0_3; + out4_3_flexoren = out3_3_flexoren; + out3_3_flexoren = out2_3_flexoren; + out2_3_flexoren = out1_3_flexoren; + out1_3_flexoren = out0_3_flexoren; + out0_3_flexoren = (NUM0_3*in0_3_flexoren + NUM1_3*in1_3_flexoren + NUM2_3*in2_3_flexoren + NUM3_3*in3_3_flexoren + NUM4_3*in4_3_flexoren - DEN1_3*out1_3_flexoren - DEN2_3*out2_3_flexoren - DEN3_3*out3_3_flexoren - DEN4_3*out4_3_flexoren ) / DEN0_3; sig_out = out0_3_flexoren; break; case 4: // signaal filteren op 15 Hz HIGHPASS - in4_extensoren = in3_extensoren; in3_extensoren = in2_extensoren; in2_extensoren = in1_extensoren; in1_extensoren = in0_extensoren; + in4_extensoren = in3_extensoren; + in3_extensoren = in2_extensoren; + in2_extensoren = in1_extensoren; + in1_extensoren = in0_extensoren; in0_extensoren = emg_extensoren.read(); - out4_extensoren = out3_extensoren; out3_extensoren = out2_extensoren; out2_extensoren = out1_extensoren; out1_extensoren = out0_extensoren; - out0_extensoren = (NUM0*in0_extensoren + NUM1*in1_extensoren + NUM2*in2_extensoren + NUM3*in3_extensoren + NUM4*in4_extensoren - DEN1*out1_extensoren - DEN2*out2_extensoren - DEN3*out3_extensoren - DEN4*out4_extensoren ) / DEN0; - + out4_extensoren = out3_extensoren; + out3_extensoren = out2_extensoren; + out2_extensoren = out1_extensoren; + out1_extensoren = out0_extensoren; + out0_extensoren = (NUM0*in0_extensoren + NUM1*in1_extensoren + NUM2*in2_extensoren + NUM3*in3_extensoren + NUM4*in4_extensoren - DEN1*out1_extensoren - DEN2*out2_extensoren - DEN3*out3_extensoren - DEN4*out4_extensoren ) / DEN0; + //signaal filteren op 40 HZ LOWPASS - in4_2_extensoren = in3_2_extensoren; in3_2_extensoren = in2_2_extensoren; in2_2_extensoren = in1_2_extensoren; in1_2_extensoren = in0_2_extensoren; + in4_2_extensoren = in3_2_extensoren; + in3_2_extensoren = in2_2_extensoren; + in2_2_extensoren = in1_2_extensoren; + in1_2_extensoren = in0_2_extensoren; in0_2_extensoren = out0_extensoren; - out4_2_extensoren = out3_2_extensoren; out3_2_extensoren = out2_2_extensoren; out2_2_extensoren = out1_2_extensoren; out1_2_extensoren = out0_2_extensoren; + out4_2_extensoren = out3_2_extensoren; + out3_2_extensoren = out2_2_extensoren; + out2_2_extensoren = out1_2_extensoren; + out1_2_extensoren = out0_2_extensoren; out0_2_extensoren = (NUM0_2*in0_2_extensoren + NUM1_2*in1_2_extensoren + NUM2_2*in2_2_extensoren + NUM3_2*in3_2_extensoren + NUM4_2*in4_2_extensoren - DEN1_2*out1_2_extensoren - DEN2_2*out2_2_extensoren - DEN3_2*out3_2_extensoren - DEN4_2*out4_2_extensoren ) / DEN0_2; - + //signaal filteren op 5Hz LOWPASS - in4_3_extensoren = in3_3_extensoren; in3_3_extensoren = in2_3_extensoren; in2_3_extensoren = in1_3_extensoren; in1_3_extensoren = in0_3_extensoren; + in4_3_extensoren = in3_3_extensoren; + in3_3_extensoren = in2_3_extensoren; + in2_3_extensoren = in1_3_extensoren; + in1_3_extensoren = in0_3_extensoren; in0_3_extensoren = abs(out0_2_extensoren); - out4_3_extensoren = out3_3_extensoren; out3_3_extensoren = out2_3_extensoren; out2_3_extensoren = out1_3_extensoren; out1_3_extensoren = out0_3_extensoren; - out0_3_extensoren = (NUM0_3*in0_3_extensoren + NUM1_3*in1_3_extensoren + NUM2_3*in2_3_extensoren + NUM3_3*in3_3_extensoren + NUM4_3*in4_3_extensoren - DEN1_3*out1_3_extensoren - DEN2_3*out2_3_extensoren - DEN3_3*out3_3_extensoren - DEN4_3*out4_3_extensoren ) / DEN0_3; + out4_3_extensoren = out3_3_extensoren; + out3_3_extensoren = out2_3_extensoren; + out2_3_extensoren = out1_3_extensoren; + out1_3_extensoren = out0_3_extensoren; + out0_3_extensoren = (NUM0_3*in0_3_extensoren + NUM1_3*in1_3_extensoren + NUM2_3*in2_3_extensoren + NUM3_3*in3_3_extensoren + NUM4_3*in4_3_extensoren - DEN1_3*out1_3_extensoren - DEN2_3*out2_3_extensoren - DEN3_3*out3_3_extensoren - DEN4_3*out4_3_extensoren ) / DEN0_3; sig_out = out0_3_extensoren; break; } @@ -242,23 +300,14 @@ int main() { //LOCAL VARIABLES - /*Potmeter input*/ - AnalogIn potmeterA(PTC2); - AnalogIn potmeterB(PTB2); - /* Encoder, using my encoder library */ - /* First pin should be PTDx or PTAx */ - /* because those pins can be used as */ - /* InterruptIn */ - Encoder motorA(PTD4,PTC8); + Encoder motorA(PTD4,PTC8); // encoder, using my encoder library Encoder motorB(PTD0,PTD2); - /* MODSERIAL to get non-blocking Serial*/ - MODSERIAL pc(USBTX,USBRX); - /* PWM control to motor */ - PwmOut pwm_motorA(PTA12); - PwmOut pwm_motorB(PTA5); - /* Direction pin */ - DigitalOut motordirA(PTD3); - DigitalOut motordirB(PTD1); + MODSERIAL pc(USBTX,USBRX); // MODSERIAL to get non-blocking Serial + PwmOut pwm_motorA(PTA12); // PWM control to motor + PwmOut pwm_motorB(PTA5); // PWM control to motor + DigitalOut motordirA(PTD3); // Direction pin + DigitalOut motordirB(PTD1); // Direction pin + /* variable to store setpoint in */ float setpointA; float setpointB; @@ -274,30 +323,23 @@ float pwm_to_motorB; float pwm_to_rechtsonder_motorA; float pwm_to_rechtsonder_motorB; - + + /* variable for PD controller*/ const float dt = 0.002; - float Kp = 0.001; //0.0113 - float Ki = 0.0759; - float Kd = 0.0004342; + float Kp = 0.001; //0.0208 + float Kd = 0.00004342; //0.0006897 float error_t0_A = 0; float error_t0_B = 0; float error_ti_A; float error_ti_B; - float error_t_1_A; - float error_t_1_B; float P_regelaar_A; float P_regelaar_B; - float I_regelaar_A; - float I_regelaar_B; float D_regelaar_A; float D_regelaar_B; float output_regelaar_A; float output_regelaar_B; - float integral_i_A; - float integral_i_B; - float integral_0_A = 0; - float integral_0_B = 0; + /* variable to store positions in*/ int32_t positionmotorA_t0; int32_t positionmotorB_t0; int32_t positionmotorA_t_1; @@ -307,223 +349,193 @@ //START OF CODE - while(1) { - while(!toggle); - { // wait while toggle == 0 - toggle_on(); + pc.baud(115200); // Set the baudrate (use this number in RealTerm too!) + + // In dit stukje code zorgen we ervoor dat de arm gaat draaien naar rechts en stopt als het tegen het frame komt. Eerst motor B botsen dan motor A botsen. + // motor B zit onder en motor A zit boven en dus op zijn kop (en dus setpoint moet - zijn). - /*Set the baudrate (use this number in RealTerm too!) */ - pc.baud(115200); - - // in dit stukje code zorgen we ervoor dat de arm gaat draaien naar rechts en stopt als het tegen het frame komt. Eerst motor B botsen dan motor A botsen. - // motor B zit onder en motor A zit boven en dus op zijn kop (en dus setpoint moet - zijn). + motordirB.write(0); + pwm_motorB.write(.08); + positionmotorB_t0 = motorB.getPosition(); + do { + wait(0.2); + positionmotorB_t_1 = positionmotorB_t0 ; + positionmotorB_t0 = motorB.getPosition(); + positiondifference_motorB = abs(positionmotorB_t0 - positionmotorB_t_1); + } while(positiondifference_motorB > 10); + motorB.setPosition(0); + pwm_motorB.write(0); - motordirB.write(0); - pwm_motorB.write(.08); - positionmotorB_t0 = motorB.getPosition(); - do { - wait(0.2); - positionmotorB_t_1 = positionmotorB_t0 ; - positionmotorB_t0 = motorB.getPosition(); - positiondifference_motorB = abs(positionmotorB_t0 - positionmotorB_t_1); - } while(positiondifference_motorB > 10); - motorB.setPosition(0); - pwm_motorB.write(0); - - wait(1); // willen nu even dat tussen ene actie en andere actie 1 seconde wacht. + wait(1); // willen nu even dat tussen ene actie en andere actie 1 seconde wacht. - motordirA.write(1); - pwm_motorA.write(.08); - positionmotorA_t0 = motorA.getPosition(); - do { - wait(0.2); - positionmotorA_t_1 = positionmotorA_t0 ; - positionmotorA_t0 = motorA.getPosition(); - positiondifference_motorA = abs(positionmotorA_t0 - positionmotorA_t_1); - } while(positiondifference_motorA > 10); - motorA.setPosition(0); - pwm_motorA.write(0); + motordirA.write(1); + pwm_motorA.write(.08); + positionmotorA_t0 = motorA.getPosition(); + do { + wait(0.2); + positionmotorA_t_1 = positionmotorA_t0 ; + positionmotorA_t0 = motorA.getPosition(); + positiondifference_motorA = abs(positionmotorA_t0 - positionmotorA_t_1); + } while(positiondifference_motorA > 10); + motorA.setPosition(0); + pwm_motorA.write(0); - wait(1); // willen nu even dat tussen ene actie en andere actie 1 seconde wacht. + wait(1); // willen nu even dat tussen ene actie en andere actie 1 seconde wacht. - // Hierna willen we de motor van zijn alleruiterste positie naar de x-as hebben. Hiervoor moet motor A eerst op de x-as worden gezet. Hiervoor moet motor A 4.11 graden (63) naar links. + // Hierna willen we de motor van zijn alleruiterste positie naar de x-as hebben. Hiervoor moet motor A eerst op de x-as worden gezet. Hiervoor moet motor A 4.11 graden (63) naar links. - motordirA.write(0); - pwm_motorA.write(.08); - do { - setpoint_beginA = -63; // x-as - pwm_to_begin_motorA = abs((setpoint_beginA + motorA.getPosition()) *.001); // + omdat men met een negatieve hoekverdraaiing werkt. - wait(0.2); - keep_in_range(&pwm_to_begin_motorA, -1, 1 ); - motordirA.write(0); - pwm_motorA.write(pwm_to_begin_motorA); - } while(pwm_to_begin_motorA <= 0); - motorA.setPosition(0); - pwm_motorA.write(0); - - wait(1); // willen nu even dat tussen ene actie en andere actie 1 seconde wacht. + motordirA.write(0); + pwm_motorA.write(.08); + do { + setpoint_beginA = -63; // x-as + pwm_to_begin_motorA = abs((setpoint_beginA + motorA.getPosition()) *.001); // + omdat men met een negatieve hoekverdraaiing werkt. + wait(0.2); + keep_in_range(&pwm_to_begin_motorA, -1, 1 ); + motordirA.write(0); + pwm_motorA.write(pwm_to_begin_motorA); + } while(pwm_to_begin_motorA <= 0); + motorA.setPosition(0); + pwm_motorA.write(0); - // hierna moet motor A naar de rechtsonder A4. Motor A 532. + wait(1); // willen nu even dat tussen ene actie en andere actie 1 seconde wacht. - motordirA.write(0); - pwm_motorA.write(0.08); - do { - setpoint_beginA = -532; // rechtsonder positie A4 - pwm_to_begin_motorA = abs((setpoint_beginA + motorA.getPosition()) *.001); - wait(0.2); - keep_in_range(&pwm_to_begin_motorA, -1, 1 ); - motordirA.write(0); - pwm_motorA.write(pwm_to_begin_motorA); - } while(pwm_to_begin_motorA <= 0); - pwm_motorA.write(0); - - wait(1); + // hierna moet motor A naar de rechtsonder A4. Motor A 532. - // Hierna moet motor B 21.6 (192) graden naar links om naar x-as te gaan. + motordirA.write(0); + pwm_motorA.write(0.08); + do { + setpoint_rechtsonderA = -532; // rechtsonder positie A4 + pwm_to_rechtsonder_motorA = abs((setpoint_rechtsonderA + motorA.getPosition()) *.001); + wait(0.2); + keep_in_range(&pwm_to_rechtsonder_motorA, -1, 1 ); + motordirA.write(0); + pwm_motorA.write(pwm_to_rechtsonder_motorA); + } while(pwm_to_rechtsonder_motorA <= 0); + pwm_motorA.write(0); - motordirB.write(1); - pwm_motorB.write(.08); - do { - setpoint_beginB = 192; // x-as - pwm_to_begin_motorB = abs((setpoint_beginB - motorB.getPosition()) *.001); - wait(0.2); - keep_in_range(&pwm_to_begin_motorB, -1, 1 ); - motordirB.write(1); - pwm_motorB.write(pwm_to_begin_motorB); - } while(pwm_to_begin_motorB <= 0); - motorB.setPosition(0); - pwm_motorB.write(0); + wait(1); + + // Hierna moet motor B 21.6 (192) graden naar links om naar x-as te gaan. - wait(1); // willen nu even dat tussen ene actie en andere actie 1 seconde wacht. - - // Hierna moet motor B van x-as naar de rechtsonder A4 positie. Motor B 460. + motordirB.write(1); + pwm_motorB.write(.08); + do { + setpoint_beginB = 192; // x-as + pwm_to_begin_motorB = abs((setpoint_beginB - motorB.getPosition()) *.001); + wait(0.2); + keep_in_range(&pwm_to_begin_motorB, -1, 1 ); + motordirB.write(1); + pwm_motorB.write(pwm_to_begin_motorB); + } while(pwm_to_begin_motorB <= 0); + motorB.setPosition(0); + pwm_motorB.write(0); - motordirB.write(1); - pwm_motorB.write(0.08); - do { - setpoint_beginB = 460; // rechtsonder positie A4 - pwm_to_begin_motorB = abs((setpoint_beginB - motorB.getPosition()) *.001); - wait(0.2); - keep_in_range(&pwm_to_begin_motorB, -1, 1 ); - motordirB.write(1); - pwm_motorB.write(pwm_to_begin_motorB); - } while(pwm_to_begin_motorB <= 0); - pwm_motorB.write(0); + wait(1); // willen nu even dat tussen ene actie en andere actie 1 seconde wacht. - wait(1); + // Hierna moet motor B van x-as naar de rechtsonder A4 positie. Motor B 460. - // Nu zijn de motoren gekalibreed en staan ze op de startpositie. - // Hierna het script dat EMG wordt omgezet in een positie verandering + motordirB.write(1); + pwm_motorB.write(0.08); + do { + setpoint_rechtsonderB = 460; // rechtsonder positie A4 + pwm_to_rechtsonder_motorB = abs((setpoint_rechtsonderB - motorB.getPosition()) *.001); + wait(0.2); + keep_in_range(&pwm_to_rechtsonder_motorB, -1, 1 ); + motordirB.write(1); + pwm_motorB.write(pwm_to_rechtsonder_motorB); + } while(pwm_to_rechtsonder_motorB <= 0); + pwm_motorB.write(0); - /*Create a ticker, and let it call the */ - /*function 'setlooptimerflag' every 0.01s */ - Ticker looptimer; - looptimer.attach(setlooptimerflag,0.01); + wait(1); + + // Nu zijn de motoren gekalibreed en staan ze op de startpositie (rechtsonderhoek van A4). + // Hierna het script dat EMG wordt omgezet in een positie verandering - //INFINITE LOOP - while(1) { + /*Create a ticker, and let it call the function 'setlooptimerflag' every 0.01s */ + Ticker looptimer; + looptimer.attach(setlooptimerflag,0.01); - while(looptimerflag != true); - looptimerflag = false; + //INFINITE LOOP + while(1) { + + while(looptimerflag != true); + looptimerflag = false; - // HIER EMG!! -float emg_value_biceps; - float emg_value_triceps; - float emg_value_flexoren; - float emg_value_extensoren; - float dy; - emg_value_biceps = ((100*(filter(1))-0.18)/0.49); - emg_value_triceps = ((100*(filter(2))-0.18)/0.35); - //emg_value_flexoren = 100*filter(3); - //emg_value_extensoren = 100*filter(4); - - if(emg_value_biceps < 0.10){ - emg_value_biceps=0; - } - else { - emg_value_biceps = emg_value_biceps; - } - if(emg_value_triceps < 0.20){ - emg_value_triceps=0; - } - else { - emg_value_triceps=emg_value_triceps; - } + // HIER EMG!! + float emg_value_biceps; + float emg_value_triceps; + float emg_value_flexoren; + float emg_value_extensoren; + float dy; // verschil tussen biceps en triceps, daarmee snelheid en richting aangeven. + + emg_value_biceps = ((100*(filter(1))-0.18)/0.49); // dit is om waarde tussen 0 en 1 te krijgen. filter(1) zegt biceps, 0.18 offset aftrekken, 0.49 maximale waarde, 100 omdat procent + emg_value_triceps = ((100*(filter(2))-0.18)/0.35); // 0.35 maximale waarde van triceps + //emg_value_flexoren = 100*filter(3); + //emg_value_extensoren = 100*filter(4); + if(emg_value_biceps < 0.10) { // lager dan 10% doe dan niks - threshold + emg_value_biceps=0; + } else { + emg_value_biceps = emg_value_biceps; + } + if(emg_value_triceps < 0.20) { // lager dan 20% doe dan niks - threshold - triceps fluctueren meer + emg_value_triceps=0; + } else { + emg_value_triceps=emg_value_triceps; + } - - dy = emg_value_biceps-emg_value_triceps; - dy=dy*10; - if(pc.rxBufferGetSize(0)-pc.rxBufferGetCount() > 30) - pc.printf("%.6f\n",dy); - - - - - //setpointA = (potmeterA.read()-0.09027)*(631); // bereik van 71 graden dit afhankelijk van waar nul punt zit en waar heel wil. Dus afh. van EMG lezen bij EMG wordt 0.5 - 0.09027 - //setpointB = (potmeterB.read())*(415); // bereik van 46.7 graden - //pc.printf("s: %f, %d ", setpointA, motorA.getPosition()); - //pc.printf("s: %f, %d ", setpointB, motorB.getPosition()); + dy = emg_value_biceps - emg_value_triceps; + dy = dy * 10; // Waarde om beter te zien, officieel zonder *10 -> straks veranderen + if(pc.rxBufferGetSize(0)-pc.rxBufferGetCount() > 30) + pc.printf("%.6f\n",dy); + + setpointB = (dy) * 415; // 451 staat voor een verdraaiing van 46.7 graden. Dit wil men met de uiteindelijke robot + setpointA = (dy) * 631; // 631 staat voor een verdraaiing van 71 graden. Dit wil men met de uiteindelijke robot. + // Zelf vraag ik mij af of de verdraaiing wordt genomen vanaf nul of niet. - setpointB = (dy); - //setpointB = (potmeterB.read() - 0.5) * (871/2); - // motor A moet de hoek altijd binnen 53.4 tot en met 124.3 graden liggen - // motor B moet de hoek altijd binnen 30.2 tot en met -16.5 graden liggen - keep_in_range(&setpointA, -1105, -474); // voor motor moet bereik zijn -1105 tot -474 - keep_in_range(&setpointB, -147, 269); // voor motor moet bereik zijn -147 tot 269 + /* motor A moet de hoek altijd binnen 53.4 tot en met 124.3 graden liggen + // motor B moet de hoek altijd binnen 30.2 tot en met -16.5 graden liggen + keep_in_range(&setpointA, -1105, -474); // voor motor moet bereik zijn -1105 tot -474 + keep_in_range(&setpointB, -147, 269); // voor motor moet bereik zijn -147 tot 269 */ + // Dit moeten we testen!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! - // PID regelaar voor motor A - //wait(dt); - //error_ti_A = setpointA - motorA.getPosition(); - //P_regelaar_A = Kp * error_ti_A; - //D_regelaar_A = Kd * ((error_ti_A - error_t0_A) / dt); - //integral_i_A = integral_0_A + (error_ti_A * dt); - //I_regelaar_A = Ki * integral_i_A; - //integral_0_A = integral_i_A; - //error_t0_A = error_ti_A; - //output_regelaar_A = P_regelaar_A; + //PD regelaar voor motor A + wait(dt); + error_ti_A = setpointA - motorA.getPosition(); + P_regelaar_A = Kp * error_ti_A; + D_regelaar_A = Kd * ((error_ti_A - error_t0_A) / dt); + error_t0_A = error_ti_A; + output_regelaar_A = P_regelaar_A + D_regelaar_A; - // PID regelaar voor motor B - //wait(dt); - //error_ti_B = setpointB - motorB.getPosition(); - //P_regelaar_B = Kp * error_ti_B; - //D_regelaar_B = Kd * ((error_ti_B - error_t0_B) / dt); - //integral_i_B = integral_0_B + (error_ti_B * dt); - //I_regelaar_B = Ki * integral_i_B; - //integral_0_B = integral_i_B; - //error_t0_B = error_ti_B; - //output_regelaar_B = P_regelaar_B; + //PD regelaar voor motor B + wait(dt); + error_ti_B = setpointB - motorB.getPosition(); + P_regelaar_B = Kp * error_ti_B; + D_regelaar_B = Kd * ((error_ti_B - error_t0_B) / dt); + error_t0_B = error_ti_B; + output_regelaar_B = P_regelaar_B + D_regelaar_B; - /* This is a PID-action! store in pwm_to_motor */ - pwm_to_motorA = (setpointA - motorA.getPosition())*.001; //output_regelaar_A; - pwm_to_motorB = (setpointB); //- motorB.getPosition())*.001; //output_regelaar_B; - - keep_in_range(&pwm_to_motorA, -1,1); - keep_in_range(&pwm_to_motorB, -1,1); + /* This is a PD-action! store in pwm_to_motor */ + pwm_to_motorA = output_regelaar_A; + pwm_to_motorB = output_regelaar_B; - if(pwm_to_motorA > 0) - motordirA.write(1); - else - motordirA.write(0); - if(pwm_to_motorB > 0) - motordirB.write(1); - else - motordirB.write(0); + keep_in_range(&pwm_to_motorA, -1,1); + keep_in_range(&pwm_to_motorB, -1,1); - pwm_motorA.write(abs(pwm_to_motorA)); - pwm_motorB.write(abs(pwm_to_motorB)); - } - } - while(toggle); - { // wait while toggle == 1 - toggle_off(); - pwm_motorA.write(0); - pwm_motorB.write(0); - } + if(pwm_to_motorA > 0) + motordirA.write(1); + else + motordirA.write(0); + if(pwm_to_motorB > 0) + motordirB.write(1); + else + motordirB.write(0); + + pwm_motorA.write(abs(pwm_to_motorA)); + pwm_motorB.write(abs(pwm_to_motorB)); } } - void keep_in_range(float * in, float min, float max) { *in > min ? *in < max? : *in = max: *in = min;