Travaux Pratiques » Polytech Ge SP4abc 2021

#include <stdio.h>

#include <stdlib.h>

#include <math.h>

#include "trame.h"

#include "sfr32c87.h"

#define PI 3.14159

/*Declaration des variables*/

typedef struct {

float latitude;

float longitude;

} Position;

typedef struct{

Position rpos;

float vitmax;

} Zone;

Zone zones[] = {

{{44.7887762, -3.012}, 50},

{{44.7891220, -3.013}, 70},

};

Position pos_prec;

//Trames de tests ? modifier si n?cessaire.

char * trames[]= {"$GPGSV,3,2,10,15,03,077,,18,04,041,42,19,85,271,,20,08,214,*7C",

"$GPGSV,3,3,10,22,39,053,50,28,15,320,*7E",

"$GPRMC,141914.00,A,4545.6424,N,00306.6036,E,0.4,99.4,010206,,*0C",

"$GPGLL,4545.6424,N,00306.6036,E,141914.00,A*0E",

"$GPGGA,141914.00,4545.0000,N,00306.6036,E,1,05,3.4,499.3,M,,M,,*7D",

"$GPGSA,A,3,,03,,22,14,,01,,18,,,,3.9,3.4,1.9*39",

"$GPVTG,99.4,T,,M,0.4,N,0.7,K*57",

"$GPZDA,141914.00,01,02,2006,00,00*69",

0};

int trame_cmp(char* trame,char* type){ //Cette fonction renvoie 1 si la trame est de la forme "GPGGA".

int i = 0;

while (type[i]!= '\0'){

if (trame[i+1] != type[i]){

return 0;

}

i+=1;

}

return 1;

}

int decode_int(char c){ //Cette fonction renvoie la valeur decimale associ?e ? un caract?re donn? en param?tre compris entre '0' et '9' et renvoie -1 sinon.

if (c<'0' || c>'9'){

return -1;

}

return c - '0';

}

int decode_nombre(char * ch,int n){ //Renvoie la valeur d?cimale des n premiers caract?res de la cha?ne ch.

int res = 0;

int i;

for (i = 0;i < n;i++){

res *= 10;

res += decode_int(ch[i]);

}

return res;

}

float conv_lat(char * latitude){ //La fonction convertie la latitude en sexa en degr?.

float val = decode_nombre(latitude,4); //* bof

float res = 0;

int i;

for (i = 5; i < 9;i++){

val += (decode_int(latitude[i]))*pow(10,-i+4); //On replace chaque chiffre au bonne endroit (unit?, dizaine ...).

}

res = (int)(val/100); //****

res += (val - res*100)/60;

return res;

}

float conv_long(char * longitude){ //La fonction convertie la longitude en sexa en degr?.

float val = decode_nombre(longitude,5); //**

float res = 0;

int i;

for (i = 6; i < 10;i++){

val += (decode_int(longitude[i]))*pow(10,-i+5);

}

res = (int)(val/100);

res += (val - res*100)/60;

return res;}

float conversion(char * chaine){ //Cette fonction convertie soit une latitude soit une longitude en degr?.

int cpt = 0;

while (chaine[cpt]!='\0'){ //***, .

cpt++;

}

if (cpt == 9){

return conv_lat(chaine);

}

else if(cpt==10){

return conv_long(chaine);

}

else{

return 1000.0;

}

}

int decode_trame(char * trame, Position *p){

char longi[10]={'\0'}; //*** trame -> @ lati -> @

char lati[9];

int cpt = 0;

int i = 0;

int j = 0;

int y = 0;

if (trame_cmp(trame,"GPGGA")==1){

while (trame[i] != '\0'){

if (trame[i] == ','){ //On compte le nombre de virgule.

cpt += 1;

i += 1;

}

if (cpt == 2) { //On regarde si le nombre de virgule correspond ? 2, ie la latitude

lati[j] = trame[i]; //On stocke caract?re par caract?re pour isoler la latitude afin de la convertir.

j ++;

}

if (cpt == 4) { //On regarde si le nombre de virgule correspond ? 4, ie la longitude

longi[y] = trame[i];

y ++;

}

i += 1;

}

p->latitude = conv_lat(lati); //***************** vous avez de la chance &trame[15], trame+15

p->longitude = conv_long(longi); //On stocke la latitude et la longitude convertie dans structure position.

return 1;

}

return 0;

}

float calcule_distance(Position p_1, Position p_2){

float rayon = 6378.14; //Rayon de la Terre (km).

float DistAng = 0;

DistAng = (PI/180)*acos((sin(p_1.latitude)*sin(p_2.latitude)+cos(p_1.latitude)*cos(p_2.latitude)*cos(p_2.longitude-p_1.longitude))); //Calcul distance angulaire

return DistAng*rayon; //On renvoi le calcul de la distance : D (km) = DistanceAngulaire*R.

}

float calcule_vitesse(Position p_1, Position p_2){

return calcule_distance(p_1,p_2)*3600; //On calcul la vitesse (en km) avec la formule v=d/t avec t=1s.

}

int distance_a_la_plus_poche_zone(Position p, Zone r[], int nb_zones, float *d){

int index = 0;

int i;

*d = calcule_distance(p, r[0].rpos); //On initialise une distance pour pouvoir ensuite la comparer et chercher le minimum.

for (i = 1; i<nb_zones; i++){ //On parcours la table de zone dangereuse.

if (calcule_distance(p,r[i].rpos) < *d){

index = i;

*d = calcule_distance(p,r[i].rpos);

}

}

return index;

}

void traitement(char * trame){

Position pos;

int alarme;

float vitesse;

float seuil = 0.5;

int index;

float distance;

if (trame_cmp(trame,"GPGGA")==1){

if (decode_trame(trame,&pos)==1){

vitesse = calcule_vitesse(pos,pos_prec);

index = distance_a_la_plus_poche_zone(pos,zones,2,&distance);

if ((distance <= seuil) && (vitesse > zones[index].vitmax)){

alarme = 1;

}

alarme = 0;

}

pos_prec.latitude = pos.latitude;

pos_prec.longitude = pos.longitude;

}

}

int t;

t = decode_nombre("3549846578254",9);

t = decode_nombre("35498465745615",9);