Good Jacobian and code Not been tested

Dependencies:   MODSERIAL biquadFilter mbed

Fork of Kinematics by Ramon Waninge

main.cpp

Committer:
Ramonwaninge
Date:
2018-10-29
Revision:
4:49dfbfcd3577
Parent:
3:de8d3ca44a3e
Child:
5:d78ed3a3e66a

File content as of revision 4:49dfbfcd3577:

#include "mbed.h"
#include <math.h>
#include <cmath>
#include "MODSERIAL.h"
#define PI 3.14159265

MODSERIAL pc(USBTX, USBRX);  // connecting to pc
DigitalIn button1(SW3);      
DigitalIn button2(SW2);     
DigitalOut LED(LED1);

//Joe dit zijn de inputsignalen
double theta1 = PI*0.5; double theta4 = PI*0.5; double emg1; double emg2; double emg3; 
// Joe dit zijn de outputsignalen
double x; double y; 

//Joe dit zijn de constantes 
double ll = 200.0; double lu = 170.0; double lb = 10.0; double le = 79.0; double xbase = 340;
//Joe dit zijn de Tickers

//forward kinematics, Check mathematica! Omdat mbed in paniek raakt met meerdere wortels, hebben we de vergelijking opgedeeld in 3 stukken
//First define the position equation of x
double xendsum = lb + xbase +ll*(cos(theta1) - cos(theta4));
double xendsqrt1 = 2*sqrt(-xbase*xbase/4 + lu*lu + ll*(xbase*(cos(theta1)+cos(theta4))/2) -ll*(1+ cos(theta1+theta4)))*(-sin(theta1)+sin(theta4));
double xendsqrt2 = sqrt(pow((-xbase/ll+cos(theta1)+cos(theta4)),2)+ pow(sin(theta1) - sin(theta4),2));
double xend = (xendsum + xendsqrt1/xendsqrt2)/2;
//Now define the pos. eq. of y
double yendsum = -le + ll/2*(sin(theta1)+sin(theta4));
double yendsqrt1 = (-xbase/ll + cos(theta1)+cos(theta4))*sqrt(-xbase*xbase/4 + lu*lu + ll/2*(xbase*(cos(theta1)+cos(theta4))- ll*(1+cos(theta1+theta4)))); 
double yendsqrt2 = sqrt(pow((-xbase/ll + cos(theta1)+ cos(theta4)),2)+ pow((sin(theta1)-sin(theta4)),2));
double yend = (yendsum + yendsqrt1/yendsqrt2);


//Hier definieren we de functies
Ticker emgcheck;


//Joe, hieronder staan de functies die door de tickers aangeroepen worden
void xcor(){
    
    if (button1 == 0){                          //als emg1==voorbij treshold, 
        theta1 = PI*(theta1/PI + 0.1);          //double theta1-> plus een paar counts (emg*richting)
                                                //double theta4-> plus een paar counts (emg*richting)
    //default = als x = xbase/2... break, okee dit moet hier niet
    }
    else {theta1 = theta1;}
    xendsum = lb + xbase +ll*(cos(theta1) - cos(theta4));
    xendsqrt1 = 2*sqrt(-xbase*xbase/4 + lu*lu + ll*(xbase*(cos(theta1)+cos(theta4))/2) -ll*(1+ cos(theta1+theta4)))*(-sin(theta1)+sin(theta4));
    xendsqrt2 = sqrt(pow((-xbase/ll+cos(theta1)+cos(theta4)),2)+ pow(sin(theta1) - sin(theta4),2));
    xend = (xendsum + xendsqrt1/xendsqrt2)/2;
}
void ycor(){
    //als emg2 == voorbij treshold,
    //double theta1 -> plus counts (emg*richting)
    //double theta4 -> plus counts (emg*richting)
    //reken y door
    //default = als y = default... break
    //end
    }
void flip(){}

int main()
{
    
    pc.baud(115200);
    //default = theta1 = theta4 = pi/2 
    emgcheck.attach(xcor, 0.1);
    while(true){
       pc.printf("%f %f \n", xend,yend);
       wait(1);
        }
    
    
    
}