Gewone differentiaalvergelijkingen: Lijnelementveld en oplossingskrommen met MATLAB
Tekenen van een lijnelemententveld met integraalkrommen [MATLAB uitwerking]
We bekijken de differentiaalvergelijking \[\frac{\dd y}{\dd t}=\frac{4-2t}{3y^2-5}\]
Opdrachten
- Bepaal de algemene oplossing van de differentiaalvergelijking in impliciete vorm via de methode van scheiden van variabelen.
- Teken een lijnelementenveld en teken daarin m.b.v de
contour
functie enkele integraalkrommen.
Uitwerking
- We schrijven de differentiaalvergelijking in differentiaalvorm, waarbij we gelijk scheiding van variabelen toepassen: \[(3y^2-5)\,\dd y=(4-2t)\,\dd t\] Als we aan beide kanten integreren krijgen we \[\int (3y^2-5)\,\dd y=\int (4-2t)\,\dd t\] Dus: \[y^3 -5y=4t-t^2+C\] voor zekere constante \(C\). Met andere woorden \[F(t,y)=C\] voor \[F(t,y)=y^3-5y+t^2-4t\]
- Er wordt gevraagd om een lijnelementenveld en enkelen integraalkrommen te tekenen in één diagram. Het laatste aspect betreft het tekenen van contourgrafiek van de functie \[F(t,y)=y^3-5y+t^2-4t\] De volgende code levert onderstaand diagram op:
>> clear
>> [t,y] = meshgrid(-3:0.5:7, -4:0.5:4); >> dy = (4-2*t)./(3*y.^2-5); >> dt = ones(size(dy)); >> L = sqrt(dt.^2 + dy.^2); >> dyu = dy./L; >> dtu = dt./L; >> figure
>> % lijnelementenveld tekenen >> quiver(t, y, dtu, dyu, 0.3, 'r'), axis tight
>> xlabel 't', ylabel 'y'; >> hold on % ga door met tekenen in hetzelfde diagram
>> % integraalkrommen tekenen in blauwe kleur >> [t,y] = meshgrid(-3:0.1:7, -4:0.1:4); >> contour(t,y, y.^3 - 5*y + t.^2 - 4*t, [-8,-4, 0, 4, 8], 'b');
De integraalkrommen beschrijven meer dan één oplossingskromme. De buitenste integraalkromme doet dit voor 3 oplossingskrommen die afhangen van een beginwaarde. We onderscheiden als oplossingskrommen
- het bovenste deel met \(y\)-waarden groter dan \(\frac{1}{3}\sqrt{15}\);
- het middelste deel met \(y\)-waarden tussen \(-\frac{1}{3}\sqrt{15}\) en \(\frac{1}{3}\sqrt{15}\);
- het onderste deel kleiner dan \(-\frac{1}{3}\sqrt{15}\)
Ontgrendel volledige toegang