用 Python 代码来模拟用核弹射击小行星会，最终能否拯救地球？

你有没有想过，如果一颗小行星正冲向地球，而人类向它发射了一颗核武器，会发生什么？它会影响小行星的轨迹吗？我们会活下来吗？瓦砾会怎样？您可能听起来很熟悉这种情况。Armageddon 和 Deep Impact 等电影就是好莱坞演绎这一壮举的例子。有一些现实项目正在努力减轻和消除小行星的威胁。美国国家航空航天局甚至有一个类似的任务称为 DART（双小行星重定向测试），该任务将利用航天器的动能使小行星偏离其地球终点路径。别担心，他们正在一颗没有威胁的小行星上进行这项研究，但这项研究将改善地球抵御宇宙威胁的防御系统。

Toutatis 是对地球生命构成潜在危害的最大小行星之一。一颗跨越惊人5.4公里、质量约50.5万亿公斤的小行星！这是一颗巨大的小行星，可能会终结我们所知的生命。这就是为什么 NASA 的 DART 任务如此重要。截至目前，地球几乎没有希望抵御如此巨大的宇宙敌人。在本文中，我们将模拟尝试防御 Toutatis 时可能发生的情况。我们将使用 Python 和一些基本的轨道力学（简要地）来模拟一个真实的结果（这是一个近似值，但是是一个不错的结果）。让我们开始编码吧！

导入包

我总是通过导入我需要的包来开始我的 Python 代码。此模拟所需的软件包包括：

NumPy：用于基本的数学函数（定义为np以便于调用）
random：用于随机化模拟以在每次运行时提供新结果
solve_ivp：用于对二体运动方程进行数值积分（我们会介绍这个，别担心）
pyplot：用于可视化结果（定义为plt以便于调用）

# Importing Pac