arctan2

Die mathematische Funktion arctan2, auch atan2, ist eine Erweiterung der inversen Winkelfunktion Arkustangens und wie diese eine Umkehrfunktion der Winkelfunktion Tangens.

Sie nimmt zwei reelle Zahlen als Argumente, im Gegensatz zum normalen Arkustangens, welcher nur eine reelle Zahl zum Argument hat. Damit hat sie genügend Information, um den Funktionswert in einem Wertebereich von $360^{\circ }$ (also allen vier Quadranten) ausgeben zu können, und muss sich nicht (wie der normale Arkustangens) auf zwei Quadranten beschränken.

Der volle Wertebereich wird häufig benötigt, beispielsweise bei der Umrechnung ebener kartesischer Koordinaten in Polarkoordinaten: wenn der Funktion Vorlage:Nowrap die beiden kartesischen Koordinaten $(x,y)$ als Argumente gegeben werden, erhält man den Polarwinkel $\phi$, der sich im richtigen Quadranten befindet, d. h. der die Beziehungen

$x=r\cos \phi$ und

$y=r\sin \phi$  mit

$r:=\sqrt{x^2+y^2}$

erfüllt.

Ein mathematisch nützlicher Zusatzeffekt ist, dass Winkel, bei denen der Tangens eine Polstelle hat, nämlich die Winkel $\pm 90^{\circ }=\pm \pi /2,$ durch ganz normale reelle Koordinaten spezifiziert werden können, nämlich durch $\pm \pi /2=\mathrm{arctan2}(0,\pm 1)$ anstatt $\arctan (\pm \mathrm{\infty }).$

Das kommt von der Definitionsmenge $E:={ℝ}^2\setminus \{(0,0)\}$ der Funktion $\mathrm{arctan2},$ der „gelochten“ Ebene, welche mit einer Gruppenstruktur versehen werden kann, die isomorph ist zur multiplikativen Gruppe ${ℂ}^{\times }$ der komplexen Zahlen ohne die Null. Diese Gruppen sind direktes Produkt der Kreisgruppe $𝕊$ der Drehungen und der Gruppe der Streckungen um einen Faktor größer null, der multiplikativen Gruppe ${ℝ}^{+}.$ Erstere Gruppe lässt sich durch den Polarwinkel $\phi$ parametrisieren, zweitere durch den (positiven) Betrag $r.$

Lösung: Zwei Argumente

Zwei vom Ursprung $𝖮(0,0)$ verschiedene Punkte ${𝖯}_{\mathrm{𝟣}}(x_1,y_1)$ und ${𝖯}_{\mathrm{𝟤}}(x_2,y_2)$ spezifizieren denselben Polarwinkel, wenn sie auf demselben Strahl durch $𝖮$ liegen. Dann sind sie bezüglich der durchVorlage:Anker

	$x_1\ne 0\vee y_1\ne 0:$
	$(x_1,y_1)\sim (x_2,y_2)$	$:⟺\mathrm{\exists }\lambda >0:x_2=\lambda x_1\wedge y_2=\lambda y_1$

definierten Relation äquivalent. Dagegen ist der Tangenswert von Polarwinkeln auch dann derselbe, wenn der Strahl um $\pm 180^{\circ }$ oder Vorlage:Nowrap also genau in den Gegenstrahl, weitergedreht ist. Informationstheoretisch betrachtet lässt der Tangens die Vorzeicheninformation von $\lambda$ (rot in den Formeln) unter den Tisch fallen:

	$x_1\ne 0\wedge x_2\ne 0:$
	$\tan \left(\frac{y_1}{x_1}\right)=\tan \left(\frac{y_2}{x_2}\right)$	$⟺\mathrm{\exists }\lambda \gtrless 0:x_2=\lambda x_1\wedge y_2=\lambda y_1,$

man nehme nur Vorlage:Nowrap

Abb. 1: Graph der Arkustangensfunktion $\phi =\arctan (t)$ über $t\in (-\mathrm{\infty },+\mathrm{\infty })=ℝ$ mit gestrichelten Asymptoten. In blauen römischen Ziffern die zu $t=\frac{y}{x}$ passenden Quadranten in der Vorlage:Nowrap.

Da der Tangens mit $\pi$ periodisch ist und der Funktionsbegriff Rechtseindeutigkeit verlangt, muss für seine Umkehrung (Spiegelung an der 1. Winkelhalbierenden) sein Definitionsbereich mindestens auf die Periodenlänge eingeschränkt werden – in diesem Artikel auf das Intervall $]-\pi /2,\pi /2[$ (s. Abb. 1). Das hat zur Folge, dass die Umkehrfunktion Arkustangens kein größeres Bild als $]-\pi /2,\pi /2[$ haben kann. Dabei ist die ganze reelle Achse $ℝ$ als Definitionsbereich des Arkustangens zulässig, weil das Bild des Tangens unter $]-\pi /2,\pi /2[$ gerade $ℝ$ ist.

Um zu einem vollwertigen Polarwinkel zu kommen, gibt es in vielen Programmiersprachen und Tabellenkalkulationen eine erweiterte Funktion, die mit den beiden kartesischen Koordinaten beschickt wird und damit genügend Information hat, um den Polarwinkel modulo $2\pi \mathbb{Z}$ (bspw. im Intervall $]-\pi ,\pi ]$ wie der Abb. 3) und in allen vier Quadranten zurückgeben zu können.

Implementierungen

Die erste Implementierung war nicht später als im Jahr 1966 in der Programmiersprache Fortran. Heute ist die Funktion auch in anderen Programmiersprachen vorhanden.

Die Funktion hat häufig den Namen Vorlage:Nowrap so bei den Programmiersprachen Vorlage:Nowrap C, C++, Java, Python, Matlab, R, Vorlage:Nowrap Vorlage:Nowrap Maxima. In vielen dieser Programmiersprachen (nicht bspw. bei LibreOffice Calc) ist die Reihenfolge der Argumente umgekehrt, also die Vorlage:Nowrap das erste Argument – und das, obwohl es hier auf die Polarachse, die üblicherweise mit der Vorlage:Nowrap identifiziert wird, in ganz besonderem Maße ankommt. Denn es hat $\mathrm{arctan2}(x,y)$ gleich null zu sein für genau die Punkte $(x|y)$ auf dieser Achse. Deshalb sollte es, wenn es um die übliche Vorlage:Nowrap geht, bei der Erstnennung der Vorlage:Nowrap bleiben; die Vorlage:Nowrap ergänzt dabei nur noch die Richtung, in welche der Polarwinkel zunimmt. In Common Lisp, wo optionale Argumente existieren, erlaubt die Vorlage:Nowrap die Vorlage:Nowrap als optionales zweites Argument zu übergeben, wobei die Standardannahme $x=1$ ist.

Ein weiterer vorkommender Name ist $\mathrm{arctan2}(x,y)$, so bei den Tabellenkalkulationen Excel und OpenOffice Calc.

In Mathematica ist eine Funktion $\mathrm{ArcTan}(x,y)$ definiert, bei der das erste Argument $x$ weggelassen werden kann.

Seit der Einführung des FPU-Befehlssatzes x87 von Intel im Jahr 1980 kennen die x87-Prozessoren einen Befehl $\mathrm{FPATAN}$ mit im Prinzip derselben Reihenfolge der Operanden und demselben Ergebnis.

Zur Beachtung

In diesem Artikel wird die Reihenfolge $(x,y)$ und der Name $\mathrm{arctan2}(x,y)$ verwendet.

Formel

Die sechs Fälle der Funktionsdefinition

$\mathrm{arctan2}:$	$E$	$\to$	$]-\pi ,+\pi ]$   oder   $[-\pi ,+\pi [$			$(\text{K}\to {\text{P}}_2)$
	$(x,y)$	$\mapsto$	$\{\begin{array}{c}\\ \\ \\ \\ \\ \\ \\ \\ \end{array}$	Vorlage:Anker$\arctan \left(\frac{y}{x}\right)$	für $x>0$	(Quadranten $\mathrm{I}$ und $\mathrm{I}\mathrm{V}$)
				$\arctan \left(\frac{y}{x}\right)+\pi$	für $x<0,y>0$	(Quadrant $\mathrm{I}\mathrm{I}$)
				Vorlage:Anker$\pm \pi$	für $x<0,y=0$	(oberer⁄unterer Rand der Bildmenge)
				$\arctan \left(\frac{y}{x}\right)-\pi$	für $x<0,y<0$	(Quadrant $\mathrm{I}\mathrm{I}\mathrm{I}$)
				$+\frac{\pi }{2}$	für $x=0,y>0$
				$-\frac{\pi }{2}$	für $x=0,y<0$

Abb. 2: 5 Punkte ● (Halbgeraden, Strahlen) in der Vorlage:Nowrap und ihr $\mathrm{arctan2}$-Hauptwert

mit $E:={ℝ}^2\setminus \{(0,0)\}$ als der „gelochten“ Ebene lassen sich zur Formel

$\begin{array}{cc}& \mathrm{arctan2}(x,y):=\arctan \left(\frac{y}{x}\right)\mathrm{sgn}(x{)}^2\\ & +\frac{1-\mathrm{sgn}(x)}{2}(\mp \mathrm{sgn}(y{)}^2+\mathrm{sgn}(y)\pm 1)\pi \end{array}$

vereinigen. Die Funktion ist bis auf den Fall $x<0,y=0⟺(x,y)\sim (-1,0)$ (die Sprungstelle, s. u.) punktsymmetrisch am Ursprung, in Formeln:

$\mathrm{arctan2}(x,-y)=-\mathrm{arctan2}(x,y)$.

Dem Argument $(x,y)=(0,0)$ wird manchmal der Funktionswert $0$ zugeordnet, wie auch andere Sonderfälle, bspw. Not a Number, unterschiedlich behandelt werden.

Der Genauigkeitsverlust der Division $\frac{y}{x}$ wegen lässt sich für $0<|x|\ll |y|$ bspw. durch die Umformung

$\arctan \left(\frac{y}{x}\right)=-\arctan \left(\frac{x}{y}\right)+\mathrm{sgn}(y)\frac{\pi }{2}$

verringern (s. jedoch auch den Abschnitt Vorlage:Nowrap

Sprungstelle und kontinuierliche Drehung des Polarwinkels

Abb. 3: Graph von $\phi =\mathrm{arctan2}(x,y)=\mathrm{atan2}(y,x)$ über $x^{-1}y=y/x$ für $|x|⪆\frac{1}{4}|y|$ mit gestrichelten Asymptoten. In blauen römischen Ziffern die entsprechenden Quadranten der Vorlage:Nowrap.

Bei zunehmendem Polarwinkel $\phi$, das heißt bei einer Drehung im mathematischen Sinn (und entgegen dem Uhrzeigersinn), so bei der Wanderung vom Quadranten $\mathrm{I}\mathrm{I}\mathrm{I}$ (entspricht $\phi \in [-\pi ,-\pi /2]$) über die Quadranten $\mathrm{I}\mathrm{V}$ (entspricht $\phi \in [-\pi /2,0]$) und $\mathrm{I}$ (entspricht $\phi \in [0,+\pi /2]$) zum Quadranten $\mathrm{I}\mathrm{I}$ (entspricht $\phi \in [+\pi /2,+\pi ]$), beginnt eine Periode in der Abb. 3 unten am (Strahl durch $𝖮$ und den) Punkt

$\phi$	$x$	$y$	$\frac{y}{x}$
$-\pi$	$-1$	$-0$	$0$	$\mathrm{I}\mathrm{I}\mathrm{I}$	,

von wo es auf dem roten Graphen von links unten nach rechts oben (immer in „ONO-Richtung“) weitergeht. Wie üblich soll in der Nähe der Null $-0$ infinitesimal unterhalb und $+0$ infinitesimal oberhalb bedeuten. Die Drehung führt weiter in den Quadranten $\mathrm{I}\mathrm{I}\mathrm{I}$ über den (auf der 1. Winkelhalbierenden liegenden und in der Abb. 3 durch eine kleine rote Kreisfläche markierten) Punkt

$-\frac{3\pi }{4}$

$-1$

$-1$

$1$

$\mathrm{I}\mathrm{I}\mathrm{I}$

,

zum Punkt

	$x_{\bullet }$	$y_{\bullet }$	$\frac{y_{\bullet }}{x_{\bullet }}$	$\phi$	$x_{𝗈}$	$y_{𝗈}$	$\frac{y_{𝗈}}{x_{𝗈}}$
$\mathrm{I}\mathrm{V}$	$+0$	$-1$	$-\mathrm{\infty }$	$-\frac{\pi }{2}$	$-0$	$-1$	$+\mathrm{\infty }$	$\mathrm{I}\mathrm{I}\mathrm{I}$	,
				$⟵$

der einer der Polstellen $-\frac{\pi }{2}$ des Tangens entspricht und deshalb für den Arkustangens $\arctan$ ein unendlich ferner Punkt ist. Der Vorlage:Nowrap wechselt (auf der waagrechten und strich-punktierten Asymptote) von $+\mathrm{\infty }$ nach Vorlage:Nowrap ein Sachverhalt, den auch der rote Kringel rechts im Quadranten $\mathrm{I}\mathrm{I}\mathrm{I}$ als Vorlage:Nowrap und der rote Knubbel links im Quadranten $\mathrm{I}\mathrm{V}$ als Vorlage:Nowrap ausdrücken soll. Aus Sicht der Funktion $\mathrm{arctan2}(x,y)$ geschieht aber nichts weiter, als dass der Vorlage:Nowrap sich von $-0$ zu $+0$ ändert.

Die weitere Drehung führt durch den Quadranten $\mathrm{I}\mathrm{V}$ über den markierten Punkt

$\mathrm{I}\mathrm{V}$

$+1$

$-1$

$-1$

$-\frac{\pi }{4}$

zum Punkt

$\mathrm{I}\mathrm{V}$	$+1$	$-0$	$0$	$0$	$+1$	$+0$	$0$	$\mathrm{I}$	,
				$⟶$

der im Koordinatenursprung $𝖮$ liegt, und von dort durch den Quadranten $\mathrm{I}$ über den markierten Punkt

$\frac{\pi }{4}$

$+1$

$+1$

$1$

$\mathrm{I}$

,

zum Punkt

$\mathrm{I}\mathrm{I}$	$-0$	$+1$	$-\mathrm{\infty }$	$\frac{\pi }{2}$	$+0$	$+1$	$+\mathrm{\infty }$	$\mathrm{I}$	.
				$⟵$

Dieser Punkt entspricht der anderen Polstelle $+\frac{\pi }{2}$ des Tangens. Bei ihm findet dasselbe Zusammenfallen der Senke im Quadranten $\mathrm{I}$ mit der Quelle im Quadranten $\mathrm{I}\mathrm{I}$ statt wie oben beim Argument Vorlage:Nowrap Die weitere Drehung durch den Quadranten $\mathrm{I}\mathrm{I}$ führt über den markierten Punkt

$\mathrm{I}\mathrm{I}$

$-1$

$+1$

$-1$

$\frac{3\pi }{4}$

schließlich zur Sprungstelle

$\mathrm{I}\mathrm{I}$

$-1$

$+0$

$0$

$\pi$

.

Dieser Fall kann durch leichte Abwandlung der Bedingungen in der Formel $(\text{K}\to {\text{P}}_2)$ entweder dem Fall in der Zeile darüber oder dem darunter zugeschlagen werden, wonach das Intervall der Bildmenge an seinem oberen Ende abgeschlossen und am unteren Ende offen ist, also Vorlage:Nowrap oder eben umgekehrt Vorlage:Nowrap

Hat die Berechnung des Polarwinkels eine kontinuierliche Drehung zu begleiten, dann kann die Funktion so angepasst oder erweitert werden, dass

• die Sprungstelle an einem beliebigen Punkt (einem beliebigen Strahl) des Definitionsbereichs ${ℝ}^2\setminus \{(0|0)\}$ liegt;
• auch bei einer Drehung über die Periodenlänge $2\pi$ hinaus der Polarwinkel kontinuierlich zu- bzw. abnimmt. Hier kommt die Umlaufzahl ins Spiel.

Beispielsweise können in Anwendungen, bei denen es auf die Stetigkeit innerhalb einer Halbebene ankommt, folgende Formeln nützlich sein:

${\mathrm{arctan2}}_{-3}(x,y):=-$	$\arctan \left(\frac{x}{y}\right)-\frac{3\pi }{2}$	$=\mathrm{arctan2}(x,y)-2\pi$	für $y>0$	(Quadranten $\mathrm{I}$ und $\mathrm{I}\mathrm{I}$)
${\mathrm{arctan2}}_{-2}(x,y):=$	$\arctan \left(\frac{y}{x}\right)-\pi$	$\equiv \mathrm{arctan2}(x,y)\text{ mod }2\pi \mathbb{Z}$	für $x<0$	(Quadranten $\mathrm{I}\mathrm{I}$ und $\mathrm{I}\mathrm{I}\mathrm{I}$)
${\mathrm{arctan2}}_{-1}(x,y):=-$	$\arctan \left(\frac{x}{y}\right)-\frac{\pi }{2}$	$\begin{array}{c}\\ \\ \\ \end{array}\}=\mathrm{arctan2}(x,y)\{\begin{array}{c}\\ \\ \\ \end{array}$	für $y<0$	(Quadranten $\mathrm{I}\mathrm{I}\mathrm{I}$ und $\mathrm{I}\mathrm{V}$)
${\mathrm{arctan2}}_0(x,y):=$	$\arctan \left(\frac{y}{x}\right)$		für $x>0$	(Quadranten $\mathrm{I}\mathrm{V}$ und $\mathrm{I}$)
${\mathrm{arctan2}}_{+1}(x,y):=-$	$\arctan \left(\frac{x}{y}\right)+\frac{\pi }{2}$		für $y>0$	(Quadranten $\mathrm{I}$ und $\mathrm{I}\mathrm{I}$)
${\mathrm{arctan2}}_{+2}(x,y)=$	$\arctan \left(\frac{y}{x}\right)+\pi$	$\equiv \mathrm{arctan2}(x,y)\text{ mod }2\pi \mathbb{Z}$	für $x<0$	(Quadranten $\mathrm{I}\mathrm{I}$ und $\mathrm{I}\mathrm{I}\mathrm{I}$)
${\mathrm{arctan2}}_{+3}(x,y):=-$	$\arctan \left(\frac{x}{y}\right)+\frac{3\pi }{2}$	$=\mathrm{arctan2}(x,y)+2\pi$	für $y<0$	(Quadranten $\mathrm{I}\mathrm{I}\mathrm{I}$ und $\mathrm{I}\mathrm{V}$)

Die 7 Zeilen sind so angeordnet, dass die Werte zweier über einander stehender Funktionen im gemeinsamen Definitionsgebiet übereinstimmen. Wegen der Sprungstelle von $\mathrm{arctan2}(x,y)$ beim Strahl $\sim (-1|0)$ ist Gleichheit mit einer der stetigen Funktionen ${\mathrm{arctan2}}_{\mp 2}(x,y)$ nicht möglich.

Isomorphie zur Kreisgruppe

Die präzise mathematische Darstellung der folgenden wohlbekannten Abbildung bedarf sowohl auf der Urbild- wie auf der Bild-Seite zusätzlicher Hilfsabbildungen.

Auf der Definitionsmenge $E:={ℝ}^2\setminus \{(0,0)\}$ von $\mathrm{arctan2}$ kann man (in Analogie zur Definition der Addition in den rationalen Zahlen) die Verknüpfung

$(x,y)\widetilde{+}(x^{\prime },y^{\prime }):=(x{x}^{\prime }-y{y}^{\prime },x{y}^{\prime }+y{x}^{\prime })$

definieren. Sie bleibt wohldefiniert unter der obigen Äquivalenzrelation Vorlage:Nowrap und die Faktormenge

$F:=E/\sim$

erweist sich als kommutative Gruppe mit dem neutralen Element $\overline{(1,0)}$ und der Inversenbildung

$\widetilde{-}\overline{(x,y)}:=\overline{(x,-y)}$.

Genauso wohldefiniert ist die induzierte Abbildung

$\begin{array}{cccc}\mathrm{arctan2}:& F& \to & ]-\pi ,+\pi ]\\ & \overline{(x,y)}& \mapsto & \mathrm{arctan2}(x,y)\end{array}$

weil äquivalente Repräsentanten denselben Vorlage:Nowrap liefern. Aus der Summenformel des Arkustangens folgt

$\begin{array}{cc}\mathrm{arctan2}((x,y)\widetilde{+}(x^{\prime },y^{\prime }))& \equiv \mathrm{arctan2}(x{x}^{\prime }-y{y}^{\prime },x{y}^{\prime }+y{x}^{\prime })\\ & \equiv \mathrm{arctan2}(x,y)+\mathrm{arctan2}(x^{\prime },y^{\prime })mod2\pi \mathbb{Z}.\end{array}$

Wendet man auf die Funktion $\mathrm{arctan2}$ die Funktion

$\begin{array}{cccc}\mathrm{mod}2\pi :& ]-\pi ,+\pi ]& \to & ℝ/(2\pi \mathbb{Z})\\ & \phi & \mapsto & \phi +2\pi \mathbb{Z}\end{array}$

an, dann ergibt sich ein Homomorphismus $\mu :=\mathrm{mod}2\pi \circ \mathrm{arctan2}$ mit

$\begin{array}{cccccc}\mu :& F& \to & ℝ/(2\pi \mathbb{Z})& & \cong 𝕊\\ & \overline{(x,y)}\widetilde{+}\overline{(x^{\prime },y^{\prime })}& \mapsto & \mu (x,y)+\mu (x^{\prime },y^{\prime })\end{array}$

auf die Kreisgruppe Vorlage:Nowrap dessen Kern das neutrale Element $\overline{(1,0)}$ mit $\mu \overline{(1,0)}=2\pi \mathbb{Z}$ ist. Sind $F$ und $ℝ/(2\pi \mathbb{Z})$ mit der natürlichen Topologie ausgestattet, dann ist $\mu$ in beiden Richtungen stetig, mithin ein Homöomorphismus.

Genauigkeitskontrolle

Mit einer kleinen Vorbereitung und mit nur einem Vergleich mehr als in den Fallunterscheidungen der Formel $(\text{K}\to {\text{P}}_2)$ lässt sich das Konvergenzverhalten der Taylorreihe (des Arkustangens) kontrollieren und ggf. verbessern.

Der Winkel von $45^{\circ }$ zeichnet sich dadurch aus, dass er ein ganzzahliger Bruchteil, nämlich ein Achtel, des vollen Winkels von $360^{\circ }$ ist und gleichzeitig sein Strahl durch ganzzahlige Koordinaten geht. Quadranten lassen sich in der Koordinatenebene so ausrichten, dass ihre Begrenzungen (die definitionsgemäß stets Strahlen sind) parallel zu den Koordinatenachsen zu liegen kommen. Bei Oktanten kommen noch die Winkelhalbierenden als Begrenzungen hinzu. Die Feststellung, zu welchem der acht Oktanten ein Punkt $(x,y)\in E$ gehört, ist bei einer derartigen Ausrichtung besonders einfach.

Schreibweise

In diesem § werden in den Beziehungen zwischen Strahlen und Winkeln die gewohnten Operatoren $+-=<>$ mit der darübergeschriebenen Tilde $\widetilde{+}\widetilde{-}\widetilde{=}\widetilde{<}\widetilde{>}$ verwendet, um auszudrücken, dass ein Strahl eine Äquivalenzklasse $\text{mod }\sim$ ist. Und bei den Vergleichsoperatoren $\widetilde{<}\widetilde{>}$ wird der Strahl stets mit dem ihm $\text{mod }2\pi \mathbb{Z}$ am nächsten liegenden Winkel verglichen.

Um Verwechslungen mit Koordinaten $(x,y)$ zu vermeiden, wird in den Dezimaldarstellungen statt des Kommas der Dezimalpunkt verwendet.

Im Folgenden wird versucht, einen beliebigen Strahl $(x,y)$ resp. Winkel $\phi$ mit einfachen und umkehrbaren Drehungen in das an der Polarachse symmetrische Winkelintervall $-22,5^{\circ }\le {\phi }_2\le 22,5^{\circ }$ zu drehen. Dann ist nämlich der Absolutbetrag des Arguments $\frac{y_2}{x_2}$ in der Taylorreihe des Arkustangens Vorlage:Nowrap

In einer ersten Drehung wird der Strahl $(x,y)\in E$ um $-22,5^{\circ }$ gedreht, d. h. der Strahl

$(x_0,y_0):=(x,y)\widetilde{-}(1,\sqrt{2}-1)=(x,y)\widetilde{-}22,5^{\circ }$

gebildet. Der Oktant, in den dieser Strahl fällt, sei der Vorlage:Nowrap und die Nummerierung der Oktanten sei so gewählt, dass der erste das Winkelintervall $]0^{\circ },45^{\circ }]$ abdeckt:

Nummer $k$ des Oktanten seine untere und obere Begrenzung

$=$

$-180\stackrel{-3}{\overbrace{{\circ }^{}-13}}5\stackrel{-2}{\overbrace{{\circ }^{}-9}}0\stackrel{-1}{\overbrace{{\circ }^{}-4}}5\stackrel{0}{\overbrace{{\circ }^{}}}0\stackrel{1}{\overbrace{{\circ }^{}4}}5\stackrel{2}{\overbrace{{\circ }^{}9}}0\stackrel{3}{\overbrace{{\circ }^{}13}}5\stackrel{4}{\overbrace{{\circ }^{}18}}{0}^{\circ }$

die obere als Strahl $(x_1,y_1)=\left[\begin{array}{c}{x}_1\\ y_1\end{array}\right]$

$\in$

$\{$

$\left[\begin{array}{c}-1\\ -1\end{array}\right]$

$,\left[\begin{array}{c}0\\ -1\end{array}\right]$

$,\left[\begin{array}{c}+1\\ -1\end{array}\right]$

$,\left[\begin{array}{c}+1\\ 0\end{array}\right]$

$,\left[\begin{array}{c}+1\\ +1\end{array}\right]$

$,\left[\begin{array}{c}0\\ +1\end{array}\right]$

$,\left[\begin{array}{c}-1\\ +1\end{array}\right]$

$,\left[\begin{array}{c}-1\\ 0\end{array}\right]$

$\}$

die obere als Polarwinkel ${\phi }_1$

$=$

$-\frac{3\pi }{4}$

$-\frac{\pi }{2}$

$-\frac{\pi }{4}$

$0$

$\frac{\pi }{4}$

$\frac{\pi }{2}$

$\frac{3\pi }{4}$

$\pm \pi$

Vom derart bestimmten Oktanten wird der obere begrenzende Strahl genommen, der durch einen Punkt $(x_1,y_1)\in E$ aus der in der Tabelle gezeigten Menge von Strahlen charakterisiert werden kann. (Alle diese Punkte haben ganzzahlige Koordinaten, und beim ersten Oktanten ist bspw. Vorlage:Nowrap Dann ist $(x_1,y_1)\widetilde{-}45^{\circ }\tilde{\le }(x_0,y_0)\tilde{\le }(x_1,y_1)$ oder

$(x_1,y_1)\widetilde{-}22,5^{\circ }\tilde{\le }(x_0,y_0)\widetilde{+}22,5^{\circ }=(x,y)\tilde{\le }(x_1,y_1)\widetilde{+}22,5^{\circ }$.

Es folgt eine Drehung von Vorlage:Nowrap die zweite, jetzt um Vorlage:Nowrap so dass

$-22,5^{\circ }\tilde{\le }(x_2,y_2):=(x,y)\widetilde{-}(x_1,y_1)\tilde{\le }+22,5^{\circ }$

im gewünschten Winkelintervall ist. Um diese zweite Drehung von ${\phi }_1:=\mathrm{arctan2}(x_1,y_1)$ $=k\cdot 45^{\circ }=k\cdot \frac{\pi }{4}$ mit ganzzahligem $k$ muss das Ergebnis, wenn der Arkustangens berechnet ist, korrigiert werden.

Die erste Drehung muss nur ungefähr $22,5^{\circ }$ betragen. Wenn sie davon etwas abweicht, etwa $(5,2)\widetilde{=}21,8^{\circ }$ oder auch $(70,29)\widetilde{=}22,5035^{\circ }$ beträgt, dann wird der Strahl $(x_2,y_2)$ möglicherweise nicht so gut in das an der Polarachse symmetrische Winkelintervall eingepasst. Das Konvergenzverhalten verschlechtert sich aber wegen $(1,1)\widetilde{+}(5,-2)\widetilde{=}(7,3)$ nur geringfügig auf Vorlage:Nowrap

Nach der zweiten Drehung kann die Taylorreihe (an der Entwicklungsstelle $0$)

${\phi }_2:=\mathrm{arctan2}(x_2,y_2)=\arctan \left(\frac{y_2}{x_2}\right)$

entwickelt und die abschließende Korrektur

$\phi :={\phi }_2+{\phi }_1=\mathrm{arctan2}(x,y)$

vorgenommen werden.

Beispiele

1. Der Ausgangsstrahl sei Vorlage:Nowrap was einem Winkel von ca. $100^{\circ }$ entspricht. Durch die Vorlage:Nowrap von $(5,-2)$ kommen wir auf Vorlage:Nowrap also in den Vorlage:Nowrap Oktanten. Dessen obere Begrenzung liegt bei Vorlage:Nowrap Wir bilden die Differenz $(x_2,y_2):=(x,y)\widetilde{-}(x_1,y_1)=(-5,27)\widetilde{-}(0,1)=(27,5)$ und berechnen ${\phi }_2=\mathrm{arctan2}(27,5)=\arctan \left(\frac{5}{27}\right)\approx 10^{\circ }$ mit $\frac{5}{27}\approx 0,185$ und korrigieren mit ${\phi }_1=+\frac{\pi }{2}$ zum Endergebnis Vorlage:Nowrap
2. Der Ausgangsstrahl sei Vorlage:Nowrap was einem Winkel von ca. $-168^{\circ }$ entspricht. Durch die Vorlage:Nowrap von $(5,-2)$ kommen wir auf Vorlage:Nowrap also in den Vorlage:Nowrap Oktanten. Das obere Ende des Oktanten liegt bei Vorlage:Nowrap Da dieser Oktant die Sprungstelle $(-1,0)$ enthält, setzen wir bei diesem Vorlage:Nowrap Oktanten im Fall $y<0$ den Korrekturwinkel auf Vorlage:Nowrap Wir bilden die Differenz $(-19,-4)\widetilde{-}(-1,0)=(19,4)$ und berechnen $\mathrm{arctan2}(19,4)=\arctan \left(\frac{4}{19}\right)$ mit $\frac{4}{19}\approx 0,211$ und korrigieren mit Vorlage:Nowrap

Verbindung zum komplexen Logarithmus

Man kann die Funktion $\mathrm{arctan2}(x,y)$ für $(x,y)\in E\cong {ℂ}^{\times }$ auch über den Hauptwert $\ln$ des komplexen Logarithmus definieren als

$\mathrm{arctan2}(x,y):=\arg (x+\mathrm{i}y)=\mathrm{Im}\ln (x+\mathrm{i}y)=\frac{1}{\mathrm{i}}\ln \frac{x+\mathrm{i}y}{\sqrt{x^2+y^2}}$

mit der Argument-Funktion $\arg$. Diese Funktion wird zum Beispiel in der inversen Kinematik benutzt, um Gelenkeinstellungen korrekt zu beschreiben. Dies ist allerdings nur eine andere formale Darstellung, denn zur Berechnung muss man $\ln z$ mit $|z|=1$ bestimmen und dazu die gegebene kartesische Darstellung von $z=:x+\mathrm{i}y$ in die Polarform überführen, wobei man im Endeffekt wieder auf die oben definierte Vorlage:Nowrap mit reellen Argumenten zurückgreift.

Ableitungen

Die Funktion $\mathrm{arctan2}$ hängt von zwei Variablen ab und ist (außer im Ursprung) stetig differenzierbar, hat also zwei partielle Ableitungen. Für die Bedingung des ersten Falls (Quadranten $\mathrm{I}$ und $\mathrm{I}\mathrm{V}$) und dessen Zuordnung ergibt sich

$\frac{\partial }{\partial x}\mathrm{arctan2}(x,y)$	$=\frac{\partial }{\partial x}\arctan \left(\frac{y}{x}\right)=\frac{-y}{x^2+y^2}$,
$\frac{\partial }{\partial y}\mathrm{arctan2}(x,y)$	$=\frac{\partial }{\partial y}\arctan \left(\frac{y}{x}\right)=\frac{x}{x^2+y^2}$

Die Einschränkung auf den ersten Fall kann nachträglich fallen gelassen werden, so dass die Gleichungen für alle $(x,y)\in E$ gelten.

Damit ist

$\mathrm{grad}\mathrm{arctan2}(x,y)=\mathrm{\nabla }\mathrm{arctan2}(x,y)=(\frac{-y}{x^2+y^2},\frac{x}{x^2+y^2})$

der Gradient der Funktion Vorlage:Nowrap und seine Richtung ist an jedem Punkt $(x,y)$ senkrecht zum Radiusvektor in mathematisch positiver Drehrichtung. Das passt zu der Tatsache, dass der Funktionswert von Vorlage:Nowrap der Polarwinkel, in dieser Richtung zunimmt.

Des Weiteren folgt für das totale Differential

$\mathrm{d}\mathrm{arctan2}=\frac{\partial \mathrm{arctan2}}{\partial x}\mathrm{d}x+\frac{\partial \mathrm{arctan2}}{\partial y}\mathrm{d}y=\frac{-y}{x^2+y^2}\mathrm{d}x+\frac{x}{x^2+y^2}\mathrm{d}y.$

Eine Integration dieses Differentials entlang eines Weges ergibt die Änderung des (Polar)winkels über den Weg. Ist der Weg geschlossen, so erhält man die Umlaufzahl (in Bezug auf den Ursprung Vorlage:Nowrap

Siehe auch

• Formelsammlung Trigonometrie
• Trigonometrische Funktionen

Weblinks

• Math. Java 1.6 SE JavaDoc
• C++ Programmer’s Reference
• MATLAB Function Reference
• Mathematica Function Reference
• atan2 auf Everything2