Fachbereich Mathematik und Informatik WS 2007/2008 der Philipps-Universität Marburg
Stephan Dahlke Manuel Werner
2. Übungsblatt zur Vorlesung Computer Aided Geometric Design
Abgabe: Donnerstag, 08.11.2007, 14.00 Uhr, Raum D 6411, Lahnberge Aufgabe 2: Baryzentrische Koordinaten
Es seien a
0, a
1, a
2affin unabhängige Punkte in E
n. Es bezeichne vol
2([a
0, a
1, a
2]) das Volumen des durch die drei Punkte aufgespannten Zwei-Simplex.
00000000000000 00000000000000 00000000000000 00000000000000 00000000000000 00000000000000 00000000000000 00000000000000 00000000000000 00000000000000 00000000000000 00000000000000
11111111111111 11111111111111 11111111111111 11111111111111 11111111111111 11111111111111 11111111111111 11111111111111 11111111111111 11111111111111 11111111111111 11111111111111
a
2a
0a
1x
Man beweise, dass für jeden Punkt x, der in der konvexen Hülle von a
0, a
1, a
2liegt, eine eindeutige Darstellung
x =
2
X
i=0
λ
ia
i,
2
X
i=0
λ
i= 1
existiert, so dass
λ
0= vol
2([a
1, a
2, x])
vol
2([a
0, a
1, a
2]) , λ
1= vol
2([a
0, a
2, x])
vol
2([a
0, a
1, a
2]) und λ
2= vol
2([a
0, a
1, x]) vol
2([a
0, a
1, a
2]) .
(2 Punkte) Aufgabe 3: Parametrisierung nach der Bogenlänge
Sei x : [a, b] → E
3eine reguläre parametrische Kurve. Weiter sei s : [a, b] → [0, R
ba
k x(θ)k ˙
2dθ], s := s(u) = R
ua
k x(θ)k ˙
2dθ die Parametrisierung nach der Bogenlänge.
Zeige
a) s ist regulär, b) L := R
ba
k x(θ)k ˙
2dθ ist unabhängig von der Parametrisierung,
c) kx
0(s)k
2= 1, wobei x
0(s) =
dsdx(s) die Ableitung nach dem Bogenlängenparameter bezeichnet.
(3 Punkte)
Aufgabe 4: Implementierung der Bézier-Bernstein-Polynome
a) Implementieren Sie eine Funktion bernstein_basis_function(n, x), welche die Bernstein-Basis-Polynome vom Grad n ∈ N
0,
B
jn(t) = n
j
t
j(1 − t)
n−j, j = 0, . . . , n
an den im Zeilenvektor x übergebenen Stellen auswertet. (3 Punkte) b) Schreiben Sie eine Funktion plot_bernstein_basis(n), welche alle
Bernstein-Basis-Polynome B
jnvom Grad n, j = 0, . . . , n, zusammen in ein
Koordinatensystem zeichnet. (2 Punkte)
c) Entwickeln Sie eine Funktion bernstein_polynom(n, c, x), die das Bézier-Bernstein-Polynom
B
n:=
n
X
j=0