From c91281669e451d03d02290ecf39fb1d44a06abc0 Mon Sep 17 00:00:00 2001 From: march3 Date: Sun, 2 Apr 2023 17:39:05 +0800 Subject: [PATCH] =?UTF-8?q?Python=E8=B6=85=E4=BA=BA-=E5=AE=87=E5=AE=99?= =?UTF-8?q?=E6=A8=A1=E6=8B=9F=E5=99=A8?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- bodies/asteroid.py | 2 -- bodies/earth.py | 2 -- bodies/jupiter.py | 2 -- bodies/mars.py | 2 -- bodies/mercury.py | 2 -- bodies/moon.py | 2 -- bodies/neptune.py | 2 -- bodies/pluto.py | 2 -- bodies/saturn.py | 2 -- bodies/uranus.py | 2 -- bodies/venus.py | 2 -- sim_scenes/fiction/fixed_stars.py | 1 - sim_scenes/interest/text_bodies.py | 3 +-- sim_scenes/interest/the_eye_of_god.py | 3 +-- sim_scenes/science/gravity_slingshot_1.py | 2 -- sim_scenes/science/gravity_slingshot_2.py | 1 - sim_scenes/solar_system/earth_moon.py | 4 ---- sim_scenes/solar_system/solar_system_3.py | 16 ++++++++-------- sim_scenes/solar_system/sun_earth.py | 2 +- sim_scenes/solar_system/sun_earth_moon.py | 2 +- sim_scenes/tri_bodies/tri_bodies_perfect.py | 9 ++++----- 21 files changed, 16 insertions(+), 49 deletions(-) diff --git a/bodies/asteroid.py b/bodies/asteroid.py index da81d3a..a8fb36a 100644 --- a/bodies/asteroid.py +++ b/bodies/asteroid.py @@ -27,8 +27,6 @@ class Asteroid(Body): """ def __init__(self, name="小行星", mass=4.1e10, - # init_position=[1.6 * AU, 0, 0], - # init_velocity=[0, 25.37, 0], init_position=[0, 0, 1.6 * AU], init_velocity=[-25.37, 0, 0], texture="", size_scale=1.0, distance_scale=1.0): diff --git a/bodies/earth.py b/bodies/earth.py index 01fc076..eaf0b43 100644 --- a/bodies/earth.py +++ b/bodies/earth.py @@ -23,8 +23,6 @@ class Earth(Body): 　平均密度: 5507.85 kg/m³ """ def __init__(self, name="地球", mass=5.97237e24,rotate_angle=0, - # init_position=[1.12 * AU, 0, 0], - # init_velocity=[0, 29.79, 0], init_position=[0, 0, 1.12 * AU], init_velocity=[-29.79, 0, 0], texture="earth1.jpg", size_scale=1.0, distance_scale=1.0, diff --git a/bodies/jupiter.py b/bodies/jupiter.py index 5ce8fdb..e2ff6e8 100644 --- a/bodies/jupiter.py +++ b/bodies/jupiter.py @@ -23,8 +23,6 @@ class Jupiter(Body): 　平均密度: 1.326 g/cm³ -> -> 1.326✕10³ kg/m³ """ def __init__(self, name="木星", mass=1.8982e27, - # init_position=[5.2 * AU, 0, 0], - # init_velocity=[0, 13.06, 0], init_position=[0, 0, 5.2 * AU], init_velocity=[-13.06, 0, 0], texture="jupiter1.jpg", size_scale=1.0, distance_scale=1.0, diff --git a/bodies/mars.py b/bodies/mars.py index 9a541dd..abd1970 100644 --- a/bodies/mars.py +++ b/bodies/mars.py @@ -23,8 +23,6 @@ class Mars(Body): 　平均密度: 3.9335 g/cm³ -> 3.9335✕10³ kg/m³ """ def __init__(self, name="火星", mass=6.4171e23, - # init_position=[1.5 * AU, 0, 0], - # init_velocity=[0, 24.13, 0], init_position=[0, 0, 1.5 * AU], init_velocity=[-24.13, 0, 0], texture="mars.jpg", size_scale=1.0, distance_scale=1.0, diff --git a/bodies/mercury.py b/bodies/mercury.py index fcf328e..a463473 100644 --- a/bodies/mercury.py +++ b/bodies/mercury.py @@ -24,8 +24,6 @@ class Mercury(Body): """ def __init__(self, name="水星", mass=3.3011e23, - # init_position=[0.4 * AU, 0, 0], - # init_velocity=[0, 47.87, 0], init_position=[0, 0, 0.4 * AU], init_velocity=[-47.87, 0, 0], texture="mercury.jpg", size_scale=1.0, distance_scale=1.0, diff --git a/bodies/moon.py b/bodies/moon.py index 6d05d92..1a8e2f4 100644 --- a/bodies/moon.py +++ b/bodies/moon.py @@ -23,8 +23,6 @@ class Moon(Body): """ def __init__(self, name="月球", mass=7.342e22, - # init_position=[363104 + 1.12 * AU, 0, 0], - # init_velocity=[0, 29.79 + 1.023, 0], init_position=[0, 0, 363104 + 1.12 * AU], init_velocity=[-(29.79 + 1.023), 0, 0], texture="moon.jpg", size_scale=1.0, distance_scale=1.0, diff --git a/bodies/neptune.py b/bodies/neptune.py index 354ffbe..13879e7 100644 --- a/bodies/neptune.py +++ b/bodies/neptune.py @@ -24,8 +24,6 @@ class Neptune(Body): """ def __init__(self, name="海王星", mass=1.0241e26, - # init_position=[30 * AU, 0, 0], - # init_velocity=[0, 5.43, 0], init_position=[0, 0, 30 * AU], init_velocity=[-5.43, 0, 0], texture="neptune.png", size_scale=1.0, distance_scale=1.0, diff --git a/bodies/pluto.py b/bodies/pluto.py index 794096c..734171a 100644 --- a/bodies/pluto.py +++ b/bodies/pluto.py @@ -24,8 +24,6 @@ class Pluto(Body): """ def __init__(self, name="冥王星", mass=1.303e22, - # init_position=[40 * AU, 0, 0], - # init_velocity=[0, 4.7, 0], init_position=[0, 0, 40 * AU], init_velocity=[-4.7, 0, 0], texture="pluto.png", size_scale=1.0, distance_scale=1.0, diff --git a/bodies/saturn.py b/bodies/saturn.py index 973ced2..457836e 100644 --- a/bodies/saturn.py +++ b/bodies/saturn.py @@ -24,8 +24,6 @@ class Saturn(Body): """ def __init__(self, name="土星", mass=5.6834e26, - # init_position=[10 * AU, 0, 0], - # init_velocity=[0, 9.64, 0], init_position=[0, 0, 10 * AU], init_velocity=[-9.64, 0, 0], texture="saturn.jpg", size_scale=1.0, distance_scale=1.0, diff --git a/bodies/uranus.py b/bodies/uranus.py index 8122a8f..f443d3e 100644 --- a/bodies/uranus.py +++ b/bodies/uranus.py @@ -24,8 +24,6 @@ class Uranus(Body): """ def __init__(self, name="天王星", mass=8.681e25, - # init_position=[19 * AU, 0, 0], - # init_velocity=[0, 6.81, 0], init_position=[0, 0, 19 * AU], init_velocity=[-6.81, 0, 0], texture="uranus.png", size_scale=1.0, distance_scale=1.0, diff --git a/bodies/venus.py b/bodies/venus.py index 1be74ac..234feb4 100644 --- a/bodies/venus.py +++ b/bodies/venus.py @@ -24,8 +24,6 @@ class Venus(Body): """ def __init__(self, name="金星", mass=4.8675e24, - # init_position=[0.72 * AU, 0, 0], - # init_velocity=[0, 35, 0], init_position=[0, 0, 0.72 * AU], init_velocity=[-35, 0, 0], texture="venus.jpg", size_scale=1.0, distance_scale=1.0, diff --git a/sim_scenes/fiction/fixed_stars.py b/sim_scenes/fiction/fixed_stars.py index 59afafc..07c4319 100644 --- a/sim_scenes/fiction/fixed_stars.py +++ b/sim_scenes/fiction/fixed_stars.py @@ -48,7 +48,6 @@ if __name__ == '__main__': d = pow((body.raduis + bodies[idx - 1].raduis) * SIZE_SCALE, 1.0) * 1.1 # d = (body.diameter + bodies[idx - 1].diameter) * SIZE_SCALE * 1.1 + D body.init_velocity = [0, 0, 0] - # body.init_position = [body.raduis * SIZE_SCALE, (distance_sum + d), AU] body.init_position = [-(distance_sum + d), AU, body.raduis * SIZE_SCALE] distance_sum += d diff --git a/sim_scenes/interest/text_bodies.py b/sim_scenes/interest/text_bodies.py index 7779550..17f1b36 100644 --- a/sim_scenes/interest/text_bodies.py +++ b/sim_scenes/interest/text_bodies.py @@ -26,8 +26,7 @@ def show_text_bodies(): params={"camera_pos": camera_pos}) # 放一个恒星作为背景 bg = FixedStar(name="bg", texture="fixed_star.png", mass=2e32, color=(0xff, 0xf8, 0xd4), - # init_position=[6000 * D, 450 * D, 100 * D], # [ 远+近- , 左+右- , 上+下-] - init_position=[-450 * D, 100 * D, 6000 * D], # [ 远+近- , 左+右- , 上+下-] + init_position=[-450 * D, 100 * D, 6000 * D], # [ 左-右+, 上+下-, 远+近- ] ignore_mass=True) bodies.append(bg) diff --git a/sim_scenes/interest/the_eye_of_god.py b/sim_scenes/interest/the_eye_of_god.py index 9cebee8..f41ebbc 100644 --- a/sim_scenes/interest/the_eye_of_god.py +++ b/sim_scenes/interest/the_eye_of_god.py @@ -25,8 +25,7 @@ def show_eye_of_god(): params={"camera_pos": camera_pos}) # 放一个恒星作为背景 bg = FixedStar(name="bg", texture="fixed_star.png", mass=2e32, color=(0xff, 0xf8, 0xd4), - # init_position=[6000 * D, 400 * D, 100 * D], # [ 远+近- , 左+右- , 上+下-] - init_position=[-400 * D, 100 * D, 6000 * D], # [ 远+近- , 左+右- , 上+下-] + init_position=[-400 * D, 100 * D, 6000 * D], # [ 左-右+, 上+下-, 远+近- ] ignore_mass=True) bodies.append(bg) diff --git a/sim_scenes/science/gravity_slingshot_1.py b/sim_scenes/science/gravity_slingshot_1.py index 9a17fcf..304380c 100644 --- a/sim_scenes/science/gravity_slingshot_1.py +++ b/sim_scenes/science/gravity_slingshot_1.py @@ -18,11 +18,9 @@ if __name__ == '__main__': bodies = [ Sun(size_scale=2e1), # 太阳放大 20 倍 Earth(size_scale=1e3, # 地球放大 1000 倍 - # init_position=[0, -3 * AU, 0], # 地球距离太阳 3 个天文单位 init_position=[3 * AU, 0, 0], # 地球距离太阳 3 个天文单位 # TODO: 尝试调整朝向太阳的速度，取值 33、38、50 或者其他 # init_velocity=[0, 33, -1], - # init_velocity=[0, 38, -1], # 朝向太阳的速度为 38km/s，-1 km/s 是为了防止地球正面对着太阳冲去 init_velocity=[-38, -1, 0], # 朝向太阳的速度为 38km/s，-1 km/s 是为了防止地球正面对着太阳冲去 # init_velocity=[0, 50, -1], ), diff --git a/sim_scenes/science/gravity_slingshot_2.py b/sim_scenes/science/gravity_slingshot_2.py index 8d2f8d6..f089e0b 100644 --- a/sim_scenes/science/gravity_slingshot_2.py +++ b/sim_scenes/science/gravity_slingshot_2.py @@ -21,7 +21,6 @@ if __name__ == '__main__': Earth(size_scale=3e2, # 地球放大 300 倍 init_position=[0, 0, 0], # # init_velocity=[0, 33, -1], - # init_velocity=[-1, 10, 0], # 朝向木星的速度为 38km/s，-1 km/s 是为了防止地球正面对着木星冲去 init_velocity=[-10, 0, -1], # 朝向木星的速度为 38km/s，-1 km/s 是为了防止地球正面对着木星冲去 # init_velocity=[0, 50, -1], ), diff --git a/sim_scenes/solar_system/earth_moon.py b/sim_scenes/solar_system/earth_moon.py index 64f2457..c8b81d8 100644 --- a/sim_scenes/solar_system/earth_moon.py +++ b/sim_scenes/solar_system/earth_moon.py @@ -14,10 +14,6 @@ from bodies.body import AU if __name__ == '__main__': """ 地球、月球 -init_position=[x, y, z], -init_position=[-y, z, x], -init_velocity=[x, y, z], -init_velocity=[-y, z, x], """ # 地球的Y方向初始速度 EARTH_INIT_VELOCITY = 0 diff --git a/sim_scenes/solar_system/solar_system_3.py b/sim_scenes/solar_system/solar_system_3.py index e99fdd4..cbb5463 100644 --- a/sim_scenes/solar_system/solar_system_3.py +++ b/sim_scenes/solar_system/solar_system_3.py @@ -27,20 +27,20 @@ if __name__ == '__main__': # region 1.构建4个小行星 ------------- # asteroids = [ # Asteroid(size_scale=3e9, # 小行星放大 1000000000 倍，距离保持不变 - # init_position=[1.6 * AU, 0, 0], - # init_velocity=[0, 25.37, 0], + # init_position=[0, 0, 1.6 * AU], + # init_velocity=[-25.37, 0, 0], # distance_scale=1), # Asteroid(size_scale=3e9, # 小行星放大 1000000000 倍，距离保持不变 - # init_position=[-1.6 * AU, 0, 0], - # init_velocity=[0, -25.37, 0], + # init_position=[0, 0, -1.6 * AU], + # init_velocity=[25.37, 0, 0], # distance_scale=1), # Asteroid(size_scale=3e9, # 小行星放大 1000000000 倍，距离保持不变 - # init_position=[0, 1.6 * AU, 0], - # init_velocity=[-25.37, 0, 0], + # init_position=[-1.6 * AU, 0, 0], + # init_velocity=[0, 0, -25.37], # distance_scale=1), # Asteroid(size_scale=3e9, # 小行星放大 1000000000 倍，距离保持不变 - # init_position=[0, -1.6 * AU, 0], - # init_velocity=[25.37, 0, 0], + # init_position=[1.6 * AU, 0, 0], + # init_velocity=[0, 0, 25.37], # distance_scale=1), # ] # endregion 1 -------------------------- diff --git a/sim_scenes/solar_system/sun_earth.py b/sim_scenes/solar_system/sun_earth.py index d78263f..96638c4 100644 --- a/sim_scenes/solar_system/sun_earth.py +++ b/sim_scenes/solar_system/sun_earth.py @@ -12,7 +12,7 @@ from sim_scenes.func import mayavi_run, ursina_run if __name__ == '__main__': """ - 太阳、地球 Sirius(size_scale=0.2, init_position=[0, 0, 0]), + 太阳、地球 """ bodies = [ Sun(size_scale=5e1), # 太阳放大 50 倍 diff --git a/sim_scenes/solar_system/sun_earth_moon.py b/sim_scenes/solar_system/sun_earth_moon.py index e5e6141..3395ccd 100644 --- a/sim_scenes/solar_system/sun_earth_moon.py +++ b/sim_scenes/solar_system/sun_earth_moon.py @@ -18,7 +18,7 @@ if __name__ == '__main__': # 地球的Y方向初始速度 EARTH_INIT_VELOCITY = 0 # 0km/s sun = Sun(init_position=[0, 0, AU], size_scale=2e1) # 太阳放大 20 倍 - # 忽略质量的引力 [x, y, z]->[-y, z, x] + # 忽略质量的引力 sun.ignore_mass = True # 观看月相变化的过程：分别是 新月、蛾眉月、上弦月、盈凸、满月、亏凸、下弦月、残月 diff --git a/sim_scenes/tri_bodies/tri_bodies_perfect.py b/sim_scenes/tri_bodies/tri_bodies_perfect.py index 68c417b..b92c829 100644 --- a/sim_scenes/tri_bodies/tri_bodies_perfect.py +++ b/sim_scenes/tri_bodies/tri_bodies_perfect.py @@ -23,11 +23,10 @@ if __name__ == '__main__': mass = 2e30 r = 2 * AU # p = 12 # TODO: 三体转圆形花 - p = 14.88 # TODO: 三体转圈近似圆形 + p = 14.95 # TODO: 三体转圈近似圆形 # p = 16 # TODO: 三体转圆形花 # p = 18 # TODO: 三体转圆形花 # p = 19 # TODO: 三体转圆形花 - # [x, y, z]->[-y, z, x] bodies = [ Sun(name="红轨太阳A", mass=mass, init_position=[0, math.sqrt(3) * r, 0], @@ -44,9 +43,9 @@ if __name__ == '__main__': init_velocity=[1 / 2 * p, math.sqrt(3) / 2 * p, 0], trail_color=(0, 0, 255), size_scale=5e1, texture="sun2.jpg"), # 太阳放大 100 倍 - # Earth(init_position=[0, -349597870.700, 0], - # init_velocity=[15.50, 0, 0], - # size_scale=4e3, distance_scale=1), # 地球放大 4000 倍，距离保持不变 + Earth(init_position=[0, 2 * r / math.sqrt(12), r * 3], + init_velocity=[0, 0, 0], + size_scale=4e3, distance_scale=1), # 地球放大 4000 倍，距离保持不变 ] # 使用 mayavi 查看的运行效果 # mayavi_run(bodies, SECONDS_PER_WEEK, view_azimuth=0) -- GitLab