Riemannian Geometry in Robotics
Riemannian Geometry in Robotics
Harald Löwe
General Assumptions
∂ ∂x −1
and .
∂xi p ∂ui x(p)
The Problem
Consider a surface M ⊆ R3 .
warped torus in R3
Question
What does a flat mathematician living on M knows about M?
Definition
I is called first fundamental form on M.
The intrinsic geometry of M is the geometry of I .
The angle between these curves at the intersection point γ(0) = δ(0) is
given by
I (γ̇(0), δ̇(0))
cos (∠(γ, δ)) = p q .
I (γ̇(0), γ̇(0)) · I (δ̇(0), δ̇(0))
Löwe (TU Braunschweig) Part III: Riemannian Geometry July 2012 5 / 61
Part III: Riemannian Geometry The First Fundamental Form
Thus, the scalar product is given by the matrix Gxp = (gij (p))i,j :
I (v , w ) = v T · Gxp · w .
Löwe (TU Braunschweig) Part III: Riemannian Geometry July 2012 6 / 61
Part III: Riemannian Geometry The First Fundamental Form
I (v , w ) = vyT Gyp wy = (Jvx )T Gyp (Jwx ) = vxT (J T Gyp J)wx = vxT Gxp wx ,
whence !T !
∂yi ∂yi
Gxp = · Gyp · .
∂xj p ∂xj p
Therefore, we have
p1
!T p1
!
−p2 −p2
|p| 1 |p|
Gyp = J T Gxp J = p2 · · p2
|p| p1 1 |p| p1
p1
!
p1 p2
−p2
|p| |p| |p| 1
= · p2 =
−p2 p1 |p| p1 |p|2
Riemannian Metrics
Gp (v , w ) := v T Gp w .
Remarks:
Example 3: Ellipsoids
x2 y2 z2
M: + 2 + 2 = 1,
a2 b c
with its first fundamental form IM , i.e. a diffeomorphism f : S2 → M with
Example 3 (cont.)
a2
(Jp F )T · (Jp F ) = Gp = b2 ,
c2
F
←→
G I
For p ∈ S2 put
1 + x2
Gp = 1 + y2
1+ z2
Then (Jp F )T (Jp F ) = Gp for every p, i.e. (S2 , G) looks like (F (S2 ), I ).
Example 4 (cont.)
Isometries
Example 5: Isometries of R2
Claim
If F : (R2 , I ) → (R2 , I ) is an isometry, and if det Jp F > 0 holds for at least
one p, then
Example 5 (cont.)
With F = (F 1 , F 2 ) we have
1
Fx Fy1
cos(α(x, y )) − sin(α(x, y ))
=
Fx2 Fy2 sin(α(x, y )) cos(α(x, y ))
G in local coordinates
G(s,r ) = m · 1.
Ekin (v ) = G(v , v ).
Example 6 (cont.)
Kinetic energy:
∂ϕ ∂ϕ ∂ϕ
v = q̇1 + q̇2 + q̇3 ⇒ Gq (v , v ) = m (l 2 q̇12 + q̇22 + q̇32 )
∂q1 q ∂q2 q ∂q3 q
Lenghts of Curves
The infimum of all lengths of all curves joining p and q is called distance
d (p, q) between p and q.
(M, d ) is a metric space.
In local coordinates:
k k
!
X ∂Y1 X ∂Yn
DP k ∂
X i ∂x
Y (p) = X i (p) , ..., X i (p)
i =1 i ∂xi p ∂xi p
i=1 i=1
Covariant Derivatives
Definition
∇X Y is called the covariant derivative of Y in direction of X .
Computation Of ∇
Let x be a chart of M:
k
X ∂
X = Xi : operator on C ∞ (M)
∂xi
i=1
k k
X ∂ j
X ∂x −1
Y = Y : function p 7→ Y j (p)
∂xj ∂uj x(p)
j=1 j=1
k k k X k
i ∂ j ∂ i ∂ j ∂
X X X
DX Y = X Y = X Y
∂xi ∂xj ∂xi ∂xj
i=1 j=1 i=1 j=1 | {z }
interesting part
Computation Of ∇ (cont.)
!
∂x −1
∂ j ∂ ∂ j
Y (p0 ) = p 7→ Y (p)
∂xi ∂xj ∂xi p0 ∂uj x(p)
∂x −1
∂ j −1
= u 7→ Y (x (u))
∂ui x(p0 ) ∂uj u
∂(Y j ◦ x −1 ) ∂x −1 ∂ 2 x −1
= + Y j (p0 )
∂ui x(p0 ) ∂uj x(p0 ) ∂ui ∂uj x(p0 )
∂Y j ∂ ∂
= +Y j (p0 ) D ∂ (p0 )
∂xi p0 ∂xj
p0
∂xi ∂xj
| {z } | {z }
tangent to M interesting part
Christoffel Symbols
We have that
k k
∂Y j ∂
X
i
X
j ∂
DX Y = X · +Y ·D ∂
∂xi ∂xj ∂xi ∂xj
i=1 j=1
k k
∂Y j ∂
X
i
X
j ∂
∇X Y = X · +Y ·∇ ∂
∂xi ∂xj ∂xi ∂xj
i=1 j=1
Definition
Let
k
∂ X ∂
∇ ∂ = Γrij .
∂xi ∂xj ∂xr
r =1
∂ ∂
Recall that ∇ ∂
∂xj is the tangent part of D ∂
∂xj , i.e.
∂xi ∂xi
∂ ∂ ∂
I ,D∂ −∇ ∂ = 0 for all s.
∂xs ∂xi ∂xj ∂xi ∂xj
∂ Pk r ∂
With ∇ ∂
∂xj = r =1 Γij ∂xr we get
∂xi
k
∂ ∂ X ∂ ∂
I ,D∂ = Γrij ·I , .
∂xs ∂xi ∂xj ∂xs ∂xr
r =1 | {z }
=grs
We compute (ϕ := x −1 ):
!
∂ϕ T ∂ϕ
∂gβγ ∂
= ·
∂uα ∂uα ∂uβ ∂uγ
T
∂ϕ T
2
∂2ϕ
∂ ϕ ∂ϕ
= · + ·
∂uα ∂uβ ∂uγ ∂uβ ∂uα ∂uγ
Putting
(g rs )rs := (grs )−1
rs ,
Compute
Γ111 = Γ211 = 0
∂2ϕ
= 0 · ϕ1 + 0 · ϕ2 − ν
∂u12
Γ112 = 0, Γ212 = − tan(u1 )
∂2ϕ
= 0 · ϕ1 − tan(u1 ) · ϕ2 + 0 · ν
∂u1 ∂u2
Γ122 = cos(u1 ) sin(u1 ), Γ222 = 0
∂2ϕ
= cos(u1 ) sin(u1 ) · ϕ1 + 0 · ϕ2 − cos2 (u1 ) · ν
∂u22
k
X
r
Ẏ + Γrij ρ̇i Y j = 0 for r = 1, . . . , k.
i,j=1
d
1 If γ is a geodesic, then dt Gγ(t) (γ̇(t), γ̇(t)) ≡ 0.
2 For p ∈ M and v ∈ Tp M, there exists a geodesic γ with γ(0) = p and
γ̇(0) = v . Locally at t = 0, γ and every other geodesic γ ∗ having
these initial values coincide.
Take Rn with its usual scalar product and the obvious chart.
Observe that Γrij ≡ 0 for all i, j, r ,
Thus, Y is parallel along γ iff Ẏ r ≡ 0, i.e. iff Y is a “constant” vector field.
This fact should explain the name “parallel”.
0 = ü + cos(u) sin(u)v̇ 2
0 = v̈ − 2 tan(u)u̇ v̇ .
Difficult? Yes!
0 = ∇γ̇ γ̇ = γ̈ >
= γ̈ − I (γ̈, γ) · γ.
γ̈ = I (γ̈, γ) · γ.
Gq1 ,q2 ,q3 ((q̇1 , q̇2 , q̇3 ), (q̇1 , q̇2 , q̇3 )) = m (l 2 q̇12 + q̇22 + q̇32 )
Observe that all Γrij vanish, whence the geodesics are lines.
Zb
1
E (γ) := G(γ̇(t), γ̇(t)) dt.
2
a
d
E (t 7→ Φ(t, s)) = 0.
ds s=0
Orientable Manifolds
Theorem
A hypersurface is orientable iff it admits a smooth normal field,
i.e. there is a smooth map N : M → Sn ⊆ Rn+1 such that N(p) is
perpendicular to M at p.
N is called the Gauss map of M.
Remarks
∂x −1
= (−(3 + cos(v )) sin(u), (3 + cos(v )) cos(u), 0)
∂u
∂x −1
= (− sin(v ) cos(u), − sin(v ) sin(u), cos(v ))
∂v
∂(N ◦ x −1 ) cos(v ) ∂x −1
= (− cos(v ) sin(u), cos(v ) cos(u), 0) =
∂u 3 + cos(v ) ∂u
∂(N ◦ x −1 ) ∂x −1
= (− sin(v ) cos(u), − sin(v ) sin(u), cos(v )) =
∂v ∂v
Example 9 (cont.)
∂ ∂
Thus, with respect to the basis ∂xu p , ∂xv p of Tp T, the Weingarten map
Lp is given by
! !
cos(v ) 3
− 3+cos(v ) p12 +p22
−1
Lp = =
−1 −1
IIp (v , w ) := hLp (v ) | w i
Normal Curvature
If v ∈ Tp M with I (v , v ) = 1,
then we call II (v , v ) the normal curvature of M in direction of v .
1 k1 , . . . , kn are called principal curvatures of M at p.
2 det Lp = k1 · . . . · kn is the Gaussian curvature of M in p.
1
3
2 tr Lp = 12 (k1 + . . . + kn ) is the mean curvature of M in p.
∂(x −1 )
∂ −1
0 = N(x (u))
∂ui x(p) ∂uj u
* + * +
∂(N ◦ x −1 ) ∂(x −1 ) ∂ 2 x −1
= + N(x −1 (x(p)))
∂ui x(p) ∂uj x(p) ∂ui ∂uj x(p)
* ! + * +
∂ ∂ ∂ 2 x −1
= −Lp + N(p)
∂xi p ∂xj p ∂ui ∂uj x(p)
Therefore
! * +
∂ ∂ ∂ 2 x −1
hij := II , = N(p)
∂xi p ∂xj p ∂ui ∂uj x(p)
(3 + cos(u2 ))2
∂ϕ ∂ϕ
(gij )ij = =
∂ui ∂uj ij
1
∂2ϕ
− cos(u2 ) (3 + cos(u2 ))
(hij )ij = N = .
∂ui ∂uj ij
−1
n
∂hij ∂hik X r
Γij hrk − Γrik hrj
0 = − +
∂uk ∂uj
r =1
(Condition of Codazzi–Mainardi)
Theorem Of O. Bonnet
A Strange Corollary
Two Theorems
Question
Are there abstractly defined manifolds / Riemannian manifolds which are
not submanifolds of some Rn ?
2
ϕ(R v ) := p p T − 1 ∈ sym(n + 1)
|p|2
1
G(A, B) := tr(AT B) for A, B ∈ sym(n + 1).
8
Neither the geometric structure nor the differential geometry of RP n is
sufficiently clear in this picture.