|
|
|
\require{AMSmath}
Reageren...
Re: Machten
Hoe voer je een parameter functie zoals
x'(t)=-1/3(x(t))3+x(t)+y(t)
en y'(t)=-ex(t) epsilon is een zeer kleine constante in matlab?
Antwoord
dag Piet,
Het gaat hier niet om een parameter functie, maar om een stelsel van twee differentiaalvergelijkingen.
Je kunt deze oplossen met de symbolic toolbox.
Hiervoor zou je de volgende commando's moeten invoeren:
clear
syms x0 y0 e
dv1='Dx=-1/3*x^3+x+y'
dv2='Dy=-e*x'
[x,y]=dsolve(dv1,dv2,'x(0)=x0,y(0)=y0')
Helaas geeft Matlab hier niet thuis:
Warning: Explicit solution could not be found.
Een poging om dit stelsel analytisch op te lossen in Maple loopt ook spaak.
Waarschijnlijk moet je hier een numerieke methode gebruiken, en daarvoor moet je de waarde van e en x0 en y0 kennen.
groet,
Gebruik dit formulier alleen om te reageren op de inhoud van de vraag en/of het
antwoord hierboven. Voor het stellen van nieuwe vragen kan je gebruik maken
van een vraag stellen in het menu aan de linker kant. Alvast bedankt!
|
