Python 调用Vpython编写三体运动模型（不含碰撞）

最新推荐文章于 2023-03-23 17:22:29 发布

寒螭镜子

最新推荐文章于 2023-03-23 17:22:29 发布

阅读量1.8k

点赞数 1

文章标签： python 力学

本文链接：https://blog.csdn.net/weixin_42387066/article/details/116771443

版权

其中的一些参数可能需要找物理大佬调整（我也好想多看到点好看的混沌运动形态），以下，欢迎指出不足。（特别是3维碰撞那块，真的不太会，管理生路过）

import time
import numpy as np
from vpython import *

scene.height=880
scene.width=1890
m1=random()*1000
m2=random()*1000
m3=random()*1000
b1 = sphere(pos=vector(random(),random(),random()),radius = (3/40/np.pi*m1)**(1/3),color=color.red,make_trail=True,trail_radius=1)  #定义球的大小和颜色
b2 = sphere(pos=vector(random(),random(),random()), radius = (3/40/np.pi*m2)**(1/3),color=color.blue,make_trail=True,trail_radius=1)  #定义球的大小和颜色
b3 = sphere(pos=vector(random(),random(),random()), radius = (3/40/np.pi*m3)**(1/3),color=color.green,make_trail=True,trail_radius=1)  #定义球的大小和颜色
b1v=[random(),random(),random()]
b2v=[random(