How to plot parameter phase plane?
Afficher commentaires plus anciens
Hi, I am having the following code for plotting the phase plane between two variables x and y, however, I would like to draw for parameters, say c and m
function my_phase()
IC = [0.1 0.1];
hold on
for ii = 1:length(IC(:,1))
%vectorfield(@EOM,0:.1:1.8,0:0.1:1);
[t,X] = ode45(@EOM,[0 100],IC(ii,:));
u = X(:,1);
w = X(:,2);
plot(u,w,'b', 'linewidth',3)
% plot(t,u, 'b', 'linewidth',3)
% plot(t,w, 'r', 'linewidth',3)
end
fsize=15;
xlabel('\itx','FontSize',fsize)
ylabel('\ity','FontSize',fsize)
grid
%quiver(X(:,1),X(:,2),u,w,'r')
end
function dX = EOM(~, x)
dX = zeros(2,1);
u = x(1);
w = x(2);
k = 3;m = 3;c=1;
dX = [u*(1-(u/k)-m*w/(1+u));...
w*(-c+m*u/(1+u))];
end
Any hint will be appreciated. Thank you very much, Sanaa Moussa
Réponses (0)
Catégories
En savoir plus sur Annotations dans Centre d'aide et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!