2017-11-25 528 views
-1

我正在编写一个函数,它可以在3D中动画随机游走,但不幸的是代码无法正常工作。情节在哪里,没有错误发生,但没有发生。我正在使用%matplotlib tk
有我的代码:3D动画随机游走[Python]

import numpy as np 
import matplotlib.pyplot as plt 
import mpl_toolkits.mplot3d.axes3d as p3 
import matplotlib.animation as animation 

def path_generator(steps, step): 
    path = np.empty((3, steps)) 
    for i in range(1, steps): 
     x_ran, y_ran, z_ran = np.random.rand(3) 
     sgnX = (x_ran - 0.5)/abs(x_ran - 0.5) 
     sgnY = (y_ran - 0.5)/abs(y_ran - 0.5) 
     sgnZ = (z_ran - 0.5)/abs(z_ran - 0.5) 
     dis = np.array([step*sgnX, step*sgnY, step*sgnZ]) 
     path[:, i] = path[:, i - 1] + dis 

    return path 

def animate(i): 
    global particles, trajectories 
    for trajectory, particle in zip(trajectories, particles): 

    trajectory.set_data(particle[0:2, :i]) 

    trajectory.set_3d_properties(particle[2, :i]) 

return trajectories 

def random_walk_3D_animated(n, traj = 1): 

    fig = plt.figure() 
    ax = p3.Axes3D(fig) 

    particles = [path_generator(n, 1) for i in range(traj)] 

    trajectories = [ax.plot(particle[0, 0:1], particle[1, 0:1], particle[2, 
        0:1])[0] for particle in particles] 

    ax.set_xlim3d([-100, 100]) 
    ax.set_ylim3d([-100, 100]) 
    ax.set_zlim3d([-100, 100]) 

    animacion = animation.FuncAnimation(fig, animate, 1000, interval=50, 
             blit=False) 

    plt.show() 

什么奇怪的是,当没有功能random_walk_3D_animated(n, traj = 1)ntraj给出的数值的代码做的工作。有时代码不会从(0,0,0)开始随机游走。我想知道为什么。

+0

我已经评论说,你需要返回动画参考你的最后一个问题如下。你会不会永远问“随机游走”问题? – ImportanceOfBeingErnest

+0

我知道我一直在问这个问题一段时间,我真的很抱歉。自从上一次以来,我取得了很大进步:)不幸的是,我不知道如何返回对动画的引用 – Hendrra

回答

1
  1. 开始位置将是emty数组的内容。这可能是任何值,所以它在这里并不真正有用。而是用零初始化path
  2. 您需要返回对动画的引用。从animation documentation:“[..]保持对实例对象的引用至关重要。”

完整的示例:

import numpy as np 
import matplotlib.pyplot as plt 
import mpl_toolkits.mplot3d.axes3d as p3 
import matplotlib.animation as animation 

def path_generator(steps, step): 
    path = np.zeros((3, steps)) 
    for i in range(1, steps): 
     x_ran, y_ran, z_ran = np.random.rand(3) 
     sgnX = (x_ran - 0.5)/abs(x_ran - 0.5) 
     sgnY = (y_ran - 0.5)/abs(y_ran - 0.5) 
     sgnZ = (z_ran - 0.5)/abs(z_ran - 0.5) 
     dis = np.array([step*sgnX, step*sgnY, step*sgnZ]) 
     path[:, i] = path[:, i - 1] + dis 

    return path 

def animate(i): 
    global particles, trajectories 
    for trajectory, particle in zip(trajectories, particles): 
     trajectory.set_data(particle[0:2, :i]) 
     trajectory.set_3d_properties(particle[2, :i]) 

def random_walk_3D_animated(n, traj = 1): 
    global particles, trajectories 
    fig = plt.figure() 
    ax = p3.Axes3D(fig) 

    particles = [path_generator(n, 1) for i in range(traj)] 
    trajectories = [ax.plot(particle[0, 0:1], particle[1, 0:1], particle[2, 
        0:1])[0] for particle in particles] 
    ax.set_xlim3d([-100, 100]) 
    ax.set_ylim3d([-100, 100]) 
    ax.set_zlim3d([-100, 100]) 

    animacion = animation.FuncAnimation(fig, animate, 1000, interval=50, 
             blit=False) 
    return animacion 

ani = random_walk_3D_animated(100, traj = 1) 
plt.show() 
+0

非常感谢!现在我完成了。再一次 - 我真的很抱歉,但这是我第一次尝试动画剧情。 – Hendrra