Python-ode上的奇怪绘图点

作为我硕士论文的一部分，我目前正在为一辆汽车开发一个动态模型

我使用python附带的ode解算器，得到的结果看起来不错。 但我注意到一些情节有问题

我的部分代码如下：

 f = [#0 - x_ctr     
            dx_ctrdt,
        #1 - dx_ctrdt
            (flx_car1_A + flx_car1_B + flx_car2_A + flx_car2_B)/M_car + (x_ctr)*(V_tr**2)/R_curve/R_ctr - F_drag/M_car,
        #2 - y_ctr     
            dy_ctrdt,
        #3 - dy_ctrdt
            (fly_car1_A + fly_car1_B + fly_car2_A + fly_car2_B)/ M_car + (R_curve + y_ctr)*(V_tr**2)/R_curve/R_ctr - 9.8*np.sin(beta),
        #4 - omega_1
            -(flx_car1_A + flx_car1_B)*r/Inertia_1 - 0.5*F_rr*r/Inertia_1 + F_motor_1*r/Inertia_1,
        #5 - omega_2
            -(flx_car2_A + flx_car2_B)*r/Inertia_2 - 0.5*F_rr*r/Inertia_2 + F_motor_2*r/Inertia_2,
        #6 - alpha
            dalphadt,  
        #7 - dalphadt
            (mz_car_1  + mz_car_2 ) / Inertia_car,                            
        #8 - x_1A
            v_1A*cos_theta_1 - (V_tr/R_curve)*R_1A*cos_psi_1A,
        #9 - y_1A
            v_1A*sin_theta_1 - (V_tr/R_curve)*R_1A*sin_psi_1A,
        #10 - x_1B
            v_1B*cos_theta_1 - (V_tr/R_curve)*R_1B*cos_psi_1B,
        #11 - y_1B
            v_1B*sin_theta_1 - (V_tr/R_curve)*R_1B*sin_psi_1B,
        #12 - x_2A
            v_2A*cos_theta_2 - (V_tr/R_curve)*R_2A*cos_psi_2A,
        #13 - y_2A
            v_2A*sin_theta_2 - (V_tr/R_curve)*R_2A*sin_psi_2A,
        #14 - x_2B
            v_2B*cos_theta_2 - (V_tr/R_curve)*R_2B*cos_psi_2B,
        #15- y_2B
            v_2B*sin_theta_2 - (V_tr/R_curve)*R_2B*sin_psi_2B,
        ]

当我绘制f的解的图时，我得到了一个平滑的图

dif_var_initial = [0.000, 0.000,    y_ctr_0,  -0.000, 140.0027, 140.0027, 0.000 ,  0.000,  x_1A_0, y_1A_0,  x_1B_0,  y_1B_0,  x_2A_0,  y_2A_0, x_2B_0, y_2B_0]

def dif_eqts(dif_var, t, kx1, ky1, cy1, M_car):
    global v_ax_1, v_ax_2,R_ctr,psi_ctr,Lamb_tk, A_tk, Time_max, delta_t, y_track_1_vector ,y_track_2_vector,R_curve, dRdt,z,j,v_1A,v_1B,v_2A,v_2B,P

    x_ctr, dx_ctrdt, y_ctr, dy_ctrdt, omega_1, omega_2, alpha, dalphadt, x_1A, y_1A, x_1B, y_1B, x_2A, y_2A, x_2B, y_2B = dif_var 

    #   Track disturbances
    eps_ax_1 = A_tk * np.sin(2.*np.pi*(V_tr/Lamb_tk)*t)
    eps_ax_1_tracker.append(eps_ax_1)
    .
    .
    .

但是，当我使用例如eps_ax_1_tracker绘制一个图形时，我会得到不平滑的奇怪绘图（该程序显然比这更复杂，但我不想让你读太多代码行而感到厌烦）

例如：

更接近：

这种行为正常吗？ 我搜索过类似的问题，但没有找到任何可以帮助我的

我希望你能提出一些建议，并感谢你的帮助

编辑（1）

我得到x_ctr的输出，如下所示：

对于以下部队：

振荡行为与轨道中的不规则性有关。

---

我希望这有帮助

数值解算器在接近但不在解轨迹上的状态下使用“探测”值。粗略地说，这会对轨迹周围的方向场进行采样，以便下一步的计算值可能更准确。在

scipy.integrate.ode

中，通过

sol_out

或类似方式使用密集输出，使用记录的路径获取输出，请参阅，对于

scipy.integrate.odeint

使用时间列表输入参数定义输出中报告的采样点。

您能告诉我们您的数据是什么样的吗？我又放了两个输出图，我不知道您是否需要其他类型的数据。非常感谢您的回答。我一直在探索，我相信你是对的。现在我对正在发生的事情更加肯定了。我会落实你的建议。