nia.m

Matlab program to plot the data acquired from the USB (and placed into a TXT file) for the NIA devic

%cd "D:/julio"
clear all;
offset_TXT=13;

  fid = fopen('sample_nia2.txt','r');
 %fid = fopen('UsbSnoop1_48HZ_julio.log','r');

pasada=1;
 tline = fgetl(fid);
while ischar(tline)   
   %disp(tline);    

  % data = fscanf(fid,"%s" ,2 );
  cadena_texto=tline;
%texto=tline (8:9)    
  %cosa1=(hex2dec(  tline (8:9) )) * 1%256;
 % cosa2= (hex2dec(  tline (11:12) )) * 1%256;
 %cosa3= (hex2dec(  tline (14:15) )) *1%65535;
 %dec2hex ( cosa1);
% dec2hex ( cosa2);
% dec2hex ( cosa3)  ;
  puntero=1+offset_TXT;  
  paquetes_validos = hex2dec(tline(55));
  tiempo_de_paquete= hex2dec(tline(54)) * 256 + hex2dec(tline(53)) - paquetes_validos;
  
  % Valor m醲imo de los 24 bits q se transmite es (0xFFFF00). Ahora,  0xFFFF00/2=  0x7FFF80 (que en decimal es 8388480)

    canal1_tot=hex2dec(tline (puntero:puntero+1) ) * 1+  (hex2dec(  tline (puntero+3:puntero+4) )) * 2^8+ (hex2dec(  tline (puntero+6:puntero+7) )) *2^16;%-8388480;
     puntero=puntero+9;
    canal2_tot=hex2dec(tline (puntero:puntero+1) ) * 1+  (hex2dec(  tline (puntero+3:puntero+4) )) * 2^8+ (hex2dec(  tline (puntero+6:puntero+7) )) *2^16;%-8388480; 
    puntero=puntero+9;
    canal3_tot=hex2dec(tline (puntero:puntero+1) ) * 1+  (hex2dec(  tline (puntero+3:puntero+4) )) * 2^8+ (hex2dec(  tline (puntero+6:puntero+7) )) *2^16;%-8388480; 
    puntero=puntero+9;
    canal4_tot=hex2dec(tline (puntero:puntero+1) ) * 1+  (hex2dec(  tline (puntero+3:puntero+4) )) * 2^8+ (hex2dec(  tline (puntero+6:puntero+7) )) *2^16;%-8388480; 
    puntero=puntero+9;
    canal5_tot=hex2dec(tline (puntero:puntero+1) ) * 1+  (hex2dec(  tline (puntero+3:puntero+4) )) * 2^8+ (hex2dec(  tline (puntero+6:puntero+7) )) *2^16;%-8388480;
     puntero=puntero+9;
    canal6_tot=hex2dec(tline (puntero:puntero+1) ) * 1+  (hex2dec(  tline (puntero+3:puntero+4) )) * 2^8+ (hex2dec(  tline (puntero+6:puntero+7) )) *2^16;%-8388480; 
    puntero=puntero+9;
    canal7_tot=hex2dec(tline (puntero:puntero+1) ) * 1+  (hex2dec(  tline (puntero+3:puntero+4) )) * 2^8+ (hex2dec(  tline (puntero+6:puntero+7) )) *2^16;%-8388480;
     puntero=puntero+9;
    canal8_tot=hex2dec(tline (puntero:puntero+1) ) * 1+  (hex2dec(  tline (puntero+3:puntero+4) )) * 2^8+ (hex2dec(  tline (puntero+6:puntero+7) )) *2^16;%-8388480; 
    puntero=puntero+9;

%    hexa_final1=dec2hex (canal1_tot) 
%    hexa_final2=dec2hex (canal2_tot) 
%    hexa_final3=dec2hex (canal3_tot) 
%   hexa_final4=dec2hex (canal4_tot) 
%    hexa_final5=dec2hex (canal5_tot) 
%   hexa_final6=dec2hex (canal6_tot)    
%    hexa_final7=dec2hex (canal7_tot) 
%   hexa_final8=dec2hex (canal8_tot)
 
   canal1_data(pasada)=canal1_tot;
   canal2_data(pasada)=canal2_tot;
   canal3_data(pasada)=canal3_tot;
   canal4_data(pasada)=canal4_tot;
   canal5_data(pasada)=canal5_tot;
   canal6_data(pasada)=canal6_tot;
   canal7_data(pasada)=canal7_tot;
   canal8_data(pasada)=canal8_tot;
     
    %double(canal1);
    pasada=pasada+1
   % pause(0.1);
   tline = fgetl(fid);
end	
 pasada=pasada-1;
t=[0:1:pasada-1];
subplot (3, 3, 1); plot (t,canal1_data );
subplot (3, 3, 2); plot (t,canal2_data );
subplot (3, 3, 3); plot (t,canal3_data );
subplot (3, 3, 4); plot (t,canal4_data );
subplot (3, 3, 5); plot (t,canal5_data );
subplot (3, 3, 6); plot (t,canal6_data );
subplot (3, 3, 7); plot (t,canal7_data );
subplot (3, 3, 8); plot (t,canal8_data );

plot (t,canal1_data )

fclose(fid);