如何在Python上制作3D图形动画

Question

我正在为我的图形制作动画（如下），由于我没有动画经验，所以我不确定从哪里开始或如何开始。我不确定它的工作方式或代码的结构，因此，如果有人可以提供伪代码或算法，我将不胜感激。我还提供了用于绘制下面的图的代码。

enter code here
from scipy.integrate import odeint
import scipy as sci
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
import matplotlib.animation as ani

# Universal Gravitational Const.

G = 6.674e-11

# Defining Mass

m1 = 2
m2 = 3.5
m3 = 2.3



# Init positions in graph (array)

pos1 = [-5,0,1]
pos2 = [5,0,10]
pos3 = [0,1,3]

p01 = np.array(pos1)
p02 = np.array(pos2)
p03 = np.array(pos3)


# Init velocities (array)

vi1 = [1,0.01,0]
vi2 = [-5,0,1]
vi3 = [0,-1,0]


v01 = np.array(vi1)
v02 = np.array(vi2)
v03 = np.array(vi3)

#Function
def derivs_func(y,t,G,m1,m2,m3):
    d1 = y[:3]
    d2 = y[3:6]   
    d3 = y[6:9]
    v1 = y[9:12]
    v2 = y[12:15]
    v3 = y[15:18]


    dist12 = np.sqrt((pos2[0]-pos1[0])**2 + (pos2[1]-pos1[1])**2 + (pos2[2]-pos1[2])**2)
    dist13 = np.sqrt((pos3[0]-pos1[0])**2 + (pos3[1]-pos1[1])**2 + (pos3[2]-pos1[2])**2)
    dist23 = np.sqrt((pos3[0]-pos2[0])**2 + (pos3[1]-pos2[1])**2 + (pos3[2]-pos2[2])**2)

    dv1dt = m2 * (d2-d1)/dist12**3 + m3 * (d3-d1)/dist13**3 
    dv2dt = m1 * (d1-d2)/dist12**3 + m3 * (d3-d2)/dist23**3 
    dv3dt = m1 * (d1-d3)/dist13**3 + m2 * (d2-d3)/dist23**3 
    dd1dt = v1 
    dd2dt = v2
    dd3dt = v3

    derivs = np.array([dd1dt,dd2dt,dd3dt,dv1dt,dv2dt,dv3dt]) 
    derivs3 = derivs.flatten() 

    return derivs3 

yo = np.array([p01, p02, p03, v01, v02, v03]) 
y0 = yo.flatten()  

time = np.linspace(0,200,500) 

sol = odeint(derivs_func, y0, time, args = (G,m1,m2,m3)) 

x1 = sol[:,:3]
x2 = sol[:,3:6]
x3 = sol[:,6:9]


fig = plt.figure(figsize = (15,15))
ax = fig.add_subplot(111,projection = '3d')

ax.plot(x1[:,0],x1[:,1],x1[:,2],color = 'b')
ax.plot(x2[:,0],x2[:,1],x2[:,2],color = 'm')
ax.plot(x3[:,0],x3[:,1],x3[:,2],color = 'g')

ax.scatter(x1[-1,0],x1[-1,1],x1[-1,2],color = 'b', marker = 'o', s=30, label = 'Mass 1')
ax.scatter(x2[-1,0],x2[-1,1],x2[-1,2],color = 'm', marker = 'o',s=90, label = 'Mass 2')
ax.scatter(x3[-1,0],x3[-1,1],x3[-1,2],color = 'g', marker = 'o',s=60, label = 'Mass 3')
ax.legend()