script zover

Dependencies:   HIDScope MODSERIAL mbed

main.cpp

Committer:
wiesdat
Date:
2014-10-20
Revision:
6:6f0bc2e465b0
Parent:
5:9415e5be8235
Child:
7:5f1c0bcd3f69

File content as of revision 6:6f0bc2e465b0:

#include "mbed.h"
#include "MODSERIAL.h"
#include "HIDScope.h"
#include <iostream>

#define A1LP1    0.018180963222803
#define A0LP1     0.016544013176248
#define B1LP1    -1.718913340044714
#define B0LP1     0.753638316443765

#define A1HP1  -1.999801878951505
#define A0HP1   0.999801878951505
#define B1HP1   -1.971717601075000
#define B0HP1  0.972111984032897

#define A0N 0.99436777112
#define A1N -1.89139989664
#define A2N 0.99436777112
#define B1N 1.89070035439
#define B2N -0.988036

#define TSAMP 0.001

AnalogIn    emg0(PTB1);
MODSERIAL pc(USBTX,USBRX);
HIDScope scope(5);
Ticker timer;
float emg_value, ylp1, yhp1, yn;
volatile bool looptimerflag;
float ysum = 0, yave=0 ;


float readEMG()
{

    emg_value=emg0.read();
    return emg_value;
}

float hpfilter1(float emg_value)
{
    static float x1=0,y1=0,x2=0, y2=0,x;
    x = emg_value;
    yhp1 = x + A1HP1*x1 + A0HP1*x2 - B1HP1*y1 + B0HP1*y2;
    x2 = x1;
    x1 = x;
    y2 = y1;
    y1 = yhp1;
    return yhp1;
}

float lpfilter1(float yhp1)
{
    static float x1=0,y1=0,x2=0, y2=0,x;
    x = yhp1;
    ylp1 = A1LP1*x1-B1LP1*y1+A0LP1*x2-B0LP1*y2;
    x2 = x1;
    x1 = x;
    y2 = y1;
    y1 = ylp1;
    return ylp1;
}

float notchfilter(float ylp1)
{
    static float yn,x1=0,x2=0,y1=0,y2=0,x;
    x = ylp1;
    yn = A0N*x + A1N*x1+A2N*x2+B1N*y1+B2N*y2;
    x2 = x1;
    x1 = x;
    y2 = y1;
    y1 = yn;
    return yn;
}

void viewer()
{
    scope.set(0,emg_value);
    scope.set(1,yhp1);
    scope.set(2,ylp1);
    scope.set(3,yn);
    scope.set(4,yave);
    
    scope.send();
}

void setlooptimerflag(void)
{
    looptimerflag = true;
}

int main()
{
    pc.baud(115200);
    timer.attach(setlooptimerflag,TSAMP);
    int n = 0;
    while(1) {

        while(!looptimerflag) /* do nothing */;
        looptimerflag = false;
        emg_value =  readEMG();
        yhp1 = hpfilter1(emg_value);  //function hpfilter
        ylp1 = lpfilter1(yhp1);       //function filter
        yn = notchfilter(ylp1);
        ysum = ysum+yn;
        n++;
       
        if(n==100) {
            yave = ysum/100;
            ysum = 0;
            n = 0;
            }

         viewer();                  //functie hidscope
        

    }
}