Keywords: MDKQ 3D.svg de Fit einer 2D-Funktion own Chrischi 2011-04-07 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 from matplotlib font_manager import FontProperties from mpl_toolkits mplot3d import axes3d import matplotlib pyplot as plt import numpy numpy random seed 11 x0 numpy random rand 10 x0 x0-min x0 ;x0 10/max x0 x1 numpy random rand 10 10 x1 x1-min x1 ;x1 10/max x1 y 10- x0-5 2+ x1-5 2 0 05+numpy random rand 10 10 x zip x0 x1 ordnung 2 2 ab hier gueltig fuer beliebige Dimension faktor i+1 for i in ordnung+1 for i in reversed range 1 len faktor faktori-1 faktori-1 faktori base numpy asarray i/faktorj+1 ordnungj+1 for j in range len ordnung for i in range faktor0 def basefunc b x return numpy prod xi bi for xi bi in zip x b N len base A numpy zeros N N for i bivec in zip range N base for j bjvec in zip range N base Ai j sum basefunc bivec xvec basefunc bjvec xvec for xvec in x b numpy zeros N for i bivec in zip range N base bi sum basefunc bivec xvec yi for xvec yi in zip x y c numpy linalg solve A b bis hier gueltig fuer beliebige Dimension xneu numpy meshgrid numpy linspace 0 10 num 20 numpy linspace 0 10 num 20 xneuPairs zip xneu0 ravel xneu1 ravel yneu numpy array sum ci basefunc b xvec for b ci in zip base c for xvec in xneuPairs yneu shape xneu0 shape fig plt figure ax2 fig add_subplot 111 ax fig add_subplot 111 projection '3d' X Y Z xneu0 xneu1 yneu p0 ax2 scatter 0 0 p1 ax scatter x0 x1 y zdir 'z' p2 ax plot_wireframe X Y Z color 'r' rstride 1 cstride 1 lx0 ax set_xlabel ' x1 ' lx1 ax set_ylabel ' x2 ' ly ax set_zlabel ' y ' font FontProperties font set_size 'medium' leg plt legend p0 p2 'Messpunkte' 'Modellfunktion' frameon True loc 'lower right' prop font leg draggable plt show </source> Regression analysis Numerical analysis Images with Python source code |