nap=0.715; %napeti tlak pred clonkou vs atmosfera [V]
adr=[73,71];
n=10000; % otacky dmychadla
k=0;
v_tun=0.71; %rychlost v tunelu
rh=53.7; %vlhkost
t=24.5; %teplota
po=750; % tlak atmosfericky v torr
savefile='C:\Users\Robin\ownCloud\škola\cvut FS\6.semestr- letní\bakalaøka\charakteristiky ventilatoru\Matlab\mereni_2_10000\2_10000_071.mat';
if s.bytesavailable>0,fread(s,s.bytesavailable);end
clear hodproud I c;
while (1)
    try
    k=k+1;
    fprintf('%03i> ',k);
    for i=1:1:length(adr)
        fwrite(s,[adr(i) 0 7]);tmp=fread(s,6);
        c(i)=tmp(5)+tmp(6)*2^8;
        I(i)=(c(i)*3.3)/(2^16*165);
        fprintf('%.3f\t ',I(i)*1e3);  
    end
    tlak_1(i)=625000*I(1)-2500;
    tlak_2(i)=625000*I(2)-2500;
    tlak_3(i)=2222.22*nap-1111.11;
    rychlostzt(i)=v_tun*31.085; %rychlost v potrubí vypoèítaná pøes rychlost v tunelu
    objprutzt(i)=rychlostzt(i)*0.01327;%objemový prùtok v potrubí vypoèítaný pøes rychlost v tunelu
    [objprutzc(i), v(i)]=vypocet_objemoveho_prutoku(tlak_3(i),tlak_2(i));%objemový prùtok v potrubí vypoèítaný z clony
    fprintf('%.3f\t ',tlak_1(i));
    fprintf('%.3f\t ',tlak_2(i));
    fprintf('%.3f\t ',objprutzc(i));
    fprintf('%.3f\t ',objprutzt(i));
    fprintf('%.3f\t ',v(i));
   
    fprintf('\n');
    hodproud(k,:)=[c(1) c(2) I(1) I(2) nap rh t po n v_tun tlak_1 tlak_2 tlak_3 objprutzc objprutzt];% hodnoty proudu tlakovych prrvodniku
    save(savefile, 'hodproud');
    pause(0.5)
    end
end

