8 files changed, 256 insertions, 72 deletions
diff --git a/buch/papers/erdbeben/Teil_Fabio.tex b/buch/papers/erdbeben/Teil_Fabio.tex
index afb2110..2ab6052 100644
--- a/buch/papers/erdbeben/Teil_Fabio.tex
+++ b/buch/papers/erdbeben/Teil_Fabio.tex
@@ -1,73 +1,51 @@
 \section{Anwendung des Kalman-Filters}
 \subsection{Ziel}
 Bis jetzt haben wir gelesen, was das Kalman-Filter bewirkt und wie es funktioniert.
-Nun möchten wir mit einem Beispiel herausfinden, ob das Filter unsere gesuchte Grösse $f(t)$ bestimmen kann.
+Nun möchten wir mit einem Beispiel herausfinden,
+ob das Filter unsere gesuchte Grösse $f(t)$ bestimmen kann.
 
 \subsection{Künstliche Erdbebendaten}
-Da wir keine Rohdaten über vergangene Erdbeben zur Hand haben, müssen wir mittels Matlab künstliche Daten  erzeugen und sie dann in das Filter eingeben.
+Da wir keine Rohdaten über vergangene Erdbeben zur Hand haben, müssen wir mittels Matlab künstliche Daten erzeugen und sie dann in das Filter eingeben.
 Diese Vorgehensweise erlaubt uns das Erdbeben beliebig zu gestalten
-und weil es digital simuliert wird, haben wir keine Bauschäden zu beklagen.
+und weil es digital simuliert wird,
+haben wir keine Bauschäden zu beklagen.
 
 \subsection{Wahl der Schwingung}
 Wir müssen uns überlegen, mit welcher Schwingung wir ein realitätsnahes Beben erzeugen können.
-
 Mit einer ungedämpften harmonischen Schwingung können wir zwar die meisten Vorgänge in der Physik erklären.
 Da aber unser Erdbeben irgendwann abklingen muss, wählen wir die gedämpfte harmonische Schwingung.
 Die dazugehörige Schwingungsgleichung lautet
-
 \begin{equation}
-	y = A e^{-\lambda t} \sin(\omega t)
+	y = A e^{-\lambda t} \sin(\omega t).
 \end{equation}
-
-Für die Variablen der harmonisch gedämpften Schwingung setzen wir die Werte
-
-\begin{equation}
-A = 5
-\end{equation}
-
-ein.
-
-$A$ ist die Amplitude der Schwingung, die uns die Heftigkeit des Erdebebens beschreibt.
+Dabei ist $A=5$ die anfängliche Amplitude der Schwingung,
+die uns die Heftigkeit des Erdebebens beschreibt.
 Sie ist vergleichbar mit der Magnitude.
-
-$\omega$ definiert sich durch
-
+$\lambda$ bezeichnet die Bodendämpfung, für die wir $0.2$ wählen.
+Sie ist dafür verantwortlich, dass unser Erdbeben abklingt
+und kreiert die bei gedämpften Schwingungen typische Hüllkurve.
+Wir nehmen an, dass $\lambda$ ein Materialparameter von geologischen Böden ist.
+Die Kreisfrequenz $\omega$ ist durch
 \begin{equation}
 	\omega = 2 \pi f
 \end{equation}
-
-wobei die Frequenz $f$ mit
-
+gegeben, 
+wobei die Momentanfrequenz $f = \mathcal N(\mu_f, \sigma_f) $ einer Normalverteilung mit
 \begin{equation}
-	f = E(\mathrm{Frequenz}) + \sigma^2(\mathrm{Frequenz})
+  \mu_f = \SI{15}{\hertz}
+  \qquad \text{und} \qquad
+  \sigma_f = \SI{10}{\hertz}
 \end{equation}
+folgt.
 
-erzeugt wird.
-
-Zusätzlich haben wir $f$ mit dem Savitzky-Golay-Filter gefiltert.
+Zusätzlich haben wir $f$ mit einem Savitzky-Golay-Filter gefiltert.
 Das Savitzky-Golay-Filter schaut sich immer eine definierte Anzahl von Datenpunkte an
-und bildet ein Polynom $n$-ter Ordnung.
-In unserer Anwendung schaut sich das Filter, im Sinne eines verschieblichen Fensters,
-jeweils zehn aufeinanderfolgende Datenpunkte an und bildet ein Polynom $0$-ter Ordnung.
-Da wir den Grad $0$ gewählt haben, erhalten wir pro zehn Punkte eine Gerade mit der Steigung $0$.
-Diese Art von der Filterung nennt sich gleitender Mittelwert.
-
-Für den Erwartungswert und die Standardabweichung setzen wir die Zahlen
-
-\begin{equation}
-E(f) = \SI{15}{\hertz}
-\end{equation}
-
-und
-\begin{equation}
-\sigma^2 = \SI{10}{\hertz}
-\end{equation}
-
-ein.
-
-$\lambda$ ist die Bodendämpfung, für die wir $0.2$ wählen.
-Sie ist dafür verantwortlich, dass unser Erdbeben abklingen wird und kreiert bei der gedämpften Schwingung die typische Hüllkurve der Amplitude.
-Wir nehmen an, dass $\lambda$ ein Materialparameter von geologischen Böden ist.
+und bildet  darüber ein Polynom $n$-ter Ordnung.
+In unserer Anwendung schaut sich das Filter, im Sinne eines verschiebbaren Fensters,
+jeweils elf aufeinanderfolgende Datenpunkte an
+und approximiert diese mit ein Polynom $0$-ter Ordnung,
+also einer Konstanten.
+Somit erhalten wir mit Matlab-Standardfunktionen einen gleitenden Mittelwert.
 
 \subsection{Versuch im Standardfall}
 Im nächsten Schritt müssen wir sinnvolle Systemparameter für unseren Seismographen definieren.
@@ -90,34 +68,37 @@ Wir nehmen an, dass
 		{0.00001}^2& 0& 0 \\
 		0 & {0.00001}^2& 0\\
 		0 & 0& {1 }^2\\
-	\end{array}\right)
+	\end{array}\right).
 \end{equation}
 
 Auch für die Messung setzen wir ein Rauschen voraus und definieren
 
 \begin{equation}
 R= ({\sigma_x}^2)=
-({0.00001}^2)
+({0.00001}^2).
 \end{equation}
 
-Sind nun die benötigten Systemparameter und das Rauschen definiert, erzeugen wir das Erdbeben und schauen, wie gut das Kalman-Filter die äussere Beschleunigung schätzen kann.
-
-\subsection*{Ergebnis}
+Sind nun die benötigten Systemparameter und Varianzen definiert,
+erzeugen wir ein Erdbeben mittels Simulation und schauen,
+wie gut das Kalman-Filter die äussere Beschleunigung schätzen kann.
 
-Wie wir in Abbildung~\ref{erdbeben:fig:standard-alles} im Positions-Zeit-Diagramm sehen, erzeugen unsere vorher gewählten Parameter eine realistische Erdbebenaufzeichnung.
-Leiten wir die Position einmal ab, erhalten wir die Geschwindigkeit.
-Die zweite Ableitung ergibt uns die Kraft, die in unserer Aufgabenstellung gesucht ist.
+\subsubsection{Ergebnis}
 
-Zoomen wir näher ran, erkennen wir in der Abbildung~\ref{erdbeben:fig:standard-zoom} im Positions-Diagramm eine Überlagerung der Massen-Eigenschwingung mit der Erdbebenschwingung.
-Die Masse schwingt mit einer tiefer Frequenz und hoher Amplitude, hingegen das Erdbeben mit einer hohen Frequenz und tiefer Amplitude.
+Abbildung~\ref{erdbeben:fig:standard-alles} zeigt zuoberst unsere Messwerte,
+also die Position der Masse relativ zum Seismografen.
+Wir sehen, dass unsere vorher gewählten Parameter eine realistische Erdbebenaufzeichnung erzeugen.
+Leiten wir die Position einmal ab, erhalten wir die Geschwindigkeit,
+und die zweite Ableitung ergibt uns die Kraft, welche für unsere Aufgabenstellung relevant ist.
+Sehr gut ersichtlich ist die typische Hüllkurve, wie wir sie bei einer gedämpften Schwingung erwarten.
 
-Vergleichen wir nun die Position mit der Kraft, stellen wir fest, dass das Kalman-Filter eine Schätzung wiedergibt, die auch eine Frequenz von \SI{15}{\hertz} hat.
-Das Filter war imstande die Eigenfrequenz zu eliminieren und die tatsächliche Kraft des Erdbebens zu wiedergeben.
+Die blaue Kurve ist die geschätzte äussere Kraft des Kalman-Filters.
+Erst wenn wir näher zoomen, erkennen wir in Abbildung~\ref{erdbeben:fig:standard-zoom},
+wie nahe die Schätzung an der idealen Schwingung liegt.
 
 \begin{figure}
 	\begin{center}
 		\includegraphics[width=\linewidth,keepaspectratio]{papers/erdbeben/Standard_alles.PNG}
-		\caption{Das Position-Zeit-Diagramm zeigt uns die typische Aufzeichnung eines Seismographen während eines Erdbebens. Um die Geschwindigkeit zu erhalten müssen wir die Position einmal ableiten. Ein weiteres Ableiten erzeugt uns die Beschleunigung, respektive die Kraft. Sehr gut ersichtlich ist die Hüllkurve der Amplitude, wie wir sie bei einer gedämpften Schwingung erwarten.}
+		\caption{Das Position-Zeit-Diagramm zeigt uns die typische Aufzeichnung eines Seismographen während eines Erdbebens. Um die Geschwindigkeit zu erhalten müssen wir die Position einmal ableiten. Ein weiteres Ableiten erzeugt uns die Beschleunigung, respektive die Kraft.}
     \label{erdbeben:fig:standard-alles}
 	\end{center}
 \end{figure}
@@ -131,21 +112,25 @@ Das Filter war imstande die Eigenfrequenz zu eliminieren und die tatsächliche K
 \end{figure}
 
 \subsection{Veränderung der Systemparameter}
-Was wir nun austesten möchten, sind die Auswirkungen wenn z.B. der Seismograph andere Systemparameter aufweist.
-Wir nehmen an, dass sich im Vergleich zum Standardfall die Masse erhöht, die Federkonstante schwächer und die Bodendämpfung doppelt so stark wirkt.
+Was wir nun testen möchten, sind die Auswirkungen wenn zum Beispiel der Seismograph andere Systemparameter aufweist.
+Wir nehmen an, dass sich im Vergleich zum Standardfall die Masse erhöht, die Federkonstante schwächer und die Federdämpfung doppelt so stark wirkt.
 Somit gilt neu
 \[
-m = 0.05
-\qquad \qquad
+m = 0.05,
+\qquad
 D = 0.5
 \qquad \text{und} \qquad
 k = 0.02.
 \]
 
-Da wir mit dieser Anpassung die Trägheit des Seismogrammes erhöht haben, erwarten wir sicher eine langsamere Bewegung der Masse, das heisst die Frequenz wird sich reduzieren.
+Da wir mit dieser Anpassung die Trägheit des Seismogrammes erhöht haben,
+erwarten wir eine langsamere Bewegung der Masse,
+das heisst die Frequenz wird kleiner.
 
-Betrachten wir die Abbildung~\ref{erdbeben:fig:systemparameter-geaendert} können wir diese Erwartung bestätigen.
-Nebst dem bemerken wir eine grössere Auslenkung der Position, die wir auf die höhere Energie der Masse und geringeren Rücklenkkraft der Feder begründen können.
+Betrachten wir Abbildung~\ref{erdbeben:fig:systemparameter-geaendert},
+können wir diese Erwartung bestätigen.
+Zudem bemerken wir eine grössere Auslenkung der Position,
+die wir mir durch die höhere Energie der Masse und geringeren Rücklenkkraft der Feder erklären können.
 
 \begin{figure}
 	\begin{center}
@@ -181,7 +166,7 @@ Die Auswertung in Abbildung~\ref{erdbeben:fig:prozessrauschen-geaendert} zeigt a
 \subsection{Verstärkung des Messrauschens}
 Als letztes verstärken wir das Messrauschen um den Faktor $100$ und belassen wieder den Rest wie im Standardfall.
 Wie man eigentlich schon erwarten kann, zeigt uns die Abbildung~\ref{erdbeben:fig:messrauschen-geaendert}, dass das Signal des Messsensors vom Messrauschen gestört wird.
-Weil die Messung somit ungenau wird, kann das Kalman-Filter nicht mehr genau arbeiten und produziert einen ungenauen Output.
+Weil die Messung somit ungenau wird, kann das Kalman-Filter nicht mehr genau arbeiten und produziert eine ungenaues Resultat.
 
 \begin{figure}
 	\begin{center}
@@ -207,8 +192,6 @@ Der Seismograph war fähig die Position der Masse während der Einwirkung des Er
 $a(t)$ kann zwar nicht mit Sensoren gemessen werden, jedoch erhalten wir $a(t)$ durch zweifaches Ableiten.
 Da wir so aber die innere Beschleunigung erhalten, mussten wir das Kalman-Filter anwenden.
 Das Kalman-Filter half uns die äussere Beschleunigung zu schätzen und lieferte erstaunlich genaue Werte.
-Ausserdem hat es das Filter geschafft, die Eigenfrequenz der Masse und die Erdbebenfrequenz zu separieren.
-Folglich erhielten wir eine Schätzung, die nur das Erdbeben betraf.
 
-Zuletzt haben wir aufgezeigt, das Veränderungen an den System- und Rauschparametern die Genauigkeit und Zuverlässigkeit des Kalman-Filters beeinträchtigen können.
+Schlussendlich haben wir aufgezeigt, das Veränderungen an den System- und Rauschparametern die Genauigkeit und Zuverlässigkeit des Kalman-Filters beeinträchtigen.
 
diff --git a/buch/papers/erdbeben/images/Makefile b/buch/papers/erdbeben/images/Makefile
new file mode 100644
index 0000000..3c82cb7
--- /dev/null
+++ b/buch/papers/erdbeben/images/Makefile
@@ -0,0 +1,15 @@
+#
+# Makefile
+#
+# (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule
+#
+all: standard.pdf messrauschen_geaendert.pdf prozessrauschen_geaendert.pdf
+
+standard.pdf: standard.tex
+	pdflatex standard.tex
+
+messrauschen_geaendert.pdf: messrauschen_geaendert.tex
+	pdflatex messrauschen_geaendert.tex
+
+prozessrauschen_geaendert.pdf: prozessrauschen_geaendert.tex
+	pdflatex prozessrauschen_geaendert.tex
diff --git a/buch/papers/erdbeben/images/messrauschen_geaendert.pdf b/buch/papers/erdbeben/images/messrauschen_geaendert.pdf
new file mode 100644
index 0000000..5056dd0
--- /dev/null
+++ b/buch/papers/erdbeben/images/messrauschen_geaendert.pdf
diff --git a/buch/papers/erdbeben/images/messrauschen_geaendert.tex b/buch/papers/erdbeben/images/messrauschen_geaendert.tex
new file mode 100644
index 0000000..a5a7509
--- /dev/null
+++ b/buch/papers/erdbeben/images/messrauschen_geaendert.tex
@@ -0,0 +1,61 @@
+%
+% messrauschen_geaendert.tex -- Kombination der Messrauschen-Bilder
+%
+% (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule
+%
+\documentclass[tikz]{standalone}
+\usepackage{amsmath}
+\usepackage{times}
+\usepackage{txfonts}
+\usepackage{pgfplots}
+\usepackage{csvsimple}
+\usetikzlibrary{arrows,intersections,math}
+\begin{document}
+\newboolean{showgrid}
+\setboolean{showgrid}{false}
+\def\skala{1}
+\begin{tikzpicture}[>=latex,thick,scale=\skala]
+
+\node at (0,0) {\includegraphics[width=14cm]{../Messrauschen_geaendert.PNG}};
+
+\def\yten{-2.06}
+\def\yfifteen{0.16}
+
+\def\links{12.2}
+\def\rechts{14.2}
+
+\pgfmathparse{(\yfifteen-(\yten))/5}
+\xdef\m{\pgfmathresult}
+\xdef\b{\yten}
+
+\pgfmathparse{\m*(\links-10)+(\b)}
+\xdef\Links{\pgfmathresult}
+
+\pgfmathparse{\m*(\rechts-10)+(\b)}
+\xdef\Rechts{\pgfmathresult}
+
+\begin{scope}[yshift=-9cm]
+\node at (0,0) {\includegraphics[width=14cm]{../Messrauschen_geaendert_zoom.PNG}};
+\end{scope}
+
+\def\breite{7}
+\def\hoehe{2}
+
+% Gitter
+\ifthenelse{\boolean{showgrid}}{
+\draw[step=0.1,line width=0.1pt] (-\breite,-\hoehe) grid (\breite, \hoehe);
+\draw[step=0.5,line width=0.4pt] (-\breite,-\hoehe) grid (\breite, \hoehe);
+\draw                            (-\breite,-\hoehe) grid (\breite, \hoehe);
+\fill (0,0) circle[radius=0.05];
+}{}
+
+\def\unten{-4.15}
+
+\draw[line width=0.7pt] (\Links,\unten) rectangle (\Rechts,4.15);
+\draw[line width=0.7pt] (\Links,\unten) -- (\Links,{\unten-0.05})
+	-- (-6.45,-4.6) -- (-6.45,-4.9);
+\draw[line width=0.7pt] (\Rechts,\unten) -- (\Rechts,{\unten-0.05})
+	-- (6.86,-4.6) -- (6.86,-4.9);
+\end{tikzpicture}
+\end{document}
+
diff --git a/buch/papers/erdbeben/images/prozessrauschen_geaendert.pdf b/buch/papers/erdbeben/images/prozessrauschen_geaendert.pdf
new file mode 100644
index 0000000..e0bf605
--- /dev/null
+++ b/buch/papers/erdbeben/images/prozessrauschen_geaendert.pdf
diff --git a/buch/papers/erdbeben/images/prozessrauschen_geaendert.tex b/buch/papers/erdbeben/images/prozessrauschen_geaendert.tex
new file mode 100644
index 0000000..ad4dcf9
--- /dev/null
+++ b/buch/papers/erdbeben/images/prozessrauschen_geaendert.tex
@@ -0,0 +1,61 @@
+%
+% prozessrauschen_geaendert.tex -- Kombination der Prozessrauschen-Bilder
+%
+% (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule
+%
+\documentclass[tikz]{standalone}
+\usepackage{amsmath}
+\usepackage{times}
+\usepackage{txfonts}
+\usepackage{pgfplots}
+\usepackage{csvsimple}
+\usetikzlibrary{arrows,intersections,math}
+\begin{document}
+\newboolean{showgrid}
+\setboolean{showgrid}{false}
+\def\skala{1}
+\begin{tikzpicture}[>=latex,thick,scale=\skala]
+
+\node at (0,0) {\includegraphics[width=14cm]{../Prozessrauschen_geaendert.PNG}};
+
+\def\yten{-2.1}
+\def\yfifteen{0.12}
+
+\def\links{13.27}
+\def\rechts{14.2}
+
+\pgfmathparse{(\yfifteen-(\yten))/5}
+\xdef\m{\pgfmathresult}
+\xdef\b{\yten}
+
+\pgfmathparse{\m*(\links-10)+(\b)}
+\xdef\Links{\pgfmathresult}
+
+\pgfmathparse{\m*(\rechts-10)+(\b)}
+\xdef\Rechts{\pgfmathresult}
+
+\begin{scope}[yshift=-9cm]
+\node at (0,0) {\includegraphics[width=14cm]{../Prozessrauschen_geaendert_zoom.PNG}};
+\end{scope}
+
+% Gitter
+\def\breite{7}
+\def\hoehe{2}
+\ifthenelse{\boolean{showgrid}}{
+\draw[step=0.1,line width=0.1pt] (-\breite,-\hoehe) grid (\breite, \hoehe);
+\draw[step=0.5,line width=0.4pt] (-\breite,-\hoehe) grid (\breite, \hoehe);
+\draw                            (-\breite,-\hoehe) grid (\breite, \hoehe);
+\fill (0,0) circle[radius=0.05];
+}{}
+
+\def\unten{-4.15}
+
+\draw[line width=0.7pt] (\Links,\unten) rectangle (\Rechts,4.18);
+\draw[line width=0.7pt] (\Links,\unten) -- (\Links,{\unten-0.05})
+	-- (-6.62,-4.6) -- (-6.62,-4.9);
+\draw[line width=0.7pt] (\Rechts,\unten) -- (\Rechts,{\unten-0.05})
+	-- (6.80,-4.6) -- (6.80,-4.9);
+
+\end{tikzpicture}
+\end{document}
+
diff --git a/buch/papers/erdbeben/images/standard.pdf b/buch/papers/erdbeben/images/standard.pdf
new file mode 100644
index 0000000..f2ca85e
--- /dev/null
+++ b/buch/papers/erdbeben/images/standard.pdf
diff --git a/buch/papers/erdbeben/images/standard.tex b/buch/papers/erdbeben/images/standard.tex
new file mode 100644
index 0000000..74ac7a1
--- /dev/null
+++ b/buch/papers/erdbeben/images/standard.tex
@@ -0,0 +1,64 @@
+%
+% standard.tex -- Kombination der Standard-Bilder
+%
+% (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule
+%
+\documentclass[tikz]{standalone}
+\usepackage{amsmath}
+\usepackage{times}
+\usepackage{txfonts}
+\usepackage{pgfplots}
+\usepackage{csvsimple}
+\usetikzlibrary{arrows,intersections,math}
+\begin{document}
+\newboolean{showgrid}
+\setboolean{showgrid}{false}
+\def\skala{1}
+\begin{tikzpicture}[>=latex,thick,scale=\skala]
+
+\node at (0,0) {\includegraphics[width=14cm]{../Standard_alles.PNG}};
+
+\def\yten{-2.12}
+\def\yfifteen{0.1}
+
+\def\links{12.5}
+\def\rechts{13.4}
+
+\pgfmathparse{(\yfifteen-(\yten))/5}
+\xdef\m{\pgfmathresult}
+\xdef\b{\yten}
+
+%\node at (0,7) {$m=\m$}; \node at (0,8) {$b=\b$};
+
+\pgfmathparse{\m*(\links-10)+(\b)}
+\xdef\Links{\pgfmathresult}
+
+\pgfmathparse{\m*(\rechts-10)+(\b)}
+\xdef\Rechts{\pgfmathresult}
+
+\begin{scope}[yshift=-9cm]
+\node at (0,0) {\includegraphics[width=14cm]{../Standard_Zoom.PNG}};
+\end{scope}
+
+\def\breite{7}
+\def\hoehe{2}
+
+% Gitter
+\ifthenelse{\boolean{showgrid}}{
+\draw[step=0.1,line width=0.1pt] (-\breite,-\hoehe) grid (\breite, \hoehe);
+\draw[step=0.5,line width=0.4pt] (-\breite,-\hoehe) grid (\breite, \hoehe);
+\draw                            (-\breite,-\hoehe) grid (\breite, \hoehe);
+\fill (0,0) circle[radius=0.05];
+}{}
+
+\def\unten{-4.15}
+
+\draw[line width=0.7pt] (\Links,\unten) rectangle (\Rechts,4.18);
+\draw[line width=0.7pt] (\Links,\unten) -- (\Links,{\unten-0.05})
+	-- (-6.36,-4.6) -- (-6.36,-4.9);
+\draw[line width=0.7pt] (\Rechts,\unten) -- (\Rechts,{\unten-0.05})
+	-- (6.87,-4.6) -- (6.87,-4.9);
+
+\end{tikzpicture}
+\end{document}
+