本帖最后由 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
http://activity.ntsec.gov.tw/activity/race-1/40/hight/%E9%AB%9816H/%E9%AB%9816%E5%9C%96%E5%A4%BE/%E9%AB%9816.jpg
本帖最后由 gohomeman1 于 2010-5-19 15:19 编辑
这里转一些相对粗略的模拟结论，因为他们模拟时，忽略了许多条件，但基本条件是足够的。文章前面的研究目的、软件、误差分析等从略，大家可以自己去看链接。
双星系统的分析结论2
三恒星系统是极不稳定，这是其计算中的图例
其中图1，我还看见过一个大些的图
楼主的配文有问题啊，看不到全部了
从本图起，作者开始探讨小行星系统。
最后是讨论与评语。虽然只是高一的学生，但可以看出，整体文章写得很有条理，分析有相当的科学性。
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
後面ｋ行： m, x, y, z, vx, vy, vz
k代表星體個數，step代表畫點的間隔時間（建議 0.01）。Tlimit代表強制停止的時間，如果你不想使用這個功能，通常可以設定為一個很大的數字（例如 1e100）。split代表在一個 step之中要進行幾次運算，越大就越精確（建議 1000）。程式碼（見報告本文）中的Tstep其實就是step/split。
後面ｋ行的資料分別代表那ｋ個星體的初始狀態。m是質量， x y z是位置的三個分量， vx vy vz是速度的三個分量。
当年需要计算20分钟的图，现在20s都不到了。
厉害！！看不懂！！！
JAVA虚拟程序实在太慢了，比VB都慢得多。我又测试了一开始的FLASH程序，我现在怀疑这个模拟是否正确。可以看到，行星本身的初始条件就不好，而后面既使是0.03太阳质量的伴星也会造成紊乱，而0.3的虽然不会导致行星在两个恒星间跳跃，但我们可以注意到，这个行星其实早就撞毁在A星上了（此程序没考虑碰撞）。
永不重复的轨迹，当然这种情况下，行星也是早就撞毁了。
学习中，很复杂了，好多东西都忘了
太科普::42::
不错 受教了
这个很给力，应该转一下刘慈欣《三体》爱好者们常去的地儿。。。
回复 38# Eldghost
转过去不会被人喷吗？还以为我是拆台的呢。这些演示的本质就是：《三体》中的设定根本无法维持1000万年的。
刘慈欣小说里的bug并非唯一，比如超新星纪元里星等的问题。。。