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:
3:de8d3ca44a3e
Parent:
2:0a7a3c0c08d3
Child:
4:49dfbfcd3577

File content as of revision 3:de8d3ca44a3e:

#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)
    theta4 = PI*0.5;          //double theta4-> plus een paar counts (emg*richting)
    xend = (xendsum + xendsqrt1/xendsqrt2)/2;
    //default = als x = xbase/2... break
    //end
    }
}
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.001);
    emgcheck.attach(ycor, 0.001);
    emgcheck.attach(flip, 0.001);
    while(true){
        pc.printf("%f %f \n", xend,yend);
        wait(1);
        }
    
    
    
}