1
svar
87
visningar
Rita fält och fältlinjer
Har följande uppgift (min kod längst ner)
I kompendiet förklaras det inte hur man använder ode45 funktionen. Är därför osäker till allt man stoppar in i ode45, tex. [t,U]=ode45(F,[0 5],[-2;2]); vad är 0 5 här? är det något intervall? [-2;2] verkar vara området vi ritar över?
Sen så får jag en konstig fältlinje som inte är kontinuerlig, bild nedan
%% Uppgift 4
clear,clf,clc
x=linspace(-2,2,15);
[X,Y]=meshgrid(x,x);
quiver(X,Y,X.*Y,X-Y,1)
axis equal, axis([-2 2 -2 2])
F=@(t,u)[u(1).*u(2); u(1)-u(2)];
[t,U]=ode45(F,[0 5],[-2;2]);
hold on
plot(U(:,1),U(:,2),'g','linewidth',2)
[t,U]=ode45(F,[0 5],[0;0.5]);
plot(U(:,1),U(:,2),'g','linewidth',2)
hold off
Kolla på 'help ode45' så ser du vad parametrarna är. Du har ritat två fältlinjer med olika startpunkter ([-2;2] resp [0;0.5]).
De gör liknande grejor här:
http://www.math.chalmers.se/Math/Grundutb/CTH/Matlab/ArkitekturTeknik/ht13/Kurvor.pdf