01.12.2012 Views

Geodætiske kurver 1 Introduktion - DTU Matematik ∙ Institut for ...

Geodætiske kurver 1 Introduktion - DTU Matematik ∙ Institut for ...

Geodætiske kurver 1 Introduktion - DTU Matematik ∙ Institut for ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

1 <strong>Introduktion</strong><br />

1.0.1 Hvad er en geodætisk kurve?<br />

<strong>Geodætiske</strong> <strong>kurver</strong><br />

Køreplan<br />

01005 <strong>Matematik</strong> 1 - FORÅR 2006<br />

Hvad er den korteste vej fra A til B? Det er et spørgsmål der igennem tiderne har optaget<br />

dovne mennesker. I R3 er den korteste vej mellem to punkter en ret linie, og den afstand der<br />

tilbagelægges ved at følge den rette linie fra A til B er den euklidiske afstand, ℓ:<br />

�<br />

ℓ = (xA − xB) 2 +(yA − yB) 2 +(zA − zB) 2<br />

Hvis man nu ikke kan bevæge sig frit i rummet, men er tvunget til at følge en bestemt flade,<br />

f.eks <strong>for</strong>di man bor på en kugle, er det imidlertid ikke sikkert, at man kan <strong>for</strong>binde to punkter<br />

med en ret linie, så denne er indeholdt i fladen. Så er man nødt til at følge en krum kurve <strong>for</strong><br />

at komme fra A til B. Spørgsmålet er så: Hvordan vælger man den kurve, man skal følge <strong>for</strong> at<br />

komme fra A til B og tilbagelægge kortest mulig afstand? Den kurve, der løser dette problem,<br />

er en geodætisk kurve på den pågældende flade. Mere generelt siges geodætiske <strong>kurver</strong> at være<br />

lokalt længdeminimerende, men det er ikke noget vi kommer nærmere ind på her.<br />

Dette projekt handler om bestemmelse af geodætiske <strong>kurver</strong> på <strong>for</strong>skellige flader.<br />

Projektet består af tre dele:<br />

Del A Her skal I bruge Maples ’NLPSolve’-kommando til at bestemme en “polygonisk geodæt”<br />

på en flade.<br />

Del B Her betragtes en indskreven polygon til en kurve, og I skal vise teoretisk og i Maple, at<br />

polygonlængden går mod kurvelængden.<br />

Del C Her skal I udlede de differentialligninger, der beskriver geodætiske <strong>kurver</strong> generelt, samt<br />

bestemme geodæter i R 2 , på en cylinderflade, og på flader af eget valg.<br />

Før vi går videre, er der imidlertid nogle definitioner, der skal på plads.<br />

Grundlæggende definitioner:<br />

En rumkurve kan beskrives ved en parameter:<br />

En glat kurve i R 3 er en uendeligt mange gange differentiabel funktion, γ : [a;b] → R 3 , der fører<br />

en reel variabel fra et interval [a;b] ⊆ R over i R 3 .<br />

Det betyder, at en kurve, γ, kan skrives som en vektorfunktion af en variabel:<br />

γ(t) = � x(t),y(t),z(t) � , t,x,y,z ∈ R .<br />

<strong>Matematik</strong> 1 – 05/06 side 1


p 1<br />

pi<br />

p i+1<br />

p 2<br />

Figur 1: Enhver glat kurve, α kan tilnærmes med en polygonkurve. Dette gøres ved vælge et antal punkter,<br />

(p1, p2,...pn), på kurven og <strong>for</strong>binde dem med rette linier. Punkterne kaldes hjørnepunkter. Jo flere<br />

punkter, jo bedre tilnærmelse. En kurve kan tilnærmes vilkårligt godt med med en polygon.<br />

En flade i rummet kan beskrives ved to parametre:<br />

En parameteriseret flade i R 3 er en funktion, F : Ω → R 3 , der fører to reelle variable fra et<br />

område i planen Ω ⊆ R 2 over i R 3 .<br />

Det betyder, at en flade F kan skrives som en vektorfunktion af to variable:<br />

F(u,v) = � x(u,v),y(u,v),z(u,v) � u,v,x,y,z ∈ R .<br />

I dette projekt betragtes geodætiske <strong>kurver</strong> på flader som er graf <strong>for</strong> en funktion af to variable.<br />

Det vil sige flader, F , der kan beskrives ved en parameterfremstilling på <strong>for</strong>men:<br />

her er det x og y, der er parametrene<br />

F(x,y) = � x,y, f(x,y) � , (x,y) ∈ R 2 , f(x,y) ∈ R<br />

En kurve ligger i en flade, når alle punkter på kurven er indeholdt i fladen.<br />

Nu er vi klar til <strong>for</strong>mulere problemet med at bestemme en geodætisk kurve mere <strong>for</strong>melt.<br />

Find den kurve, som har den mindste kurvelængde blandt de <strong>kurver</strong>, der <strong>for</strong>binder A og B<br />

og som ligger i fladen.<br />

I første omgang vil vi prøve at tilnærme en løsning, idet vi i stedet <strong>for</strong> en kurve i fladen betragter<br />

en polygonkurve, hvis hjørnepunkter ligger i fladen.<br />

En polygonkurve, Γ, består af et antal hjørnepunkter p1, p2,..., pi,..., pn i R 3 , der er <strong>for</strong>bundet<br />

af rette liniestykker. En indskreven polygon <strong>for</strong> en glat kurve, γ, er en polygonkurve hvis<br />

hjørnepunkter ligger på γ. Det betyder, at ved at vælge nogle værdier, ti, i intervallet [a;b], så<br />

a = t1 < t2 < ... < ti < ... < tn = b ,<br />

kan man angive hjørnepunkterne pi = (xi,yi,zi) <strong>for</strong> en indskreven polygon sådan:<br />

pi = γ(ti) = � x(ti),y(ti),z(ti) �<br />

Figur 1 viser, hvordan det kan se ud.<br />

<strong>Matematik</strong> 1 – 05/06 side 2<br />

pn


x<br />

z<br />

Γ<br />

pi<br />

pi+1<br />

Figur 2: Fladen F og polygonen Γ.<br />

Når geodætproblemet skal løses <strong>for</strong> en polygonkurve, har vi altså to faste endepunkter A = p1<br />

og B = pn samt et antal frie hjørnepunkter pi, i = 1...n, der kan flyttes rundt på fladen. De frie<br />

hjørnepunkter skal så placeres sådan, at længden af polygonkurven bliver mindst mulig.<br />

Længden, ℓn, af en polygonkurve er summen af sidelængderne:<br />

2 DEL A<br />

n<br />

ℓn = ∑<br />

i=2<br />

�<br />

(xi − xi−1) 2 +(yi − yi−1) 2 +(zi − zi−1) 2<br />

2.1 Geodæter repræsenteret ved polygoner<br />

I første del af projektet betragter vi polygoner, Γ, der som tidligere nævnt bliver beskrevet ved<br />

hjørnepunkterne p1,..., pi,..., pn ∈ R 3 , på en flade F(x,y) = (x,y, f(x,y)), se figur 2.<br />

Ved at holde Γ’s endepunkter p1 og pn fast er problemet at minimere længden af Γ, og derved<br />

opnå en tilnærmelse til en geodæt på fladen F . Dette kan <strong>for</strong>muleres således:<br />

hvor δi er afstanden mellem pi og pi+1:<br />

min<br />

p j, j=2,...,n−1<br />

y<br />

F<br />

n−1<br />

∑ δi , (2.1)<br />

i=1<br />

δi =<br />

�<br />

(xi+1 − xi) 2 +(yi+1 − yi) 2 +( f(xi+1,yi+1) − f(xi,yi)) 2 , i = 1,...,n − 1 . (2.2)<br />

Opgave 1. Hvis problemet (2.1) blev minimeret uden yderligere betingelser, så ville alle punkterne<br />

blive placeret i p1 og pn. Forklar hvor<strong>for</strong> dette er tilfældet.<br />

Vi er der<strong>for</strong> nød til at tilføje en begrænsning til minimerings problemet, <strong>for</strong> eksempel at:<br />

δi = δi+1, i = 1,...,n − 2 (2.3)<br />

<strong>Matematik</strong> 1 – 05/06 side 3


Dette betyder, at at vores minimeringsproblem bliver ligning (2.1) under hensyn til begrænsningerne<br />

givet ved (2.3).<br />

Opgave 2. Hvad betyder begrænsning (2.3) <strong>for</strong> polygonen Γ?<br />

2.1.1 Løsning i Maple<br />

Det<br />

�Ö�×Ø�ÖØ�Û�Ø�ÔÐÓØ×�Û�Ø�ÇÔØ�Ñ�Þ�Ø�ÓÒ�<br />

er typisk ikke muligt at løse (2.1) og (2.3) i hånden, og det er der<strong>for</strong> nødvendigt at bruge en<br />

numerisk algoritme i Maple.<br />

MaplesÇÔØ�Ñ�Þ�Ø�ÓÒpakke indeholder mange <strong>for</strong>skellige algoritmer til løsning af minimeringsproblemer.<br />

Vores problem er ikke-lineært, der<strong>for</strong> bruger vi procedurenÆÄÈËÓÐÚ�(NLP-<br />

Solve:<br />

����ÜÝ<br />

Non-Linear Program Solver). Prøv at udføre følgende kommandoer i Maple:<br />

���ÑÑ�Ü���×�Õ ���ÑÑ�Ý���×�Õ �×�ÒÜ� È�Ò È�Ò Ó×Ý�� Ò� Ò� �� ��<br />

Funktionen f , der giver fladen, bliver defineret som en funktion af to variable, og kurven Γ gives<br />

som to vektorer indeholdende x og y koordinaterne.<br />

�Ò��ÒÓÔ×ÓÒÚ�ÖØ��ÑÑ�ÜÐ�×Ø� ���ÐØ����ÜÝ�<br />

Antallet<br />

���ÐØ���ÜÝ<br />

af punkter n defineres som antallet af elementer i gammax vektoren. Ligning (2.2), (2.1)<br />

og<br />

����ÜÝ� �Ü��℄Ý��℄<br />

(2.3) defineres som<br />

���ÐØ��ÜÝ���ÐØ��ÜÝ� ������ÐØ��ÜÝ��� �×ÕÖØÜ��℄Ü��℄� �Ü��℄Ý��℄�� Ò �Ý��℄Ý��℄�<br />

henholdsvis deltai, delta og bc:<br />

�<br />

��Ò�ÓÐ�Ú�Ð�ÇÔØ�Ñ�Þ�Ø�ÓÒ℄��� �ÇÈ��ÆÄÈËÓÐÚ���ÐØ�ÜÝß×�Õ�ÜÝ����<br />

MaplesÆÄÈËÓÐÚ�har syntaksenÆÄÈËÓÐÚ�(Funktion der skal minimeres, {begrænsninger}, initialpoint<br />

= {begyndelsebetingelser},<br />

×�ÕÝ��℄���ÑÑ�Ý��℄�� ��ÑÑ�Ü�℄Ü�Ò℄���ÑÑ�Ü�Ò℄Ý�℄���ÑÑ�Ý�℄ Ý�Ò℄���ÑÑ�Ý�Ò℄��Ò�Ø��ÐÔÓ�ÒØ�ß×�ÕÜ��℄���ÑÑ�Ü��℄�� Ò Ü�℄�<br />

ekstra parametre),<br />

Ò�Ñ�Ø�Ó��×ÕÔ�Ø�Ö�Ø�ÓÒÐ�Ñ�Ø�Ò læs mere omÆÄÈËÓÐÚ�iMaples��ÐÔ.<br />

�<br />

��ÓÔ���×�Õ×Ù�×ÇÈ�℄Ü��℄�� Ò��<br />

Da outputtet fraÆÄÈËÓÐÚ�er en liste bestående af de x og y koordinater, der minimerer længden<br />

�Ô��ÔÓ�ÒØÔÐÓØ��×�Õ��ÓÔ��℄�ÓÔ��℄��ÓÔ��℄�ÓÔ��℄℄�� ��ÓÔ���×�Õ×Ù�×ÇÈ�℄Ý��℄�� Ò��<br />

af Γ, skal kordinaterne i outputtet konverteres til to vektorer:<br />

�Ô��ÔÓ�ÒØÔÐÓØ��×�Õ���ÑÑ�Ü��℄��ÑÑ�Ý��℄���ÑÑ�Ü��℄��ÑÑ�Ý��℄℄ ÓÒÒ�Ø�ØÖÙ��Ü�×��ÓÜ��ÓÐÓÖ��Ð��Ø���Ò�××�� Ò℄<br />

Resultatet kan nu plottes:<br />

<strong>Matematik</strong> 1 – 05/06 side 4


�Ô��ÔÐÓØ������������×ØÝÐ��Ô�Ø�ÒÓ�Ö�� �� �ÐÓ××�Ò�××��Ö���� Ò℄ÓÒÒ�Ø�ØÖÙ��Ü�×��ÓÜ��ÓÐÓÖ��ÐÙ�Ø���Ò�××�� ℄�<br />

���ÑÑ�Ü���×�ÕÒ�Ò��Ò�� ����ÜÝ �×�ÒÜ Ó×Ý�� ×�ÕÒ� �Ò��Ò� ��<br />

���×ÔÐ�ÝÔÔÔ�<br />

���ÑÑ�Ý���×�ÕÒ�Ò�����<br />

Opgave 3. Prøv at ændre f og polygonen Γ i Maplen koden til:<br />

Eksperimentér også med egen valg af flader.<br />

3 DEL B<br />

3.1 Fra længden af en polygon til kurvelængden af en kurve<br />

Vi har i det <strong>for</strong>egående tilnærmet kurvelængden med længden af en indskreven polygon, hvilket<br />

intuitivt virker rimeligt. Det svarer jo til at <strong>for</strong>søge at måle længden af en krum kurve med en<br />

lineal. Denne del af projektet handler om <strong>for</strong>bindelsen mellem kurvelængde og polygonlængde.<br />

Vi skal vise, at når antallet af hjørnepunkter i Γ går mod uendelig, så går længden af polygonen<br />

mod kurvelængden.<br />

Kurvelængden, L (γ), er kendt fra bogen Matematisk Analyse 2. For en kurve γ : [a;b] → R 3 er<br />

kurvelængden givet ved:<br />

L (γ) =<br />

� b<br />

a<br />

� ˙γ(t)� dt (3.1)<br />

Opgave 4. En indskreven polygon med n hjørnepunkter har længden ℓn, givet ved:<br />

n �<br />

�x(ti) ℓn = − x(ti−1) �2 �<br />

+ y(ti) − y(ti−1) �2 �<br />

+ z(ti) − z(ti−1) �2 .<br />

∑<br />

i=2<br />

1. Vis ved hjælp af middelværdisætningen 1 anvendt på hver af funktionerne x,y,z, at ℓn<br />

også kan skrives på <strong>for</strong>men:<br />

ℓn =<br />

n �<br />

∑ ˙x<br />

i=2<br />

2 (ζi)+ ˙y 2 (ηi)+ ˙z 2 (ξi) ·(ti −ti−1) , (3.2)<br />

med passende valg af værdier ζi, ηi og ξi i intervallet ]ti−1;ti[ (<strong>for</strong> i = 2,...,n) .<br />

2. Opstil middelsummen, Sn <strong>for</strong> integralet (3.1). Hvorved adskiller dette udtryk sig fra<br />

(3.2)?<br />

3. Forklar med ord hvor<strong>for</strong> det virker rimeligt at:<br />

|ℓn − Sn| → 0 <strong>for</strong> n → ∞<br />

Hvilken konsekvens har det? (Vink: Hvad går Sn mod?)<br />

1 Se evt. hjemmmesiden: http://en.wikipedia.org/wiki/Mean_Value_Theorem#Formal_statement<br />

<strong>Matematik</strong> 1 – 05/06 side 5


Efter at have redegjort teoretisk <strong>for</strong>, at ℓn er en tilnærmelse til L (γ) vil det være interessant<br />

at undersøge, hvor god en tilnærmelse den egentligt er, og hvad der er den bedste tilnærmelse,<br />

polygonlængden eller middelsummen.<br />

Opgave 5. Denne opgave skal regnes med Maple. I vil få brug <strong>for</strong> kommandoerne:×Ô��ÙÖÚ�, ×�ÕogÔÓ�ÒØÔÐÓØ�.<br />

1. Udregn ℓn og Sn <strong>for</strong> en selvvalgt2 kurve <strong>for</strong> <strong>for</strong>skellige finheder af inddelingen (dvs.<br />

<strong>for</strong>skellige værdier af n).<br />

4 DEL C<br />

2. Plot <strong>for</strong> hver finhed den tilhørende polygonkurve sammen med kurven selv.<br />

3. Udregn kurvelængden. Hvad er den bedste tilnærmelse, ℓn eller Sn?<br />

4. Undersøg om dette gælder generelt. Prøv med nogle <strong>for</strong>skellige <strong>kurver</strong>, f.eks<br />

γ1(t) = (t 2 ,t 4 ,t 6 ), t ∈ [−1;1] ,<br />

γ2(t) = � cos(t),sin(t),t � , t ∈ [0;10] .<br />

5. Kan man ved at betragte billederne af polygon<strong>kurver</strong>ne <strong>for</strong>klare, hvor<strong>for</strong> Sn hhv. ℓn<br />

er den bedste tilnærmelse?<br />

4.1 Variationsregning<br />

I afsnit A har vi arbejdet med kurvelængden af polygon<strong>kurver</strong> og har betragtet længden som en<br />

funktion af nogle punktplaceringer. Vi skal nu se på] hvordan dette kan generaliseres til at man<br />

betragter kurvelængden som en funktion af kurven og at man minimerer over alle mulige (glatte)<br />

<strong>kurver</strong>.<br />

For en given kurve γ er kurvelængden L (γ) er et reelt tal. Vi opfatter nu kurvelængden som en<br />

afbildning L : C∞ ([a,b],R) → R , der <strong>for</strong> en glat kurve γ : [a,b] → R3 (dvs. γ ∈ C∞ ([a,b],R))<br />

giver værdien af længden af kurven:<br />

L (γ) =<br />

� b<br />

hvor L( ˙γ(t)) = � ˙γ(t)� kaldes Lagrangefunktionen.<br />

a<br />

L( ˙γ(t))dt (4.1)<br />

L er altså ikke en almindelig funktion, der har en delmængde i R n som definitionsmængde; da<br />

L (γ) er en reel værdi kalder man afbildningen <strong>for</strong> en funktional.<br />

4.1.1 Variationer<br />

For en differentiabel funktion af flere variable har vi begrebet retningsafledede og et tilsvarende<br />

begreb kan defineres <strong>for</strong> funktionaler; man udnytter her, at , C ∞ ([a,b],R) er et vektorrum (jf<br />

Lineær Algebra).<br />

2 For inspiration, se hjemmsiden: http://mathworld.wolfram.com/topics/Curves.html<br />

<strong>Matematik</strong> 1 – 05/06 side 6


Definition 4.1. For to glatte <strong>kurver</strong> γ,v ∈ C∞ ([a,b],R 3 ) med fastholdte endepunkter, γ(a) = A,<br />

γ(b) = B og v(a) = v(b) = 0 og en funktional L : C∞ ([a,b],R) → R er variationen af L <strong>for</strong><br />

kurven γ i retningen v givet ved<br />

L (γ + εv) −L (γ)<br />

δL (γ,v) = lim<br />

=<br />

ε→0 ε<br />

∂<br />

L (γ + εv)<br />

∂ε<br />

⏐<br />

ε=0<br />

Man kan tænke på δL på følgende måde: For et fastholdt γ varieres L med funktionen v ved<br />

at gå stykket ε i v’s “retning”, og så lade ε → 0; man udnytter her, at funktioner kan “lægees<br />

sammen” (vektorrumsstrukturen af C ∞ ([a,b],R)). Samtidig giver udregningen L (γ + εv) <strong>for</strong><br />

faste valg af γ og v en sædvanlig funktion af en reel variabel ε; dvs. L (γ + εv) differentieres på<br />

den "normale" måde m.h.t. ε.<br />

Opgave 6. Vis at variation <strong>for</strong> L (γ) = � b a � ˙γ(t)�dt , med γ(t) = (t,y(t)), bliver<br />

δL (γ,v) =<br />

� b<br />

a<br />

˙y<br />

� 1+ ˙y(t) ˙v(t)dt .<br />

Opgave 7. Udregn variationen <strong>for</strong> følgende <strong>kurver</strong> ved at bruge ligning (4.1) og definition 4.1.<br />

4.1.2 Minimum<br />

γ2D(t) = (x(t),y(t)) , γ3D(t) = (x(t),y(t),z(t)) . (4.2)<br />

Når vi søger en geodæt, søger vi den kurve, som har den kortest kurvelængde, altså den γ0 <strong>for</strong><br />

hvilket funktionalen L er minimeret. For et minimum L (γ0) må der gælde, at (hvor<strong>for</strong>?):<br />

L (γ0) ≤ L (γ0 + εv), <strong>for</strong> alle v og alle ε . (4.3)<br />

Opgave 8. Vis, at med den variation, der er defineret i definition 4.1, er en nødvendig betingelse<br />

<strong>for</strong> et minimum <strong>for</strong> L , at<br />

δL (γ0,v) = 0, <strong>for</strong> alle v . (4.4)<br />

4.2 Euler-Lagrange ligningerne <strong>for</strong> geodæter<br />

Vi vender nu tilbage til variationen af L ; vi vil gerne udnytte ligning (4.4) og definition 4.1 til<br />

at finde et udtryk der måske gør det muligt (nemt?) at finde minima <strong>for</strong> L .<br />

4.2.1 Euler-Lagrange ligningen <strong>for</strong> γ(t) = (t,y(t))<br />

Lad γ være en plan kurve i R 2 , givet ved grafen <strong>for</strong> en funktion t ↦→ y(t), således at γ(t)= (t,y(t)).<br />

Langragefunktionen bliver nu<br />

L( ˙γ(t)) = � ˙γ� =<br />

�<br />

˙t 2 + ˙y 2 �<br />

(t) = 1+ ˙y 2 (t) ,<br />

og den kan der<strong>for</strong> skrives som L( ˙y(t)). Vi vil nu udlede den såkaldte Euler-Lagrange ligning <strong>for</strong><br />

en Lagrangefunktion på denne <strong>for</strong>m.<br />

<strong>Matematik</strong> 1 – 05/06 side 7


Opgave 9. Vis (uden at bruge det eksplicitte udtryk <strong>for</strong> L( ˙γ(t))= � ˙γ(t)�), at med γ(t)= (t,y(t)),<br />

er δL (γ,v) givet som:<br />

� b<br />

δL (γ,v) = ˙v(t) ∂L<br />

( ˙y(t))dt .<br />

∂ ˙y<br />

(4.5)<br />

Opgave 10. Benyt delvis integration på resultatet i ligning (4.5) til at udlede, at<br />

δL (γ,v) =<br />

� b<br />

a<br />

a<br />

− d ∂L<br />

( ˙y(t))v(t)dt (4.6)<br />

dt ∂ ˙y<br />

Vink: Benyt, at værdierne af v i endepunkterne t = a og t = b er kendte.<br />

For at komme videre er det nødvendigt at bruge følgende resultat:<br />

Lemma 4.2 (Det fundamentale lemma i variationsregning). Lad f og g være glatte funktioner<br />

(altså f,g ∈ C∞ ([a,b],R)). Så gælder der, at<br />

� b<br />

f(t)g(t)dt = 0 <strong>for</strong> alle g ,<br />

hvis og kun hvis<br />

a<br />

f(t) = 0, t ∈ [a,b].<br />

Ved at bruge lemma 4.2 på ligningerne (4.4) og (4.6) kan man konkludere, at en nødvendig<br />

betingelse <strong>for</strong> et minimum af L i dette tilfælde er:<br />

d ∂L<br />

( ˙y(t)) = 0<br />

dt ∂ ˙y<br />

Dette er den såkaldte Euler-Lagrange ligning <strong>for</strong> en Lagrangefunktion L( ˙y(t)). Euler-Lagrange<br />

ligningen er en differentialligning. For at opsummere har vi at:<br />

Sætning 4.3. (Euler-Lagrange) Lad γ : [a,b] → R 2 være en glat kurve, der opfylder γ(a) = A og<br />

γ(b) = B. Desuden er γ givet ved γ(t) = (t,y(t)).<br />

Hvis γ er et minimum <strong>for</strong> L med Lagrangefunktion L( ˙y(t)), da opfylder y(t) følgende ligning:<br />

d ∂L<br />

( ˙y(t)) = 0 (4.7)<br />

dt ∂ ˙y<br />

Det er værd at understrege, at sætningen kun giver en nødvendig betingelse <strong>for</strong> at γ er et minimum.<br />

Ved at løse ligning (4.7) findes generelt også hvad der svarer til saddelpunkter og maksima.<br />

Ud fra alle løsninger til (4.7) må man så identificere geodæterne ved f.eks. at bruge ligning (4.3);<br />

dette bliver dog ikke nødvendigt i det følgende. Bemærk også, at det mange tilfælde er overordentligt<br />

vanskeligt at finde løsninger til Euler-Lagrange ligningen.<br />

Bemærkning 4.4. Ligning (4.7) kan uden problemer integreres med hensyn til den variable t.<br />

Dermed får Euler-Lagrange ligningen <strong>for</strong>men<br />

hvor C er en arbitrær konstant.<br />

∂L<br />

( ˙y(t)) = C , (4.8)<br />

∂ ˙y<br />

<strong>Matematik</strong> 1 – 05/06 side 8


4.2.2 Geodæter i R 2<br />

y(t)<br />

A<br />

γ = (t,y(t))<br />

Figur 3: .<br />

Vi vil nu prøve at finde en geodæt mellem to punkter A og B i R 2 .<br />

Vi antager, at kurven γ skal være grafen <strong>for</strong> en funktion, således at γ(t) = (t,y(t)), A = (a,yA)<br />

og B = (b,yB), se figur 3. Vi har der<strong>for</strong> fra <strong>for</strong>rige afsnit, at Lagrangefunktionen bliver L( ˙y(t)) =<br />

� 1+ ˙y 2 (t).<br />

Opgave 11. Vis, at Euler-Lagrange ligningen (4.7) <strong>for</strong> ovenstående Lagrangefunktion bliver<br />

�<br />

1<br />

(1+ ˙y 2 (t)) 3/2<br />

�<br />

¨y(t) = 0<br />

Opgave 12. Forklar at vi nødvendigvis må have at<br />

og find <strong>for</strong>skriften <strong>for</strong> kurven γ.<br />

4.2.3 <strong>Geodætiske</strong> <strong>kurver</strong> på cylinderen<br />

¨y(t) = 0<br />

En cylinder med radius 1 har parameterfremstillingen:<br />

hvor θ ∈ [0;2π[ og z ∈ R<br />

F∞(θ,z) = � cosθ,sinθ,z � ,<br />

Cylinderen ligner planen, idet man ved at skære cylinderen op og folde den ud kan få et stykke<br />

af planen. Når en geodætisk kurve i planen der<strong>for</strong> kan antages at være graf <strong>for</strong> en funktion, så<br />

må noget lignende gøre sig gældende på cylinderen.<br />

Opgave 13. Lad os antage, at den geodætiske kurve opfylder, at z er en funktion af θ. Ved<br />

at indsætte dette i cylinderens parameterfremstilling kan de <strong>kurver</strong> γ, som vi undersøger<br />

beskrives således:<br />

γ(θ) = � cosθ,sinθ,z(θ) �<br />

<strong>Matematik</strong> 1 – 05/06 side 9<br />

B<br />

t


a. Opskriv udtrykket <strong>for</strong> kurvelængden af γ(θ), og bestem z(θ) ved hjælp af resultater<br />

fra tidligere i projektet, idet randbetingelserne er givet ved:<br />

hvor θa,θB ∈ [0;2π[ og θa < θb.<br />

z(θa) = zA, z(θb) = zB,<br />

b. Plot <strong>for</strong> <strong>for</strong>skellige θa, θb, zA og zB geodætiske <strong>kurver</strong> <strong>for</strong> cylinderen. Hvad sker der<br />

når zA = zB, og hvad sker der når zA �= zB ?<br />

c. Betragt to <strong>for</strong>skellige punkter på cylinderen, der opfylder at θa = θb. Vis, at de kan<br />

<strong>for</strong>bindes af en ret linie, der er indeholdt i fladen. Er denne rette linie en geodætisk<br />

kurve? Hvis ja, hvor<strong>for</strong> er dette tilfælde så ikke dækket af spørgsmål (a)?<br />

4.2.4 Geodæter på nogle flere flader i R 3<br />

Opgave 14. Betragt fladen i R 3 med parameterfremstillingen<br />

F∈(x,y) = � x,y, 1 3 x3� , (x,y) ∈ R 2 .<br />

Antag, at vi kan finde en geodæt på fladen, der svarer til en kurve (t,y(t)) i R 2 .<br />

i. Find udtrykket <strong>for</strong> længden af en sådan kurve på fladen og opskriv Euler-Lagrange<br />

ligningen <strong>for</strong> kurven. Dette er en differentialligning. Hvad er ordenen? Er den lineær?<br />

ii. Find løsninger til denne Euler-Lagrange ligningen. Man kan her vælge at lade kurven<br />

være bestemt af værdier af y(a) og ˙y(a) (fx med a = 0) (giver dette mening?). Det er<br />

sandsynligt, at man må benytte numerisk løsning, se maples�×ÓÐÚ�ÒÙÑ�Ö�.<br />

iii. Sammenlign nogle af <strong>kurver</strong>ne yg fundet i (ii) med de tilsvarende <strong>kurver</strong>, der stammer<br />

fra at punkterne (a,yg(a)) og (b,yg(b)) <strong>for</strong>bindes med en ret linie i (x,y)-planen.<br />

iv. Er der typer af geodæter man ikke kan finde med denne fremgangsmåde?<br />

Opgave 15. Betragt paraboloiden i R 3 med parameterfremstillingen<br />

F∋(x,y) = � x,y, 1 2 (x2 + y 2 )) , (x,y) ∈ R 2 .<br />

Gennemfør samme undersøgelser som i opgave 14. Betragt <strong>for</strong> eksempel geodæter, der<br />

svarer til <strong>kurver</strong> (r(φ)cosφ,r(φ)sinφ) i R 2 (med φ ∈ [a,b]). Bemærk, at der i dette tilfælde<br />

skal opskrives en ny <strong>for</strong>m <strong>for</strong> Euler-Lagrange ligning da L nu afhænger af både r(φ) og<br />

r ′ (φ).<br />

Opgave 16. Gennemfør samme undersøgelser som i opgave 14 <strong>for</strong> en flade af eget valg.<br />

Opgave 17. Overvej hvad der skal til <strong>for</strong> at finde geodæter på <strong>for</strong>men<br />

på en flade i rummet.<br />

�<br />

�<br />

x(t),y(t),F(x(t),y(t))<br />

<strong>Matematik</strong> 1 – 05/06 side 10

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

Saved successfully!

Ooh no, something went wrong!