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