《三体3》读后感

人类一思考,上帝就发笑。

但是上帝给予了人类大脑,人类就得思考。不光思考自己应该如何思考,而且还要思考上帝怎么思考。小到不能再小的微观粒子,大到不能再大的浩瀚宇宙,无数人都渴望着了解每一样事物的来龙去脉。但对于无尽的宇宙来说,人类的存在至今也只是短暂一瞬,地球的存在也只是沧海一粒。对于我们来说,许许多多的问题也许永远都得不到答案,只有内心的猜测和幻想,一直留在大脑深处。

但我们有幻想的权利。仰望星空不去思索星星是什么,我们的远祖也无法走出非洲,人类也永远只能在树林里面调来调去。对于上帝如何思考,我们也要追根到底,找出一个合乎解释的推理。无数人这么尝试了,但是无法得到大家的共鸣。有一个深居内陆偏僻小镇的中年男人,也这么去尝试了。虽然我们都清楚,这只是千千万万幻想中的一个特例而已,但是,我们都会发自内心地感慨,宇宙原来可以是这么回事,地球原来可以是这个样子,人类文明原来可以这样理解!从此以后,我们可以堂堂正正地思考上帝的想法!一个小虫子只能看到一片树叶,对于地球人之外的文明,我们又何尝不是虫子?要想不做虫子,就得想虫子不敢想之事,抛开禁锢我们思维的一切束缚,去做一次想象力上的无比酣畅之旅吧!

这就是刘慈欣新作《三体3》给我的震撼。此作作为《地球往事》的最后一步,从情节上全面超越了《三体1》和《三体2——黑暗森林》。如同软件作品一样,科幻小说的BUG很难消除,此书也有BUG。但是瑕永远不能掩瑜,网络上铺天盖地的评论和讨论可以得以证实。此作有希望摆脱科幻书籍的小众化市场和小众化销量,和主流文学作品同上畅销书的排行榜。这是大刘的福气,还是中国科幻的福气,甚至也是我的福气。在极度网络化的时代,书评一搜就一大堆,我就不加以细评和剧透了,只在这里吼几下,发表自己的感慨。无论如何,这都是一本好书。对于不喜欢科幻的人来说,里面的3个嵌入式童话也值得一看。

三体运动是指三个物体之间的相互作用和运动,可以使用Python进行模拟。其中,可以使用scipy库中的odeint函数求解微分方程组,来模拟三体运动。 以下是一个简单的三体运动模拟代码示例: ```python import numpy as np from scipy.integrate import odeint # 定义微分方程组 def three_body_equations(w, t, G, m1, m2, m3): x1, y1, vx1, vy1, x2, y2, vx2, vy2, x3, y3, vx3, vy3 = w r12 = np.sqrt((x2 - x1)**2 + (y2 - y1)**2) r13 = np.sqrt((x3 - x1)**2 + (y3 - y1)**2) r23 = np.sqrt((x3 - x2)**2 + (y3 - y2)**2) dx1dt = vx1 dy1dt = vy1 dvx1dt = G * m2 * (x2 - x1) / r12**3 + G * m3 * (x3 - x1) / r13**3 dvy1dt = G * m2 * (y2 - y1) / r12**3 + G * m3 * (y3 - y1) / r13**3 dx2dt = vx2 dy2dt = vy2 dvx2dt = G * m1 * (x1 - x2) / r12**3 + G * m3 * (x3 - x2) / r23**3 dvy2dt = G * m1 * (y1 - y2) / r12**3 + G * m3 * (y3 - y2) / r23**3 dx3dt = vx3 dy3dt = vy3 dvx3dt = G * m1 * (x1 - x3) / r13**3 + G * m2 * (x2 - x3) / r23**3 dvy3dt = G * m1 * (y1 - y3) / r13**3 + G * m2 * (y2 - y3) / r23**3 return dx1dt, dy1dt, dvx1dt, dvy1dt, dx2dt, dy2dt, dvx2dt, dvy2dt, dx3dt, dy3dt, dvx3dt, dvy3dt # 定义初始状态和参数 w0 = [1, 0, 0, 6, -1, 0, 0, -6, 0, 0, 0, 0] t = np.linspace(0, 10, 1000) G = 1 m1 = 1 m2 = 1 m3 = 1 # 求解微分方程组 wsol = odeint(three_body_equations, w0, t, args=(G, m1, m2, m3)) # 绘制轨迹图 import matplotlib.pyplot as plt plt.plot(wsol[:,0], wsol[:,1], label='Body 1') plt.plot(wsol[:,4], wsol[:,5], label='Body 2') plt.plot(wsol[:,8], wsol[:,9], label='Body 3') plt.legend() plt.show() ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值