##// END OF EJS Templates
ABS Monitoring changes
imanay -
r145:146
parent child
Show More
@@ -1,1106 +1,1140
1 1 /*
2 2 * Servidor.c
3 3 *
4 4 * Created on: Nov 3, 2009
5 5 * Author: Jose Francisco Quenta
6 6 *
7 7 * Se implementa:
8 8 * -Carga en memoria los apuntes contenidos en un archivo de experimentos: apunte0 -> GPIO
9 9 * -Cambio de apunte.
10 10 * -Lectura del estado actual del apunte y grabado del mismo en un archivo
11 11 *
12 12 * Modified by Iván Manay since Nov 2012
13 13 * -From UDP to TCP.
14 14 * -Use of a frame for TCP communications with the central control module.
15 15 */
16 16
17 17 #include <stdio.h>
18 18 #include <stdlib.h>
19 19 #include <string.h>
20 20 #include <unistd.h>
21 21 #include <errno.h>
22 22
23 23 #include <sys/types.h>
24 24 #include <sys/socket.h>
25 25 #include <netinet/in.h>
26 26 #include <arpa/inet.h>
27 27 #include <netdb.h>
28 28 #include <time.h>
29 29 #include <math.h>
30 30
31 31 #include "./Librerias/at91gpio.h"
32 32 //#include "./Librerias/Mensajes.h"
33 33 #include "./Librerias/at91adc.h"
34 34 //clock
35 35 #include "./Librerias/at91sysclock.h"
36 36
37 37 #define PUERTO_SERVIDOR 5500
38 38 #define TAM_BUFFER 1024
39 39
40 40 #define maskc_out_beam PC30+PC28+PC26+PC24+PC22+PC20 //MSB-UP-LSB MSB-DOWN-LSB //APUNTE
41 41 #define maskc_out_cntl PC4+PC5+PC6 //MONITORING CONTROL
42 42
43 43 #define mask_sel_canal PC4 //MONITORING CONTROL
44 44 #define mask_sel_atenuacion PC5
45 45 #define mask_sel_calibracion PC6
46 46
47 47 #define maskb_in PB16+PB18+PB20+PB30+PB24+PB22 //MSB-UP-LSB MSB-DOWN-LSB //VERIFICACION
48 48
49 49 #define maskc_in_adc PC0+PC1
50 50
51 51 #define bit_up_2 0x00010000 //Mascara de cada bit a revisar: bit_up_2 es MSB
52 52 #define bit_up_1 0x00040000
53 53 #define bit_up_0 0x00100000
54 54 #define bit_dow_2 0x40000000
55 55 #define bit_dow_1 0x01000000
56 56 #define bit_dow_0 0x00400000
57 57
58 58 #define MyID 11
59 59 #define MAXPENDING 5 /* Maximum outstanding connection requests */
60 60
61 61 //parameters for the name of the output file
62 62 #define FPRE "AD" //prefix for the output file name
63 63 #define FEXT ".out" //file extension for the output file
64 64 #define FNAMELEN 40
65 65
66 66 //ADC parameters
67 67 #define REP 1 //defines how many times the data acquisation loop is repeated
68 68 #define NSAMPLES 10 //defines how many samples are taken in one data acqu-
69 // isation loop
69 // isation loop
70
71 #define NLOOPS 200 // Number of loops for the transmisor checking
70 72 #define CNVTIME 14.3 //defines how long it takes to get one sample. Value
71 73 // is only needed for the output file, doesn't change
72 74 // any ADC configurations
73 75 #define UREF 3.3 //Reference Voltage of ADC (max. ADC Voltage)
74 76 #define ADCRES 1023 //Resolution of ADC (10bit=1023)
75 77
76 78
77 79 char *buff_experimento= NULL;
78 80
79 81 AT91S_PIO *pioc;
80 82 AT91S_PIO *piob;
81 83 AT91S_PIO *pio_adc;
82 84 AT91S_ADC *padc;
83 85 AT91S_ADC *padd;
84 86
85 87 struct control_module_parameters {
86 88 char ID[20];
87 89 char param2[20];
88 90 char param3[20];
89 91 char param4[20];
90 92 };
91 93
92 94 typedef struct control_module_parameters cmp;
93 95
94 96 char *header = NULL;
95 97 char *TypeOfInstrument = NULL;
96 98 char *iDSource = NULL;
97 99 char *iDDestino = NULL;
98 100 char *rx_len = NULL;
99 101 char *cmd = NULL;
100 102 char *rx_data = NULL;
101 103 char *crc = NULL;
102 104
103 105 cmp parameters;
104 106 /*
105 107 * Zona de declaracion de cabeceras.
106 108 */
107 109 cmp inicializa_modulo(cmp p);
108 110 int inicializa_ethernet();
109 111 int rxData(int, char*);
110 112 void txData(int, char*);
111 113 void inicializa_gpio();
112 114 void inicializa_adc();
113 115 void procesa_peticion(char *rx_buffer, char *tx_buffer);
114 116 int cambia_apuntamiento(char *puntero_char);
115 117 int carga_experimento(char *nombre_archivo);
116 118 char *chequeo_sistema(char *filename, char *numero_muestras);
117 119 void recibe_experimento(char *data, char filename[]);
118 120 void SplitFrame(char *frame);
119 121 void intToStr( int number, char* str );
120 122
121 123 //ABS monitoring
122 124 //int ABS_monitoreo(int sel_atenuador, int sel_calibracion, float umbral, int pulsewidth);
123 125 char* ABS_monitoreo(int sel_atenuador, int sel_calibracion, float umbral, int pulsewidth);
124 126 AT91S_ADC * configADC1(void);
125 127 AT91S_ADC * configADC2(void);
126 128
127 129 FILE * create_Output(char*, time_t);
128 130
129 131 void writeOutput(float resultado, FILE * output);
130 132
131 133 int checkTx(long int results1[],long int results2[], float umbral, int pulsewidth);
132 134
133 135 double mediana(long int *results, unsigned int cuenta);
134 136 float getPhase(long int results1[], long int results2[]);
135 137
136 138 int fExists(char *);
137 139 int configCLK();
138 140 //
139 141
140 142 /*
141 143 *
142 144 */
143 145 int main(){
144 146
145 147 int servSocket;
146 148 int clntSocket;
147 149
148 150
149 151 char *rx_buffer = (char *) malloc(TAM_BUFFER);
150 152 char *tx_buffer = (char *) malloc(TAM_BUFFER);
151 153 /* Inicializa parametros del modulo*/
152 154 parameters = inicializa_modulo(parameters);
153 155 printf("%s\n%s\n%s\n%s\n",parameters.ID, parameters.param2, parameters.param3, parameters.param4);
154 156 /* Inicializa red*/
155 157 servSocket = inicializa_ethernet();
156 158 /* Inicializamos el puerto GPIO del sistema embebido GSBC-9260S */
157 159 inicializa_gpio();
158 160 inicializa_adc();
159 161
160 162 while(1){
161 163 // Recepción TCP de petición
162 164 clntSocket = rxData(servSocket, rx_buffer);
163 165 //testpoint
164 166 printf("rx:%s\n",rx_buffer);
165 167 // Procesamiento de la petición
166 168 procesa_peticion(rx_buffer, tx_buffer);
167 169 //testpoint
168 170 printf("tx:%s\n",tx_buffer);
169 171 // Respuesta del modulo de control
170 172 txData(clntSocket, tx_buffer);
171 173
172 174 }
173 175 }
174 176
175 177
176 178 int inicializa_ethernet(){
177 179
178 180 struct sockaddr_in inf_servidor;
179 181
180 182 int servSocket;
181 183
182 184 int resultado;
183 185
184 186 /* Haciendo la estructura local*/
185 187 memset(&inf_servidor, 0, sizeof(inf_servidor));
186 188 inf_servidor.sin_family= AF_INET;
187 189 inf_servidor.sin_port= htons(PUERTO_SERVIDOR);
188 190 inf_servidor.sin_addr.s_addr= INADDR_ANY;
189 191
190 192 /* Se establece el socket */
191 193 servSocket = socket(AF_INET,SOCK_STREAM, IPPROTO_TCP);
192 194 if (servSocket == -1){
193 195 printf("No se establecio correctamente el socket: socket()\n");
194 196 //ERROR_FATAL("No se establecio correctamente el socket: socket()\n");
195 197 exit(-1);
196 198 }
197 199
198 200 /* Se asocia el socket a un puerto y una IP */
199 201 resultado = bind(servSocket,(struct sockaddr *)&inf_servidor,sizeof(inf_servidor));
200 202 if (resultado== -1){
201 203 printf("No se establecio correctamente el socket: bind()\n");
202 204 //ERROR_FATAL("No se establecio correctamente el socket: bind()\n");
203 205 exit(-1);
204 206 }
205 207
206 208 if (listen(servSocket, MAXPENDING) < 0){
207 209 printf("listen() failed\n");
208 210 exit(-1);
209 211 }
210 212
211 213 return servSocket;
212 214
213 215 }
214 216
215 217 int rxData(int servSocket, char* rx_buffer){
216 218
217 219 int clntSocket;
218 220 struct sockaddr_in inf_cliente;
219 221 int numbytes_recibidos;
220 222 unsigned int inf_client_Len;
221 223
222 224 printf("\nEsperando solicitud de cliente...\n");
223 225
224 226 /* Set the size of the in-out parameter */
225 227 inf_client_Len = sizeof(inf_cliente);
226 228 /* Se espera hasta que un cliente se conecte */
227 229 if ((clntSocket = accept(servSocket, (struct sockaddr *) &inf_cliente,
228 230 &inf_client_Len)) < 0)
229 231 printf("accept() failed\n");
230 232
231 233 if ((numbytes_recibidos = recv(clntSocket, rx_buffer, TAM_BUFFER, 0)) < 0)
232 234 printf("recv() failed\n");
233 235
234 236 /* Se procede a procesar los datos recibidos */
235 237 rx_buffer[numbytes_recibidos]= '\0';
236 238
237 239 return clntSocket;
238 240 }
239 241
240 242 void txData(int clntSocket, char* data){
241 243
242 244 /* Echo message back to client */
243 245 if (send(clntSocket, data, strlen(data), 0) != strlen(data))
244 246 printf("send() failed\n");
245 247
246 248 close(clntSocket); /* Close client socket */
247 249 }
248 250
249 251 /*
250 252 * Esta funcion incializa el puerto GPIO
251 253 */
252 254 void inicializa_gpio(){
253 255
254 256 int desplazamiento_beam[6]={30,28,26,24,22,20}; //Desplazamientos para los pines de apunte
255 257 int desplazamiento_cntl[3]={4,5,6}; //Desplazamientos para los pines de control de monitoreo
256 258 unsigned long acumulado_ceros=0;
257 259 int i;
258 260
259 261 // Configuracion de los pines de APUNTE y control de monitoreo
260 262 pioc = pio_map(PIOC_BASE);
261 263 pio_enable(pioc, maskc_out_beam + maskc_out_cntl);
262 264 pio_disable_irq(pioc, maskc_out_beam + maskc_out_cntl);
263 265 pio_disable_multiple_driver(pioc, maskc_out_beam + maskc_out_cntl);
264 266 pio_disable_pull_ups(pioc, maskc_out_beam + maskc_out_cntl);
265 267 pio_synchronous_data_output(pioc, maskc_out_beam + maskc_out_cntl);
266 268 pio_output_enable(pioc, maskc_out_beam + maskc_out_cntl);
267 269
268 270 // Configuracion de los pines de VERIFICACION
269 271 piob = pio_map(PIOB_BASE);
270 272 pio_enable(piob, maskb_in);
271 273 pio_disable_irq(piob, maskb_in);
272 274 pio_disable_multiple_driver(piob, maskb_in);
273 275 pio_disable_pull_ups(piob, maskb_in);
274 276 pio_input_enable(piob, maskb_in);
275 277
276 278 //Pînes de ADC
277 279 pio_adc = pio_map(PIOC_BASE);
278 280 pin_adc_enable(pio_adc,maskc_in_adc); //Habilitamos PC0 para usar con ADC0 y 1
279 281 pio_disable_irq(pio_adc, maskc_in_adc);
280 282 pio_disable_multiple_driver(pio_adc, maskc_in_adc);
281 283 pio_disable_pull_ups(pio_adc, maskc_in_adc);
282 284 pio_input_enable(pio_adc, maskc_in_adc);
283 285
284 286 //Inicializando a cero los pines de apunte
285 287 for(i=0;i<6;i++)
286 288 acumulado_ceros= acumulado_ceros + (1 << desplazamiento_beam[i]);
287 289 pio_out(pioc, maskc_out_beam, acumulado_ceros, 0);
288 290
289 291 //Inicializando a cero los pines de control de monitoreo: atenuacion, calibracion y canal
290 292 for(i=0;i<3;i++)
291 293 acumulado_ceros= acumulado_ceros + (1 << desplazamiento_cntl[i]);
292 294 pio_out(pioc, maskc_out_cntl, acumulado_ceros, 0);
293 295
294 296 }
295 297
296 298 void inicializa_adc(){
297 299
298 300 if (configCLK() == 1)
299 301 printf("clock ADC enable.\n");
300 302 //configure ADC Settings
301 303 padc=configADC1();
302 304 padd=configADC2();
303 305 }
304 306
305 307
306 308 /*
307 309 * Divide rx frame into the frame components
308 310 */
309 311 void SplitFrame(char *frame){
310 312 header = malloc(4);
311 313 *header = *frame;
312 314 *(header + 1) = *(frame + 1);
313 315 *(header + 2) = *(frame + 2);
314 316 *(header + 3) = '\0';
315 317 TypeOfInstrument = malloc(4);
316 318 *TypeOfInstrument = *(frame + 3);
317 319 *(TypeOfInstrument + 1) = *(frame + 4);
318 320 *(TypeOfInstrument + 2) = *(frame + 5);
319 321 *(TypeOfInstrument + 3) = '\0';
320 322 iDSource = malloc(8);
321 323 *iDSource = *(frame + 6);
322 324 *(iDSource + 1) = *(frame + 7);
323 325 *(iDSource + 2) = *(frame + 8);
324 326 *(iDSource + 3) = *(frame + 9);
325 327 *(iDSource + 4) = *(frame + 10);
326 328 *(iDSource + 5) = *(frame + 11);
327 329 *(iDSource + 6) = *(frame + 12);
328 330 *(iDSource + 7) = '\0';
329 331 iDDestino = malloc(8);
330 332 *iDDestino = *(frame + 13);
331 333 *(iDDestino + 1) = *(frame + 14);
332 334 *(iDDestino + 2) = *(frame + 15);
333 335 *(iDDestino + 3) = *(frame + 16);
334 336 *(iDDestino + 4) = *(frame + 17);
335 337 *(iDDestino + 5) = *(frame + 18);
336 338 *(iDDestino + 6) = *(frame + 19);
337 339 *(iDDestino + 7) = '\0';
338 340 rx_len = malloc(7);
339 341 *rx_len = *(frame + 20);
340 342 *(rx_len + 1) = *(frame + 21);
341 343 *(rx_len + 2) = *(frame + 22);
342 344 *(rx_len + 3) = *(frame + 23);
343 345 *(rx_len + 4) = *(frame + 24);
344 346 *(rx_len + 5) = *(frame + 25);
345 347 *(rx_len + 6) = '\0';
346 348 cmd = malloc(5);
347 349 *cmd = *(frame + 26);
348 350 *(cmd + 1) = *(frame + 27);
349 351 *(cmd + 2) = *(frame + 28);
350 352 *(cmd + 3) = *(frame + 29);
351 353 *(cmd + 4) = '\0';
352 354
353 355 int l = atoi(rx_len);
354 356 rx_data = malloc(l + 1);
355 357 int i;
356 358 for (i = 30; i < 30 + l; i++)
357 359 *(rx_data + (i-30)) = *(frame + i);
358 360 *(rx_data + l) = '\0';
359 361 crc = malloc(2);
360 362 *crc = *(frame + 30 + l);
361 363 *(crc + 1) = '\0';
362 364 }
363 365
364 366 /*
365 367 * Esta funcion procesa el mensaje de peticion y genera respuesta
366 368 */
367 369 void procesa_peticion(char *rx_buffer, char *tx_buffer){
368 370 // int n = 0;
369 371 char filename1[50];
370 372 char filename2[] = "verificacion.txt";
371 373 char *tx_data = NULL;
372 374 char *tx_len = NULL;
373 375 SplitFrame(rx_buffer);
374 376
375 377 if ((cmd == NULL) || (rx_data == NULL)){
376 378 printf("procesarPeticion: formato de mensaje incorrecto");
377 379 //ERROR("procesarPeticion: formato de mensaje incorrecto");
378 380
379 381 }
380 382 else{
381 383 if(strcmp(cmd,"SNDF") == 0){
382 384 recibe_experimento(rx_data,filename1);
383 385 carga_experimento(filename1);
384 386 cambia_apuntamiento("0");
385 387 tx_data = (char*)malloc(3);
386 388 tx_data = "OK";
387 389 }
388 390 else if(strcmp(cmd,"CHGB") == 0){
389 391 cambia_apuntamiento(rx_data);
390 392 tx_data = (char*)malloc(3);
391 393 tx_data = "OK";
392 394 }
393 395 else if(strcmp(cmd,"ANST") == 0){
394 396 tx_data = chequeo_sistema(filename2,rx_data);
395 397 printf("%s\n",tx_data);
396 398 }
397 else if(strcmp(cmd,"ANPH") == 0){
398 tx_data = ABS_monitoreo(1, 1, 50, 10);
399 else if(strcmp(cmd,"BGPH") == 0){
400 tx_data = ABS_monitoreo(1, 0, 50, 10);
401 // tx_data = "Not implemented\n";
402 printf("%s\n",tx_data);
403 }
404 else if(strcmp(cmd,"LWPH") == 0){
405 tx_data = ABS_monitoreo(0, 0, 50, 10);
399 406 // tx_data = "Not implemented\n";
400 407 printf("%s\n",tx_data);
401 408 }
402 409 else if(strcmp(cmd,"NTST") == 0){
403 410 tx_data = malloc(strlen(parameters.ID) + 1);
404 411 strcpy(tx_data,parameters.ID);
405 412 printf("%s\n",tx_data);
406 413 }
407 414 else{
408 415 tx_data = (char*)malloc(6);
409 416 tx_data = "Error";
410 417 printf("procesa_peticion: comando no reconocido");
411 418 //ERROR("procesa_peticion: comando no reconocido");
412 419 }
413 420
414 421 tx_len = malloc(7);
415 422 int number = strlen(tx_data);
416 423 intToStr(number, tx_len );
417 424
418 425 strcpy(tx_buffer,header); //3
419 426 strcat(tx_buffer,TypeOfInstrument); //3
420 427 strcat(tx_buffer,parameters.ID); //7
421 428 strcat(tx_buffer,iDSource); //7
422 429 strcat(tx_buffer,tx_len); //6
423 430 strcat(tx_buffer,cmd); //4
424 431 strcat(tx_buffer,tx_data); //?
425 432 strcat(tx_buffer,crc); //1
426 433
427 434 }
428 435
429 436 }
430 437
431 438 /*
432 439 * Esta función genera el archivo de experimento a partir de la trama TCP recibida
433 440 */
434 441 void recibe_experimento(char *p_data, char filename[]){
435 442 FILE *fd;
436 443 int i = 0;
437 444
438 445 while (*p_data != '\n'){
439 446 filename[i] = *p_data;
440 447 i++;
441 448 p_data++;
442 449 }
443 450 filename[i] = '\0';
444 451 p_data = p_data - i;
445 452 fd = fopen(filename,"w");
446 453 fprintf(fd, p_data);
447 454 fclose(fd);
448 455 }
449 456
450 457 /*
451 458 * Esta funcion carga un archivo en un buffer que esta ubicado en memoria, luego
452 459 * este buffer es usado en la funcion "cambia_apuntamiento" para obtener el dato
453 460 * que sera usado en el cambio de apuntamiento.
454 461 */
455 462 int carga_experimento(char *nombre_archivo){
456 463
457 464 FILE *Archivo_Fd;
458 465
459 466 char *cadena = (char *) malloc(25);
460 467
461 468 int longitud_cadena;
462 469 int num_bytes= 0;
463 470 int num_filas= 0;
464 471
465 472 Archivo_Fd = fopen(nombre_archivo,"r"); // Se procede a abrir el archivo, segun la ruta especificada
466 473 if(!Archivo_Fd){
467 474 printf("carga_archivo: No se pudo abrir el archivo!!! --> fopen()\n");
468 475 //ERROR("carga_archivo: No se pudo abrir el archivo!!! --> fopen()\n");
469 476 return -1;
470 477 }else{
471 478
472 479 while(!feof(Archivo_Fd)){ // Se procede a calcular la longitud del archivo para separar memoria
473 480 fgets(cadena,20,Archivo_Fd);
474 481 longitud_cadena= strlen(cadena);
475 482 cadena[longitud_cadena-1] = '\0';
476 483 num_bytes = num_bytes + longitud_cadena;
477 484 num_filas++;
478 485 }
479 486
480 487 rewind(Archivo_Fd); // Se reinicia el puntero del archivo
481 488
482 489 char *buffer_temporal = (char *) malloc(num_bytes+1); // Se separa espacio de memoria segun
483 490 // la longitud del archivo
484 491 fread(buffer_temporal, sizeof(char), num_bytes, Archivo_Fd);
485 492
486 493 char *puntero= strstr(buffer_temporal,".ab1"); // Se procede a eliminar la cabecera del archivo
487 494 puntero= puntero + 12;
488 495
489 496 buff_experimento = (char *) malloc(7*(num_filas-3)); // num_bytes_fila*(num_filas-3);
490 497 strncpy(buff_experimento,puntero,7*(num_filas-3)); // Se carga en memoria la informacion del archivo
491 498
492 499 fclose(Archivo_Fd);
493 500
494 501 return 1;
495 502 }
496 503 }
497 504
498 505 /*
499 506 * Esta funcion recibe un numero en formato char, el dato se transforma a su equivalente en
500 507 * un numero entero, que sera usado para sacar un dato del buffer "buff_experimento", esta
501 508 * dato es el valor que se enviara al sistema de conmutacion RF para el cambio de apunte a
502 509 * traves del puerto GPIO.
503 510 */
504 511 int cambia_apuntamiento(char *puntero_char){
505 512
506 513 /*MSB-UP-LSB MSB-DOWN-LSB*/
507 514 int desplazamiento[6]={30,28,26,24,22,20}; // Defino los dezplazamientos que se aplicara
508 515 // al dato que ingresa para formar el número
509 516 // entero que se le pasara al puerto GPIO
510 517 // Estos números son los pines del puerto GPIO
511 518 // que se estan usando para el control
512 519
513 520 int puntero= atoi(puntero_char); // Se convierte a entero la direccion del puntero
514 521
515 522 int base= 7*puntero; // base= cantidad_bytes del dato x puntero
516 523 // cantidad de bytes es el numero de bytes que
517 524 printf("%s\n",puntero_char); // contiene cada dato, para este caso es 7
518 525 // porque es 6 bits de datos + 1 bit del cambio
519 526 // de linea.
520 527 char valor_char;
521 528 unsigned long valor;
522 529 unsigned long acumulado_ceros=0;
523 530 unsigned long acumulado_unos=0;
524 531
525 532 int offset; // Defino offset para el desplazamiento a traves
526 533 for(offset=0;offset<6;offset++){ // de cada dato que se obtiene del "buff_experimento"
527 534
528 535 valor_char= buff_experimento[base+offset]; // Obtengo el dato
529 536
530 537 if (valor_char == '0'){ // Obtengo el número acumulado segun sea un cero o un uno
531 538 valor= 0;
532 539 acumulado_ceros= acumulado_ceros + (1 << desplazamiento[offset]);
533 540 }else{
534 541 valor= 1;
535 542 acumulado_unos= acumulado_unos + (1 << desplazamiento[offset]);
536 543 }
537 544 }
538 545 pio_out(pioc, maskc_out_beam, acumulado_unos, 1);
539 546 pio_out(pioc, maskc_out_beam, acumulado_ceros, 0);
540 547
541 548 return 1;
542 549
543 550 }
544 551
545 552 /*
546 553 * Esta funcion lee "n" veces el estado del APUNTE actual y reporta
547 554 * una cadena de Verificacion.
548 555 */
549 556 char* chequeo_sistema(char *filename, char *numero_muestras){
550 557
551 558 int i;
552 559 int cnt = 0;
553 560 unsigned int entradac= 0;
554 561
555 562 char page0[250];
556 563
557 564 /*strcpy(page0,"Verificacion\n");
558 565 strcat(page0,parameters.ID);*/
559 566 strcpy(page0,parameters.ID);
560 567 strcat(page0,"\n-------\n");
561 568
562 569 char page1[8];
563 570
564 571 do{
565 572 //Inicializando arreglo
566 573 for(i=0;i<6;i++)
567 574 page1[i]='0';
568 575 page1[6] = '\n';
569 576 page1[7] = '\0';
570 577 //Lectura de puerto
571 578 entradac= pio_in(piob,maskb_in);
572 579 //Dandole formato al dato
573 580 if (!(entradac & bit_up_2))
574 581 page1[0] = '1';
575 582 if (!(entradac & bit_up_1))
576 583 page1[1] = '1';
577 584 if (!(entradac & bit_up_0))
578 585 page1[2] = '1';
579 586 if (!(entradac & bit_dow_2))
580 587 page1[3] = '1';
581 588 if (!(entradac & bit_dow_1))
582 589 page1[4] = '1';
583 590 if (!(entradac & bit_dow_0))
584 591 page1[5] = '1';
585 592
586 593 strcat(page0, page1);
587 594 cnt=cnt+1;
588 595 usleep(1*1000*1000);
589 596
590 597 }while(cnt < atoi(numero_muestras));
591 598
592 599 page0[strlen(page0)] = '\0';
593 600
594 601 char *all_pages = malloc(strlen(page0)+1);
595 602 strcpy(all_pages, page0);
596 603 return all_pages;
597 604 }
598 605
599 606 /*
600 607 *
601 608 */
602 609 cmp inicializa_modulo(cmp p){
603 610 FILE *fd = fopen("configuration.txt","r");
604 611 fgets(p.ID,20,fd);
605 612 p.ID[7]='\0';
606 613 fgets(p.param2,20,fd);
607 614 p.param2[10]='\0';
608 615 fgets(p.param3,20,fd);
609 616 p.param3[10]='\0';
610 617 fgets(p.param4,20,fd);
611 618 p.param4[10]='\0';
612 619 fclose(fd);
613 620 return p;
614 621 }
615 622
616 623 /*
617 624 *
618 625 */
619 626 void intToStr( int number, char* str )
620 627
621 628 {
622 629 int index = 0;
623 630
624 631 while( number > 0 )
625 632 {
626 633 int digit = number%10;
627 634 str[index++] = digit + '0';
628 635 number /= 10;
629 636 }
630 637 str[index] = '\0';
631 638 //Adding zero to the left
632 639 int n= strlen(str);
633 640 if (n == 1) {
634 641 strcat(str,"00000");
635 642 index = index + 5;
636 643 }else if(n == 2){
637 644 strcat(str,"0000");
638 645 index = index + 4;
639 646 }else if(n == 3){
640 647 strcat(str,"000");
641 648 index = index + 3;
642 649 }else if(n == 4){
643 650 strcat(str,"00");
644 651 index = index + 2;
645 652 }else if(n == 5){
646 653 strcat(str,"0");
647 654 index = index + 1;
648 655 }
649 656 //Now reverse the numbers in the string.
650 657 int position;
651 658 for( position = 0; position <= (index-1)/2; ++position )
652 659 {
653 660 char tmp = str[position];
654 661 str[position] = str[(index-1)-position];
655 662 str[(index-1)-position] = tmp;
656 663 }
657 664 }
658 665
659 666
660 667 //*****************************************************************
661 668 //ABS_monitoreo es la funci�n principal del proyecto ABS_Monitoreo.
662 669 //Esta funci�n es la que se debe agregar en otros c�digos.
663 670 //*****************************************************************
664 671 char* ABS_monitoreo(int sel_atenuador, int sel_calibracion, float umbral, int pulsewidth){
665 672
666 673 //local variables
667 674 /* AT91S_PIO *pioc;
668 675 pioc = pio_map(PIOC_BASE);
669 676 unsigned int mask_sel_canal =PC4; //Aqu� se indican los pines que se desean usar como salidas. Las constantes PCx est�n defiidas en el header at91gpio.h
670 677 unsigned int mask_sel_atenuacion =PC5;
671 678 unsigned int mask_sel_calibracion =PC6;*/
672 679
673 680 unsigned long acumulado_ceros=0;
674 681 unsigned long acumulado_unos=0;
675 682 /* AT91S_ADC *padc;
676 683 AT91S_ADC *padd;*/
677 684 FILE *fp;
678 685 long int results1[NSAMPLES], results2[NSAMPLES], results3[NSAMPLES], results4[NSAMPLES];
679 686 unsigned int i=0;
680 687 char fname[FNAMELEN];
681 688 int j=0;
682 689 time_t now;
683 690 FILE *archivo;
684 691 float phase1;
685 692 float phase2;
686 693 char page0[30];
687 694 char page1[20];
695
696 int cnt = 0;
688 697 //system("./map_clock");
689 698 /*
690 699 if (configCLK() == 1)
691 700 printf("clock ADC enable.\n");
692 701 */
693 702 /*
694 703 //configurar tres pines como salida usando als m�scaras mask_sel_canal, mask_sel_atenuacion y mask_sel_calibracion. En este caso corresponden a los pines pc4, pc5 y pc6.
695 704 pio_enable(pioc, mask_sel_canal);
696 705 pio_enable(pioc, mask_sel_atenuacion);
697 706 pio_enable(pioc, mask_sel_calibracion);
698 707 pio_output_enable(pioc, mask_sel_canal); //configurar pc4 como salida
699 708 pio_output_enable(pioc, mask_sel_atenuacion); //configurar pc5 como salida
700 709 pio_output_enable(pioc, mask_sel_calibracion); //configurar pc6 como salida
701 710 */
702 711
703 712 //Se modifican las salidas correspondientes a la selecci�n del atenuador y calibraci�n, de acuerdo a los par�metros ingresados en la funci�n ABS_monitoreo.
704 713 /*if ( sel_atenuador == 1)
705 714 pio_out(pioc, mask_sel_atenuacion, sel_atenuador,1);
706 715 else
707 716 pio_out(pioc, mask_sel_atenuacion, sel_atenuador,0);
708 717 if ( sel_calibracion == 1)
709 718 pio_out(pioc, mask_sel_calibracion, sel_calibracion,1);
710 719 else
711 720 pio_out(pioc, mask_sel_calibracion, sel_calibracion,0);*/
712 721
713 722 if ( sel_atenuador == 1)
714 723 acumulado_unos = acumulado_unos + (1 << 5);
715 724 else
716 725 acumulado_ceros = acumulado_ceros + (1 << 5);
717 726 if ( sel_calibracion == 1)
718 727 acumulado_unos = acumulado_unos + (1 << 6);
719 728 else
720 729 acumulado_ceros = acumulado_ceros + (1 << 6);
721 730
722 731
723 732 strcpy (fname, "/mnt/sd/absmonitoreo.txt"); //Direcci�n y nombre del archivo donde se desea guardar los datos.
724 733
725 734 if (fExists(fname)==0){ //si el archivo no existe, crea uno y le asigna el titulo
726 735 archivo = fopen(fname,"a+");
727 736 fprintf(archivo,"%s"," Registro de datos del ABS Control \n");
728 737 fprintf(archivo,"%s"," Fecha y hora Fase UP Fase DOWN\n");
729 738 fclose(archivo);
730 739 }
731 740
732 741 /*
733 742 //configure ADC Settings
734 743 padc=configADC1();
735 744 padd=configADC2();
736 745 */
737 746 //while (1){
738 747
739 748 ENABLE_CHANNEL(padc, ADC_CH0+ADC_CH1);
740 749 printf("\nAdquiriendo datos...\n"); //Indica en el terminal que se est�n adquiriendo datos (muestreando la se�al).
741 750
742 751
743 752 now = time(0); //Get current Time for File Name
744 753
745 754 //Se pone la salida de selecci�n de canal para seleccionar el canal 1 del detector de fase
746 755 acumulado_ceros = acumulado_ceros + (1 << 4);
747 756 pio_out(pioc, maskc_out_cntl, acumulado_ceros, 0);
748 757 pio_out(pioc, maskc_out_cntl, acumulado_unos, 1);
749 758 sleep(1);
750 759
751 760 //Se toman muestras para el canal 1 del detector de fase
752 //while(1){
761 /*while(1){
753 762 for(i=0; i < NSAMPLES; i++){
754 763
755 764 ADC_INIT(padc);
756 765 results1[i] = GET_ADC0(padc);
757 766 results2[i] = GET_ADC1(padd);
758 767 }
759 /*
768
760 769 if (checkTx(results1, results2, umbral, pulsewidth)==1){ //Se verifica que las muestras tomadas del canal 1 del datector de fase //correspondan a un pulso.
761 770 break;
762 771 }
763 */
764 //}
765
772
773 }*/
774
775 do{
776 for(i=0; i < NSAMPLES; i++){
777
778 ADC_INIT(padc);
779 results1[i] = GET_ADC0(padc);
780 results2[i] = GET_ADC1(padd);
781 }
782
783 if (checkTx(results1, results2, umbral, pulsewidth)==1){ //Se verifica que las muestras tomadas del canal 1 del datector de fase //correspondan a un pulso.
784 break;
785 }
786 cnt += 1;
787 }while (cnt < NLOOPS);
766 788
767 789 //Se pone la salida de selecci�n de canal para seleccionar el canal 2 del detector de fase
768 790 acumulado_ceros = acumulado_ceros - (1 << 4);
769 791 acumulado_unos = acumulado_unos + (1 << 4);
770 792 pio_out(pioc, maskc_out_cntl, acumulado_ceros, 0);
771 793 pio_out(pioc, maskc_out_cntl, acumulado_unos, 1);
772 794 sleep(1);
773 795
774 796
775 797 //Setoman muestras para el canal 2 del detector de fase
776 // while(1){
798 /* while(1){
777 799 for(i=0; i < NSAMPLES; i++){
778 800
779 801 ADC_INIT(padc);
780 802 results3[i] = GET_ADC0(padc);
781 803 results4[i] = GET_ADC1(padd);
782 804 }
783 /*
805
784 806 if (checkTx(results3, results4, umbral, pulsewidth)==1){ //Se verifica que las muestras tomadas del canal 2 del detector de fase //correspondan a un pulso.
785 807 break;
786 808 }
787 */
788 // }
789
809 }*/
810 cnt = 0;
811 do{
812 for(i=0; i < NSAMPLES; i++){
813
814 ADC_INIT(padc);
815 results3[i] = GET_ADC0(padc);
816 results4[i] = GET_ADC1(padd);
817 }
818
819 if (checkTx(results3, results4, umbral, pulsewidth)==1){ //Se verifica que las muestras tomadas del canal 2 del detector de fase //correspondan a un pulso.
820 break;
821 }
822 cnt += 1;
823 }while(cnt < NLOOPS);
790 824
791 825 //Una vez que se ha encontrado un pulso en cada canal, se calcula la fase de ambos.
792 826
793 827 phase1 = getPhase(results1, results2); //Calcular la fase del canal 1 del detector de fase.
794 828 phase2 = getPhase(results3, results4); //Calcular la fase del canal 2 del detector de fase.
795 829 //create Output File
796 830
797 831 strcpy (fname, "/mnt/sd/absmonitoreo.txt");
798 832 printf("\nTerminada la prueba # %d \n", j++);
799 833 fp=create_Output(fname, now); //Coloca la fecha y la hora en el archivo de texto
800 834 printf("mediana ch1 = %1.2f\n", phase1); //muestra resultado en terminal
801 835 printf("mediana ch2 = %1.2f\n", phase2); //muestra resultado en terminal
802 836 writeOutput(phase1, fp); //graba el resultado en el archivo de texto
803 837 writeOutput(phase2, fp); //graba el resultado en el archivo de texto
804 838 fprintf(fp, "\n"); //Pasa a la siguiente l�nea del archivo de texto
805 839 fclose(fp);
806 840 printf("Resultado guardado en %s \n", fname);
807 841
808 842 sleep(1);
809 843
810 844 strcpy(page0,parameters.ID);
811 845 strcat(page0,"\n-------\n");
812 846
813 847 sprintf(page1,"UP:%1.2f DW:%1.2f\n",phase1, phase2);
814 848 strcat(page0,page1);
815 849 char *all_pages = malloc(strlen(page0)+1);
816 850 strcpy(all_pages, page0);
817 851 return all_pages;
818 852
819 853 // }
820 854 // return 0;
821 855 }
822 856 /*=============================================================================
823 857 Function definitions
824 858 =============================================================================*/
825 859
826 860 // Configures ADC registers in order to get a sample every 10us
827 861 AT91S_ADC * configADC1(void){
828 862 //Variables a usar:
829 863 /*
830 864 unsigned int maskc_adc =PC0; //Usamos ADC0 y ADC1
831 865
832 866 //configuro pin:
833 867 AT91S_PIO *pioc;
834 868 pioc = pio_map(PIOC_BASE);
835 869 pin_adc_enable(pioc,maskc_adc); //Habilitamos PC0 para usar con ADC0 y 1
836 870 pio_disable_irq(pioc, maskc_adc);
837 871 pio_disable_multiple_driver(pioc, maskc_adc);
838 872 pio_disable_pull_ups(pioc, maskc_adc);
839 873 pio_input_enable(pioc, maskc_adc);
840 874 */
841 875
842 876 //Configuro el ADC:
843 877 //AT91S_ADC *padc;
844 878
845 879 padc = adc_map1(ADC_BASE);
846 880
847 881 //clock ADC = 1MHz
848 882 //time startup = 8us
849 883 //time sample and hold = 2us
850 884 // hold
851 885 // ___________
852 886 // start ___________| |___________
853 887 //
854 888 // | --1.2us-- | --0.15us-- |
855 889 //ADC_RESET(padc);
856 890 CONFIG_ADC(padc,ADC_TRGEN_DIS | ADC_RES_10BIT | ADC_SLEEP_NORMAL_MODE | ADC_PRESCAL | ADC_STARTUP | ADC_SHTIM);
857 891 ENABLE_CHANNEL(padc,ADC_CH0); //habilito canal 0
858 892 return padc;
859 893 }
860 894
861 895 AT91S_ADC * configADC2(void){
862 896 //Variables a usar:
863 897 /*
864 898 unsigned int maskc_adc =PC1; //Usamos ADC0 y ADC1
865 899
866 900 //configuro pin:
867 901 AT91S_PIO *piod;
868 902 piod = pio_map(PIOC_BASE);
869 903 pin_adc_enable(piod,maskc_adc); //Habilitamos PC1 para usar con ADC0 y 1
870 904 pio_disable_irq(piod, maskc_adc);
871 905 pio_disable_multiple_driver(piod, maskc_adc);
872 906 pio_disable_pull_ups(piod, maskc_adc);
873 907 pio_input_enable(piod, maskc_adc);
874 908 */
875 909 //Configuro el ADC:
876 910 //AT91S_ADC *padd;
877 911
878 912 padd = adc_map1(ADC_BASE);
879 913
880 914 //clock ADC = 1MHz
881 915 //time startup = 8us
882 916 //time sample and hold = 2us
883 917 // hold
884 918 // ___________
885 919 // start ___________| |___________
886 920 //
887 921 // | --1.2us-- | --0.15us-- |
888 922 //ADC_RESET(padc);
889 923 CONFIG_ADC(padd,ADC_TRGEN_DIS | ADC_RES_10BIT | ADC_SLEEP_NORMAL_MODE | ADC_PRESCAL | ADC_STARTUP | ADC_SHTIM);
890 924 ENABLE_CHANNEL(padd,ADC_CH1); //habilito canal 1
891 925 return padd;
892 926 }
893 927
894 928
895 929 //++++++++++++++++++++
896 930
897 931 //creats the output file with a timestamp in the name
898 932 FILE * create_Output(char *fname, time_t rawtime){
899 933 FILE *file;
900 934 char timestamp[80];//, counter[5]="dcv";
901 935 //char str[4];
902 936 struct tm * timeinfo;
903 937
904 938 //format time
905 939 timeinfo = localtime ( &rawtime );
906 940 strftime (timestamp,sizeof(timestamp),"%a %y-%m-%d %H:%M:%S %Z",timeinfo);
907 941
908 942
909 943 //Creates the file name out of the #define parameters
910 944
911 945 strcpy (fname, "/mnt/sd/absmonitoreo.txt");
912 946 file = fopen(fname,"a+");
913 947 fprintf(file,"%s", timestamp);
914 948 //printf("\nTerminada la prueba # %d. Guardando resultado en %s\n",r, fname);
915 949 //printf("\nTerminada la prueba # %d/%d. Writing data to the file %s\n",r+1 , REP, fname);
916 950 //printf("\nAAAAAAAAAA %d...%s\n", counter[1], fname);
917 951 // return file pointer
918 952 return file;
919 953 }
920 954
921 955 //++++++++++++++++++++
922 956
923 957 //tests if a file already exists. returns 1 if it exists and 0 if it doesn't
924 958
925 959
926 960
927 961 //Funci�n checkTx verifica que la se�al muestreada corresponda a un pulso.
928 962 //results1 y results2 son los arreglos que contienen los datos muestreados por ambos canales del ADC del embebido.
929 963 //umbral indica qu� valor debe superar una muestra para considerarla un posible pulso o pico.
930 964 //pulsewidth indica cu�ntas muestras consecutivas deben superar el umbral para que se considere que se ha detectado un pulso.
931 965 int checkTx(long int results1[],long int results2[], float umbral, int pulsewidth){
932 966
933 967 int i, cont;
934 968 float z[NSAMPLES], sum, avg;
935 969 int isSignal, pulse;
936 970
937 971 for(i=0;i<NSAMPLES;i++){
938 972
939 973 z[i] =sqrt(1.0*results1[i]*results1[i]+1.0*results2[i]*results2[i]);
940 974 }
941 975
942 976 pulse = 0;
943 977 isSignal = 0;
944 978 cont =0;
945 979
946 980 sum = 0;
947 981 for(i=0;i<NSAMPLES;i++){
948 982
949 983 sum += z[i];
950 984 avg = sum/(i+1);
951 985 if ((z[i] - avg) > umbral){
952 986 if (isSignal == 1){
953 987 cont += 1;
954 988 }
955 989 if (cont == pulsewidth){
956 990 pulse = 1;
957 991 break;
958 992 }
959 993 isSignal = 1;
960 994 continue;
961 995 isSignal = 0;
962 996 cont = 0;
963 997 }
964 998 }
965 999
966 1000 return pulse; //devuelve un entero: 1 si se ha detectado pulso, de lo contrario, 0.
967 1001 }
968 1002
969 1003
970 1004 int fExists(char * fname){
971 1005 FILE * file;
972 1006
973 1007 file = fopen (fname, "r");
974 1008 if (file == NULL)
975 1009 {
976 1010 return 0;
977 1011 }
978 1012 fclose(file);
979 1013 return 1;
980 1014 }
981 1015
982 1016
983 1017 //Funci�n que calcula la mediana de un conjunto de muestras
984 1018 double mediana(long int *results,unsigned int cuenta){
985 1019 unsigned int i=0,j=0,aux=0;
986 1020
987 1021 double median=0;
988 1022 /*Calculo mediana */
989 1023
990 1024 for(i=0;i<cuenta-1;i++){
991 1025 for (j=i+1;j<cuenta;j++){
992 1026 if(results[i]>results[j] ){
993 1027
994 1028 aux=results[i];
995 1029 results[i]=results[j];
996 1030 results[j]=aux;
997 1031
998 1032 }
999 1033 }
1000 1034
1001 1035 }
1002 1036 median=results[cuenta/2];
1003 1037 return median;
1004 1038 }
1005 1039
1006 1040
1007 1041
1008 1042 //Funci�n que halla la fase de la se�al.
1009 1043 //Tiene como entradas las muestras correspondientes a la parte real e imaginaria de la se�al.
1010 1044 float getPhase(long int results1[],long int results2[]){
1011 1045
1012 1046 unsigned int count=0, i=0,umbral=1000;
1013 1047 //long int results1[];
1014 1048 //long int results2[];
1015 1049 long int power[NSAMPLES];
1016 1050 long int sumI=0,sumQ=0,I[NSAMPLES], Q[NSAMPLES],II[NSAMPLES], QQ[NSAMPLES];
1017 1051 double median1=0,median2=0;
1018 1052 long int promedioI=0,promedioQ=0;/*Calculo mediana 1*/
1019 1053 float resultado=0;
1020 1054
1021 1055 for(i=0;i<NSAMPLES;i++){
1022 1056
1023 1057 I[i] =results1[i];
1024 1058 Q[i] =results2[i];
1025 1059 }
1026 1060
1027 1061 /*Calculo mediana 1*/
1028 1062 median1=mediana(I,NSAMPLES);
1029 1063
1030 1064 /*Calculo mediana 2*/
1031 1065 median2=mediana(Q,NSAMPLES);
1032 1066
1033 1067
1034 1068
1035 1069
1036 1070
1037 1071
1038 1072 for(i=0;i<NSAMPLES;i++){
1039 1073
1040 1074 I[i] =results1[i];
1041 1075 Q[i] =results2[i];
1042 1076
1043 1077 }
1044 1078
1045 1079
1046 1080
1047 1081 for(i = 0; i < NSAMPLES ; i++){
1048 1082
1049 1083 I[i]=(I[i]-median1);
1050 1084 Q[i]=(Q[i]-median2);
1051 1085
1052 1086 }
1053 1087
1054 1088 for(i = 0; i < NSAMPLES ; i++){
1055 1089
1056 1090 power[i]=I[i]*I[i]+Q[i]*Q[i];
1057 1091
1058 1092 if(power[i] > umbral)
1059 1093 {
1060 1094
1061 1095 II[count]=I[i];
1062 1096 QQ[count]=Q[i];
1063 1097 count=count+1;
1064 1098
1065 1099 }
1066 1100
1067 1101 }
1068 1102
1069 1103 for(i = 0; i < count ; i++){
1070 1104
1071 1105 sumI=sumI+II[i];
1072 1106 sumQ=sumQ+QQ[i];
1073 1107
1074 1108 }
1075 1109
1076 1110 promedioI=sumI;
1077 1111 promedioQ=sumQ;
1078 1112
1079 1113 resultado = atan2(1.0*promedioI,1.0*promedioQ)*180/3.1416+62-44;
1080 1114
1081 1115
1082 1116 return resultado;
1083 1117
1084 1118 }
1085 1119
1086 1120
1087 1121
1088 1122 //Funci�n que muestra la fase detectada en el terminal y tambi�n la graba en el archivo de texto.
1089 1123 void writeOutput(float resultado, FILE * output){
1090 1124
1091 1125
1092 1126 //
1093 1127
1094 1128 fprintf(output," %1.2f ",resultado); //graba resultado en archivo .txt
1095 1129 //
1096 1130
1097 1131 }
1098 1132
1099 1133 int configCLK(){
1100 1134 //configuro pin:
1101 1135 AT91S_PMC *sys_clock;
1102 1136 sys_clock = clock_map(CLOCK_BASE);
1103 1137 enable_clock_adc(sys_clock);
1104 1138 //printf("clock ADC enable.\n");
1105 1139 return 1;
1106 1140 }
General Comments 0
You need to be logged in to leave comments. Login now