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.

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

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

Saved successfully!

Ooh no, something went wrong!