Differential-algebraische Gleichung

aus Wikipedia, der freien Enzyklopädie
Wechseln zu: Navigation, Suche

In einer Differential-algebraischen Gleichung (auch differentiell-algebraische Gleichung, Algebro-Differentialgleichung oder Deskriptor-System) sind gewöhnliche Differentialgleichungen und algebraische (d. h. hier: ableitungsfreie) Nebenbedingungen gekoppelt und werden als eine Gleichung bzw. Gleichungssystem aufgefasst. In einigen Fällen ist diese Struktur schon in der Form des Gleichungssystems angelegt, z. B. in

\begin{align}
\dot x(t)&= f(t,x(t),y(t))\\
0&=g(t,x(t),y(t))
\end{align}

Diese Form ergibt sich regelmäßig bei Problemen aus der Mechanik von Körpern unter Zwangsbedingungen, als instruktives Beispiel wird oft das Pendel gewählt.

Die allgemeinste Form einer differentiell-algebraischen Gleichung ist eine implizite Differentialgleichung in der Form

 F(\dot{x}(t),x(t),t)=0 , F:\R^n\times\R^n\times \R\to\R^m

für eine vektorwertige Funktion x:I\rightarrow\R^m mit I\subset\R. Eine Gleichung in dieser impliziten Form ist (lokal) nach \dot{x} auflösbar, wenn die partielle Ableitung F_{\dot{x}} regulär ist. Dies folgt aus dem klassischen Satz über implizite Funktionen. In diesem speziellen Fall kann man die implizite Gleichung umschreiben in die Form

 \dot{x}=g(x,t)

und hat damit wieder eine explizite gewöhnliche Differential-Gleichung.

Eine echte differentiell-algebraische Gleichung liegt dann vor, wenn die partielle Ableitung F_{\dot{x}} singulär ist. Dann zerfällt die implizite Differentialgleichung lokal in eine inhärente Differentialgleichung und eine algebraische Nebenbedingung. Dies entspricht praktisch einer Differentialgleichung, die auf einer Mannigfaltigkeit betrachtet wird. Das praktische Problem bei der impliziten Differentialgleichung ist jedoch, dass diese Mannigfaltigkeit zunächst nicht explizit bekannt ist.

Im Gegensatz zu gewöhnlichen Differentialgleichungen, deren Lösung durch Integration bestimmt wird, ergeben sich Teile der Lösung einer differentiell-algebraischen Gleichung durch Differentiation. Dies stellt weitere Anforderungen an die Systemfunktion F. Muss diese bei gewöhnlichen Differentialgleichungen nur stetig bzw. stetig differenzierbar sein, um die Lösbarkeit zu garantieren, so werden nun auch höhere Ableitungen für die Lösung benötigt. Die genaue Ordnung der benötigten Ableitungen hängt vom gewählten Lösungsansatz ab und wird allgemein als Index der differentiell-algebraischen Gleichung bezeichnet.

Durch die im Lösungsprozess hinzuzuziehenden Ableitungen von Komponenten des Gleichungssystems entsteht ein überbestimmtes System. Eine Folge davon ist, dass die Lösungen auch noch einer Anzahl expliziter oder impliziter algebraischer Nebenbedingungen genügen müssen. Insbesondere gilt dies für Anfangswerte von Anfangswertproblemen. Die Suche nach konsistenten Anfangswerten, z. B. in der Nähe vorgegebener inkonsistenter Anfangswerte, ist ein nichttriviales erstes Problem bei der praktischen Lösung von differentiell-algebraischen Gleichungen.

Typen differentiell-algebraischer Gleichungen[Bearbeiten]

Semi-explizite differentiell-algebraische Gleichung[Bearbeiten]

Ein spezieller Fall für eine differentiell-algebraische Gleichung ist ein System in der Form

 \dot{x}_1=f_1(x_1,x_2,t) \quad,\quad 0=f_2(x_1,x_2,t).

Durch Differenzieren der zweiten Differentialgleichung und Einsetzen der ersten erhält man als weitere Bedingung an eine Lösung

0=\partial_1f_2(x_1,x_2,t)f_1(x_1,x_2,t)+\partial_2f_2(x_1,x_2,t)\dot x_2+\partial_tf_2(x_1,x_2,t).

Ist der Faktor vor  \dot x_2 von Null verschieden, so ergibt sich ein explizites System gewöhnlicher Differentialgleichungen. Anfangswerte für dieses System müssen aber auch die undifferenzierte zweite Gleichung erfüllen, so dass nur ein Parameter frei gewählt werden kann.

Lineare differentiell-algebraische Gleichung[Bearbeiten]

Sehr häufig treten differentiell-algebraische Gleichungen auf in der Form

 E\dot{x}+Cx=q \,,\qquad x(t)\in\R^m \,,\quad t\in I\subset\R

mit stetigen Matrix-Koeffizienten

 E(t)\in L(\R^m,\R^k) \,,\quad C(t)\in L(\R^m,\R^k) \,,\quad t\in I

und einer Funktion

 q:I\rightarrow\R^k.

Eine echte differentiell-algebraische Gleichung liegt hier dann vor, wenn die Matrix-Funktion  E auf  I einen nichttrivialen Kern hat. Ein besonders einfacher Fall tritt ein, wenn die Matrizen quadratisch mit konstanten Einträgen sind.

Lineare differentiell-algebraische Gleichung mit proper formuliertem Hauptterm[Bearbeiten]

Eine andere Schreibweise für lineare differentiell-algebraische Gleichungen ist die Form

 A(Bx)'+Cx=q \,,\qquad x(t)\in\R^m \,,\quad t\in I\subset\R

mit (wenigstens) stetigen Matrix-Koeffizienten


  A(t)\in L(\R^n,\R^k)
    \,,\quad
  B(t)\in L(\R^m,\R^n)
    \,,\quad
  C(t)\in L(\R^m,\R^k)
    \,,\quad
  t\in I

und einer Funktion

 q:I\to\R^k.

In dieser Schreibweise wird der Tatsache Rechnung getragen, dass bei einer differentiell-algebraischen Gleichung nur ein Teil des Variablenvektors  x differenziert wird. Tatsächlich wird hier nur die Komponente  Bx differenziert und nicht der gesamte Variablenvektor  x . Als klassische Lösungen dieser Gleichung werden Funktionen aus dem Raum

 C_B^1(I,\R^m):=\left\{x\in C(I,\R^m)\mid Bx\in C^1(I,\R^n)\right\}

betrachtet, also dem Raum der stetigen Funktionen  x , für die die Komponente  Bx stetig differenzierbar ist.

Die beiden Matrix-Funktionen  A und  B bilden den Hauptterm der Gleichung und dieser heißt proper formuliert, wenn zwei Eigenschaften erfüllt sind:

1. Es gilt

 \ker A(t)\oplus\mbox{im}\,B(t)=\R^n \quad,\quad t\in I .

2. Es existiert eine stetig differenzierbare Projektor-Funktion

 R(t)\in L(\R^n,\R^n) \quad,\quad R^2(t)=R(t) \quad,\quad t\in I

mit der Eigenschaft

 \ker R(t)=\ker A(t) \quad,\quad \mbox{im}\,R(t)=\mbox{im}\,B(t) \quad,\quad t\in I .

Hier stellt die erste Bedingung sicher, dass zwischen den beiden Matrix-Funktionen A und B „nichts verloren geht“. Im Kern der Matrix  A kann nichts aus dem Bild der Matrix  B verschwinden. Die Projektor-Funktion  R realisiert genau die durch die Matrix-Funktionen  A und  B gegebene Zerlegung des Raumes  \R^n und ist für die Analyse der Gleichung hilfreich.

Ein einfacher Spezialfall für einen proper formulierten Hauptterm ist gegeben durch Matrix-Funktionen  A und  B mit der Eigenschaft

 \ker A(t)=\{0\} \quad,\quad \mbox{im}\,B(t)=\R^n \quad,\quad t\in I .

Für die Projektor-Funktion  R kann dann die Einheitsmatrix gewählt werden.

Indexbegriffe für DAEs[Bearbeiten]

Differentiationsindex[Bearbeiten]

Oftmals kann die Lösung eines Algebro-Differentialgleichungssystems durch (spezielle) Lösungskurven eines gewöhnlichen Differentialgleichungssystems dargestellt werden, obwohl f_{\dot x} singulär ist. Eine Schlüsselrolle nimmt hierbei der Differentiationsindex des Algebro-Differentialgleichungssystems ein.

Numerische Verfahren zur Lösung von Algebro-Differentialgleichungssystemen können meist nur Systeme integrieren, deren Differentiationsindex einen gewissen Maximalwert nicht überschreitet. So darf der Differentiationsindex des Systems beim impliziten Euler-Verfahren zum Beispiel nicht größer als eins sein.

Der Differentiationsindex eines Algebro-Differentialgleichungssystems

f(\dot{x},x,t) = 0

ist die Anzahl N\geq 0 der Zeitableitungen, die notwendig sind, um aus dem entstehenden Gleichungssystem


\begin{align}
f(\dot{x},x, t) &= 0\\
\frac{d}{dt}\left( f(\dot{x}, x, t)\right) &= 0\\
&\vdots\\
\frac{d^{N}}{dt^{N}}\left( f(\dot{x}, x, t)\right) &= 0
\end{align}

durch algebraische Umformungen ein gewöhnliches Differentialgleichungssystem

\dot x = g(x,t)

extrahieren zu können.

Beispiele[Bearbeiten]

Ein Algebro-Differentialgleichungssystem mit regulärer Matrix f_{\dot{x}}, das also algebraisch nach \dot{x} umgestellt werden kann, hat den Differentiationsindex null.

Eine rein algebraische Gleichung

F(x,t) = 0

mit regulärer Jacobi-Matrix F_{x}(x,t), die als Algebro-Differentialgleichung mit  f(\dot{x},x,t) = F(x,t) interpretiert wird, hat Differentiationsindex eins: Nach einmaligem Differenzieren erhält man die Gleichung

0=\frac{d}{dt} f(\dot{x},x,t) = F_{t}(x,t) +  F_{x}(x,t) \dot{x},

die nach \dot{x} auflösbar ist:

\dot{x} = -F_{x}(x,t)^{-1} F_{t}(x,t).

Diese Tatsache wird manchmal zur Konstruktion von Homotopieverfahren genutzt.

Die Euler-Lagrange-Gleichungen für das mathematische Pendel (mit auf eins normierter Erdbeschleunigung und Pendellänge) lauten


\begin{align}
\ddot{x}_{1} &=2x_{1}\lambda\\
\ddot{x}_{2} &=2x_{2}\lambda-1\\
0&=x_{1}^{2}+x_{2}^{2}-1
\end{align}

Dieses Algebro-Differentialgleichungssystem hat den Differentiationsindex drei: Zweifache Zeitableitung der Zwangsbedingung (dritte Gleichung) nach der Zeit liefert

2\left(x_{1}\ddot{x}_{1} + \dot{x}_{1}^{2} + x_{2}\ddot{x}_{2} + \dot{x}_{2}^{2}\right) = 0.

Mit Hilfe der zwei Differentialgleichungen in den Euler-Lagrange-Gleichungen lassen sich die zweiten zeitlichen Ableitungen \ddot{x}_1 und \ddot{x}_2 ersetzen, was

2\left( 2x_{1}^{2}\lambda + \dot{x}_{1}^{2} + 2x_{2}^{2}\lambda - x_{2}  + \dot{x}_{2}^{2} \right)=0

liefert. Mit x_{1}^{2}+x_{2}^{2}=1 erhält man daraus die Gleichung

\lambda = -\frac{1}{2}\left(\dot{x}_{1}^{2}+\dot{x}_{2}^{2}\right)+\frac{1}{2}x_{2}.

Durch Zeitableitung dieser Gleichung (das ist die dritte Zeitableitung) kommt man dann auf die fehlende Differentialgleichung für \lambda


\begin{align}
\dot\lambda &=
-\dot{x}_1\ddot{x}_1-\dot{x}_2\ddot{x}_2 +\frac{1}{2}\dot x_{2}\\
&= -2\lambda(\dot{x}_1 x_1 + \dot{x}_2 x_2 )-\frac{1}{2}\dot x_{2}\\
&= -\frac{1}{2}\dot x_{2}
\end{align},

wobei wieder die Differentialgleichungen aus den Euler-Lagrange-Gleichungen genutzt wurden, um \ddot{x}_{1} und \ddot{x}_{2} zu ersetzen und außerdem berücksichtigt wurde, dass 2(\dot{x}_1 x_1 + \dot{x}_2 x_2) = \frac{d}{dt}\left(x_{1}^{2}+x_{2}^{2}\right)=\frac{d}{dt}(1) = 0 gilt.

Geometrischer Index[Bearbeiten]

Ein mathematisch klar gefasster und geometrisch gut interpretierbarer Begriff ist der geometrische Index eines Algebro-Differentialgleichungssystems. Die Grundidee ist, dass man nach dem im Folgenden dargestellten iterativen Verfahren die maximale Zwangsmannigfaltigkeit ermittelt, auf der die Algebro-Differentialgleichung ein Vektorfeld (als Vektorfeld auf einer Mannigfaltigkeit) beschreibt. Der geometrische Index des Algebro-Differentialgleichungssystems ist dann die minimale Anzahl an Iterationsschritten, die bei diesem Verfahren benötigt wird.

Der geometrische Index ist gleich dem Differentiationsindex[1].

Gegeben sei eine autonome Algebro-Differentialgleichung

f(\dot{x}(t),x(t))=0

mit hinreichend oft differenzierbarer Funktion f:\R^n\times\R^n\rightarrow\R^n.

Im Rahmen des Algorithmus wird der \R^n als Mannigfaltigkeit M_0:=\R^n mit dem Tangentialbündel TM_0=\R^n\times\R^n interpretiert. Die Paare (x,y)\in TM_0 werden auch als Tangentialvektoren des \R^n bezeichnet.

Durch die Funktion f ist die Menge 
N:=\left\{(x,v)\in \R^n\times\R^n\mid f(v,x)=0\right\}
festgelegt, die jedem Punkt x\in\R^n alle für Lösungen des Algebro-DGL-Systems zulässigen Geschwindigkeitsvektoren v in diesem Punkt zugeordnet.

Es ist möglich, dass für ein Punkt x\in\R^n überhaupt keine Paar (x,v), genau ein solches Paar oder mehrere solcher Paare in N existieren.

Die Punkte, durch die eventuell Lösungen gehen können, erfasst man in der Menge

\left.M_1 := \mathrm{pr}_1 N,\right.

(mit der Projektion \mathrm{pr}_1 auf die erste Komponente, also \mathrm{pr}_1 N = \{x\mid (x,y)\in N\}). An dieser Stelle soll davon ausgegangen werden, dass M_1 eine differenzierbare Untermannigfaltigkeit des \R^n darstellt.

Jeder Tangentialvektor (x(t),\dot x(t)) an eine Lösung x:I\rightarrow\R^n der Algebro-Differentialgleichung muss auch im Tangentialbündel

TM_1 :=\left\{
(\bar{x}(0),\dot{\bar{x}}(0))\in\R^n\times\R^n
\mid
\bar{x}\in C^1((-\epsilon,\epsilon),M_1)
\mbox{ mit einem }\epsilon>0\right\}

von M_1 liegen (dabei bedeutet \bar{x}\in C^1((-\epsilon,\epsilon),M_1), dass \bar{x} eine auf einem Intervall (-\epsilon,\epsilon) definierte, einmal stetig differenzierbare Kurve ist, die vollständig in M_1 liegt).

Die Tangentialvektoren an Lösungen der Algebro-Differentialgleichung müssen auch in der Menge N\cap TM_1 und damit die Lösungen selber in der Menge M_2 := \mathrm{pr}_1(N\cap TM_1) liegen.

Diesen Prozess kann man (unter bestimmten Bedingungen) fortsetzen und aus der Zwangsmannigfaltigkeit M_k die Zwangsmannigfaltigkeit

M_{k+1} := \mathrm{pr}_1(N\cap TM_k)

bilden. Es ist möglich, dass ab einem k\in\{0,1,\ldots\} jedem Punkt x\in M_{k+1} in  N\cap TM_k genau ein Tangentialvektor (x,v) zugeordnet ist. Dann beschreibt N\cap TM_k ein Vektorfeld auf der Mannigfaltigkeit M_{k+1}.

Der geometrische Index der Algebro-Differentialgleichung ist gerade die minimale Zahl k\in\{0,1,\ldots\} für die N\cap TM_k ein Vektorfeld auf der Mannigfaltigkeit M_{k+1} beschreibt.

Beispiel

Die durch die Gleichung

f(\dot x,x) :=
\begin{pmatrix}
x_3\\
\dot x_1 - \cos(x_2)\\
\dot x_3 - \sin(x_2)
\end{pmatrix}

definierte Funktion und die zugehörige Algebro-Differentialgleichung dienen im folgenden Text als mitlaufendes Beispiel.

Im Beispiel gibt es für alle Punkte x\in\R^3, die nicht in der durch x_3=0 definierten Ebene liegen, keine Paare (x,v)\in N. Also verlaufen in diesem Beispiel außerhalb dieser Ebene keine Lösungen der Algebro-Differentialgleichung.

Es ergibt sich  M_1 = \{x\in\R^3\mid x_3 = 0\} und TM_1 = \{(x,v)\in\R^3\times\R^3\mid x_3 = v_3 = 0 \} und damit

N\cap TM_1 =\left\{(x,v)\in \R^3\times\R^3\mid x_3 = 0,\,v_1 = \cos(x_2),\, v_3 = \sin(x_2),\, v_3=0
\right\}.

Wie man sieht, liegt der durch N vorgegebene Tangentialvektor (x,v) (des \R^3) für Werte x_2\neq k\pi mit k\in\Z wegen v_3 = \sin(x_2) \neq 0 nicht im Tangentialraum TM_1, kann also nicht zu einer Lösung des Algebro-Differentialgleichungssystems korrespondieren. Damit ergibt sich

M_2 = \mathrm{pr}_1 (N\cap TM_1) = \left\{ x\in \R^3\mid
k\in\Z,\,
x_2 = k\pi, x_3 = 0
\right\}

Wir erhalten

TM_2 = \{(x,v)\in\R^n\times\R^n \mid
k\in\Z,\,x_2=k\pi,\,x_3=0,\,v_2=v_3=0 \}

und die Menge

N\cap TM_2 =\{
(x,v)\in\R^n\times\R^n
\mid
k\in\Z,\,x_2=k\pi,\,x_3=0,\,v_2=v_3=0,\,
v_1 = (-1)^k
\}

ordnet jedem Punkt x aus der Menge M_3 = \mathrm{pr}_1(N\cap TM_2) (die hier gerade gleich  M_2 ist) genau einen Tangentialvektor zu. Bei der Menge N\cap TM_1 ist das noch nicht der Fall, da bei Tangentialvektoren aus dieser Menge die Komponente v_2 noch nicht eingeschränkt ist.

Der geometrische Index des Algebro-Differentialgleichungssystems in diesem Beispiel ist also gleich zwei.

Ist M_1 eine Mannigfaltigkeit, so kann diese mit Hilfe einer Funktion g_1:\R^n\rightarrow\R^{m_1} in der Form

M_1 = \{x\in\R^n\mid g_1(x)=0\}

dargestellt werden. Die einschränkenden Gleichungen g_1(x)=0 in dieser Darstellung werden als Zwangsbedingungen der Algebro-Differentialgleichung bezeichnet.

Im Beispiel: g_1(x) = x_3.

Darüber hinaus kann für k=2,3,\ldots die Mannigfaltigkeit M_{k} mit Hilfe einer Funktion g_{k}:\R^n\rightarrow\R^{m_k} aus der Mannigfaltigkeit M_{k-1} ausgesondert werden: 
M_{k} = \{x \in M_{k-1}\mid g_k(x) = 0\}.
Die Gleichungen g_k(x)=0 mit k=2,3,\ldots werden auch als verdeckte Zwangsbedingungen der Algebro-Differentialgleichung bezeichnet (engl.: hidden constraints).

Im Beispiel: g_2(x) = \sin(x_2).

Bemerkungen:

  • Dass in diesem Abschnitt nur autonome Algebro-Differentialgleichungen betrachtet werden, erleichtert die geometrische Interpretation und ist nicht wirklich eine Einschränkung, da jede zeitabhängige Algebro-Differentialgleichung  f(\dot{x},x,t)=0 durch Einführen einer zusätzlichen Variable x_{n+1}:=t und einer zusätzlichen Differentialgleichung \dot x_{n+1} = 1 in eine autonome Algebro-Differentialgleichung umgeschrieben werden kann.
  • In diesem Abschnitt wurde vorausgesetzt, dass M_{k+1}:=\mathrm{pr}_1(N\cap TM_k) eine Untermannigfaltigkeit des \R^n ist. Falls dies nicht der Fall ist, ist für die betreffende Algebro-Differentialgleichung der geometrische Index nicht erklärt.
  • Es existieren auch Algebro-Differentialgleichungen bei denen der geometrische Index unendlich ist.

Konsistente Anfangswerte[Bearbeiten]

Gegeben sei wieder eine Algebro-Differentialgleichung

f(\dot{x},x,t) = 0

mit f:\R^n\times\R^n\times\R\rightarrow\R^n hinreichend oft differenzierbar.

Ein Punkt x_0\in\R^n heißt konsistenter Anfangswert zur Zeit t_0\in\R, falls es eine in einem offenen Intervall I\subset\R mit t_0\in I definierte Lösung x der Algebro-Differentialgleichung gibt, für die  x(t_0) = x_0 gilt.

Bei der Berechnung ist zu beachten, dass von konsistenten Anfangswerten außer den Zwangsbedingungen auch die verdeckten Zwangsbedingungen zu erfüllen sind (siehe Abschnitt Geometrischer Index).

Literatur[Bearbeiten]

  • Ernst Hairer und Gerhard Wanner: Solving Ordinary Differential Equations II, Stiff and Differential-Algebraic Problems. Second Revised Edition, Springer-Verlag, Berlin, 1996, ISBN 978-3-642-05220-0 (Print), ISBN 978-3-642-05221-7 (Online), DOI 10.1007/978-3-642-05221-7.
  • Uri M. Ascher und Linda R. Petzold: Computer Methods for Ordinary Differential equations and Differential-Algebraic equations. SIAM, Philadelphia, 1998, ISBN 0-89871-412-5.
  • Peter Kunkel und Volker Mehrmann: Differential-Algebraic Equations. EMS Textbooks in Mathematics, EMS Publishing House, Zürich, 2006, ISBN 3-03719-017-5, DOI 10.4171/017.
  • René Lamour, Roswitha März und Caren Tischendorf. Differential-Algebraic Equations: A Projector Based Analysis. Differential-Algebraic Equations Forum, Springer Berlin Heidelberg, 2013, ISBN 978-3-642-27554-8 (Print), ISBN 978-3-642-27555-5 (Online), DOI 10.1007/978-3-642-27555-5.

Einzelnachweise[Bearbeiten]

  1. G. Reißig: Beiträge zur Theorie und Anwendungen impliziter Differentialgleichungen. Dissertation, Dresdner Universitätsverlag, 1998.