% Beräkning av tidskonstant för Pt-100 element
% samt anpassning av kurva till mätpunkter.

hold off
hold on; grid on;
% Importera tabellerna med uppmätta värden

DATA1 = importdata('1_2Cooldown_DOT_DECIMAL.txt');
DATA2 = importdata('1_2Reference_DOT_DECIMAL.txt');
Tid_exp = DATA1(:,1);
T_exp = DATA1(:,2);
Tid_rums = DATA2(:,1);
T_rums = DATA2(:,2);

%Tid_exp = [0 2 3 6 10];
%T_exp = [5792 1027 548 32 8];
% Beräkna medelvärdet av rumstemperaturen
%Tmean = mean(T_rums);
Tmean = mean(T_rums);

% Minsta kvadrat metoden y = x*beta
y = log(T_exp - Tmean);



x = [ones(size(Tid_exp)), Tid_exp];
beta = pinv(x)*y;

% beta innehåller nu ln(C_0) och C_1

C_0 = exp(beta(1))
C_1 = beta(2)

T_anpassad = C_0.*exp(C_1 .* Tid_exp) + Tmean;
plot(Tid_exp, T_anpassad)

plot(Tid_exp, T_exp,'o')
xlabel('Tid'); ylabel('Temperatur C'); title('Passiv nedkylning av Pt-100 element');
legend('Anpassade exponentialkurvan', 'Mätpunkter'); 

dim = [.2 .5 .3 .3];
str = 'C_0 = 5,9659 C_1 = -0,087';
annotation('textbox',dim,'String',str,'FitBoxToText','on')



