try this for a change

Dependencies:   HIDScope MODSERIAL QEI biquadFilter mbed

Fork of EMG_5 by Ralph Gerlings

Revision:
23:e5db011bd410
Parent:
22:68ab712b62b2
Child:
24:26659f1de039
diff -r 68ab712b62b2 -r e5db011bd410 main.cpp
--- a/main.cpp	Wed Nov 02 12:55:41 2016 +0000
+++ b/main.cpp	Thu Nov 03 12:17:16 2016 +0000
@@ -4,68 +4,185 @@
 #include "MODSERIAL.h"
 
 //Define objects
-AnalogIn    emg1_in( A0 );
-AnalogIn    emg2_in( A1 );
+    AnalogIn    emg1_in( A0 );
+    AnalogIn    emg2_in( A1 );
+    AnalogIn    emg3_in( A2 );
+    DigitalIn   max_reader1( SW2 );
+    DigitalIn   max_reader3( SW3 );
+
 
-Ticker      sample_timer;
-HIDScope    scope( 3 );
-DigitalOut  red(LED_RED);
-DigitalOut  blue(LED_BLUE);
-DigitalOut  green(LED_GREEN);
-DigitalOut  emg1_out( D8 );
-DigitalOut  emg2_out( D9 );
-MODSERIAL   pc(USBTX, USBRX);
+    Ticker      main_timer;
+    Ticker      max_read1;
+    Ticker      max_read3;
+    HIDScope    scope( 5 );
+    DigitalOut  red(LED_RED);
+    DigitalOut  blue(LED_BLUE);
+    DigitalOut  green(LED_GREEN);
+    MODSERIAL   pc(USBTX, USBRX);
 
 
 // EMG variables
-double emg1;
-double emg1highfilter;
-double emg1notchfilter;
-double emg1abs;
-double emg1lowfilter;
-double emgpeak;
-
-double emg2;
-double emg2highfilter;
-double emg2notchfilter;
-double emg2abs;
-double emg2lowfilter;
+    //Right Biceps
+    double emg1;
+    double emg1highfilter;
+    double emg1notchfilter;
+    double emg1abs;
+    double emg1lowfilter;
+    double emg1peak;
+    // Left Biceps
+    double emg2;
+    double emg2highfilter;
+    double emg2notchfilter;
+    double emg2abs;
+    double emg2lowfilter;
+    double emg2peak;
+    double max1;
+    double maxpart1;
+    // Left Lower Arm OR Triceps
+    double emg3;
+    double emg3highfilter;
+    double emg3notchfilter;
+    double emg3abs;
+    double emg3lowfilter;
+    double emg3peak;
+    double max3;
+    double maxpart3;
 
 // Filter settings
-BiQuad filterhigh(9.704e-01,-1.9408,9.704e-01,-1.9389,9.427e-01);
-BiQuad filternotch(9.495e-01,-1.8062,9.495e-01,-1.8062,8.992e-01);
-BiQuad filterlow(1.439e-02,2.794e-02,1.397e-02,-1.7229,7.788e-01);
-BiQuad filterpeak(1.0878,-1.950,8.771e-01,-1.95032,9.5499e-01);
-
-// Filtering
-void filter() {
-    emg1 = emg1_in.read();
-    emg1highfilter = filterhigh.step(emg1);
-    emg1notchfilter = filternotch.step(emg1highfilter);
-    emg1abs = fabs(emg1notchfilter);
-    emg1lowfilter = filterlow.step(emg1abs);
-    emgpeak = filterpeak.step(emg1lowfilter);
+    // Right Biceps
+    BiQuad filterhigh1(9.704e-01,-1.9408,9.704e-01,-1.9389,9.427e-01); /* Filter at 10 Hz */
+    BiQuad filternotch1(9.495e-01,-1.8062,9.495e-01,-1.8062,8.992e-01); /* Filter at 50 Hz */
+    BiQuad filterlow1(2.119e-03,4.238e-03,2.119e-03,-1.9016,9.101e-01); /* Filter at 15 Hz */
+    BiQuad filterpeak1(1.0123,-1.983,9.761e-01,-1.9838,9.8855e-01); /* Gain peak at 11 Hz */
+    // Left Biceps
+    BiQuad filterhigh2(9.704e-01,-1.9408,9.704e-01,-1.9389,9.427e-01);
+    BiQuad filternotch2(9.495e-01,-1.8062,9.495e-01,-1.8062,8.992e-01);
+    BiQuad filterlow3(2.119e-03,4.238e-03,2.119e-03,-1.9016,9.101e-01);
+    BiQuad filterpeak2(1.0123,-1.983,9.761e-01,-1.9838,9.8855e-01);
+    // Left Lower Arm OR Triceps
+    BiQuad filterhigh3(9.704e-01,-1.9408,9.704e-01,-1.9389,9.427e-01);
+    BiQuad filternotch3(9.495e-01,-1.8062,9.495e-01,-1.8062,8.992e-01);
+    BiQuad filterlow3(2.119e-03,4.238e-03,2.119e-03,-1.9016,9.101e-01);
+    BiQuad filterpeak3(1.0123,-1.983,9.761e-01,-1.9838,9.8855e-01);
+    //
 
-    emg2 = emg2_in.read();
-    emg2highfilter = filterhigh.step(emg2);
-    emg2notchfilter = filternotch.step(emg2highfilter);
-    emg2abs = fabs(emg2notchfilter);
-    emg2lowfilter = filterlow.step(emg2abs);
+// Finding max values for correct motor switch
+void get_max1(){
+    if (max_reader1==0){
+            green = !green;
+            red = 1;
+            blue = 1;
+            for(int n=0;n<2000;n++){
+                        
+            emg1 = emg1_in.read();
+            emg1highfilter = filterhigh1.step(emg1);
+            emg1notchfilter = filternotch1.step(emg1highfilter);
+            emg1abs = fabs(emg1notchfilter);
+            emg1lowfilter = filterlow1.step(emg1abs);
+            emg1peak = filterpeak1.step(emg1lowfilter);
+            
+            if (max1<emg1peak){
+                max1 = emg1peak;
+            }
+            wait(0.001f);    
+            }
+            wait(0.2f);
+            green = 1;
+    }
+    maxpart1 = 0.35*max1;
+}
 
-    /* Set the sampled emg values in channel 0 (the first channel) and 1 (the second channel) in the 'HIDScope' instance named 'scope' */
-    scope.set(0, emg1_in.read() );
-    scope.set(1, emg1lowfilter );
-    scope.set(2, emgpeak );
-    /*  Finally, send all channels to the PC at once */
-    scope.send();
-    /* To indicate that the function is working, the LED is toggled */
-    green = !green;
+void get_max3(){
+    if (max_reader3==0){
+            green = 1;
+            blue = 1;
+            red = !red;
+            for(int n=0;n<2000;n++){
+                        
+            emg3 = emg3_in.read();
+            emg3highfilter = filterhigh3.step(emg3);
+            emg3notchfilter = filternotch3.step(emg3highfilter);
+            emg3abs = fabs(emg3notchfilter);
+            emg3lowfilter = filterlow3.step(emg3abs);
+            emg3peak = filterpeak3.step(emg3lowfilter);
+            
+            if (max3<emg3peak){
+                max3 = emg3peak;
+            }
+            wait(0.001f);    
+            }
+            wait(0.2f);
+            red = 1;
+    }
+    maxpart3 = 0.3*max3;
 }
 
-int main()
-{   
+// Filtering & Scope
+void filter() {
+    // Right Biceps
+    emg1 = emg1_in.read();
+    emg1highfilter = filterhigh1.step(emg1);
+    emg1notchfilter = filternotch1.step(emg1highfilter);
+    emg1abs = fabs(emg1notchfilter);
+    emg1lowfilter = filterlow1.step(emg1abs);
+    emg1peak = filterpeak1.step(emg1lowfilter);
+    // Left Biceps
+    emg2 = emg2_in.read();
+    emg2highfilter = filterhigh2.step(emg2);
+    emg2notchfilter = filternotch2.step(emg2highfilter);
+    emg2abs = fabs(emg2notchfilter);
+    emg2lowfilter = filterlow2.step(emg2abs);
+    emg2peak = filterpeak2.step(emg2lowfilter);
+    // Left Lower Arm OR Triceps
+    emg3 = emg3_in.read();
+    emg3highfilter = filterhigh3.step(emg3);
+    emg3notchfilter = filternotch3.step(emg3highfilter);
+    emg3abs = fabs(emg3notchfilter);
+    emg3lowfilter = filterlow3.step(emg3abs);
+    emg3peak = filterpeak3.step(emg3lowfilter);
+    
+    
+    /* Compare measurement to the calibrated value to decide actions */
+    if (maxpart1<emg1peak){
+        red = 0;
+        blue = 1;
+        green = 1;
+        }
+    else {
+        if (maxpart1<emg2peak){
+            red = 1;
+            blue = 0;
+            green = 1;
+            }
+        
+    else {
+        if (maxpart3<emg2peak){
+            red = 1;
+            blue = 1;
+            green = 0;
+            }
+        
+    else {
+        red = 1;
+        blue = 1;
+        green = 1;
+        }
+        }
+        }
+    /* Set the sampled emg values in channel 0 (the first channel) and 1 (the second channel) in the 'HIDScope' instance named 'scope' */
+    scope.set(0, emg1peak );
+    scope.set(1, emg2peak );
+    scope.set(2, maxpart1 );
+    scope.set(3, emg3peak );
+    scope.set(4, maxpart3 );    
+ 
+    scope.send();
+}
 
-    sample_timer.attach(&filter, 0.001);
+int main(){   
 
+    main_timer.attach(&filter, 0.001);
+    max_read1.attach(&get_max1, 2);
+    max_read3.attach(&get_max3, 2);
     while(1) {}
 }
\ No newline at end of file