13 svar
211 visningar
Moffen är nöjd med hjälpen
Moffen 1873
Postad: 3 apr 2020 15:25

Ordinär differentialekvation - Liouvilles ekvation

Hej!

Jag sitter fast på följande ode:

V''(r)+1r·V'(r)=eV(r).

Vänsterledet härleddes från V(x2+y2) där vi ersätter en tvåvariabels funktion med vår funktion V(r)=V(x2+y2), för att hitta en lösning som bara beror av avståndet från origo. Det är egentligen en enklare version av Liouvilles ekvation, U(x,y)=eU(x,y), som jag vill börja med att lösa. Jag började med att försöka hitta en (reell) lösning som bara beror på avståndet från origo, men då kommer jag fram till den här ode:n och fastnar. 

Vet inte riktigt vart jag ska börja, har testat att skriva om den som V''(r)·e-V(r)+e-V(r)·V'(r)r=1 men det verkar inte leda någonstans...?

Har även testat att ansätta en potensserie som lösning, för att på nåt sätt kunna identifiera koefficienter eftersom vi vet vad ex har för potensserie. Men där fastnade jag också ganska snabbt...

All hjälp är uppskattad.

Moffen 1873
Postad: 5 apr 2020 16:48

Jag vet att det finns komplexa lösningar till Liouvilles ekvation, men just här är väl frågan om det finns några reella lösningar, som geometriskt kan tolkas genom att bara veta avståndet från origo (bump).

AlvinB 4014
Postad: 5 apr 2020 19:07 Redigerad: 5 apr 2020 19:11

Wolfram Alpha ger två gräsligt fula, men reella, lösningar:

V1r=ln4e22-A·B·A-22·r22-A-2A-2r22-A+e22-A·B2V_{1}\left(r\right)=\ln\left(\frac{4e^{\sqrt{2}\sqrt{2-A}\cdot B}\cdot\left(A-2\right)^{2}\cdot r^{\sqrt{2}\sqrt{2-A}-2}}{\left(\left(A-2\right)r^{\sqrt{2}\sqrt{2-A}}+e^{\sqrt{2}\sqrt{2-A}\cdot B}\right)^{2}}\right)

V2r=ln4e22-A·B·A-22·r22-A-2r22-A+A-2e22-A·B2V_{2}\left(r\right)=\ln\left(\frac{4e^{\sqrt{2}\sqrt{2-A}\cdot B}\cdot\left(A-2\right)^{2}\cdot r^{\sqrt{2}\sqrt{2-A}-2}}{\left(r^{\sqrt{2}\sqrt{2-A}}+\left(A-2\right)e^{\sqrt{2}\sqrt{2-A}\cdot B}\right)^{2}}\right)

Hur man tar fram dessa framgår tyvärr inte...

Qetsiyah 6503 – Livehjälpare
Postad: 5 apr 2020 19:11 Redigerad: 5 apr 2020 19:12

Är det meningen att du ska lösa den algebraiskt? wow

AlvinB: ja, väldigt fula. WA brukar va bra på att förenkla dessutom

Moffen 1873
Postad: 6 apr 2020 13:29

Ah, tack för svar. Nej jag kan inte riktigt heller se hur man ska komma fram till dom där lösningarna... Inte speciellt fina direkt, och inget jag känner att jag kan koppla till något geometriskt resonemang heller.

AlvinB 4014
Postad: 6 apr 2020 22:45

Satt och klurade lite till på detta idag, men utan att komma fram till hur man kan få fram dessa lösningar. Jag lyckades dock snygga till svaret lite grann. Det visar sig att båda lösningarna ovan ger samma mängd av lösningar när man låter AA och BB variera, så vi kan skriva ihop dem till en. Om man dessutom döper om 22-A\sqrt{2}\sqrt{2-A} till AA får man den snyggare lösningen:

Vr=lneABA4rA-2rA-A22eAB2V\left(r\right)=\ln\left(\frac{e^{AB}A^{4}r^{A-2}}{\left(r^{A}-\frac{A^{2}}{2}e^{AB}\right)^{2}}\right)

Moffen 1873
Postad: 6 apr 2020 23:09 Redigerad: 6 apr 2020 23:09
AlvinB skrev:

Satt och klurade lite till på detta idag, men utan att komma fram till hur man kan få fram dessa lösningar. Jag lyckades dock snygga till svaret lite grann. Det visar sig att båda lösningarna ovan ger samma mängd av lösningar när man låter AA och BB variera, så vi kan skriva ihop dem till en. Om man dessutom döper om 22-A\sqrt{2}\sqrt{2-A} till AA får man den snyggare lösningen:

Vr=lneABA4rA-2rA-A22eAB2V\left(r\right)=\ln\left(\frac{e^{AB}A^{4}r^{A-2}}{\left(r^{A}-\frac{A^{2}}{2}e^{AB}\right)^{2}}\right)

Ja det blev ju lite finare i alla fall. Det får nog bli komplex analys som får lösa problemet helt enkelt. Men tack för försöket :)

AlvinB 4014
Postad: 11 apr 2020 00:52 Redigerad: 11 apr 2020 00:54

Jag tror jag har spenderat lite väl mycket tid på det här nu, men jag har faktiskt lyckats lösa detta (med en knuff i rätt riktning tack vare MathStackExchange). Det gick ut på att transformera ekvationen med hjälp av variabelbytet V(x)=ln(u(x))V(x)=\ln(u(x)) (jag använder xx som oberoende variabel istället för rr) till

u''x-u'(x)2u(x)+u'(x)x=ux2u''\left(x\right)-\dfrac{u'(x)^2}{u(x)}+\dfrac{u'(x)}{x}=u\left(x\right)^2

Att lösa denna ekvation är ganska invecklat. Jag har skrivit upp resonemanget i sin helhet på MathSE, men det gick ut på att göra några kluriga transformationer (bland annat att låta ux=1x2v(x)u\left(x\right)=\frac{1}{x^2v(x)} och t=ln(x)t=\ln(x)) som sedan kunde lösas.

Jag lyckades dessutom få lösningen på ännu enklare form. Den är

Vr=ln(2A2BrA-2(BrA-1)2)V\left(r\right)=\ln(\dfrac{2A^2Br^{A-2}}{(Br^A-1)^2})

för godtyckliga konstanter AA och BB.

Moffen 1873
Postad: 12 apr 2020 14:01

Wow, tack så hemskt mycket AlvinB!

Det var precis något sånt där jag letade efter. Hade bara hoppats att jag kunde kommit på det själv haha.

Det är några väldigt speciella substitutioner man gör, och hur ni kom på dom är helt magiskt tycker jag.

Moffen 1873
Postad: 14 apr 2020 14:02

Jag vet inte om det här borde vara en separat tråd eller inte men det handlar om uppgiften i alla fall:

Hur kommer det sig att 

är en "autonomous equation"? Jag är väldigt ny till det här, och hänger inte riktigt med hur vi ser (eller visar) att det där är en "autonomous equation". 

Min tanke är att vi måste kunna skriva om det på nåt bra sätt till ett system av ode'er där varje ekvation ger en derivata som en godtycklig funktion av funktionen (läs det en gång till). Men jag vet inte riktigt hur man går tillväga i det här fallet.

dr_lund 1177 – Fd. Medlem
Postad: 14 apr 2020 14:15 Redigerad: 14 apr 2020 14:20

En autonom ekvation saknar oberoende variabel i explicit form.

Enkelt exempel: y'=f(y)y'=f(y).

I  ditt exempel: F(v,v',v'')=0F(v,v',v'')=0, där v'=dvdtv'=\dfrac{dv}{dt} osv.

Moffen 1873
Postad: 14 apr 2020 15:00
dr_lund skrev:

En autonom ekvation saknar oberoende variabel i explicit form.

Enkelt exempel: y'=f(y)y'=f(y).

I  ditt exempel: F(v,v',v'')=0F(v,v',v'')=0, där v'=dvdtv'=\dfrac{dv}{dt} osv.

Jaha, tack så mycket dr_lund! Då hade jag helt enkelt missförstått vad jag läste på Wikipedia. Så i ditt/mitt exempel får inte t förekomma i ekvationen?

Så vi kan helt enkelt ha typ vilken funktion som helst av "funktionen", så länge inte dess oberoende variabel förekommer i ekvationen? Typ:

F(y,y',y'') =y-y'y''+1y32=0 skulle alltså vara en autonom ekvation om exempelvis y=y(x)?

Fråga: Får det förekomma några konstanter i vänsterledet (i funktionen F)?

AlvinB 4014
Postad: 14 apr 2020 15:53 Redigerad: 14 apr 2020 15:53
Moffen skrev:
dr_lund skrev:

En autonom ekvation saknar oberoende variabel i explicit form.

Enkelt exempel: y'=f(y)y'=f(y).

I  ditt exempel: F(v,v',v'')=0F(v,v',v'')=0, där v'=dvdtv'=\dfrac{dv}{dt} osv.

Jaha, tack så mycket dr_lund! Då hade jag helt enkelt missförstått vad jag läste på Wikipedia. Så i ditt/mitt exempel får inte t förekomma i ekvationen?

Så vi kan helt enkelt ha typ vilken funktion som helst av "funktionen", så länge inte dess oberoende variabel förekommer i ekvationen? Typ:

F(y,y',y'') =y-y'y''+1y32=0 skulle alltså vara en autonom ekvation om exempelvis y=y(x)?

Fråga: Får det förekomma några konstanter i vänsterledet (i funktionen F)?

Ja. Det får förekomma konstanter, och vad som helst som inte beror av den oberoende variabeln (xx eller tt).

Den där substitutionen u=dvdtu=\frac{dv}{dt} är lite utav en klassiker när man har en autonom ekvation. Det gör att man kan omvandla en ekvation av andra ordningen till en ekvation av första ordningen.

Moffen 1873
Postad: 14 apr 2020 22:18

Då är jag med, tack så mycket! :)

Svara Avbryt
Close