所以我正在编写一个斯诺克游戏，我决定弄清楚如何让球相互碰撞的最好方法是在一个单独的程序中这样做，然后将它复制进去。我是一个非常称职的数学家，所以我'已经坐下来，绘制了事件并计算了实际发生的事情。
我的方法是将每个球的初始速度分解为 xr 和 yr 分量，在参考系中 xr 分量与通过每个球中心的矢量一致，并且 yr 分量垂直于这个。然后我对球的 xr 分量进行简单切换，并保持 yr 分量不变，然后计算速度的 x 和 y 分量回到标准参考系中的值。
出于某种原因，无论是通过数学还是编程错误，我似乎都无法让它工作。以下是我到目前为止所拥有的内容，并且我已经查看了互联网上我能找到的几乎所有相关页面以及在该网站上提出的所有类似问题。我也不是一个精通程序员。
from visual import *
dt = 0.01
r = 5
red = sphere(pos=(-25,25,0),radius = r,color=color.red)
green = sphere(pos=(25,-25,0),radius = r,color=color.green)
red.velocity = vector(10,-10,0)
green.velocity = vector(-10,10,0)
def posupdate(ball):
ball.pos = ball.pos + ball.velocity*dt
def ballhit(ball1,ball2):
v1 = ball1.velocity
v1x = ball1.velocity.x
v1y = ball1.velocity.y
v2 = ball2.velocity
v2x = ball2.velocity.x
v2y = ball2.velocity.y
xaxis = vector(1,0,0)
btb = ball2.pos - ball1.pos
nbtb = btb/abs(btb)
if abs(btb) < 2*r:
phi = acos(dot(nbtb,xaxis)/abs(nbtb)*abs(xaxis))
ang1 = acos(dot(v1,xaxis)/abs(v1)*abs(xaxis))
ang2 = acos(dot(v2,xaxis)/abs(v2)*abs(xaxis))
v1xr = abs(v1)*cos((ang1-phi))
v1yr = abs(v1)*sin((ang1-phi))
v2xr = abs(v2)*cos((ang2-phi))
v2yr = abs(v2)*sin((ang2-phi))
v1fxr = v2xr
v2fxr = v1xr
v1fyr = v1yr
v2fyr = v2yr
v1fx = cos(phi)*v1fxr+cos(phi+pi/2)*v1fyr
v1fy = sin(phi)*v1fxr+sin(phi+pi/2)*v1fyr
v2fx = cos(phi)*v2fxr+cos(phi+pi/2)*v2fyr
v2fy = sin(phi)*v2fxr+sin(phi+pi/2)*v2fyr
ball1.velocity.x = v1fx
ball1.velocity.y = v1fy
ball2.velocity.x = v2fx
ball2.velocity.y = v2fy
return ball1.velocity, ball2.velocity
while 1==1:
rate(100)
posupdate(red)
posupdate(green)
ballhit(red,green)
提前感谢您提供的任何帮助。
编辑：碰撞检测没有问题，只是碰撞后球的速度矢量的计算。抱歉，我应该说得更清楚。