Gewone differentiaalvergelijkingen: Lijnelementveld en oplossingskrommen met MATLAB
Tekenen van een lijnelemententveld met integraalkrommen [MATLAB uitwerking]
We bekijken de differentiaalvergelijking
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: Als we aan beide kanten integreren krijgen weDus:voor zekere constante . Met andere woordenvoor
- 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 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 -waarden groter dan ;
- het middelste deel met -waarden tussen en ;
- het onderste deel kleiner dan
Ontgrendel volledige toegang