Matlab上缺少数据

Question

我有一个问题-为什么图中的数据部分丢失了？ 请帮助我，我被困住了，不知道=（

这是ode的代码：

function xDot = parallelDegreesW4(t,x,~,w_af,g_af,gcurr)
xDot  = [x(2); ... %//= theta'

    sin(x(1)) * cos(x(1)) * ((x(4))^2) + ...
    (w_af^2) * ((sin(x(1)))^3 * cos(x(1)) * (cos(x(3)))^4 + (sin(x(1)))^3 * cos(x(1)) * (sin(x(3)))^4 - (cos(x(1)))^3 * sin(x(1))  ) + ...
    - 2 * g_af * x(2) + ...
    + 2 * gcurr * sin(x(1)) * cos(x(1)) * sin(x(3)) * cos(x(3)); ...

    x(4); ... %//= phi'

    - x(2) * x(4) * cot(x(1)) + ...
    (w_af^2) * (sin(x(1)))^2 * cos(x(3)) * sin(x(3)) *(-(cos(x(3)))^2 + (sin(x(3)))^2) + ...
    - 2 * g_af * x(4) + ...
    + gcurr * ((cos(x(3)))^2 - (sin(x(3)))^2)];

这是我运行的代码：

function f = ps(phi_vel,gcurr)
    w_af = 0.1;
    g_af = 0.01;
    %initial positions
    initial_theta = pi/2;
    initial_phi = (3*pi/4+pi/2)/2;
    %initial velocities
    initial_velocity_theta = 0;
    initial_velocity_phi = phi_vel;
    %finaly x0
    x0 = [initial_theta,initial_velocity_theta,initial_phi,initial_velocity_phi];

    time = 1000;
    tspan = 0:0.01:time;
    options = odeset('RelTol',1e-8,'AbsTol',[1e-8 1e-8 1e-8 1e-8]);

    [t,x]=ode45('parallelDegreesW4',tspan,x0, options,w_af, g_af, gcurr);
    xx1 = x(:,1);%theta
    xx2 = x(:,3);%phi
    xx3 = x(:,2);%V_theta
    xx4 = x(:,4);%V_phi
    %

    %
    %%
    ini_1=xx1(1);
    ini_2=xx2(1);
    ini_3=xx3(1);
    ini_4=xx4(1);
    last_1 = xx1(length(xx1)-1);
    last_2 = xx2(length(xx2)-1);
    last_3 = xx3(length(xx3)-1);
    last_4 = xx4(length(xx4)-1);


    figure(1)
  plot(initial_phi,initial_velocity_phi,'ok','MarkerSize',5,'MarkerFaceColor','g');
        hold on;
        plot(last_3,last_4,'ok','MarkerSize',5,'MarkerFaceColor', 'r');
        hold on;
        plot(xx3,xx4,'--','LineWidth' ,4);
        hold on;    

之后，我得到了这样的图片。 （绿色表示-开始，红色-变量phi在相空间中的运动结束）

所以我的问题是，所有必须显示从绿色到红色的运动的数据在哪里？ 这很奇怪，因为实际上是动态的，下面我在其中包含了变量phi与时间有关的图片。 感谢你们对我的帮助！

Answer 1

您正在V-theta和V_phi轴上绘制theta和phi-太多初始条件，这意味着Matlab选择显示theta_0和phi_0的轴被缩小了太多，无法显示折线图。

尝试简单地注释掉第一个绘图命令

%plot(initial_phi,initial_velocity_phi,'ok','MarkerSize',5,'MarkerFaceColor','g');

得到V_phi vs V_theta的图并运行ps(2,1)获得速度的相空间。 或者尝试注释掉figure(1)之后的所有内容并替换为

    plot(t,xx2,'--','LineWidth' ,4);

并运行ps(.2,.3)获得类似您正在寻找的示例的东西。
同样出于测试目的，为了使其运行更快，我将误差范围增加到

options = odeset('RelTol',1e-3,'AbsTol',[1e-3 1e-3 1e-3 1e-3]);

减少的time time = 400; %1000; time = 400; %1000; tspan至tspan = 0:0.05:time; %0:0.01:time; tspan = 0:0.05:time; %0:0.01:time; -如果您已经可以猜测正在寻找的解决方案呈指数衰减，那么继续在非常大的时域内进行求解，尤其是使用固定时间步长的ODE求解器，将学到的东西很少。 我还需要像function xDot = parallelDegreesW4(t,x,w_af,g_af,gcurr)这样删除函数标头中的~ ，以使其在八度中工作。