专栏首页Crossin的编程教室用python模拟地球飞越火星

用python模拟地球飞越火星

本文涉及影片内容,但算不上剧透。想看还未看的同学请谨慎阅读。

春节档的几部电影里,《流浪地球》成为一匹黑马,刷爆微博和朋友圈,从最初不高的排片率,一跃上升到票房榜首。

看过影片的人可能会有一个疑问,太阳系那么空旷,为什么地球在“流浪”的过程中,非要作死往木星上靠?

这就牵涉到一个时常在影视作品中拿来做文章,实际在太空探测中很常见的一个概念:

引力弹弓效应

当一个飞行器飞过一个行星时,如果距离较近,就会受到行星的引力影响:靠近时会逐渐加速,远离时会逐渐减速。由于能量守恒,飞行器在进入和离开行星轨道时,与行星的相对速度大小没有变化,但方向却发生了变化。再加上行星本身的速度,从旁观者角度来看,飞行器的速度就发生了变化。举个更通俗的例子:你把一个网球以速度 v 砸向一面墙,球会以同样的速度 v 反弹回来;而如果你以速度 v 砸向一个以速度 u 迎面开来卡车头上,网球则会以 v + 2u 的速度反弹回来(相对速度 v + u 加上卡车速度 u)

利用这个办法,可以在不额外消耗燃料的情况下,使太空探测器得到加速或者减速,取决于与行星接触的角度)。这看起来似乎不符合能量守恒定律和动量守恒定律,实则是因为行星的质量相较飞行器来说太大了,其速度的变化可以忽略不计。

目前离地球最远人造物体:旅行者1号,当年就连续利用了木星土星进行了加速。

而在影片中,也就是为了获取更大的速度以脱离太阳系,人类才选择借助木星的引力弹弓效应

出于好奇,我用 python 写了一个简化的仿真,模拟地球飞越木星时可能发生的情况:

利用引力弹弓加速

速度太快或距离太远,加速效果不明显

速度太慢或距离太近,坠入木星

利用引力弹弓减速

视频版

通过几次测试,我发现这“弹弓”可不是那么好用的,离得远了没啥效果,离得太近又可能被引力捕获,撞在木星上(事实在撞上之前,地球就会因达到洛希极限而被撕裂)。

至于剧中发生了什么,我这里也就不多说了。

我这个模拟代码中,为了视觉效果,地球和木星的大小,以及他们的相对距离是不成比例的,否则会更小。而且为了让速度变化更明显,对运行速度也做了放大。所以从绝对数值上来说是不准确的,但相对变化都是通过基本物理定律计算得出的

核心部分代码:

# 地木坐标差
delta_x = (jupiter[0] - earth[0]) * k
delta_y = (jupiter[1] - earth[1]) * k
# 地木距离平方
r2 = delta_x ** 2 + delta_y ** 2
# 地木间引力,万有引力定律
F = G * m * M / r2
# 地木夹角
theta = math.acos(delta_x / r2 ** 0.5)
# x、y 轴引力分量
fx = abs(F * math.cos(theta)) * sign(delta_x)
fy = abs(F * math.sin(theta)) * sign(delta_y)
# x、y 轴加速度,牛顿第二定律 F = ma
ax = fx / m
ay = fy / m
# 速度变化,vt = v0 + at
vel_x += ax * t
vel_y += ay * t
# 位移变化,st = s0 + vt
pos_x += vel_x * t / k
pos_y += vel_y * t / k

相关参数:

k = 1e7            # 距离缩放参数
m = 5.9742e24      # 地球质量
M = 1898.7e27      # 木星质量
G = 6.67259e-17    # 万有引力常量
t = 1e5            # 时间缩放参数
pos_x= 0           # 地球坐标
pos_y= 550earth = pos_x, pos_y
vel_x= 300         # 地球速度
vel_y= 0jupiter = 700, 150 # 木星坐标
v_j = 3            # 木星速度

想要测试不同效果,可以调整坐标、速度等参数。

“引力弹弓”的概念在几年前的《星际穿越》和《火星救援》中都有提及,我一直想写个模拟,直到这次趁《流浪地球》热映的机会才动手实现。这部电影我从去年7月尚未跳票之时比较期待,也看过原著,本来还担心会不会排片太少,没想到会火成这样。这里面当然有很多偶然的商业因素,影片本身也有很多不足,但我作为一个“科学爱好者”,乐见如此。之前流产的《三体》电影版或许也因此有机会重新投胎。

希望以后会有更多精彩的科幻作品出现,也希望更多的人在低头忙碌之余能抬头仰望星空。

本文分享自微信公众号 - Crossin的编程教室(crossincode)

原文出处及转载信息见文内详细说明,如有侵权,请联系 yunjia_community@tencent.com 删除。

原始发表时间:2019-02-13

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

我来说两句

0 条评论
登录 后参与评论

相关文章

  • 10个对Web开发者最有用的Python包

    Python最近成为了开发人员最喜欢的语言之一。无论你是专业的,业余的,还是一个初学者,你都可以从Python语言及其程序包中受益。Python已经被证明是当今...

    Crossin先生
  • 【Python 第0课】Why Python?

    为什么用Python作为编程入门语言? 原因很简单。 嗯。。。原因就是,很简单。。。 每种语言都会有它的支持者和反对者。去Google一下“why python...

    Crossin先生
  • 【Python 第29课】 连接list

    今天要说的方法是join。它和昨天说的split正好相反:split是把一个字符串分割成很多字符串组成的list,而join则是把一个list中的所有字符串连接...

    Crossin先生
  • 解决lenove e450c安装centos7无线网卡(RTL8723BE)无法使用

    tanmx
  • TW洞见 | 我和敏捷团队的五个约定

    我——作为一名测试人员——有一个与众不同的习惯:每当要加入一个新项目的时候,我总会找到项目中的同伴,真诚而亲切地说:“为了更好地合作,我有5个约定,希望大家能尽...

    ThoughtWorks
  • TW洞见 | 我和敏捷团队的五个约定

    今日洞见 文章作者来自ThoughtWorks: 覃其慧。 本文转载自InfoQ:http://www.infoq.com/cn/articles/though...

    ThoughtWorks
  • 软件测试方法课程笔记(1)

    举某些例子, 软件测试方法有黑盒测试, 白盒测试 按阶段来区分的话有单元测试, 集成测试, 系统测试 按目的来分有性能测试等

    Mezereon
  • 直播软件搭建,现阶段搭建后期注意事项

    很多人在直播软件搭建完成后就觉得万事大吉,只需要上架就可,其实不然,直播软件源码为什么具备二次开发功能,就是为了在后期不断对系统进行升级,对运营商来说,注重直播...

    yunbaokeji柯基
  • Python让你自己做一个软件,自己开个聊天室,厉害吧!

    内容 知识点 asyncore 、asynchat模块运用 环境 python 3.5 功能描绘 在本实验中,我们将实现一个简略的图形界面谈天系统。我们可以经...

    企鹅号小编
  • Spring Boot使用CORS解决跨域问题

    Web开发经常会遇到跨域问题,解决方案有:jsonp,iframe,CORS等等。 CORS 与 JSONP 相比: 1、 JSONP 只能实现 GET 请...

    itlemon

扫码关注云+社区

领取腾讯云代金券