8 svar
189 visningar
jjnte 3 – Fd. Medlem
Postad: 25 sep 2018 14:14 Redigerad: 25 sep 2018 14:18

Variationskalkyl: Eulers ekvation-problem

Uppgiften består i att bestämma den allmänna lösningen till Eulers ekvation för funktionalen

     J(y)=abF(x,y(x),y'(x))dx=abf(x)1+(y')2dx

     f(x)=y(x)

F beror inte av x så då utgår jag ifrån

     F-y'Fy'=C=konstant

     y1+(y')2-y(y')21+(y')2=C

Jag förenklar och löser ut y'

      dydx=y2-C2C2dx=Cy2-C2dy

Hur går jag vidare efter det här och är jag på rätt spår?

Svaret ska bli

     y=cy2-C2dx

Albiki 5096 – Fd. Medlem
Postad: 25 sep 2018 22:33

Välkommen till Pluggakuten!

Funktionalen

    J(y)=abF(x,y,y')dx\displaystyle J(y) = \int_{a}^{b}F(x,y,y')dx

är stationär när integranden FF uppfyller Euler-Lagranges ekvation

    Fy-ddxFy'=0.\frac{\partial F}{\partial y} - \frac{d}{dx}\frac{\partial F}{\partial y'} = 0.

Integranden F(x,y,y')=y1+(y')2F(x,y,y') = y\sqrt{1+(y')^2} har de partiella derivatorna

    Fy=1+(y')2\frac{\partial F}{\partial y} = \sqrt{1+(y')^2} och Fy'=y'y1+(y')2\frac{\partial F}{\partial y'} = \frac{y'y}{\sqrt{1+(y')^2}}

vilket ger E-L-ekvationen

    1+(y')2-ddxy'y1+(y')2=0.\sqrt{1+(y')^2} - \frac{d}{dx}\frac{y'y}{\sqrt{1+(y')^2}} = 0.

Albiki 5096 – Fd. Medlem
Postad: 25 sep 2018 22:36 Redigerad: 25 sep 2018 22:37

Derivering av den andra termen ger

    ddxy'y1+(y')2=11+(y')2·({y''y+(y')2}1+(y')2-y(y')21+(y')2)\frac{d}{dx}\frac{y'y}{\sqrt{1+(y')^2}} = \frac{1}{1+(y')^2}\cdot(\{y''y+(y')^2\}\sqrt{1+(y')^2}-\frac{y(y')^2}{\sqrt{1+(y')^2}})

Albiki 5096 – Fd. Medlem
Postad: 25 sep 2018 22:47 Redigerad: 25 sep 2018 22:49

Förenkling av denna derivata ger 

    11+(y')2(y''y+(y')2-y(y')21+(y')2)\frac{1}{\sqrt{1+(y')^2}}(y''y+(y')^2 - \frac{y(y')^2}{1+(y')^2})

så att E-L-ekvationen blir 

    1-y''y+y(y')21+(y')2=0y(y')21+(y')2=y''y-1y(y')2=y''y-1+y''(y')2y-(y')21- y''y+\frac{y(y')^2}{1+(y')^2} = 0 \iff \frac{y(y')^2}{1+(y')^2}=y''y-1 \iff y(y')^2=y''y-1+y''(y')^2y-(y')^2

jjnte 3 – Fd. Medlem
Postad: 25 sep 2018 22:57

Tack för dina fina svar, uppskattas verkligen. Det kan hända att jag beskriver problemet lite felaktigt då min kunskap inom detta område är tämligen begränsat. Men då integranden F ej beror av x (F(x,y,y')=F(y,y')), används då inte ett specialfall:

     Fy'Fy'=C=konstant

Albiki 5096 – Fd. Medlem
Postad: 25 sep 2018 23:05
jjnte skrev:

Tack för dina fina svar, uppskattas verkligen. Det kan hända att jag beskriver problemet lite felaktigt då min kunskap inom detta område är tämligen begränsat. Men då integranden F ej beror av x (F(x,y,y')=F(y,y')), används då inte ett specialfall:

     Fy'Fy'=C=konstant

 Jag såg att du skrev det i ditt första inlägg, men jag ser inte hur det kommer från E-L-ekvationen. Att integranden FF ej beror på xx har ingen direkt påverkan på de två derivatorna som ingår i E-L-ekvationen, som du ser. Har du själv kommit fram till "specialfallet" eller har du blivit tillsagd att det är på det sättet av någon annan person?

jjnte 3 – Fd. Medlem
Postad: 25 sep 2018 23:15
Albiki skrev:

 Jag såg att du skrev det i ditt första inlägg, men jag ser inte hur det kommer från E-L-ekvationen. Att integranden FF ej beror på xx har ingen direkt påverkan på de två derivatorna som ingår i E-L-ekvationen, som du ser. Har du själv kommit fram till "specialfallet" eller har du blivit tillsagd att det är på det sättet av någon annan person?

Jag har en pdf med okänt ursprung som används i kursen jag deltar i. Jag delar en snapshot där specialfallen tas upp, där villkor b är det jag hänvisar till.

Albiki 5096 – Fd. Medlem
Postad: 25 sep 2018 23:40

Okej, så man får trixa litet för att få resultatet som du nämner. Men man måste ändå bearbeta resultatet för att få en diffekvation för funktionen yy, så jag fortsätter med mina beräkningar.

Notera att om y(x)=konstanty(x)=konstant så ger beräkningarna ekvationen 0=0-1+0-00=0-1+0-0, vilket visar att funktionen (y')2 (y')^2 aldrig antar värdet noll; det går därför bra att dividera diffekvationen med denna funktion för att få följande. 

    y=y''y'yy'+y''y-1-1(y')2.y=\frac{y''}{y'}\frac{y}{y'}+y''y-1-\frac{1}{(y')^2}.

Notera att (logy')'=y''/y'(\log y')' = y''/y' och (logy)'=y'/y

emmynoether 663 – Fd. Medlem
Postad: 28 sep 2018 15:08

Specialfallet är det som kallas för Beltramis identitet (Beltrami Identity).

Svara
Close