2024第十三届认证杯数学中国数学建模国际赛（小美赛） - 数学中国数学建模国际赛(CAMCM) A题 木星：保护者还是威胁者？ 解题思路和代码

目录

任务一：小行星带威胁分析

代码实现

1. 初始化参数

2. 定义动力学方程

3. 轨道模拟

4. 分析威胁和保护事件

5. 综合模拟与结果统计

解题思路

任务二：假设威胁仅源于奥尔特云

思路

代码

任务三：木星质量和轨道变化对“晚期重轰炸效应”的影响

关键因素：

模型分析：

模型假设：

数值模拟步骤：

代码实现框架：

---

任务一：小行星带威胁分析

代码实现

以下代码使用Python结合NumPy和SciPy实现数值模拟，主要模拟小行星在木星引力作用下的轨道演化。

1. 初始化参数

import numpy as np
from scipy.integrate import solve_ivp

# 常量定义
G = 6.67430e-11  # 万有引力常数 (m^3/kg/s^2)
M_sun = 1.989e30  # 太阳质量 (kg)
M_jupiter = 1.898e27  # 木星质量 (kg)
M_earth = 5.972e24  # 地球质量 (kg)

# 天体轨道参数 (AU -> m)
AU = 1.496e11  # 1天文单位 (m)
R_jupiter = 5.2 * AU  # 木星轨道半径
R_earth = 1.0 * AU  # 地球轨道半径
R_asteroid_min = 2.1 * AU  # 小行星带内侧
R_asteroid_max = 3.3 * AU  # 小行星带外侧