Carsten Svaneborg
De to forgående artikler
har handlet om at løse
optimeringsproblemer ved hjælp af algoritmer, hvor en funktion er
kendt som funktion af en mængde af kontinuerte og/eller diskrete
variable, men hvor de afledte af funktionen ikke er kendt eller ikke
er defineret. I det følgende vil vi antage, at der er et kontinuum af
ukendte variable givet ved en differentiabel funktion f(x), samt en
funktional. Vi leder efter den funktion, der minimerer/maksimerer
funktionalen.
Der er mange eksempler på problemer, der kan formuleres på denne måde, for eksempel Fermats princip, der siger, at lyset fra en lyskilde i punktet A, der bliver detekteret i B, vil følge den kurve fra A til B, hvor lyset når hurtigst (eller langsomst) frem, og af Fermats Princip følger geometrisk optik. Af samme skuffe er en geodæt, der er den korteste rute mellem to punkter. En klassiker er brachistochronproblemet (gr. kortest tid), der var det første variationsproblem og blev stillet af Jakob Bernoulli i 1696. Populært formuleret lyder det som følger: Hvilken profil skal en bakke have, for at en bold, der starter fra toppen, hurtigst vil nå bunden (når man ser bort fra rotation og friktion)?
Givet et linært normeret rum
af passende differentiable
funktioner (her
) og en funktional
, kan partiel differentation generaliseres til
Hvis ovenstående eksisterer for alle
, kaldes
udtrykket for gâteaux-variationen af funktionalen;
er
mængen af funktioner, man varierer over, og antages at være en lineær
delmængde af
, idet man oftest vil have bestemte kriterier
for h. I det følgende er
, svarende til at endepunkterne i variationen
fastholdes.
Et interessant eksempel (folk, der har haft analytisk mekanik, vil
genkende eksemplet og notationen): Hvis
, hvor L(x,y,z) er en differentiabel
funktion, og
, så er gâteaux-variationen:
Fra den 3. til den 4. ligning er udtrykket for differention af produktfunktioner brugt, og i den sidste ligning er brugt, at grænseledet forsvinder da h(a)=h(b)=0.
For en funktion
, der er et ekstremum (relativt til
og normen på dette rum) for funktionalen I[f], gælder
denne ligning kaldes Euler-Lagrange (EL) ligningen, og nu er det
praktisk at variationsregningens fundementallemma giver (for
), at
dvs., at man kan slippe af med h(t) og integralet på letteste vis.
Med funktionalen fra eksemplet indsat, giver EL og lemmaet følgende:
Denne ligning kaldes også for Euler-Lagrange-ligningen for funktionalen S[q(t)].
Ligningen kan umiddelbart generaliseres til det tilfælde, hvor
, dvs., funktionalen
er afhængig af N funktioner i stedet
for kun en funktion, og resultatet er N EL-ligninger, med samme
form som den forrige EL-ligning, bare krydret med j'er:
Eksemplet er ikke valgt tilfældigt. Beskrives et fysisk system ved N
koordinater
i et vilkårligt koordinatsystem (fx sphæriske
koordinater), og identificeres
(Lagrange funktionen) med forskellen i kinetisk og potentiel
energi (til tiden t i punktet
) udtrykt i det valgte
koordinatsystem, da kaldes
for aktion-funktionalen (og det er
ikke tilfældigt, at den har samme dimension som
(energi gange
tid, red.)), og sættet af Euler-Lagrange-ligningerne er en
generalisering af Newtons bevægelsesligning til det valgte
koordinatsystem.
Hvis et koordinatsystem har en metrisk tensor
, dvs.,
at pythagoras's læresætning i et punkt
for et differentielt
skridt
er
, da er
Lagrange funktionen for en fri partikel
, og
, men
funktionalen
har
samme extrema, da EL-ligningen vil blive multipliceret med
(hvis T>0), som derfor kan divideres ud.
Men
er længden af den kurve, partiklen vil følge.
Extrema for denne funktional er faktisk geodæter, dvs., at en fri
partikel vil følge geodæter. Dette udsagn er korrekt både i klassisk
mekanik såvel som generel relativitetsteori, og kan fx observeres ved
en
afbøjning af lys fra stjerner set tæt ved solens
overflade.
For en partikel beskrevet af 1 (kartesisk) dimension (dvs.
q(t)=x(t)), og i et potential U(x) er
, indsat i EL-ligningen:
Dvs., EL-ligningerne reducerer til Newtons anden lov i kartesiske koordinater. Det er ofte praktisk at kunne anvende et smart koordinatsystem, i stedet for kartesiske koordinater. Fx et pendul kræver 2 koordinater, og i kartesiske koordinater er Newton kedelig at løse, men i polære koordinater er problemet trivielt.
Nu er de forgående eksempler mekaniske, men EL-ligningen gælder for
enhver funktional, der kun indeholder
og
, og det er
let at generalisere EL til højere afledede. For problemerne nævnt i
indledningen er (i 2D kartesiske koordinater) Fermats princip
,
hvor T er passagetiden (lyshastigheden c=1) og n(x,y) er
brydningsindekset, og både geodæt- og brachistochronfunktionalerne er
specialtilfælde med n(x,y(x))=1 hhv.
.
Der findes mange andre problemer, der kan formuleres som variationsproblemer udover al klassisk mekanik (inkl. Maxwells ligninger og generel relativitetsteori). En anden gruppe er kontrolproblemer, hvor man ønsker en kontrolfunktion, der maksimerer ens kontrol over et kendt system, fx en laserpulsform, der selektivt eksiterer atomer til en ønsket tilstand, eller bryder en bestemt binding i et molekyle. Et andet og meget visuelt eksempel er sæbebobler på en ståltrådsfigur (hvor eksperimentet nok er mere gennemsigtigt end matematikken) hvor overfladeenergien skal minimeres, mens volume skal være konstant, og der skal være passende grænsebetingelser svarende til den ståltrådsfigur, som boblen sidder på.