50 likes | 292 Views
Übung zur Numerische Geologie. Entwicklung eines Summenkurven-Programms in Matlab . Import der Excel-Dateien über die Import-Option im Menu 'File‚. who. 2. Berechnung der Wertesumme und der Länge des Datensatzes. Summe = sum (data) Laenge = length (data).
E N D
Übung zur Numerische Geologie Entwicklung eines Summenkurven-Programms in Matlab • Import der Excel-Dateien über die Import-Option im Menu 'File‚ who 2. Berechnung der Wertesumme und der Länge des Datensatzes Summe = sum (data) Laenge = length (data) 3. Zuordnung der einzelnen Prozentwerte für jede Klasse (3-Satz-Rechnung zur Ermittlung der Einzelwerte pro Klasse) Prozentwert(1) = (data(1) * 100) / Summe Prozentwert(2) = (data(2) * 100) / Summe Bessere Lösung ist eine Programmierung als allgemeingültige while-Schleife i=1; while i <= Laenge i = i+1; end Prozentwert(i) = (data(i) * 100) / Summe Klemens Seelos 13-17.05.2005
Übung zur Numerische Geologie 4. Kurze Grafikausgabe um vorläufiges Programm zu testen: figure; bar (Prozentwert) 5. Erstellung einer Summenkurve aus den errechneten Prozentwerten (Einzelschritte) Summenkurve(1) = Prozentwert(1) Summenkurve(2) = Prozentwert(1) + Prozentwert(2) oder einfacher und eleganter Summenkurve(2) = sum (Prozentwert(1:2)) bessere Lösung ist eine Programmierung als allgemeingültige while-Schleife i=1; while i <= Laenge i = i+1; end Summenkurve(i) = sum (Prozentwert(1:i) ) Klemens Seelos 13-17.05.2005
Übung zur Numerische Geologie 4. Kurze Grafikausgabe um vorläufiges Programm zu testen: figure; bar (Summenkurve) 5. Eliminierung des Startfehlers ( Kurve muss bei Null beginnen) Summenkurve(1) = 0; i = 1; while i <= Laenge Summenkurve(1+i) = sum (Prozentwert(1:i) ) i = 1+i; end LaengeSummenkurve = length (Summenkurve) 6. Kurze Grafikausgabe um vorläufiges Programm zu testen: figure; bar (Summenkurve) Klemens Seelos 13-17.05.2005
Übung zur Numerische Geologie 7. Interpolation der Summenkurve, um glatten Verlauf zu gewährleisten und um später Perzentile rechnen zu können pchip-Interpolation der Kurve 7a. Festlegung des Start- und Endpunktes der Interpolation Startendpunkt = [ 0.01 : LaengeSummenkurve ] 7b. Festlegung des Interpolationsintervalls Intervall = 0.01 : 0.1 : LaengeSummenkurve - 0.9 7c. Durchführung der Interpolation InterpolKurve = pchip (Startendpunkt,Summenkurve,Intervall) 8. Kurze Grafikausgabe um vorläufiges Programm zu testen: figure; subplot(2,1,1); bar (Prozentwert) hold on plot (Intervall, InterpolKurve, 'r') subplot(2,1,2); plot (Intervall, InterpolKurve, 'r') hold on plot (Summenkurve) Klemens Seelos 13-17.05.2005
Übung zur Numerische Geologie 9. Automatische Berechnung der Verteilungsbreite UeberNull = find (data > 0) Verteilungsbreite = length (UeberNull) 10. Berechnung der Perzentile Perz25 = find (InterpolKurve >= 24.5 & 25.5 >= InterpolKurve) Perzentil25 = mean(Intervall(Perz25)) Perzentil25 = 2^(4+Perzentil25) Dieser Algorithmus muss nun für alle übrigen Perzentile angewendet werden: Perz10, Perz50, Perz60, Perz75, Perz95 text(0.05,0.92‚‘Perz10: ','FontSize',7,'Color','k') text(0.72,0.92,[num2str(Perz10)],'FontSize',7,'Color',[0.2 0.2 0.4]) Klemens Seelos 13-17.05.2005