有人能写出这个二体系统的稳定一些的解吗

歌者文明清理员

代码如下

1. # -*- coding:utf-8 -*-
   
2. """
   
3. @Author: 歌者文明清理员1410382
   
4. @Datetime: 2023-03-04 16:26:16
   
5. @Name: 双星系统.py
   
6. @OriginLink:
   
7. """
   
8. from calc import *
   
9. 
   
10. a = obj(1, 0.5, vec(0.2, 0, 0), color.yellow, True, vec(0, 0, -10))
    
11. b = obj(1, 0.5, vec(-0.2, 0, 0), color.cyan, True, vec(0, 0, 10))
    
12. earth = obj(0.000000000000000000000000000000000001, 0.2, vec(0.2, 0.295, 0.01), color.blue, False, vec(0, 15, 0))
    
13. mainloop(1)
    

复制代码

calc.py内容如下

1. # -*- coding:utf-8 -*-
   
2. """
   
3. @Author: 歌者文明清理员1410382
   
4. @Datetime: 2023-03-04 16:21:03
   
5. @Name: calc.py
   
6. @OriginLink:
   
7. """
   
8. from vpython import *
   
9. 
   
10. 
    
11. def obj(mass, radius, momentum, color, emissive, pos, lock=False):
    
12.     a = sphere(color=color, make_trail=True, retain=150, trail_width=radius / 10)
    
13.     a.mass = mass
    
14.     a.radius = radius
    
15.     a.momentum = momentum
    
16.     a.trail_color = color
    
17.     a.emissive = emissive
    
18.     a.pos = pos
    
19.     a.locked = lock
    
20.     if a.emissive:
    
21.         attach_light(a)
    
22.     stars.append(a)
    
23.     return a
    
24. 
    
25. 
    
26. def move(g=1):
    
27.     for sa in stars:
    
28.         sa.pos += sa.momentum * (not sa.locked)
    
29.         for sb in stars:
    
30.             if sa == sb:
    
31.                 continue  # samestar
    
32.             distance = mag2(sa.pos-sb.pos)**1.5
    
33.             if distance < limit:
    
34.                 distance = limit
    
35.             force = g * sa.mass / \
    
36.                 distance * (sa.pos-sb.pos)
    
37.             sb.momentum += force
    
38. 
    
39. 
    
40. def mainloop(dt=1, g=1):
    
41.     scene.width, scene.height = 1900, 1050
    
42.     scene.range = 100
    
43. 
    
44.     while True:
    
45.         rate(dt * 100)
    
46.         move(g)
    
47. 
    
48. 
    
49. stars = []
    
50. limit = 0.001
    

复制代码

每次都是甩出去，而且两个恒星慢慢的总是会相撞弹开
我知道弹开是因为，它们撞在一起时，距离的平方突然上升导致速度上升
准备工作

1. pip install vpython -i https://mirrors.aliyun.com/pypi/simple

复制代码

歌者文明清理员

只需要改参数就行了，10-12行第三个参数（速度）最后一个参数（初始位置）
因为模拟中的二体系统中的行星轨迹极不规律，到处乱飞，有三体星系那么难以生存
改一点点参数（0.001）整个轨迹就完全不一样了

1. obj(质量, 半径, 速度, 颜色, 是否发光, 位置[, 位置是否锁定])

复制代码

wyhpylc

歌者文明清理员 发表于 2023-3-4 18:44
只需要改参数就行了，10-12行第三个参数（速度）最后一个参数（初始位置）
因为模拟中的二体系统中的行星 ...

这难道就是芝士的力量吗

歌者文明清理员

wyhpylc 发表于 2023-3-5 10:14
这难道就是芝士的力量吗

你说哪个，公式我是从 glowscript 的例子上复制的（已做出太阳系，正在完善，有轨道倾角，离心率也有）

wyhpylc

歌者文明清理员 发表于 2023-3-5 10:17
你说哪个，公式我是从 glowscript 的例子上复制的（已做出太阳系，正在完善，有轨道倾角，离心率也有）

az，粘贴的也厉害啊，我甚至看不懂这个公式]

歌者文明清理员

wyhpylc 发表于 2023-3-5 10:27
az，粘贴的也厉害啊，我甚至看不懂这个公式]

我也看不懂

wyhpylc

歌者文明清理员 发表于 2023-3-5 10:31
我也看不懂

az