VPython - example - 模拟斜上抛运动（X-Y- -Z）

本作品采用知识共享署名-非商业性使用-相同方式共享 3.0 Unported许可协议进行许可。允许非商业转载，但应注明作者及出处。

作者：liuyuan_jq

2011-04-10

Fromhttp://wiki.showmedo.com/index.php/PythonThompsonVPythonSeries

---

from visual import *

scene.width = 400
scene.height = 300

scene.autoscale = 0
scene.range = (100,100,100)
scene.center = (50,40,0)

ball = sphere(pos=(0,2,0),radius=2, color=color.green)
ground = box(pos=(50,-1,0),size=(100,2,50))


gravity = 9.8  # m/s**2
velocity = 25   # m/s
angle = 45     # degrees
angle = angle * (pi/180) # converted to radians

# sin = opp / hyp
# cos = adj / hyp

# therefore

# opp = hyp * sin
# adj = hyp * cos

VelocityY = velocity * sin(angle)
VelocityX = velocity * cos(angle)
VelocityZ = 0

VelocityThrown = vector(VelocityX,VelocityY,0)

VelocityWind = vector(2,1,-5)

TotalVelocity = VelocityThrown + VelocityWind

seconds = 0
dt = .01

finished = False
while not finished:
    rate(100) # go thru the loop no more than 100 times/s
    seconds += dt

    # position equation: y(t) = y0 + v0*t + .5 * a * t**2
    ballY = 2 + TotalVelocity.y * seconds - .5 * gravity * seconds**2
    ballX = TotalVelocity.x * seconds
    ballZ = TotalVelocity.z * seconds

    ball.pos = vector(ballX,ballY,ballZ)

    if ballY - 2 <= 0:
        finished = True
        print "initial velocity: " + str(velocity)
        print "angle thrown: " + str(angle)
        print "seconds in flight: " + str(seconds)
        print "distance in the x direction: " + str(ballX)
        mylabel = label(pos=(50,60,0),text="seconds in flight: " + str(seconds),height=10)