..
Dependencies: HIDScope MODSERIAL biquadFilter mbed
main.cpp@10:39ec51206c8b, 2018-10-19 (annotated)
- Committer:
- MarijkeZondag
- Date:
- Fri Oct 19 14:13:59 2018 +0000
- Revision:
- 10:39ec51206c8b
- Parent:
- 9:c722418997b5
- Child:
- 11:b95b0e9e1b89
EMG filter erin geprobeerd te bouwen
Who changed what in which revision?
User | Revision | Line number | New contents of line |
---|---|---|---|
vsluiter | 0:c8f15874531b | 1 | #include "mbed.h" |
vsluiter | 0:c8f15874531b | 2 | #include "MODSERIAL.h" |
MarijkeZondag | 10:39ec51206c8b | 3 | #include "BiQuad.h" |
MarijkeZondag | 10:39ec51206c8b | 4 | #include "HIDScope.h" |
MarijkeZondag | 10:39ec51206c8b | 5 | #include <math.h> |
vsluiter | 0:c8f15874531b | 6 | |
MarijkeZondag | 10:39ec51206c8b | 7 | AnalogIn emg0_in (A0); |
MarijkeZondag | 10:39ec51206c8b | 8 | AnalogIn emg1_in (A1); |
MarijkeZondag | 10:39ec51206c8b | 9 | AnalogIn emg2_in (A2); |
MarijkeZondag | 10:39ec51206c8b | 10 | |
MarijkeZondag | 9:c722418997b5 | 11 | DigitalIn button2 (D10); //Let op, is deze niet bezet? |
MarijkeZondag | 9:c722418997b5 | 12 | InterruptIn encoderA (D9); |
MarijkeZondag | 9:c722418997b5 | 13 | InterruptIn encoderB (D8); |
MarijkeZondag | 6:f4bbb73f3989 | 14 | |
MarijkeZondag | 9:c722418997b5 | 15 | DigitalOut directionpin1 (D4); |
MarijkeZondag | 9:c722418997b5 | 16 | DigitalOut directionpin2 (D7); |
MarijkeZondag | 9:c722418997b5 | 17 | PwmOut pwmpin1 (D5); |
MarijkeZondag | 9:c722418997b5 | 18 | PwmOut pwmpin2 (D6); |
MarijkeZondag | 9:c722418997b5 | 19 | |
MarijkeZondag | 6:f4bbb73f3989 | 20 | |
vsluiter | 0:c8f15874531b | 21 | MODSERIAL pc(USBTX, USBRX); |
vsluiter | 0:c8f15874531b | 22 | |
MarijkeZondag | 9:c722418997b5 | 23 | |
MarijkeZondag | 9:c722418997b5 | 24 | //Global variables |
MarijkeZondag | 8:895d941a5910 | 25 | int encoder = 0; |
MarijkeZondag | 10:39ec51206c8b | 26 | |
MarijkeZondag | 10:39ec51206c8b | 27 | //Biquad |
MarijkeZondag | 10:39ec51206c8b | 28 | BiQuadChain emg0band; |
MarijkeZondag | 10:39ec51206c8b | 29 | BiQuad emg0band1( 7.29441e-01, -1.89276e-08, -7.29450e-01, -1.64507e-01, -7.26543e-01 ); |
MarijkeZondag | 10:39ec51206c8b | 30 | BiQuad emg0band2( 1.00000e+00, 1.99999e+00, 9.99994e-01, 1.72349e+00, 7.79616e-01 ); |
MarijkeZondag | 10:39ec51206c8b | 31 | BiQuad emg0band3( 1.00000e+00, -1.99999e+00, 9.99994e-01, -1.93552e+00, 9.39358e-01 ); |
MarijkeZondag | 10:39ec51206c8b | 32 | |
MarijkeZondag | 10:39ec51206c8b | 33 | BiQuadChain emg1band; |
MarijkeZondag | 10:39ec51206c8b | 34 | BiQuad emg1band1( 7.29441e-01, -1.89276e-08, -7.29450e-01, -1.64507e-01, -7.26543e-01 ); |
MarijkeZondag | 10:39ec51206c8b | 35 | BiQuad emg1band2( 1.00000e+00, 1.99999e+00, 9.99994e-01, 1.72349e+00, 7.79616e-01 ); |
MarijkeZondag | 10:39ec51206c8b | 36 | BiQuad emg1band3( 1.00000e+00, -1.99999e+00, 9.99994e-01, -1.93552e+00, 9.39358e-01 ); |
MarijkeZondag | 10:39ec51206c8b | 37 | |
MarijkeZondag | 10:39ec51206c8b | 38 | BiQuadChain emg2band; |
MarijkeZondag | 10:39ec51206c8b | 39 | BiQuad emg2band1( 7.29441e-01, -1.89276e-08, -7.29450e-01, -1.64507e-01, -7.26543e-01 ); |
MarijkeZondag | 10:39ec51206c8b | 40 | BiQuad emg2band2( 1.00000e+00, 1.99999e+00, 9.99994e-01, 1.72349e+00, 7.79616e-01 ); |
MarijkeZondag | 10:39ec51206c8b | 41 | BiQuad emg2band3( 1.00000e+00, -1.99999e+00, 9.99994e-01, -1.93552e+00, 9.39358e-01 ); |
MarijkeZondag | 10:39ec51206c8b | 42 | |
MarijkeZondag | 10:39ec51206c8b | 43 | BiQuad notch1( 9.91104e-01, -1.60364e+00, 9.91104e-01, -1.60364e+00, 9.82207e-01 ); //Notch filter |
MarijkeZondag | 10:39ec51206c8b | 44 | |
MarijkeZondag | 10:39ec51206c8b | 45 | |
MarijkeZondag | 10:39ec51206c8b | 46 | //Tickers |
MarijkeZondag | 10:39ec51206c8b | 47 | Ticker emg_filter_tick; |
MarijkeZondag | 10:39ec51206c8b | 48 | Ticker Mov_av_tick; |
MarijkeZondag | 9:c722418997b5 | 49 | |
MarijkeZondag | 9:c722418997b5 | 50 | //Functions |
MarijkeZondag | 10:39ec51206c8b | 51 | void EMGFilter0() |
MarijkeZondag | 10:39ec51206c8b | 52 | { |
MarijkeZondag | 10:39ec51206c8b | 53 | double emg0 = emg0_in.read(); |
MarijkeZondag | 10:39ec51206c8b | 54 | double bandpass = emg0band.step(emg0); |
MarijkeZondag | 10:39ec51206c8b | 55 | double absolute = fabs(bandpass); |
MarijkeZondag | 10:39ec51206c8b | 56 | double notch0 = notch1.step(absolute); |
MarijkeZondag | 10:39ec51206c8b | 57 | } |
MarijkeZondag | 10:39ec51206c8b | 58 | |
MarijkeZondag | 10:39ec51206c8b | 59 | void EMGFilter1() |
MarijkeZondag | 10:39ec51206c8b | 60 | { |
MarijkeZondag | 10:39ec51206c8b | 61 | double emg1 = emg1_in.read(); |
MarijkeZondag | 10:39ec51206c8b | 62 | double bandpass = emg1band.step(emg1); |
MarijkeZondag | 10:39ec51206c8b | 63 | double absolute = fabs(bandpass); |
MarijkeZondag | 10:39ec51206c8b | 64 | double notch1 = notch1.step(absolute); |
MarijkeZondag | 10:39ec51206c8b | 65 | } |
MarijkeZondag | 10:39ec51206c8b | 66 | |
MarijkeZondag | 10:39ec51206c8b | 67 | void EMGFilter2() |
MarijkeZondag | 10:39ec51206c8b | 68 | { |
MarijkeZondag | 10:39ec51206c8b | 69 | double emg2 = emg2_in.read(); |
MarijkeZondag | 10:39ec51206c8b | 70 | double bandpass = emg2band.step(emg2); |
MarijkeZondag | 10:39ec51206c8b | 71 | double absolute = fabs(bandpass); |
MarijkeZondag | 10:39ec51206c8b | 72 | double notch2 = notch1.step(absolute); |
MarijkeZondag | 10:39ec51206c8b | 73 | } |
MarijkeZondag | 10:39ec51206c8b | 74 | |
MarijkeZondag | 8:895d941a5910 | 75 | void encoderA_rise() |
MarijkeZondag | 8:895d941a5910 | 76 | { |
MarijkeZondag | 8:895d941a5910 | 77 | if(encoderB==false) |
MarijkeZondag | 8:895d941a5910 | 78 | { |
MarijkeZondag | 8:895d941a5910 | 79 | encoder++; |
MarijkeZondag | 8:895d941a5910 | 80 | } |
MarijkeZondag | 8:895d941a5910 | 81 | else |
MarijkeZondag | 8:895d941a5910 | 82 | { |
MarijkeZondag | 8:895d941a5910 | 83 | encoder--; |
MarijkeZondag | 8:895d941a5910 | 84 | } |
MarijkeZondag | 8:895d941a5910 | 85 | } |
MarijkeZondag | 8:895d941a5910 | 86 | |
MarijkeZondag | 8:895d941a5910 | 87 | void encoderA_fall() |
MarijkeZondag | 8:895d941a5910 | 88 | { |
MarijkeZondag | 8:895d941a5910 | 89 | if(encoderB==true) |
MarijkeZondag | 8:895d941a5910 | 90 | { |
MarijkeZondag | 8:895d941a5910 | 91 | encoder++; |
MarijkeZondag | 8:895d941a5910 | 92 | } |
MarijkeZondag | 8:895d941a5910 | 93 | else |
MarijkeZondag | 8:895d941a5910 | 94 | { |
MarijkeZondag | 8:895d941a5910 | 95 | encoder--; |
MarijkeZondag | 8:895d941a5910 | 96 | } |
MarijkeZondag | 8:895d941a5910 | 97 | } |
MarijkeZondag | 8:895d941a5910 | 98 | |
MarijkeZondag | 8:895d941a5910 | 99 | void encoderB_rise() |
MarijkeZondag | 8:895d941a5910 | 100 | { |
MarijkeZondag | 8:895d941a5910 | 101 | if(encoderA==true) |
MarijkeZondag | 8:895d941a5910 | 102 | { |
MarijkeZondag | 8:895d941a5910 | 103 | encoder++; |
MarijkeZondag | 8:895d941a5910 | 104 | } |
MarijkeZondag | 8:895d941a5910 | 105 | else |
MarijkeZondag | 8:895d941a5910 | 106 | { |
MarijkeZondag | 8:895d941a5910 | 107 | encoder--; |
MarijkeZondag | 8:895d941a5910 | 108 | } |
MarijkeZondag | 8:895d941a5910 | 109 | } |
MarijkeZondag | 8:895d941a5910 | 110 | |
MarijkeZondag | 8:895d941a5910 | 111 | void encoderB_fall() |
MarijkeZondag | 8:895d941a5910 | 112 | { |
MarijkeZondag | 8:895d941a5910 | 113 | if(encoderA==false) |
MarijkeZondag | 8:895d941a5910 | 114 | { |
MarijkeZondag | 8:895d941a5910 | 115 | encoder++; |
MarijkeZondag | 8:895d941a5910 | 116 | } |
MarijkeZondag | 8:895d941a5910 | 117 | else |
MarijkeZondag | 8:895d941a5910 | 118 | { |
MarijkeZondag | 8:895d941a5910 | 119 | encoder--; |
MarijkeZondag | 8:895d941a5910 | 120 | } |
MarijkeZondag | 8:895d941a5910 | 121 | } |
MarijkeZondag | 8:895d941a5910 | 122 | |
MarijkeZondag | 8:895d941a5910 | 123 | |
MarijkeZondag | 9:c722418997b5 | 124 | // Main function start. |
MarijkeZondag | 8:895d941a5910 | 125 | |
vsluiter | 0:c8f15874531b | 126 | int main() |
MarijkeZondag | 3:a3ad58758696 | 127 | { |
vsluiter | 0:c8f15874531b | 128 | pc.baud(115200); |
MarijkeZondag | 6:f4bbb73f3989 | 129 | pc.printf("hello\n\r"); |
MarijkeZondag | 10:39ec51206c8b | 130 | |
MarijkeZondag | 10:39ec51206c8b | 131 | //EMG signaal filteren |
MarijkeZondag | 10:39ec51206c8b | 132 | |
MarijkeZondag | 10:39ec51206c8b | 133 | // Deze? of die hierondbqc.add( &bq1 ).add( &bq2 ).add( &bq3 ); |
MarijkeZondag | 10:39ec51206c8b | 134 | bqc = bq1 * bq2 * bq3; |
MarijkeZondag | 10:39ec51206c8b | 135 | emgSampleTicker.attach(&emgSample,0.01); //100 Hz |
MarijkeZondag | 10:39ec51206c8b | 136 | while(true){/*do not return from main()*/} |
MarijkeZondag | 10:39ec51206c8b | 137 | |
MarijkeZondag | 10:39ec51206c8b | 138 | |
MarijkeZondag | 6:f4bbb73f3989 | 139 | pwmpin1.period_us(60); //60 microseconds PWM period, 16.7 kHz |
MarijkeZondag | 6:f4bbb73f3989 | 140 | |
MarijkeZondag | 8:895d941a5910 | 141 | encoderA.rise(&encoderA_rise); |
MarijkeZondag | 8:895d941a5910 | 142 | encoderA.fall(&encoderA_fall); |
MarijkeZondag | 8:895d941a5910 | 143 | encoderB.rise(&encoderB_rise); |
MarijkeZondag | 8:895d941a5910 | 144 | encoderB.fall(&encoderB_fall); |
vsluiter | 0:c8f15874531b | 145 | |
MarijkeZondag | 2:2f4444f1504d | 146 | while (true) |
MarijkeZondag | 2:2f4444f1504d | 147 | { |
MarijkeZondag | 10:39ec51206c8b | 148 | //Motor aansturen en encoder uitlezen |
MarijkeZondag | 6:f4bbb73f3989 | 149 | float u1 = potmetervalue1; |
MarijkeZondag | 6:f4bbb73f3989 | 150 | float u2 = potmetervalue2; |
MarijkeZondag | 6:f4bbb73f3989 | 151 | |
MarijkeZondag | 6:f4bbb73f3989 | 152 | float m1 = ((u1*2.0f)-1.0f); |
MarijkeZondag | 6:f4bbb73f3989 | 153 | float m2 = ((u2*2.0f)-1.0f); |
MarijkeZondag | 7:f32005d13749 | 154 | |
MarijkeZondag | 7:f32005d13749 | 155 | pwmpin1 = fabs(m1*0.6f)+0.4f; //pwm duty cycle can only be positive, floating, 0.4f is "inefficiënt", dit tellen we erbij op, en keer 0.6 om te corrigeren voor de helling. |
MarijkeZondag | 7:f32005d13749 | 156 | directionpin1.write(m1>0); //Indien waar, motor draait rechtsom. Indien niet waar, motor draait linksom. |
MarijkeZondag | 7:f32005d13749 | 157 | wait(0.01f); //zodat de code niet oneindig doorgaat. |
MarijkeZondag | 7:f32005d13749 | 158 | pwmpin2 = fabs(m2*0.6f)+0.4f; |
MarijkeZondag | 8:895d941a5910 | 159 | directionpin2.write(m2>0); |
MarijkeZondag | 10:39ec51206c8b | 160 | |
MarijkeZondag | 10:39ec51206c8b | 161 | float encoderDegrees = float(encoder)*(360.0/8400.0); |
MarijkeZondag | 9:c722418997b5 | 162 | |
MarijkeZondag | 10:39ec51206c8b | 163 | pc.printf("Encoder count: %f \n\r",encoderDegrees); |
MarijkeZondag | 10:39ec51206c8b | 164 | |
vsluiter | 0:c8f15874531b | 165 | } |
vsluiter | 0:c8f15874531b | 166 | } |
MarijkeZondag | 4:fbea05c8fe2b | 167 |