Springen naar inhoud

Plaatjes gezocht van oplossingen DV


  • Log in om te kunnen reageren

#1

*_gast_Bartjes_*

  • Gast

Geplaatst op 17 juli 2013 - 18:34

In een ander topic ben ik vastgelopen op het tekenen van de oplossingen van een differentiaalvergelijking voor een aantal waarden van een daarin voorkomende constante c.

Dit is het topic:

http://www.wetenscha...-chain-fenomeen

En dit zijn mijn relevante berichtjes:

De oplossing van de hieronderstaande zuiver wiskundige (want dimensieloze) differentiaalvergelijking (*) die als beginpunt (u,z) = (0,0) heeft, die met een positieve afgeleide dz/du start en als positieve constante c bevat schrijven we als z = boog(u;c) :
LaTeX



Ik heb niet de middelen om het gevraagde plaatje (voor bijvoorbeeld c = 1/8, 1/4, 1/2, 1, 2, 4 en 8) te maken. Ook op het internet vond ik geen makkelijke tool om dat voor elkaar te krijgen. Wie weet raad?


Je kan die differentiaalvergelijking splitsen in twee wortels voor respectievelijk de opgaande en neergaande flank van de grafiek. Met een sign-functie kunnen die dan weer aan elkaar geplakt worden.

Hoe dan ook, de numerieke programmatjes die ik op internet vond lopen steeds vast in de top.

Om eens wat goeds in huis te halen heb ik vervolgens onderstaande programma geïnstalleerd:

http://nl.wikipedia....er_math_toolbox

Maar dat is een studie op zich, zodat het mij daarmee ook nog niet gelukt is.

Wie weet een eenvoudige oplossing?

Veranderd door Bartjes, 17 juli 2013 - 18:43


Dit forum kan gratis blijven vanwege banners als deze. Door te registeren zal de onderstaande banner overigens verdwijnen.

#2

Axioma91

    Axioma91


  • >250 berichten
  • 264 berichten
  • Ervaren gebruiker

Geplaatst op 17 juli 2013 - 20:51

Als je de oplossing hebt, dan kun je die in Matlab stoppen. Je kunt er ook voor kiezen om de DV numeriek op te lossen, maar dat is misschien niet nodig. Als je geen Matlab hebt, dan kun je de functie die je geplot wilt hebben hier posten, dan kijk ik of het me lukt (mits gewenst).

Veranderd door Axioma91, 17 juli 2013 - 20:51


#3

Marko

    Marko


  • >5k berichten
  • 8933 berichten
  • VIP

Geplaatst op 17 juli 2013 - 20:52

Wolfram Alpha? Dit is de oplossing voor dz/du < 0:

LaTeX


De oplossing voor dz/du > 0 vind je op analoge wijze

Cetero censeo Senseo non esse bibendum


#4

*_gast_Bartjes_*

  • Gast

Geplaatst op 17 juli 2013 - 21:16

Met het afleiden van een formule voor de grafiek ben ik gestopt omdat die zo ingewikkeld is dat je daar niets meer aan ziet. Best mogelijk dat die uitkomst van WolframAlpha klopt. Verder moet de integratieconstante daarin nog bepaald worden.

Wellicht moeten we van die formule dan maar een plaatje maken?

Veranderd door Bartjes, 17 juli 2013 - 21:54


#5

*_gast_Bartjes_*

  • Gast

Geplaatst op 18 juli 2013 - 22:54

dvboog1.GIF

dvboog2.GIF

dvboog3.GIF

Ziehier mijn probeersels met Euler Math Toolbox.

Om te beginnen heb ik c=1 gekozen. De absolute waarde onder de wortel is toegevoegd om foutmeldigen te voorkomen. Eerst heb ik de opwaartse flank laten tekenen. Daaruit heb ik xtop geschat, dat komt op ca. 0,6. En vervolgens heb ik het plaatje van de boog-functie met behulp van de sign-functie compleet gemaakt.

Lijkt dit ergens naar?

#6

*_gast_Bartjes_*

  • Gast

Geplaatst op 19 juli 2013 - 22:21

Met onderstaande opdrachten voor Euler Math Toolbox kom ik een heel eind:

c:=C; p:=P; vectorfield("sign(p-x)*sqrt( abs( ( c^2 +1)*( 1 - y)^2 - 1 ) )",0,3,-1.5,1.5):
t=0:0.01:3; c:=C;  p:=P; s=runge("sign(p-x)*sqrt( abs( ( c^2 +1)*( 1 - y)^2 - 1 ) )",t,0); plot2d(t,s,add=1);  plot2d(0,0,>points,>add):

Kies daarin C en bepaal vervolgens P zodat de opgaande en neergaande flank mooi op elkaar aansluiten. Voorbeeldje van een resultaat:

c1.png


Nog mooier zou zijn wanneer we p (=xtop) in c konden uitdrukken. De plaats van de top is uit het plaatje immers niet zo nauwkeurig vast te stellen.

#7

*_gast_Bartjes_*

  • Gast

Geplaatst op 21 juli 2013 - 17:31

Gelukt:

http://www.wetenscha...post__p__967053

De formule om p (=xtop) in c uit te drukken is:

LaTeX .


Bedenk hierbij dat het programma x en y i.p.v. u en z gebruikt. De opdrachten die we nu voor Euler Math Toolbox kunnen gebruiken zijn:

c:=C; p:=( ln( c + sqrt( c^2 + 1) ))/(sqrt( c^2 + 1) ) ; vectorfield("sign(p-x)*sqrt( abs( ( c^2 +1)*( 1 - y)^2 - 1 ) )",0,3,-1.5,1.5):

t=0:0.01:3; c:=C; p:=( ln( c + sqrt( c^2 + 1) ))/(sqrt( c^2 + 1) ) ; s=runge("sign(p-x)*sqrt( abs( ( c^2 +1)*( 1 - y)^2 - 1 ) )",t,0); plot2d(t,s,add=1); plot2d(0,0,>points,>add):

Daarmee kunnen we dan grafiekjes van de boogfunctie voor gekozen waarden C van c tekenen.


Voor c tussen ca. 10-7 en 150 komen er fraaie plaatjes uit.

Veranderd door Bartjes, 21 juli 2013 - 17:58


#8

*_gast_Bartjes_*

  • Gast

Geplaatst op 21 juli 2013 - 20:33

Nog mooier worden de plaatjes voor een iets andere schaal:

c:=C; p:=( ln( c + sqrt( c^2 + 1) ))/(sqrt( c^2 + 1) ) ; vectorfield("sign(p-x)*sqrt( abs( ( c^2 +1)*( 1 - y)^2 - 1 ) )",0,2,-1,1):

t=0:0.01:3; c:=C; p:=( ln( c + sqrt( c^2 + 1) ))/(sqrt( c^2 + 1) ) ; s=runge("sign(p-x)*sqrt( abs( ( c^2 +1)*( 1 - y)^2 - 1 ) )",t,0); plot2d(t,s,add=1); plot2d(0,0,>points,>add):


Zo vinden we dan:


c = 0,001
c-is-0.001.png


c = 0,01
c-is-0.01.png


c = 0.1
c-is-0.1.png


c = 1
c-is-1.png


c = 5
c-is-5.png


c = 10
c-is-10.png


c = 100
c-is-100.png

Veranderd door Bartjes, 21 juli 2013 - 20:54






0 gebruiker(s) lezen dit onderwerp

0 leden, 0 bezoekers, 0 anonieme gebruikers

Ook adverteren op onze website? Lees hier meer!

Gesponsorde vacatures

Vacatures