how to draw a trajectory

Question

mostafa 2014-4-22

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/126638-how-to-draw-a-trajectory

评论： mostafa 2014-4-22

采纳的回答： Mischa Kim

Hi

i have a differential equation like this:

x'' + x + x^2 = 0

its a chaos problem and i want to display the trajectories.it has two fix point, a saddle at (-1,0) and a repellor at (0,0).

we can transform the equation to ODE by this:

x'= y

y'= -x -x^2

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

mostafa 2014-4-22

在 MATLAB Online 中打开

if true
  x0 = input('Enter the initial position [x y] - ');
  tspan=[0,10*pi];
  [t,x]=ode45(@simple,tspan,x0,[]);
  plot(x(:,1),x(:,2));
  title('Lorenz Model Trajectory XZ')
  xlabel('X')
  ylabel('Z')
end

but it give me error. i write equations in simple.m .

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Mischa Kim 2014-4-22

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/126638-how-to-draw-a-trajectory#answer_134128

在 MATLAB Online 中打开

Mostafa, try something like

 function my_phase()
 IC = rand(10,2);
 hold on
 for ii = 1:length(IC(:,1))
    [~,X] = ode45(@EOM,[0 2],IC(ii,:));
    x = X(:,1);
    y = X(:,2);
    plot(x,y,'r')
 end
 xlabel('x')
 ylabel('y')
 grid
 end
 function dZ = EOM(t, z)
 dZ = zeros(2,1);
 x  = z(1);
 y  = z(2);
 dZ = [  y;...
       - x - x^2];
 end