EMG filter step by step
Dependencies: HIDScope MODSERIAL mbed-dsp mbed
Fork of emg_filter_step_by_step by
Revision 0:e24c2f32e23e, committed 2014-10-20
- Comitter:
- Tanja2211
- Date:
- Mon Oct 20 13:06:38 2014 +0000
- Commit message:
- d
Changed in this revision
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/EMGfilter.cpp Mon Oct 20 13:06:38 2014 +0000 @@ -0,0 +1,260 @@ +#include "mbed.h" +#include "HIDScope.h" +#include "MODSERIAL.h" +#include "arm_math.h" + + +// ****** emg filter shizzle ****** + +//Define objects +AnalogIn emgB(PTB0); //Analog input bicep +AnalogIn emgT(PTB1); //Analog input tricep + +float filtered_emgB; +float filtered_emgT; +float MeanT2, MeanB1, MeanB2, MeanB3; +float filtered_emgT10_C, filtered_emgT11_C, filtered_emgT12_C, filtered_emgT13_C, filtered_emgT14_C, filtered_emgT15_C, filtered_emgT16_C, filtered_emgT17_C, filtered_emgT18_C, filtered_emgT19_C, filtered_emgT110_C, filtered_emgT111_C, filtered_emgT112_C, filtered_emgT113_C, filtered_emgT114_C, filtered_emgT115_C, filtered_emgT116_C, filtered_emgT117_C, filtered_emgT118_C, filtered_emgT119_C; +float B0, B1, B2, B3, B4, B5, B6, B7, B8, B9, MOVAVG_B; +float T0, T1, T2, T3, T4, T5, T6, T7, T8, T9, MOVAVG_T; +float drempelwaardeB1, drempelwaardeB2, drempelwaardeB3, drempelwaardeT1, drempelwaardeT2; + +MODSERIAL pc(USBTX,USBRX); + +HIDScope scope(4);//uitgang scherm + +arm_biquad_casd_df1_inst_f32 lowpass; +//constants for 50Hz lowpass +float lowpass_const[] = {0.2928920553, 0.5857841107, 0.2928920554, -0, -0.17156822136};//{a0 a1 a2 -b1 -b2} van online calculator +//state values +float lowpass_states[4]; + +arm_biquad_casd_df1_inst_f32 highpass; +//constants for 10Hz highpass +float highpass_const[] = {0.8005910267, -1.6011820533, 0.8005910267, 1.5610153913, -0.6413487154};//{a0 a1 a2 -b1 -b2} +//state values +float highpass_states[4]; + + +/** Looper function +* functions used for Ticker and Timeout should be of type void <name>(void) +* i.e. no input arguments, no output arguments. +* if you want to change a variable that you use in other places (for example in main) +* you will have to make that variable global in order to be able to reach it both from +* the function called at interrupt time, and in the main function. +* To make a variable global, define it under the includes. +* variables that are changed in the interrupt routine (written to) should be made +* 'volatile' to let the compiler know that those values may change outside the current context. +* i.e.: "volatile uint16_t emg_value;" instead of "uint16_t emg_value" +* in the example below, the variable is not re-used in the main function, and is thus declared +* local in the looper function only. +**/ + +// Calibratie + +pc.printf("Calibratie drempelwaarde Triceps stand 1\n"); +wait(0.5); +{ + int i; + int j=19; + + for (i = 0, i<=j; i ++) { + /*variable to store value in*/ + uint16_t emg_valueT1i_C; + + float emg_value_f32T1i_C; + /*put raw emg value both in red and in emg_value*/ + emg_valueT1i_C = emgT1i_C.read_u16(); // read direct ADC result, converted to 16 bit integer (0..2^16 = 0..65536 = 0..3.3V) + emg_value_f32T1i_C = emgT1i_C.read(); + + //process emg + arm_biquad_cascade_df1_f32(&highpass, &emg_value_f32T1i_C, &filtered_emgT1i_C, 1 ); + filtered_emgT1i_C = fabs(filtered_emgT1i_C); + arm_biquad_cascade_df1_f32(&lowpass, &filtered_emgT1i_C, &filtered_emgT1i_C, 1 ); + } +} +// Mean Triceps stand 1 + + +float MeanT1=filtered_emgT10_C*0.05+filtered_emgT11_C*0.05+filtered_emgT12_C*0.05+filtered_emgT13_C*0.05+filtered_emgT14_C*0.05+filtered_emgT15_C*0.05+filtered_emgT16_C*0.05+filtered_emgT17_C*0.05+filtered_emgT18_C*0.05+filtered_emgT19_C*0.05+filtered_emgT110_C*0.05+filtered_emgT111_C*0.05+filtered_emgT112_C*0.05+filtered_emgT113_C*0.05+filtered_emgT114_C*0.05+filtered_emgT115_C*0.05+filtered_emgT116_C*0.05+filtered_emgT117_C*0.05+filtered_emgT118_C*0.05+filtered_emgT119_C*0.05; + + + + + +//BICEP EMG LEZEN +void looperB() +{ + /*variable to store value in*/ + uint16_t emg_valueB; + + float emg_value_f32B; + /*put raw emg value both in red and in emg_value*/ + emg_valueB = emgB.read_u16(); // read direct ADC result, converted to 16 bit integer (0..2^16 = 0..65536 = 0..3.3V) + emg_value_f32B = emgB.read(); + + //process emg + arm_biquad_cascade_df1_f32(&highpass, &emg_value_f32B, &filtered_emgB, 1 ); + filtered_emgB = fabs(filtered_emgB); + arm_biquad_cascade_df1_f32(&lowpass, &filtered_emgB, &filtered_emgB, 1 ); + + /*send value to PC. */ + scope.set(0,emg_valueB); //uint value + scope.set(1,filtered_emgB); //processed float + scope.send(); + + // Moving Average Filter Biceps + + { + B0=filtered_emgB; + MOVAVG_B=B0*0.1+B1*0.1+B2*0.1+B3*0.1+B4*0.1+B5*0.1+B6*0.1+B7*0.1+B8*0.1+B9*0.1; + B9=B8; + B8=B7; + B7=B6; + B6=B5; + B5=B4; + B4=B3; + B3=B2; + B2=B1; + B1=B0; + } +} +// Triceps EMG lezen +void looperT() +{ + /*variable to store value in*/ + uint16_t emg_valueT; + + float emg_value_f32T; + /*put raw emg value both in red and in emg_value*/ + emg_valueT = emgT.read_u16(); // read direct ADC result, converted to 16 bit integer (0..2^16 = 0..65536 = 0..3.3V) + emg_value_f32T = emgT.read(); + + //process emg + arm_biquad_cascade_df1_f32(&highpass, &emg_value_f32T, &filtered_emgT, 1 ); + filtered_emgT = fabs(filtered_emgT); + arm_biquad_cascade_df1_f32(&lowpass, &filtered_emgT, &filtered_emgT, 1 ); + + /*send value to PC. */ + scope.set(2,emg_valueT); //uint value + scope.set(3,filtered_emgT); //processed float + scope.send(); + + // Moving Average Filter Triceps + { + T0=filtered_emgT; + MOVAVG_T=T0*0.1+T1*0.1+T2*0.1+T3*0.1+T4*0.1+T5*0.1+T6*0.1+T7*0.1+T8*0.1+T9*0.1; + + T9=T8; + T8=T7; + T7=T6; + T6=T5; + T5=T4; + T4=T3; + T3=T2; + T2=T1; + T1=T0; + } +} + +int main() +{ + Ticker log_timer; + //set up filters. Use external array for constants + arm_biquad_cascade_df1_init_f32(&lowpass,1 , lowpass_const, lowpass_states); + arm_biquad_cascade_df1_init_f32(&highpass,1 ,highpass_const, highpass_states); + + /**Here you attach the 'void looper(void)' function to the Ticker object + * The looper() function will be called every 0.01 seconds. + * Please mind that the parentheses after looper are omitted when using attach. + */ + log_timer.attach(looperB, 0.005);//?? + log_timer.attach(looperT, 0.005);//?? + while(1) { //Loop + /*Empty!*/ + /*Everything is handled by the interrupt routine now!*/ + } +} + +//filtered_emgB +//filtered_emgT + +void AntwoordT() +{ + drempelwaardeT1=4.99; + drempelwaardeT2=7; + int yT1, yT2; + + if (MOVAVG_T > drempelwaardeT1) { + yT1=1; + if (MOVAVG_T > drempelwaardeT2) { + yT2=1; + } else { + yT2=0; + } + } else { + yT1=0; + } + + int positie; + + positie=yT1+yT2; + if (positie==0) { + pc.printf("Motor 2 beweegt niet\n"); + } else { + pc.printf("Motor 2 gaat beweegen\n"); + } + if (positie==1) { + pc.printf("Motor 2 beweegt naar positie 1\n"); + } else { + pc.printf("Motor 1 beweegt niet naar positie 1\n"); + } + if (positie==2) { + pc.printf("Motor 1 beweegt naar positie 2\n"); + } else { + pc.printf("Motor 1 beweegt niet naar positie 2\n"); + } +} + + void AntwoordB() { + drempelwaardeB1=4.99; + drempelwaardeB2=6; + drempelwaardeB3=10; + int yB1; + int yB2; + int yB3; + + if (MOVAVG_B > drempelwaardeB1) { + yB1=1; + if (MOVAVG_B > drempelwaardeB2) { + yB2=1; + if (MOVAVG_B > drempelwaardeB3) { + yB3=1; + } else { + yB3=0; + } + } else { + yB2=0; + } + } else { + yB1=0; + } + + int snelheidsstand; + + snelheidsstand=yB1+yB2+yB3; + if (snelheidsstand==1) { + pc.printf("Motor 1 beweegt met snelheid 1\n"); + } else { + pc.printf("Motor 1 beweegt niet met snelheid 1\n"); + } + if (snelheidsstand==2) { + pc.printf("Motor 1 beweegt met snelheid 2\n"); + } else { + pc.printf("Motor 1 beweegt niet met snelheid 2\n"); + } + if (snelheidsstand==3) { + pc.printf("Motor 1 beweegt met snelheid 3\n"); + } else { + pc.printf("Motor 1 beweegt niet met snelheid 3\n"); + } + }
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/HIDScope.lib Mon Oct 20 13:06:38 2014 +0000 @@ -0,0 +1,1 @@ +http://mbed.org/users/tomlankhorst/code/HIDScope/#e44574634162
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/MAF.cpp Mon Oct 20 13:06:38 2014 +0000 @@ -0,0 +1,59 @@ +/* Copyright (c) <2012> MIT License + * + * Permission is hereby granted, free of charge, to any person obtaining a copy of this software + * and associated documentation files (the "Software"), to deal in the Software without restriction, + * including without limitation the rights to use, copy, modify, merge, publish, distribute, + * sublicense, and/or sell copies of the Software, and to permit persons to whom the Software is + * furnished to do so, subject to the following conditions: + * + * The above copyright notice and this permission notice shall be included in all copies or + * substantial portions of the Software. + * + * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING + * BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND + * NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, + * DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, + * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE. + * + * @section DESCRIPTION + * + * 10-Point Moving Average Filter + */ + +#include "MAF.h" + + MAF::MAF(){} + + float MAF::update(float data){ + + _k[0] = data; + _result = _k[0]*0.1 + _k[1]*0.1 + _k[2]*0.1 + _k[3]*0.1 + _k[4]*0.1 + _k[5]*0.1 + _k[6]*0.1 + _k[7]*0.1 + _k[8]*0.1 + _k[9]*0.1 ; + _k[9] = _k[8]; + _k[8] = _k[7]; + _k[7] = _k[6]; + _k[6] = _k[5]; + _k[5] = _k[4]; + _k[4] = _k[3]; + _k[3] = _k[2]; + _k[2] = _k[1]; + _k[1] = _k[0]; + + return _result; + } + + int MAF::updateInt(int data){ + + _a[0] = data; + _result = _a[0]*0.1 + _a[1]*0.1 + _a[2]*0.1 + _a[3]*0.1 + _a[4]*0.1 + _a[5]*0.1 + _a[6]*0.1 + _a[7]*0.1 + _a[8]*0.1 + _a[9]*0.1 ; + _a[9] = _a[8]; + _a[8] = _a[7]; + _a[7] = _a[6]; + _a[6] = _a[5]; + _a[5] = _a[4]; + _a[4] = _a[3]; + _a[3] = _a[2]; + _a[2] = _a[1]; + _a[1] = _a[0]; + + return _resultInt; + } \ No newline at end of file
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/MAF.h Mon Oct 20 13:06:38 2014 +0000 @@ -0,0 +1,85 @@ +/* Copyright (c) <2012> MIT License + * + * Permission is hereby granted, free of charge, to any person obtaining a copy of this software + * and associated documentation files (the "Software"), to deal in the Software without restriction, + * including without limitation the rights to use, copy, modify, merge, publish, distribute, + * sublicense, and/or sell copies of the Software, and to permit persons to whom the Software is + * furnished to do so, subject to the following conditions: + * + * The above copyright notice and this permission notice shall be included in all copies or + * substantial portions of the Software. + * + * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING + * BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND + * NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, + * DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, + * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE. + * + * @section DESCRIPTION + * + * Moving Average Filter + * + * + */ + +#ifndef MAF_H +#define MAF_H + +#include "mbed.h" + +/** 10-Point Moving Average filter + * + * Example: + * @code + * // Update the filter with gyroscope readings + * #include "mbed.h" + * #include "ITG3200.h" + * #include "MAF.h" + * + * Serial pc(USBTX, USBRX); + * ITG3200 gyro(p9, p10); + * MAF Filter(); + * int Average; + * + * int main() { + * + * gyro.setLpBandwidth(LPFBW_42HZ); + * + * while(1){ + * + * Average = Filter.updateInt(gyro.getGyroX()); + * pc.printf("%i\n",Average); + * wait_ms(50); + * } + * } + * @endcode + */ +class MAF { + + public: + + /* Creates a 10-Point Moving Average Filter object + * + */ + MAF(void); + + /* Updates the filter + * @param data new data input for the filter,returns the float average + */ + float update(float data); + + /* Updates the filter + * @param data new data input for the filter, returns the integer average + */ + int updateInt(int data); + + private : + + float _k[10]; + int _a[10]; + + float _result; + int _resultInt; +}; + +#endif \ No newline at end of file
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/MODSERIAL.lib Mon Oct 20 13:06:38 2014 +0000 @@ -0,0 +1,1 @@ +http://mbed.org/users/Sissors/code/MODSERIAL/#180e968a751e
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/mbed-dsp.lib Mon Oct 20 13:06:38 2014 +0000 @@ -0,0 +1,1 @@ +http://mbed.org/teams/mbed-official/code/mbed-dsp/#7a284390b0ce
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/mbed.bld Mon Oct 20 13:06:38 2014 +0000 @@ -0,0 +1,1 @@ +http://mbed.org/users/mbed_official/code/mbed/builds/9327015d4013 \ No newline at end of file