01.12.2012 Views

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 ...

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!