How to plot parameter phase plane?
Show older comments
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
Answers (0)
Categories
Find more on Annotations in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!