Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- function [m, b, sm, sb] = ausgleichsgerade_halb_logarithmisch()
- %Messwerte aus der Aufgabenstellung
- x=0:0.1:5;
- T=5;
- y=7*exp(x/T)
- Daten = [x;y]
- %% Aufteilung in x- und y-Werte
- x_Werte = Daten(1,:);
- y_Werte = Daten(2,:);
- % Berechnung
- [p,s] = polyfit(x_Werte,log(y_Werte),1); m = p(1), b = p(2)
- covp = inv(s.R)*inv(s.R)'*s.normr^2/s.df;
- sm = sqrt(covp(1,1))
- sb = sqrt(covp(2,2))
- %% Grafik
- Iag = linspace(0,1.2); % 100 Werte für die x-Achse, äquidistant verteilt
- Uag = polyval(p,Iag);
- plot(x_Werte,log(y_Werte),'*',Iag,Uag,'r')
- grid
- set(gca,'xminorgrid','on')
- set(gca,'yminorgrid','on')
- xlabel('X-Achse')
- ylabel('Y-Achse')
- legend('Messwerte','Ausgleichsgerade',3)
- hold off
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement