% T8A5_Ute

load reizT8
load ratenT8

% Daten plotten

figure
plot(reizT8,ratenT8,'o')
xlabel('Stromstärke [nA]')
ylabel('Aktionspotenzial [Sp/s]')
legend('Messwerte')
axis([-3 6 -100 250])
hold on
pause

% Gerade mit polyfit

x=-3:0.5:6;
[p,s]=polyfit(reizT8,ratenT8,1);
schaetz=polyval(p,x,s);

plot(x,schaetz,'g')
legend('Messwerte','Angepasste Gerade')
pause

% Schätzwerte

x=-3:0.5:6;
schaetz=polyval(p,x,s);

plot([-2 5],[schaetz(3) schaetz(end-2)],'^r')
legend('Messwerte','Angepasste Gerade','Schätzwerte')
xlabel('Stromstärke [nA]')
ylabel('Aktionspotenzial [Sp/s]')
text(-2.2,-60,['Schätzwert bei -2 nA: ' num2str(schaetz(3))])
text(1.8,160,['Schätzwert bei 5 nA: ' num2str(schaetz(end-2))])
pause

close