How to plot parameter phase plane?

9 次查看（过去 30 天）

Sama 2018-7-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/408471-how-to-plot-parameter-phase-plane

评论： Sama 2018-7-3

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