[Tip] Owon Oscilloscoop data uitlezen

% Open het .bin databestand. "le" opent het bestand als Little-Endian
fileID = fopen("20160722_255650.bin", "r", "ieee-le");

% Lees het .bin databestand, dit bestand bestaat uit signed characters
% en deze worden in de array 'data' opgeslagen.
data = fread(fileID,"signed char");

% Sluit het bestand, de data uit het bestand staat nu in het geheugen.
fclose(fileID);

% BEGIN MET HET VERWERKEN VAN DE DATA

% Deze arrays bevatten de niveau's van de volt/div en van de t/div.
voltage_level_lookup  = [2e-3 5e-3 10e-3 20e-3 50e-3 100e-3 200e-3 500e-3 1 2 5 10];
timebase_level_lookup = [2e-9 5e-9 10e-9 20e-9 50e-9 100e-9 200e-9 500e-9 1e-6 2e-6 5e-6 10e-6 20e-6 50e-6 100e-6 200e-6 500e-6 1e-3 2e-3 5e-3 10e-3 20e-3 50e-3 100e-3 200e-3 500e-3 1 2 5 10 20 50 100];


% Byte 57 bevat het kanaal nummer, geformateerd als ASCII. Door daar 48 af 
% te trekken wordt het kanaalnummer van ASCII omgezet naar decimaal. 
channel = data(57) - 48


% Byte 58 t/m 61 bevat de data block lengte (geteld vanaf Byte 58). De 
% typecast maakt van de 4 losse bytes één 32-bit integer. Als deze 
% block lengte negatief is dan gaat het om een deep memory wave.
data_length = typecast(int8(data(58:61)),'int32');


% Byte 82 t/m 85 bevat de index voor het timebase level welke daarna
% gebruikt kan worden om de bijbehorende array "timebase_level_lookup" uit
% te lezen. Er wordt 1 bij opgeteld omdat de array index beginnen te tellen
% bij 1 en de ontvangen index begint te tellen bij 0.
timebase = timebase_level_lookup(typecast(int8(data(82:85)),'int32')+1)


% Byte 86 t/m 89 bevat de voltage offset van het kanaal. Dus op welk 
% punt de 0 volt lijn ligt (op het scopescherm).
y_offset = typecast(int8(data(86:89)),'int32');


% Byte 90 t/m 93 bevat de index voor het voltage level welke daarna
% gebruikt kan worden om de bijbehorende array "voltage_level_lookup" uit
% te lezen. Er wordt 1 bij opgeteld omdat de array index beginnen te tellen
% bij 1 en de ontvangen index begint te tellen bij 0.
voltage = voltage_level_lookup(typecast(int8(data(90:93)),'int32')+1)


% Byte 94 t/m 97 bevat de ingestelde probe verzwakking. Deze waarde wordt
% gegeven als macht van 10. (Net zoals het exponent bij wetenschappelijke
% notatie).
probe_attenuation = 10 ^ typecast(int8(data(94:97)),'int32');


% Byte 98 t/m 101 bevat de x-as interval. Gegeven in uS. Deze is nodig 
% om tijdens het plotten de x-as te maken. Moet gelezen worden als float
% en wordt daarom getypecast naar 'single'.
timebase_interval = typecast(int8(data(98:101)),'single')*1e-6;


% Byte 102 t/m 105 zouden volgens de handleiding de frequentie moeten bevatten
% in alle binfiles die ik door de scope heb laten maken staat deze op 
% nul. Ik moet nog uitzoeken wat ik hiermee kan.


% Byte 106 t/m 109 zouden volgens de handleiding de cycle moeten bevatten
% in alle binfiles die ik door de scope heb laten maken staat deze op 
% nul. Ik moet nog uitzoeken wat ik hiermee kan.


% Byte 110 t/m 113 bevat de y-as interval. Gegeven in mV (of mA als op
% het kanaal de ampere probe is ingeschakeld). Deze moet vermenigvuldigd
% worden met de wave om de exacte spanningsniveaus weer te geven.
voltage_interval = typecast(int8(data(110:113)),'single')*1e-3;


% Byte 114 t/m de laatste byte van het bestand bevat de wave. Deze moet 
% worden gelezen als allemaal signed acht bit getallen. De functie 'double'
% zorgt ervoor dat deze wordt opgeslagen als een float, zodat deze later met
% de y-as interval (welke ook een float is) kan worden vermenigvuldigd.
wave = double(data(118:end) - y_offset);


% Hier begint het plotten, maak er maar wat moois van! :)
plot([0:timebase_interval:timebase_interval*(numel(wave)-1)], wave.*voltage_interval);