Advertisement
Gigli-0neiric

Untitled

Jan 9th, 2016
338
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
MatLab 1.92 KB | None | 0 0
  1. function [] = ausgleichsgerade_doppelt_logarithmisch ()
  2. %Messwerte aus der Aufgabenstellung
  3. Daten = [
  4. 0.050 0.3
  5. ...
  6. ];
  7. %% Aufteilung in x- und y-Werte
  8. x_Werte = Daten(:,1);
  9. y_Werte = Daten(:,2);
  10. %% Grafik 1 - Plot der gegebenen Daten in figure(1) als Messpunkte ohne Verbindungslinie
  11. figure(1)
  12. h1 = plot(x_Werte,y_Werte,'*'); % '*' <- Anstatt der Linie werden * verwendet
  13. hold on
  14. %% Grafik 2 - Berechnung der Ausgleichsgerade in doppelt-log-Darstellung
  15. figure(2)
  16. logX = log10(x_Werte);
  17. logY = log10(y_Werte);
  18. % Hängt eine Messgröße – vermutet oder theoretisch begründet – von einer Potenz
  19. % einer anderen Größe ab, so ist i. d. R. der Exponent dieser Potenz die wesentliche
  20. % charakterisierende Größe. Er lässt sich mit Hilfe einer Ausgleichsgeraden bestimmen,
  21. % wenn eine doppelt-logarithmische Darstellung der Größen gewählt wird.
  22. % Durch Logarithmieren erhält man die Geradengleichung in der die Steigung die gesuchte Potenz m liefert.
  23. % b liefert der y-Abschnitt.
  24. [m b] = ausgleichsgerade(logX, logY);
  25. vag = linspace(0,1.2);
  26. %Zurück logarithmieren
  27. ag = 10.^(m*log10(vag)+b);
  28. %ag und x-,y-Werte werden logairthmiert und dargestellt.
  29. loglog(x_Werte,y_Werte,'*', vag, ag, 'r')
  30. grid
  31. set(gca,'xminorgrid','on')
  32. set(gca,'yminorgrid','on')
  33. xlabel('Beleuchtungsstärke')
  34. ylabel('Dimlevel')
  35. legend('Messwerte','Ausgleichsgerade',3)
  36. hold off
  37. %% Berechnung der Ausgleichspolynoms
  38. %Die Berechnung der Steigung für die Ausgleichkurve in einer nicht
  39. %logarithmischen Darsellung findet ebenfalls über logarithmieren der Werte
  40. %und der Berechnung der Steigung von der Ausgleichkurve statt.
  41. agp = polyfit(x_Werte,y_Werte, ceil(m));
  42. pv = polyval(agp, x_Werte);
  43. %% Grafik 1 - Plot der Ausgleichkurve in figure(1)
  44. figure(1)
  45. grid
  46. h2 = plot(x_Werte, pv, 'r');
  47. legend([h1,h2],'Messwerte','Ausgleichsgerade',3)
  48. set(gca,'xminorgrid','on')
  49. set(gca,'yminorgrid','on')
  50. xlabel('xlabel')
  51. ylabel('ylabel')
  52. grid;
  53. hold off;
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement