% Anpassning av exponentialkurva

% Importera tabellerna med uppmätta värden
DATA1 = importdata('1_2Cooldown_DOT_DECIMAL.txt');
DATA2 = importdata('1_2Reference_DOT_DECIMAL.txt');
x = DATA1(:,1)';
y = DATA1(:,2)';

%Beräkna medelvärde för rumstemperatur
Tmean = mean(DATA2(:,2));
hold on; grid on;
plot(x, y, 'o')
% Subtrahera rumstemperaturen
y = y - Tmean;

% Beräkna koefficienter med Least Square metoden 
% (här används det alternativa beräkningssättet från boken)
coeff = [size(x, 2) sum(x); sum(x) sum(x.^2)] \ [sum(log(y)); sum(x.*log(y))]

c0 = exp(coeff(1))
c1 = coeff(2)

% Skapa den anpassade exponentialkurvan och återinför
% rumstemperaturoffsetten, och plotta.
ylin = c0*exp(c1.*x) + Tmean;
plot(x, ylin)


