[ < ] [ globale Übersicht ] [ Kapitelübersicht ] [ Stichwortsuche ] [ > ]
Ziel der analytischen Modellbildung ist die Gewinnung von Funktionen
,
deren Definitionsbereich B eine zusammenhängende Teilmenge
des n
, also z.B. ein endliches oder auch unendliches Intervall, ist. Abhängig
von der Dimension n des Definitionsbereichs B unterscheidet
man
univariate Modelle mit einer unabhängigen Variable (n=1) und
multivariate Modelle mit mehreren unabhängigen Variablen (n>1).
Es wird jedoch bei der Modellbildung keine strukturelle Analogie zwischen Objekt und Modell gefordert, d.h. es werden nur Funktionsmodelle betrachtet, die
Wie in der obigen Formel zur Modellbildung ersichtlich, handelt es sich um Funktionen, die jedem Element x ihres Definitionsbereichs B genau ein Element ihres Wertebereichs zuordnen. Falls als Modell Kurven und Flächen im Raum in Betracht gezogen werden, können auch diese durch Modellfunktionen dargestellt werden.
Beispiel (Raumkurve) Eine Kurve im 3,
die durch die Punkte
Pi=(ui , vi , wi )
3,
i=1,2,3,...,k
verläuft, kann durch drei interpolierende Funktionen f u(s), f v(s) und f w(s) modelliert werden, wobei folgende Bedingung erfüllt sein muß:
Wie bereits angesprochen, wird durch die Forderung nach der Verarbeitbarkeit der Modelle auf einem Computer das Ergebnis eines kontinuierlich-funktionalen Modellierungs- bzw. Approximationsprozesses auf elementare Funktionen eingeschänkt.
Definition: Elementare Funktionen lassen sich durch eine endliche Anzahl von Parametern charakterisieren. Sie sind durch Formeln definierbar, in denen nur endlich viele Auswertungen der auf einem Computer verfügbaren Standardfunktionen (algebraisch, trigonometrisch oder logarithmisch) und allenfalls Fallunterscheidungen durch bedingte Anweisungen (z.B. bei stückweise definierten Funktionen) auftreten.
Als Modell- bzw. Approximationsfunktionen kommen zum Beispiel in Frage:
Funktionen, die diese Forderung nicht erfüllen, weil sie z.B. unendlich viele Koeffizienten einer Taylor-Reihe zu ihrer Charakterisierung benötigen, müssen durch elementare Funktionen approximiert (modelliert) werden.
Elementare Funktionen für analytische Operationen
Die Klasse der elementaren Funktionen kann für verschiedene Aufgabenstellungen noch immer zu umfassend sein. Eine besonders wichtige Modellbildungssituation, wo oft weitere Einschränkungen erforderlich sind, bilden numerische Verfahren zur Lösung mathematischer Probleme wie z.B. Integration, Differentiation etc. Diese Probleme nönne im allgemeinen nicht direkt gelöst werden. Konstruktive Lösungsmethoden erfordern daher (oft automatisch-algorithmitsch ablaufende) Modellbildungs- bzw. Approximationsvorgänge:
Dabei ist die Wahl einer geeigneten Klasse von Ersatzfunktionen stark von den Besonderheiten der jeweiligen Problemstellung abhängig. Eine universelle Regel ist:
Das Modellproblem, dessen Modellfunktionen an die Stelle der Originalfunktionen getreten sind, soll exakt lösbar sein.
Bei numerischen Verfahren am häufigsten verwendete Klasse elementarer Funktionen ist die der Polynome, da diese bei den meisten Aufgabenstellungen eine exakte und mit geringem Rechenaufwand verbundene Lösung erlauben.
Beispiel: Die Ermittlung von Extremwerten (Maxima, Minima, Nullstellen) analytischer Daten eines mathematischen Problems am Computer ist auf der Basis einer Nullstellenbestimmung der ersten Ableitung meist nicht direkt möglich. Die Funktion, deren Extrema gesucht sind, muß durch geeignete Modellfunktionen ersetzt werden, die eine einfache und exakte Extremwertbestimmung erlauben. Im univariaten Fall sind dies vorzugsweise quadratische Polynome:
Beschreibung: Eine gegebenen Funktion wird diskretisiert und anhand der Diskretisierungspunkte ein quadratisches Polynom interpoliert. Der Extremwert des Interpolationspolynoms wird bestimmt und als Näherung für das gesuchte Minimum verwendet.
In den folgenden Kapiteln wird auf diese Vorgänge (Diskretisierung, Homogenisierung) weiter eingegangen.
Im Unterschied zu den elementaren Funktionen geschieht bei der Darstellung und Implementierung von Modellen mit Hilfe von Operatorgleichungen die Auswertung konkreter Werte nicht formelartig, sondern iterativ/prozeßartig. Die numerischen Werte hängen daher unter anderem von der Konvergenzgeschwindigkeit des Iterationsverfahren an der Stelle der Auswertung, vom Abbruchkriterium (der geforderten Approximationsgenauigkeit) oder der aktuellen Maschinenarithmetik ab.
Beispiel: Quadratwurzelfunktion
Die Standardfunktion wird
nicht mit Hilfe einer Approximationsfunktion implementiert, sondern
durch einen Algorithmus, der die iterative Lösung der Gleichung
(bzw. des Nullstellenproblems)
beschreibt.
Bei diesem Algorithmus wird zunächst der Exponent e(x) halbiert
und anschließend das Newton-Verfahren zur Bestimmung der Quadratwurzel
der Mantisse mx := M(x) der Gleitpunktzahl x angewendet, d.h.
wird
iterativ gelöst:
Dies wird solange fortgesetzt, bis ein gefordertes Abbruchkriterium,
also z.B. ,eintritt,
bzw. bis imax erreicht wird, das genauso wie die Startnäherung
g0 abhängig von der Maschinenarithmetik gewählt wird.