% Redovisning: Kalibreringsdata NTC-motstånd

hold off;
% Importera tabellen med uppmätta värden
DATA = importdata('1_1Matpunkter_edit_DOT_DECIMAL.txt')
T = DATA(:, 2);
U_R = DATA(:, 3);
% Plotta U_R som funktion av T
plot(T, U_R, 'or')
xlabel('Temperatur C'); ylabel('Spänning U_R');  title('Kalibreringsdata för NTC')
grid on; hold on;

% Least Square Method för att hitta koefficienterna
% Inför kolumn av 1:or i T
T = [ones(size(T)), T]
% Beräkna pseudoinvers av T
T = pinv(T);
% Multiplicera T med U_R för att erhålla koefficienterna
koeff = T*U_R
% Plotta linjen med beräknade koefficienter
x = 22:60; plot(x, koeff(2)*x + koeff(1), '-b')

legend('Mätpunkter', 'Anpassade linjen', 'Location', 'southeast')
