aboutsummaryrefslogtreecommitdiffstats
path: root/buch/papers/kra/loesung.tex
blob: ef53adc21bdebc6d6bb75c011f877a840c138d48 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
\section{Lösungsmethoden} \label{kra:section:loesung}
\rhead{Lösungsmethoden}

\subsection{Riccatische Differentialgleichung} \label{kra:loesung:riccati}
Eine allgemeine analytische Lösung der Riccati Differentialgleichung ist nicht möglich.
Es gibt aber Spezialfälle, in denen sich die Gleichung vereinfachen lässt und so eine analytische Lösung gefunden werden kann.
Diese wollen wir im folgenden Abschnitt genauer anschauen.

\subsubsection{Fall 1: Konstante Koeffizienten}
Im Fall von konstanten Koeffizienten $f(x), g(x), h(x)$, wird die Gleichung \eqref{kra:equation:riccati} zu
\begin{equation}
    y' = fy^2 + gy + h.
\end{equation}
Durch Ausschreiben des Differentialquotienten
\begin{equation}
    \frac{dy}{dx} = fy^2 + gy + h
\end{equation}
erkennt man, dass die Differentialgleichung separierbar ist. Die Lösung findet man nun durch die Berechnung des Integrals
\begin{equation} \label{kra:equation:case1_int}
    \int \frac{dy}{fy^2 + gy + h} = \int dx.
\end{equation}

\subsubsection{Fall 2: Bekannte spezielle Lösung}
Kennt man eine spezielle Lösung $y_p$, so kann die riccatische Differentialgleichung mit Hilfe einer Substitution auf eine lineare Gleichung reduziert werden.
Wir wählen als Substitution
\begin{equation} \label{kra:equation:substitution}
    z = \frac{1}{y - y_p},
\end{equation}
durch Umstellen von \eqref{kra:equation:substitution} folgt
\begin{equation}
    y = y_p + \frac{1}{z^2} \label{kra:equation:backsubstitution}
\end{equation}
\begin{equation}
    y' = y_p' - \frac{1}{z^2}z',
\end{equation}
mit Einsetzten in die Differentialgleichung \eqref{kra:equation:riccati} resultiert
\begin{equation}
    y_p' - \frac{1}{z^2}z' = f(x)(y_p + \frac{1}{z}) + g(x)(y_p + \frac{1}{z})^2 + h(x)
\end{equation}
\begin{equation}
    -z^{2}y_p' + z' = -z^2\underbrace{(y_{p}f(x) + g(x)y_p^2 + h(x))}_{\displaystyle{y_p'}} - z(f(x) + 2y_{p}g(x)) - g(x)
\end{equation}
was uns direkt auf die lineare Differentialgleichung 1. Ordnung
\begin{equation}
    z' = -z(f(x) + 2y_{p}g(x)) - g(x)
\end{equation}
führt.
Diese kann nun mit den Methoden zur Lösung von linearen Differentialgleichungen 1. Ordnung gelöst werden.
Durch die Rücksubstitution \eqref{kra:equation:backsubstitution} erhält man dann die Lösung von \eqref{kra:equation:riccati}.

\subsection{Matrix-Riccati-Differentialgleichung} \label{kra:loesung:riccati}
Im Folgenden wollen wir uns anschauen wie die Matrix-Riccati-Differentialgleichung entsteht und wie sie gelöst werden kann.
Der Ausgangspunkt bildet die Matrix-Differentialgleichung
\begin{equation}
    \label{kra:equation:matrix-dgl}
    \begin{pmatrix}
        \dot{X}(t) \\
        \dot{Y}(t)
    \end{pmatrix}
    =
    \underbrace{
        \begin{pmatrix}
            A & B \\
            C & D
        \end{pmatrix}
    }_{\displaystyle{H}},
\end{equation}
mit den allgemeinen quadratischen Matrizen $A, B, C$ und $D$ welche zusammen die sogenannte Hamiltonsche-Matrix bilden.
Betrachten wir das Verhältniss von $Y$ zu $X$
\[
    U(t) = Y(t)X(t)^{-1}
\]
und deren Ableitung $\dot{P}(t)$, so erhalten wir die Riccati-Matrix-Differentialgleichung
\[
    \dot{U}(t) = C(t)  + DU(t) - U(t)A - U(t)BU(t).
\]

Die Lösung erhalten wir dann mit
\begin{equation}
    \label{kra:matrixriccati-solution}
    \begin{pmatrix}
        X(t) \\
        Y(t)
    \end{pmatrix}
    =
    \Phi(t_0, t)
    \begin{pmatrix}
        I(t) \\
        P_0(t)
    \end{pmatrix}
    =
    \begin{pmatrix}
        \Phi_{11}(t_0, t) & \Phi_{12}(t_0, t) \\
        \Phi_{21}(t_0, t) & \Phi_{22}(t_0, t)
    \end{pmatrix}
    \begin{pmatrix}
        I(t) \\
        P_0(t)
    \end{pmatrix}
\end{equation}
\begin{equation}
    P(t) =
    \begin{pmatrix}
        \Phi_{21}(t_0, t) + \Phi_{22}(t_0, t)
    \end{pmatrix}
    \begin{pmatrix}
        \Phi_{11}(t_0, t) + \Phi_{12}(t_0, t)
    \end{pmatrix}
    ^{-1}
\end{equation}
wobei $\Phi(t_0, t) = e^{H(t - t_0)}$ die sogenannte Zustandsübergangsmatrix von \eqref{kra:equation:matrix-dgl} ist \cite{kra:kalmanisae}.