答案很简单:matplotlib
动画不会擦除帧之间的图像。关键是您自己必须更改屏幕上对象的属性。现在,当您执行以下操作时,您可以使用一些新对象绘制新图像plt.scatter
in circle
.
我更改了代码中的几行以避免添加新对象,请参阅标记为的注释行####
。现在应该快一点了。 (尽管月球正在逃离地球的引力场。遗憾。)
import numpy as np
import matplotlib.pyplot as plt
import math
import matplotlib.animation as animation
er = 6378100*10#m #earth radius
mr = 1737400*10#m #moon radius
em = 5.97219*10**24#kg #earth mass
mm = 7.34767309*10**22#kg #moon mass
d = 384400000#m #distance earth-moon
G = 6.67384*10**(-11) #gravity constant
mv = -1023#m/s #Moon velocity
nts = 10000 #no. time steps
def circle(r, h, k, a):
x = r*math.cos(a)+h
y = r*math.sin(a)+k
#### CHANGED
moony.center = x,y
def simData():
tmax = 10000*nts
ts = 10000
x = 0.0
t = 0.0
while t < tmax:
n = 0
for i in range(120):
circle(mr, d, mv*t, n)
n = n + math.pi/60
t = t + ts
yield x, t
def simPoints(simData):
x, t = simData[0], simData[1]
time_text.set_text(time_template%(t))
line.set_data(t, x)
return line, time_text
fig = plt.figure()
ax = plt.axes(xlim=(-430000000, 430000000), ylim=(-430000000, 430000000))
#### CHANGED: a grey circle of moony dimensions to be moved around
moony = plt.Circle((0,0), mr, facecolor=(.8,.8,.8))
ax.add_artist(moony)
time_template = 'Time = %.1f s' # prints running simulation time
time_text = ax.text(0.05, 0.9, '', transform=ax.transAxes)
ani = animation.FuncAnimation(fig, simPoints, simData, blit=False,\
interval=10, repeat=True)
plt.show()
当然,您可能还想创建一个圆圈来说明地球。您不需要有任何plt.plot
如果您只想绘制两个对象,请在文件中添加命令。