牧夫天文网

找回密码
加入牧夫(请注明天文爱好者,否则无法通过审核)

QQ登录

只需一步,快速开始

发新帖
楼主: gohomeman1
查看: 12346| 回复: 40
发表于 2010-5-19 15:06 来自: 浙江省宁波市 鄞州畅联信息技术有限公司
本帖最后由 gohomeman1 于 2010-5-19 15:40 编辑

FTP中还有一段代码,是关于三体问题的数值模拟的。首先声明,那些代码不是我写的,所以不要问我关于代码的任何问题,当然也请勿用于商业用途。
关于小行星、木星、太阳间的关系,我找到一个网页,汗颜得很,人家是高一学生哦!作者姓名:陈威尹 谢俊驹
http://activity.ntsec.gov.tw/activity/race-1/40/hight/%E9%AB%9816H/h%E9%AB%9816.html


                               
登录/注册后可看大图

回复 显示全部楼层 道具 举报

发表于 2010-5-19 15:14 来自: 浙江省宁波市 鄞州畅联信息技术有限公司
本帖最后由 gohomeman1 于 2010-5-19 15:19 编辑

这里转一些相对粗略的模拟结论,因为他们模拟时,忽略了许多条件,但基本条件是足够的。文章前面的研究目的、软件、误差分析等从略,大家可以自己去看链接。
双星系统1.JPG

回复 显示全部楼层 道具 举报

发表于 2010-5-19 15:24 来自: 浙江省宁波市 鄞州畅联信息技术有限公司
双星系统的分析结论2
双星系统2.JPG

回复 显示全部楼层 道具 举报

发表于 2010-5-19 15:29 来自: 浙江省宁波市 鄞州畅联信息技术有限公司
三恒星系统是极不稳定,这是其计算中的图例
三恒星不能稳定.JPG

回复 显示全部楼层 道具 举报

发表于 2010-5-19 15:30 来自: 浙江省宁波市 鄞州畅联信息技术有限公司
其中图1,我还看见过一个大些的图
混沌三体.jpg

回复 显示全部楼层 道具 举报

发表于 2010-5-19 15:38 来自: 江西省上饶市 电信
楼主的配文有问题啊,看不到全部了

回复 显示全部楼层 道具 举报

发表于 2010-5-19 15:39 来自: 浙江省宁波市 鄞州畅联信息技术有限公司
从本图起,作者开始探讨小行星系统。
稳定的小行星1.JPG

回复 显示全部楼层 道具 举报

发表于 2010-5-19 15:45 来自: 浙江省宁波市 鄞州畅联信息技术有限公司
稳定的小行星2.JPG
稳定的小行星3.JPG

回复 显示全部楼层 道具 举报

发表于 2010-5-19 15:50 来自: 浙江省宁波市 鄞州畅联信息技术有限公司
稳定的小行星4.JPG
研究结论.JPG

回复 显示全部楼层 道具 举报

发表于 2010-5-19 15:53 来自: 浙江省宁波市 鄞州畅联信息技术有限公司
最后是讨论与评语。虽然只是高一的学生,但可以看出,整体文章写得很有条理,分析有相当的科学性。
讨论与评语.JPG

回复 显示全部楼层 道具 举报

发表于 2010-5-19 16:41 来自: 浙江省宁波市 鄞州畅联信息技术有限公司
2000年的时候,那位陈同学读高一。他写的程序现在再来运行,我都要采用兼容模式+下载DLL文件了。这是他当年找出了一个可以运行1.3亿步的三恒星系统,初始条件如下:
3, 0.01 , 1e10 ,  1000
1,-0.1,-0.2,0,0.2,-0.07,0,     0.1,1,0,0,0,0.7,0,       0.2,0,1,0,-1,0,0
第一行: k, step, Tlimit, split
後面k行: m, x, y, z, vx, vy, vz

k代表星體個數,step代表畫點的間隔時間(建議 0.01)。Tlimit代表強制停止的時間,如果你不想使用這個功能,通常可以設定為一個很大的數字(例如 1e100)。split代表在一個 step之中要進行幾次運算,越大就越精確(建議 1000)。程式碼(見報告本文)中的Tstep其實就是step/split。
後面k行的資料分別代表那k個星體的初始狀態。m是質量, x y z是位置的三個分量, vx vy vz是速度的三個分量。


当年需要计算20分钟的图,现在20s都不到了。

混沌三体2.jpg

回复 显示全部楼层 道具 举报

发表于 2010-5-19 20:28 来自: 辽宁省大连市 大连理工大学西山生活区20舍
厉害!!看不懂!!!

回复 显示全部楼层 道具 举报

发表于 2010-5-19 21:44 来自: 浙江省宁波市 鄞州畅联信息技术有限公司
JAVA虚拟程序实在太慢了,比VB都慢得多。我又测试了一开始的FLASH程序,我现在怀疑这个模拟是否正确。可以看到,行星本身的初始条件就不好,而后面既使是0.03太阳质量的伴星也会造成紊乱,而0.3的虽然不会导致行星在两个恒星间跳跃,但我们可以注意到,这个行星其实早就撞毁在A星上了(此程序没考虑碰撞)。
0.3 1情况.JPG
0.03 1也不稳定?.JPG

回复 显示全部楼层 道具 举报

发表于 2010-5-19 22:00 来自: 浙江省宁波市 鄞州畅联信息技术有限公司
永不重复的轨迹,当然这种情况下,行星也是早就撞毁了。
2 1演化1.JPG
2 1演化2.JPG

回复 显示全部楼层 道具 举报

发表于 2010-10-22 22:00 来自: 天津市 电信
学习中,很复杂了,好多东西都忘了

回复 显示全部楼层 道具 举报

发表于 2010-10-24 16:48 来自: 福建省福州市 电信
太科普  

回复 显示全部楼层 道具 举报

发表于 2010-11-29 10:26 来自: 云南省楚雄州 联通
不错 受教了

回复 显示全部楼层 道具 举报

发表于 2010-11-30 21:35 来自: 辽宁省沈阳市 沈阳广电电信合作宽带
这个很给力,应该转一下刘慈欣《三体》爱好者们常去的地儿。。。

回复 显示全部楼层 道具 举报

发表于 2010-11-30 21:43 来自: 浙江省宁波市 鄞州畅联信息技术有限公司
回复 38# Eldghost

转过去不会被人喷吗?还以为我是拆台的呢。这些演示的本质就是:《三体》中的设定根本无法维持1000万年的。

回复 显示全部楼层 道具 举报

发表于 2010-11-30 23:07 来自: 辽宁省沈阳市 沈阳广电电信合作宽带
刘慈欣小说里的bug并非唯一,比如超新星纪元里星等的问题。。。

回复 显示全部楼层 道具 举报

APP下載|爱牧夫天文淘宝店|牧夫天文网 ( 公安备案号21021102000967 )|网站地图|辽ICP备19018387号

GMT+8, 2024-7-12 19:58 , Processed in 0.031003 second(s), 8 queries , Gzip On, Redis On.

Powered by Discuz! X3.4

© 2001-2023 Discuz! Team.

快速回复 返回顶部 返回列表