簡體   English   中英

如何阻止第二個 plot 出現在 matplotlib 中?

[英]How to stop second plot from showing up in matplotlib?

我正在嘗試物理系統的 plot 和 animation。 我已經計算出方程式並繪制了圖,但是還有第二個 plot 我不希望它繼續出現,我無法讓它停止。

import numpy as np
from scipy import integrate
import matplotlib.pyplot as plt
from matplotlib import animation, rc
from IPython.display import HTML

# Input constants 
m = 10 # mass (kg)
L = 4 # length (m)
g = 9.81 # gravity (m/s^2)

dt = 0.1 # time step size (seconds)
t_max = 40 # max sim time (seconds)
num_steps = 1 + int(t_max/dt)
theta_0 = np.pi/2 # initial angle (radians)
theta_dot_0 = 0 # initial angular velocity (rad/s) 
state0 = [theta_0,theta_dot_0] 
# Get timesteps
time_index = np.arange(0, t_max + dt, dt)

def derivatives(state, time_index, L=L, m=m, g=g):
    theta_dot = state[1]
    theta_ddot = -g*np.sin(state[0])/L
    return theta_dot,theta_ddot

output = integrate.odeint(derivatives, state0, time_index)
theta = output[:,0]
theta_dot = output[:,1]

fig = plt.figure()
ax = fig.add_subplot(111, autoscale_on=True, xlim=(-2*L, 2*L), ylim=(-2*L, 2*L))
ax.set_aspect('equal')
ax.grid()

line, = ax.plot([], [], 'o-', lw=2)
time_template = 'time = %.1fs'
time_text = ax.text(0.05, 0.9, '', transform=ax.transAxes)

x =  L*np.sin(theta)
y = -L*np.cos(theta)

def init():
    # create empty object to put the points in 
    line.set_data([], [])
    time_text.set_text('')
    return line, time_text

def animate(t):
    x_t = [0, x[t]]
    y_t = [0, y[t]]

    # add the point to the line 
    line.set_data(x_t, y_t)

    # add the time text to plot 
    # time_template will be the text next to thetime
    time_text.set_text(time_template % (t*dt))

    return line, time_text

# we don't want the range of the steps to start at 0 
# start it at one 
ani = animation.FuncAnimation(fig, animate, range(1, num_steps),
                              interval=dt*1000, blit=True, init_func=init)

rc('animation', html='jshtml')

#ani.save('pendulum.mp4', fps=15)

ani

這是 output:

在此處輸入圖像描述

我想擺脫的plot是我用紅色圈起來的那個。 這是我的全部代碼,所以它應該是完全可重現的。

我嘗試了幾種修剪繪圖代碼的變體,但我無法調試它發生的原因。

我怎樣才能擺脫這第二個 plot?

在調用ani之前,一個簡單的plt.close()就可以完成這項工作。

最后幾行:

ani = animation.FuncAnimation(fig, animate, range(1, num_steps),
                              interval=dt*1000, blit=True, init_func=init)

rc('animation', html='jshtml')
#ani.save('pendulum.mp4', fps=15)
plt.close()
ani

演示:

在此處輸入圖像描述

更多信息在這個鏈接。

暫無
暫無

聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.

 
粵ICP備18138465號  © 2020-2024 STACKOOM.COM