Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- function [] = ausgleichsgerade_doppelt_logarithmisch ()
- %Messwerte aus der Aufgabenstellung
- Daten = [
- 0.050 0.3
- ...
- ];
- %% Aufteilung in x- und y-Werte
- x_Werte = Daten(:,1);
- y_Werte = Daten(:,2);
- %% Grafik 1 - Plot der gegebenen Daten in figure(1) als Messpunkte ohne Verbindungslinie
- figure(1)
- h1 = plot(x_Werte,y_Werte,'*'); % '*' <- Anstatt der Linie werden * verwendet
- hold on
- %% Grafik 2 - Berechnung der Ausgleichsgerade in doppelt-log-Darstellung
- figure(2)
- logX = log10(x_Werte);
- logY = log10(y_Werte);
- % Hängt eine Messgröße – vermutet oder theoretisch begründet – von einer Potenz
- % einer anderen Größe ab, so ist i. d. R. der Exponent dieser Potenz die wesentliche
- % charakterisierende Größe. Er lässt sich mit Hilfe einer Ausgleichsgeraden bestimmen,
- % wenn eine doppelt-logarithmische Darstellung der Größen gewählt wird.
- % Durch Logarithmieren erhält man die Geradengleichung in der die Steigung die gesuchte Potenz m liefert.
- % b liefert der y-Abschnitt.
- [m b] = ausgleichsgerade(logX, logY);
- vag = linspace(0,1.2);
- %Zurück logarithmieren
- ag = 10.^(m*log10(vag)+b);
- %ag und x-,y-Werte werden logairthmiert und dargestellt.
- loglog(x_Werte,y_Werte,'*', vag, ag, 'r')
- grid
- set(gca,'xminorgrid','on')
- set(gca,'yminorgrid','on')
- xlabel('Beleuchtungsstärke')
- ylabel('Dimlevel')
- legend('Messwerte','Ausgleichsgerade',3)
- hold off
- %% Berechnung der Ausgleichspolynoms
- %Die Berechnung der Steigung für die Ausgleichkurve in einer nicht
- %logarithmischen Darsellung findet ebenfalls über logarithmieren der Werte
- %und der Berechnung der Steigung von der Ausgleichkurve statt.
- agp = polyfit(x_Werte,y_Werte, ceil(m));
- pv = polyval(agp, x_Werte);
- %% Grafik 1 - Plot der Ausgleichkurve in figure(1)
- figure(1)
- grid
- h2 = plot(x_Werte, pv, 'r');
- legend([h1,h2],'Messwerte','Ausgleichsgerade',3)
- set(gca,'xminorgrid','on')
- set(gca,'yminorgrid','on')
- xlabel('xlabel')
- ylabel('ylabel')
- grid;
- hold off;
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement