MAKE A MEME
View Large Image
| View Original: | MDKQ1.svg (324x252) | |||
| Download: | Original | Medium | Small | Thumb |
| Courtesy of: | commons.wikimedia.org | More Like This | ||
| Keywords: MDKQ1.svg de Fit einer Logistischen Funktion an Messdaten own Neufassung von File LogistischeFunktion png vom Benutzer de User Philipendula Christian Schirm 2011-02-01 Cc-zero Quellen Skript zur Bildgenerierung Erzeugungsskript um die Grafik zu erstellen Anleitung Benötigte Open-Source-Software http //www python org/download Python Python-Paket http //sourceforge net/projects/numpy numpy Python-Paket http //matplotlib sourceforge net matplotlib Nach der Installation von Python den Quelltext in eine Datei mdkq py kopieren und starten durch Doppelklicken oder in der Konsole durch Eingabe von python mdkq py Python-Skript Created with Matplotlib <source lang python > This source code is public domain Author Christian Schirm import numpy pylab from matplotlib font_manager import FontProperties import matplotlib pyplot as plt G 1 k 1 f0 0 5 x numpy linspace -8 8 16 1 -1 numpy random seed 50 y numpy random normal G/ 1+numpy exp -k G x G/f0-1 0 1 numpy random seed int time time err 1E8 err numpy mean numpy square y-G/ 1+numpy exp -k G x G/f0-1 print err G k f0 numpy random seed 2 for i in range 5000 faktor 1+0 01 numpy random rand - 5 for ivar in 1 2 3 backup err G k f0 var backup varivar varivar faktor err G k f0 var err_neu numpy mean numpy square y-G/ 1+numpy exp -k G x G/f0-1 if err_neu<err err err_neu else varivar backupivar err G k f0 var print err G k f0 Fehlerquadrat minimiert xneu numpy linspace -8 8 50 yneu G/ 1+numpy exp -k G xneu G/f0-1 xr x yr G/ 1+numpy exp -k G xr G/f0-1 residuen for i in range len x residuen+ xi+8 xi+8 yi 10 yri 10 'g-' plt clf fig plt figure figsize 4 5 3 5 fig subplotpars bottom 0 15 y0 plt plot residuen -3 plt setp y0 color ' 60c060' linewidth 1 5 y0 plt plot residuen-3 label 'Residuum' plt setp y0 color ' 60c060' linewidth 1 5 y2 plt plot xneu+8 yneu 10 'r-' label 'Modellfunktion' plt setp y2 linewidth 1 5 y1 plt plot x+8 y 10 'o' label 'Messungen' plt xlabel 'x' plt ylabel 'y' font FontProperties font set_size 'medium' leg plt legend y1 y2 y0 'Messpunkte' 'Modellfunktion' 'Residuen' frameon True loc 'lower right' labelspacing 0 3 prop font plt grid True plt show plt savefig 'MDKQ1 svg' </source> Regression analysis Numerical analysis | ||||
