MM08 Geometri I Ugeseddel 10 - Institut for Matematik og Datalogi ...
MM08 Geometri I Ugeseddel 10 - Institut for Matematik og Datalogi ...
MM08 Geometri I Ugeseddel 10 - Institut for Matematik og Datalogi ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
<strong>Institut</strong> <strong>for</strong> <strong>Matematik</strong> <strong>og</strong> Datal<strong>og</strong>i 2. maj 2007<br />
Syddansk Universitet, Odense MSV<br />
<strong>MM08</strong> <strong>Geometri</strong> I <strong>Ugeseddel</strong> <strong>10</strong><br />
Forelæsningerne i uge 18 <strong>og</strong> 19<br />
Vi har set flere eksempler af flade flader, såsom tangentflader, cylinderflader <strong>og</strong> kegleflader. Vi har<br />
<strong>og</strong>så kigget på pseudosfæren som et eksempel af en omdrejningsflade med konstant negativ Gausskrumning.<br />
Fælles <strong>for</strong> disse flader er at de ikke er kompakte, <strong>og</strong> dette er ikke et tilfælde. Vi skal vise<br />
at det på en kompakt flade findes mindst et punkt hvor Gausskrumningen er > 0.<br />
Derefter skal vi betragte den så kaldte Gaussafbildning <strong>for</strong> et kortomegn σ(U) ⊆ S på en flade. Dette<br />
er afbildningen fra σ(U) (eller, ved sammensætning med σ, fra U) til enhedssfæren S 2 , der til hvert<br />
punkt p ∈ σ(U) giver punktet N σ(p) ∈ S 2 . Vi skal vise, at hvis R δ ⊂ U er en disk med center i<br />
(u0, v0) ∈ U <strong>og</strong> radie δ, så gælder at<br />
areal af N<br />
lim<br />
σ(σ(R δ))<br />
δ→0 areal af σ(Rδ) = |K(σ(u0, v0))|.<br />
Vi kan altså få in<strong>for</strong>mation om Gausskrumningen fra Gaussafbildningen.<br />
Nøgleord er således<br />
• Kompakte flader;<br />
• Gaussafbildningen;<br />
• areal <strong>og</strong> Gausskrumning.<br />
Øvelser i uge 19<br />
1. 7.1, 7.11, 7.18, 7.19.<br />
2. Antag at S er en flade <strong>og</strong> p ∈ S. En tangentretning i TpS siges at være en asymptotretning (“asymptotic<br />
direction“) hvis normalkrumningen er nul i den retning. En kurve γ på S siges at være<br />
en asymptotekurve (“asymptotic“) på S hvis ˙γ(t) er en asymptoteretning <strong>for</strong> alle t, dvs. normalkrumningen<br />
<strong>for</strong> γ er nul overalt. Antag at γ er en enhedsparametriseret asymptotekurve på S,<br />
hvis krumning ikke er nul.<br />
a) Vis at K ≤ 0 langs γ.<br />
b) Antag at γ(t) = σ(u(t), v(t)) <strong>for</strong> en lokal parametrisering σ af S. Lad N σ være standardnormalen<br />
af σ, <strong>og</strong> betegn ved b bi-normalen af γ. Vis at b = ±N σ ◦ γ.<br />
c) Lad t1, t2 være en ortonormal basis <strong>for</strong> TpS af hovedvektorer, p = γ(t). Betegn ved θ<br />
vinklen mellem ˙γ(t) <strong>og</strong> t1. Vis at<br />
τn = ±(κ1 cos θt1 + κ2 sin θt2)<br />
i punktet p, der n er hovednormalen til γ i p.
3. 7.14.<br />
d) Brug Eulers <strong>for</strong>mel <strong>for</strong> at vise at det langs γ gælder<br />
|τ| = √ −K (Beltrami-Ennepers Sætning).<br />
4. En kurve på en flade siges at være en krumningskurve (“line of curvature“), hvis dens tangent<br />
i enhvert punkt er en hovedretning i dette punkt. Antag at S1 <strong>og</strong> S2 skær hinanden langs en<br />
regulær kurve γ, <strong>og</strong> lad θ(p) være vinklen mellem S1 <strong>og</strong> S2 i punktet p ∈ S1 ∩ S2. Antag at γ er<br />
en krumningskurve på S1. Vis at θ er konstant hvis <strong>og</strong> kun hvis γ <strong>og</strong>så er en krumningskurve<br />
på S2 (Joachimsthals Sætning).<br />
5. Brug maplekoden neden til at vise at alle punkter på fladen z = 3 2 xy − 1 2 x3 er hyperbolske.<br />
Maple<br />
Vi kan bruge Maple til at plotte omdrejningsflader med konstant Gausskrumning, <strong>og</strong>så i tilfælderne<br />
hvor vi ikke kan integrere per hånd. Bemærk brug af Int i stedet <strong>for</strong> int. Vi begynder med<br />
with(plots):setoptions3d(scaling=constrained):<br />
Lad os så prøve med Gausskrumning 0:<br />
sigmaZero:=(u,v) -> [(a*u+b)*cos(v),(a*u+b)*sin(v),sqrt(1-a^2)*u];<br />
a:=0.5; b:=0.5; plot3d(sigmaZero(u,v), u=-1..1, v=0..2*Pi);<br />
positiv Gausskrumning:<br />
sigmaPlus:=(u,v)->[a*cos(u)*cos(v),a*cos(u)*sin(v),<br />
Int(sqrt(1-a^2*sin(s)^2),s=0..u)];<br />
a:=0.5; b:=0.5; plot3d(sigmaPlus(u,v),u=-1..1, v=0..2*Pi);<br />
<strong>og</strong> så negativ Gausskrumning:<br />
sigmaMinus:=(u,v)->[(a*exp(u)+b*exp(-u))*cos(v),<br />
(a*exp(u)+b*exp(-u))*sin(v),Int(sqrt(1-(a*exp(s)-b*exp(-s))^2),s=0..u)];<br />
a:=0.5; b:=0.5; plot3d(sigmaMinus(u,v), u=-1..0, v=0..2*Pi);<br />
Prøv <strong>og</strong>så andre parametre, <strong>og</strong> <strong>for</strong>sikre jer om at fladerne har den specificerede Gausskrumning.<br />
Det er <strong>og</strong>så nemt at skrive rutiner <strong>for</strong> at beregne Gauss- <strong>og</strong> middelkrumning <strong>for</strong> en flade. Vi begynder<br />
med at hente de pakke vi har brug <strong>for</strong>:<br />
with(plots): setoptions3d(scaling=constrained): with(LinearAlgebra):<br />
Derefter definerer vi n<strong>og</strong>le operationer:<br />
DotProd:=proc(u::list,v::list)<br />
convert(zip(‘*‘,u,v),‘+‘);<br />
end proc;<br />
CrossProd:=proc(u::list,v::list)<br />
[u[2]*v[3]-u[3]*v[2],<br />
u[3]*v[1]-u[1]*v[3],<br />
2
u[1]*v[2]-u[2]*v[1]];<br />
end proc;<br />
NormSq:=proc(u::list)<br />
sqrt(DotProd(u,u));<br />
end proc;<br />
Vi begynder med en rutine <strong>for</strong> den første fundamentale <strong>for</strong>m:<br />
F1:=proc(flade)<br />
local E, F, G;<br />
E:=simplify(DotProd(diff(flade(u,v),u),diff(flade(u,v),u)),assume=real);<br />
F:=simplify(DotProd(diff(flade(u,v),u),diff(flade(u,v),v)),assume=real);<br />
G:=simplify(DotProd(diff(flade(u,v),v),diff(flade(u,v),v)),assume=real);<br />
;<br />
end proc;<br />
<strong>og</strong> derefter den anden fundamentale <strong>for</strong>m:<br />
F2:=proc(flade)<br />
local oN, L, M, N;<br />
oN:=simplify(expand(CrossProd(diff(flade(u,v),u),diff(flade(u,v),v))<br />
/radsimp(sqrt(Determinant(F1(flade))))));<br />
L:=simplify(DotProd(diff(flade(u,v),u,u),oN));<br />
M:=simplify(DotProd(diff(flade(u,v),u,v),oN));<br />
N:=simplify(DotProd(diff(flade(u,v),v,v),oN));<br />
<br />
end proc;<br />
Gausskrumningen får vi som determinanten af Weingartenmatrixen, <strong>og</strong> middelkrumningen som<br />
sporet af densamme delt med 2:<br />
Gauss:=proc(flade)<br />
simplify(Determinant(F1(flade)^(-1).F2(flade)));<br />
end proc;<br />
Mean:=proc(flade)<br />
simplify(Trace(F1(flade)^(-1).F2(flade))/2);<br />
end proc;<br />
Så kan vi <strong>for</strong> eksempel prøve med sigmaZero:<br />
sigmaZero:=(u,v) -> [(a*u+b)*cos(v),(a*u+b)*sin(v),sqrt(1-a^2)*u];<br />
F1(sigmaZero);<br />
F2(sigmaZero);<br />
Gauss(sigmaZero);<br />
Mean(sigmaZero);<br />
Prøv <strong>og</strong>så n<strong>og</strong>le andre flader.<br />
http://www.imada.sdu.dk/Courses/<strong>MM08</strong>/<br />
3<br />
Martin Svensson