% Beräkning av data från värmepump (ver.2.3, 23.8.2004, Thomas Nilsson)
% Energipulsräkning, temperatur i varm och kall hink samlas med Pasco datainsamlingssytemet
% tillsammans med Science WorksShop-mjukvaran.

clear f_name row_1 n_cols data T_H T_C n_pts E t

[t E T_H T_C] = read_vpdata;

% Skriv meddelande till skärmen
fprintf('\nFöljande kolumnvektorer har skapats:\n');
fprintf('   t : tidvärden (i sekunder)\n');
fprintf('   E : energivärden (i joule)\n');
fprintf('   T_H : temperaturvärden i den varma hinken (i grader Celsius)\n');
fprintf('   T_C : temperaturvärden i den kalla hinken (i grader Celsius)\n');

% Polynomanpassning till experimentdata
T_H_pol = polyfit(t, T_H, 3);
T_C_pol = polyfit(t, T_C, 3);
E_pol = polyfit(t, E, 3);

figure(1)
plot(t, T_H, '.r', t, polyval(T_H_pol,t), '-k')
title('Temperaturen på den varma sidan')
xlabel('tid / s')
ylabel('T_H / grader C')
grid on

figure(2)
plot(t, T_C, '.b', t, polyval(T_C_pol,t), '-k')
title('Temperaturen på den kalla sidan')
xlabel('tid / s')
ylabel('T_C / grader C')
grid on

figure(3)
plot(t, E, '.g', t, polyval(E_pol,t), '-k')
title('Elförbrukningen')
xlabel('tid / s')
ylabel('W_e_i_n / J')
grid on

