Riccatische Differentialgleichung

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

Eine riccatische Differentialgleichung oder Riccati-Differentialgleichung ist eine spezielle Klasse nichtlinearer gewöhnlicher Differentialgleichungen erster Ordnung. Sie besitzen die Form

y'(x) = f(x)y^2(x)+ g(x)y(x) + h(x)

mit gegebenen Funktionen f, g und h.

Sie ist nach dem Mathematiker Jacopo Francesco Riccati benannt, einem italienischen Grafen (1676–1754), der sich intensiv mit der Klassifizierung von Differentialgleichungen befasste und Methoden zur Erniedrigung der Ordnung von Gleichungen entwickelte.

Eine allgemeine Lösung einer Riccati-Differentialgleichung ist im Allgemeinen nicht möglich, jedoch kann eine solche angegeben werden, falls eine spezielle Lösung bekannt ist.

Denselben Namen riccatische Differentialgleichung tragen noch zwei andere Gleichungstypen, die für verschiedene Themen von angewandter Mathematik bis zur Finanzwissenschaft von Bedeutung sind.

Transformation im Falle einer bekannten Lösung[Bearbeiten]

Angenommen, man hätte bereits eine Lösung u (etwa durch Raten) gefunden. Dann lässt sich die riccatische Differentialgleichung vollständig lösen, da das Auffinden der übrigen Lösungen sich nun auf eine bernoullische Differentialgleichung reduziert, welche leicht gelöst werden kann.[1]

Formulierung des Transformationssatzes[Bearbeiten]

Es seien x_0 \in (a, b) sowie u: (a, b) \rightarrow \mathbb{R} eine Lösung der riccatischen Differentialgleichung

\ u'(x) = f(x)u^2(x) + g(x)u(x) + h(x)

und z eine Lösung der bernoullischen Differentialgleichung

z'(x)=f(x)z(x)^2 + (2u(x)f(x) + g(x))z(x)\ .

Dann ist

\ y(x) := z(x) + u(x)

die Lösung der riccatischen Differentialgleichung

y'(x) = f(x)y^2(x) + g(x)y(x) + h(x)\ ,\ 
y(x_0) = y_0 := z(x_0) + u(x_0)\ .

Beweis[Bearbeiten]

Es gilt

\begin{array}{lll}
y'(x)&=&z'(x) + u'(x)\\
&=&f(x)z(x)^2 + (2u(x)f(x) + g(x))z(x) + f(x)u^2(x) + g(x)u(x) + h(x)\\
&=&f(x)[z(x)^2 + 2u(x)z(x) + u^2(x)] + g(x)[z(x)+u(x)] + h(x)\\
&=&f(x)y^2(x) + g(x)y(x) + h(x)\ ,\\
\end{array}

während der Anfangswert trivialerweise erfüllt ist.

Umformung auf lineare Differentialgleichung zweiter Ordnung[Bearbeiten]

Im Allgemeinen, unabhängig davon, ob man eine spezielle Lösung gefunden hat, lässt sich die riccatische Differentialgleichung auf eine lineare Differentialgleichung zweiter Ordnung mit nicht-konstanten Koeffizienten transformieren.[2] Sollten zufälligerweise die Koeffizienten konstant sein, lässt sich diese transformierte Gleichung mit Hilfe der charakteristischen Gleichung leicht vollständig lösen. Im Fall nicht-konstanter Koeffizienten kann auch die lineare Form der riccatischen Differentialgleichung nur sehr schwer lösbar sein.

Formulierung des Transformationssatzes[Bearbeiten]

Es seien x_0 \in (a, b) sowie f: (a,b) \rightarrow \mathbb{R}\setminus\{0\} stetig differenzierbar und z eine Lösung der linearen Differentialgleichung zweiter Ordnung

z''(x) - \left[g(x) + \frac{f'(x)}{f(x)}\right] \cdot z'(x) + [f(x)h(x)] \cdot z(x) = 0

mit z(x) \neq 0 für alle x \in (a,b). Dann ist

y(x) := - \frac{z'(x)}{f(x)z(x)}

die Lösung der riccatischen Differentialgleichung

y'(x) = f(x)y^2(x) + g(x)y(x) + h(x)\ ,\ 
y(x_0) = y_0 := -\frac{z'(x_0)}{f(x_0)z(x_0)}\ .

Beweis[Bearbeiten]

Der Übersicht halber werden die Argumente nicht mitgeschrieben. Nach der Quotientenregel gilt

y'=\frac{-fzz''+z'(f'z+fz')}{f^2z^2} 
=\frac{-fz[(g+\frac{f'}{f})z' - fhz] + z'(f'z+fz')}{f^2z^2}
=f\frac{z'^2}{f^2z^2} - g\frac{z'}{fz} + h=fy^2 + gy + h,

während der Anfangswert trivialerweise erfüllt ist.

Transformation auf ein lineares Differentialgleichungssystem[Bearbeiten]

Neben der Umformung der Riccati-Differentialgleichung auf eine lineare Differentialgleichung zweiter Ordnung ist auch eine Transformation auf ein lineares Differentialgleichungssystem möglich.[3] Damit eröffnen sich weitere Lösungsmöglichkeiten. So wird beispielsweise im Fall konstanter Koeffizienten mit der Matrixexponentialfunktion eine analytische Lösung der ursprünglichen Riccati-Differentialgleichung erhalten.

Die gewöhnliche Differentialgleichung zweiter Ordnung (Vgl. Formulierung des Transformationssatzes)

z''(x) - k_1 \cdot z'(x) + k_2 \cdot z(x) = 0

lässt sich mit der Substitution t_1=z und t_2=z' auf ein System gewöhnlicher Differentialgleichungen erster Ordnung zurückführen.[4] Es folgt das lineare homogene Differentialgleichungssystem

T^\prime(x)=A(x)\,T(x)

mit der Koeffizientenmatrix A(x) \in \R^{2\times 2}

A = \begin{pmatrix} 0 & 1 \\ k_2(x) & k_1(x) \end{pmatrix}

und dem Vektor T(x) \in \R^2

 T = \begin{pmatrix} t_1 \\ t_2  \end{pmatrix}.

Ist die Koeffizientenmatrix eine stetige Funktion von x (d. h. auch im Fall nicht-konstanter Koeffizienten), dann hat das zugehörige Anfangswertproblem mit den Anfangswerten (x_0,T_0) \in \R^{2+1} stets eine eindeutig bestimmte Lösung, die für alle x \in \R erklärt ist. Darüber hinaus kann die Lösung in Matrixschreibweise über die Matrixexponentialfunktion auch im Fall nicht-konstanter Koeffizienten angegeben werden.[5]

Literatur[Bearbeiten]

  • Harro Heuser: Gewöhnliche Differentialgleichungen. Teubner, 2004, ISBN 3-519-32227-7.
  • Wolfgang Walter: Gewöhnliche Differentialgleichungen. 7. Auflage. Springer, Berlin/Heidelberg 2000, ISBN 978-3-540-67642-3.

Einzelnachweise[Bearbeiten]

  1. Walter: Gewöhnliche Differentialgleichungen. 2000, S. 94.
  2. Walter: Gewöhnliche Differentialgleichungen. 2000, S. 305.
  3. Möller, T: Symbolic mathematics-based simulation of cylinder spaces for regenerative gas cycles. In: Int J Energy Environ Eng. Springer Berlin Heidelberg, Feb. 2015.
  4. Merziger, G., Wirth, Th.: Repetitorium der höheren Mathematik. Binomi Verlag, Hannover 2006.
  5. Blanes, S.; Casas, F.; Oteo, J.A.; Ros, J.: The Magnus expansion and some of its applications. Physics Reports Bd. 470, Cornell University Library, 2009.