From ded30e493c1b05f1f412f2e78636d7195ea054e0 Mon Sep 17 00:00:00 2001 From: Kuster Yanik Date: Thu, 4 Aug 2022 21:24:11 +0200 Subject: added new subsection wird das Ziel erreicht? --- buch/papers/lambertw/Bilder/Intuition.pdf | Bin 186972 -> 187016 bytes buch/papers/lambertw/Bilder/Strategie.py | 4 +- .../Bilder/lambertAbstandBauchgef\303\274hl.py" | 10 +- buch/papers/lambertw/teil0.tex | 5 +- buch/papers/lambertw/teil1.tex | 101 +++++++++++++++++---- 5 files changed, 93 insertions(+), 27 deletions(-) (limited to 'buch/papers/lambertw') diff --git a/buch/papers/lambertw/Bilder/Intuition.pdf b/buch/papers/lambertw/Bilder/Intuition.pdf index 236212a..739b02b 100644 Binary files a/buch/papers/lambertw/Bilder/Intuition.pdf and b/buch/papers/lambertw/Bilder/Intuition.pdf differ diff --git a/buch/papers/lambertw/Bilder/Strategie.py b/buch/papers/lambertw/Bilder/Strategie.py index d7d06cb..975e248 100644 --- a/buch/papers/lambertw/Bilder/Strategie.py +++ b/buch/papers/lambertw/Bilder/Strategie.py @@ -34,8 +34,8 @@ ax.quiver(X, Y, U, W, angles='xy', scale_units='xy', scale=1, headwidth=5, headl ax.plot([V[0], (VZ+V)[0]], [V[1], (VZ+V)[1]], 'k--') ax.plot(np.vstack([V, Z])[:, 0], np.vstack([V, Z])[:,1], 'bo', markersize=10) -ax.set_xlabel("x") -ax.set_ylabel("y") +ax.set_xlabel("x", size=20) +ax.set_ylabel("y", size=20) ax.text(2.5, 4.5, "Visierlinie", size=20, rotation=10) diff --git "a/buch/papers/lambertw/Bilder/lambertAbstandBauchgef\303\274hl.py" "b/buch/papers/lambertw/Bilder/lambertAbstandBauchgef\303\274hl.py" index 9031bfc..3a90afa 100644 --- "a/buch/papers/lambertw/Bilder/lambertAbstandBauchgef\303\274hl.py" +++ "b/buch/papers/lambertw/Bilder/lambertAbstandBauchgef\303\274hl.py" @@ -39,8 +39,8 @@ plt.plot(0, ymin, 'bo', markersize=10) plt.plot([0, xmin], [ymin, ymin], 'k--') #plt.xlim(-0.1, 1) #plt.ylim(1, 2) -#plt.ylabel("y") -#plt.xlabel("x") +plt.ylabel("y") +plt.xlabel("x") plt.grid(True) plt.quiver(xmin, ymin, -0.2, 0, scale=1) @@ -53,6 +53,6 @@ plt.rcParams.update({ "font.serif": ["New Century Schoolbook"], }) -plt.text(xmin-0.11, ymin-0.12, r"$\dot{v}$", size=30) -plt.text(xmin-0.02, ymin+0.05, r"$V$", size=30, c='b') -plt.text(0.02, ymin+0.05, r"$Z$", size=30, c='b') \ No newline at end of file +plt.text(xmin-0.11, ymin-0.08, r"$\dot{v}$", size=20) +plt.text(xmin-0.02, ymin+0.05, r"$V$", size=20, c='b') +plt.text(0.02, ymin+0.05, r"$Z$", size=20, c='b') \ No newline at end of file diff --git a/buch/papers/lambertw/teil0.tex b/buch/papers/lambertw/teil0.tex index 088cb7b..6632eca 100644 --- a/buch/papers/lambertw/teil0.tex +++ b/buch/papers/lambertw/teil0.tex @@ -74,7 +74,7 @@ darstellen. Der Geschwindigkeitsvektor muss auf das Ziel zeigen, woraus folgt z-v \text{.} \end{equation} -Um den Richtungsvektor zu konstruieren kann der Einheitsvektor parallel zu $z-v$ um $\dot{v}$ gestreckt werden, was zu +Um den Richtungsvektor zu konstruieren kann der Einheitsvektor parallel zu $z-v$ um $|\dot{v}|$ gestreckt werden, was zu \begin{equation} \dot{v} = @@ -86,6 +86,7 @@ führt. Dies kann noch ausgeschrieben werden zu = |\dot{v}|\cdot\frac{z-v}{|z-v|} \text{.} + \label{lambertw:richtungsvektor} \end{equation} % Aus dem Verfolgungsproblem ist auch ersichtlich, dass die Punkte $V$ und $Z$ nicht am gleichen Ort starten und so eine Division durch Null ausgeschlossen ist. @@ -105,7 +106,7 @@ was algebraisch zu 1 \end{align} umgeformt werden kann. -Die Lösungen dieser Differentialgleichung sind die gesuchten Verfolgungskurven, insofern der Verfolger die Jagdstrategie verwendet. +Die Lösungen dieser Differentialgleichung sind die gesuchten Verfolgungskurven, sofern der Verfolger die Jagdstrategie verwendet. % \subsection{Ziel \label{lambertw:subsection:Ziel}} diff --git a/buch/papers/lambertw/teil1.tex b/buch/papers/lambertw/teil1.tex index 0fd0108..e8eca2c 100644 --- a/buch/papers/lambertw/teil1.tex +++ b/buch/papers/lambertw/teil1.tex @@ -30,7 +30,7 @@ Die Parametrisierung von $z(t)$ ist im Beispiel definiert als \left( \begin{array}{c} 0 \\ t \end{array} \right)\text{.} \end{equation} % -Die Parametrisierung von $v(t)$ ist von den Startbedingungen abhängig. Deshalb wird \eqref{bedingung_treffer} jeweils für die unterschiedlichen Startbedingungen separat analysiert. +Die Parametrisierung von $v(t)$ ist von den Startbedingungen abhängig. Deshalb wird die Bedingung \eqref{bedingung_treffer} jeweils für die unterschiedlichen Startbedingungen separat analysiert. % \subsection{Anfangsbedingung im ersten Quadranten} % @@ -41,7 +41,7 @@ Wenn der Verfolger im ersten Quadranten startet, dann kann $v(t)$ mit den Gleich x_0\cdot\sqrt{\frac{1}{\chi}W\left(\chi\cdot \exp\left( \chi-\frac{4t}{r_0-y_0}\right) \right)} \\ y(t) &= - \frac{1}{4}\left(\left(y_0+r_0\right)\left(\frac{x(t)}{x_0}\right)^2+\left(r_0-y_0\right)\operatorname{ln}\left(\left(\frac{x(t)}{x_0}\right)^2\right)-r_0+3y_0\right)\\ + \frac{1}{4}\left(\left(y_0+r_0\right)\left(\frac{x(t)}{x_0}\right)^2+\left(y_0-r_0\right)\operatorname{ln}\left(\left(\frac{x(t)}{x_0}\right)^2\right)-r_0+3y_0\right)\\ \chi &= \frac{r_0+y_0}{r_0-y_0}, \quad @@ -54,7 +54,7 @@ Wenn der Verfolger im ersten Quadranten startet, dann kann $v(t)$ mit den Gleich \text{.} \end{align} % -Der Folger ist durch +Der Verfolger ist durch \begin{equation} v(t) = @@ -76,31 +76,37 @@ Daher wird das Problem in zwei einzelne Teilprobleme zerlegt, wodurch die Beding &= y(t) = - \frac{1}{4}\left(\left(y_0+r_0\right)\left(\frac{x(t)}{x_0}\right)^2+\left(r_0-y_0\right)\operatorname{ln}\left(\left(\frac{x(t)}{x_0}\right)^2\right)-r_0+3y_0\right)\text{,} + \frac{1}{4}\left(\left(y_0+r_0\right)\left(\frac{x(t)}{x_0}\right)^2+\left(y_0-r_0\right)\operatorname{ln}\left(\left(\frac{x(t)}{x_0}\right)^2\right)-r_0+3y_0\right)\text{,} \end{align} % -welche Beide gleichzeitig erfüllt sein müssen, damit das Ziel erreicht wurde. +welche beide gleichzeitig erfüllt sein müssen, damit das Ziel erreicht wurde. Zuerst wird die Bedingung der $x$-Koordinate betrachtet. -Da $x_0 \neq 0$ und $\chi \neq 0$ mit +Da $x_0 \neq 0$ und $\chi \neq 0$ kann \begin{equation} 0 = x_0\sqrt{\frac{1}{\chi}W\left(\chi\cdot \exp\left( \chi-\frac{4t}{r_0-y_0}\right)\right)} \end{equation} -ist diese Bedingung genau dann erfüllt, wenn +algebraisch zu \begin{equation} 0 = W\left(\chi\cdot \exp\left( \chi-\frac{4t}{r_0-y_0}\right)\right) - \text{.} \end{equation} -% +umgeformt werden. Es ist zu beachten, dass $W(x)$ die Lambert W-Funktion ist, welche im Kapitel \eqref{buch:section:lambertw} behandelt wurde. -Diese Gleichung entspricht genau den Nullstellen der Lambert W-Funktion. Da die Lambert W-Funktion genau eine Nullstelle bei -\begin{equation} +Diese Gleichung entspricht genau den Nullstellen der Lambert W-Funktion. Mit der einzigen Nullstelle der Lambert W-Funktion bei +\begin{equation*} W(0)=0 + \text{,} +\end{equation*} +kann die Bedingung weiter vereinfacht werden zu +\begin{equation} + 0 + = + \chi\cdot \exp\left( \chi-\frac{4t}{r_0-y_0}\right) + \text{.} \end{equation} -% Da $\chi\neq0$ und die Exponentialfunktion nie null sein kann, ist diese Bedingung unmöglich zu erfüllen. Beim Grenzwert für $t\rightarrow\infty$ geht die Exponentialfunktion gegen null. Dies nützt nicht viel, da unendlich viel Zeit vergehen müsste, damit ein Einholen möglich wäre. @@ -203,16 +209,18 @@ Durch quadrieren verschwindet die Wurzel des Betrages, womit % \begin{equation} |v-z|^2e_y\cdot v$. +Aus diesem Argument würde folgen, dass beim tiefsten Punkt der Verfolgungskurve im Beispiel den minimalen Abstand befindet. +% \begin{figure} \centering \includegraphics[scale=0.4]{./papers/lambertw/Bilder/Intuition.pdf} @@ -220,7 +228,8 @@ Es lässt sich vermuten, dass bei diesem Punkt der Abstand zum Ziel minimal sein \label{lambertw:grafic:intuition} \end{figure} % -Dies kann leicht überprüft werden, indem wir lokal alle relevanten benachbarten Punkte betrachten und das Vorzeichen der Änderung des Abstandes prüfen. + +Dieses Argument kann leicht überprüft werden, indem lokal alle relevanten benachbarten Punkte betrachtet und das Vorzeichen der Änderung des Abstandes überprüft wird. Dafür wird ein Ausdruck benötigt, der den Abstand und die benachbarten Punkte beschreibt. Der Richtungsvektor wird allgemein mit dem Winkel $\alpha \in[ 0, 2\pi)$ Die Ortsvektoren der Punkte können wiederum mit @@ -280,5 +289,61 @@ unterteilt werden. Von Interesse ist lediglich das Intervall $\alpha\in\left( \frac{\pi}{2}, \frac{3\pi}{2}\right)$, da der Verfolger sich stets in die negative $y$-Richtung bewegt. In diesem Intervall ist die Ableitung negativ, woraus folgt, dass jeglicher unmittelbar benachbarte Punkt, den der Verfolger als nächstes begehen könnte, stets näher am Ziel ist als zuvor. Dies bedeutet, dass der Scheitelpunkt der Verfolgungskurve nie ein lokales Minimum bezüglich des Abstandes sein kann. +% +\subsection{Wo ist der Abstand minimal?} +Damit der Verfolger das Ziel erreicht muss die Bedingung \eqref{lambertw:minimumAbstand} erfüllt sein. +Somit ist es ausreichend zu zeigen, dass +\begin{equation} + \operatorname{min}(|z-v|)