Sources/Lexikon/Kugelkoordinaten

Aus numpedia
Zur Navigation springen Zur Suche springen

Definition

In Kugelkoordinaten oder räumlichen Polarkoordinaten wird ein Punkt P im dreidimensionalen Raum durch seinen Abstand vom Ursprung und zwei Winkel angegeben.

Kugelkoordinaten r, φ1, φ2 eines Punktes P und kartesisches Koordinatensystem mit den Achsen x1, x2,x3.

So ist der Vektor vom Ursprung zum Punkt P

r→=r⋅e→r

mit

e→r=(sin⁡(φ1)cos⁡(φ2)sin⁡(φ1)sin⁡(φ2)cos⁡(φ1)).

Die Kugelkoordinaten kann man - ähnlich wie bei den Euler-Winkeln -zur Definition eines neuen, lokalen Koordinatensystems nutzen. Neben dem Flächen-Normalenvektor e→r spannen dabei die Tangentialvektoren e→φ,1,e→φ,2 zu φ1 und φ2 eine neue Basis e_→K auf.

Einheitsvektoren der Orthogonalbasis e_→K=[e→φ,1,e→φ,2,e→r], die in Punkt P der Kugel mit e→r die Flächennormale definieren und mit e→φ,1,e→φ,2 die Tangentialebene aufspannen.

Die Koordinatentransformation erfolgt über

(e→φ,1e→φ,2e→r)=D__12(φ1(t),φ2(t))(e→x,1e→x,2e→x,3)

mit der Transformationsmatrix

D__12(φ1(t),φ2(t))=(cos⁡(φ1(t))cos⁡(φ2(t))cos⁡(φ1(t))sin⁡(φ2(t))−sin⁡(φ1(t))−sin⁡(φ2(t))cos⁡(φ2(t))0sin⁡(φ1(t))cos⁡(φ2(t))sin⁡(φ1(t))sin⁡(φ2(t))cos⁡(φ1(t)))

Title

Oft kann man räumliche Polarkoordinaten besser interpretieren als Euler-Winkel. Wir schreiben deshalb hier eine Transformationsvorschrift an, mit der Euler-Winkel in Polarkoordinaten übersetzten kann.

Dabei geben wir von einem Punkt P in Euler-Koordinaten mit

cos⁡(θ3)=sin⁡(φ2)cos⁡(φ3)sin⁡(φ3)+sin⁡(φ1)cos⁡(φ2)cos⁡(φ3)2−sin⁡(φ1)cos⁡(φ2)sin⁡(θ3)2sin⁡(φ2)sin⁡(θ3)sin⁡(θ2)=−sin⁡(φ1)sin⁡(θ3)sin⁡(φ2)sin⁡(φ3)+sin⁡(φ1)cos⁡(φ2)cos⁡(φ3)cos⁡(θ2)=cos⁡(φ1)sin⁡(φ2)sin⁡(θ3)sin⁡(φ2)sin⁡(φ3)+sin⁡(φ1)cos⁡(φ2)cos⁡(φ3)sin⁡(θ1)=sin⁡(φ2)sin⁡(φ3)+sin⁡(φ1)cos⁡(φ2)cos⁡(φ3)sin⁡(θ3)cos⁡(θ1)=cos⁡(φ1)cos⁡(φ2)sin⁡(θ3)2=sin⁡(φ2)2sin⁡(φ3)2+2sin⁡(φ1)cos⁡(φ2)sin⁡(φ2)cos⁡(φ3)sin⁡(φ3)+(1−cos⁡(φ1)2)cos⁡(φ2)2cos⁡(φ3)21−cos⁡(φ1)2cos⁡(φ2)2

/* Euler-Winkel */
D[1](φ) := matrix([1,0,0],[0,cos(φ),sin(φ)],[0,-sin(φ),cos(φ)]);
D[2](φ) := matrix([cos(φ),0,-sin(φ)],[0,1,0],[sin(φ),0,cos(φ)]);
D[3](φ) := matrix([cos(φ),sin(φ),0],[-sin(φ),cos(φ),0],[0,0,1]);

/* Kugelkoordinanten */
D[12](θ,φ):= matrix([ cos(θ)*cos(φ), cos(θ)*sin(φ), -sin(θ)],
                    [       -sin(φ),        cos(φ),     0  ],
                    [ sin(θ)*cos(φ), sin(θ)*sin(φ),  cos(θ)]);



/* Zusammenhang Euler-Winkel - Kugelkoordinaten*/
/* Euler -> Kugel */

equ: matrix([a[x],a[y],a[z]]).D[3](θ[3]).D[12](θ[1],θ[2]) - matrix([a[x],a[y],a[z]]).D[3](φ[3]).D[1](φ[1]).D[2](φ[2]);
equ: flatten(makelist(coeff(args(equ),[a[x],a[y],a[z]][i]),i,1,3));
equ: makelist(equ[i]=0,i,1,9);
repl: solve([equ[5],equ[6],equ[7],equ[8],equ[9]],[cos(θ[3]),sin(θ[2]),cos(θ[2]),sin(θ[1]),cos(θ[1])])[1];
equ: transpose(subst(repl,equ));
repl: append(trigsimp(repl), subst([sin(φ[1])^2=1-cos(φ[1])^2],trigsimp(solve(equ[1],sin(θ[3])^2))));
equ: trigsimp(transpose(subst(repl,equ)));




Links

  1. Kugelkoordinaten auf Wikipedia
  2. Sources/Lexikon/Eulersche Winkel
  3. Quaternionen für Drehungen
  4. Geographische Koordinaten

Literature

  • ...