matlab模拟三体运动_物理与科学界求解了300多年,找到三体运行轨道解意味着什么呢?...

女作者丨黄金

三体的问题因无法求解而著名,物理学乃至科学界在1687年牛顿发表其定律后发生了巨变,其中的运动与引力方程将原本看起来不稳定的宇宙转变为了可预测的完美模型,如果已知目前太阳系中各星体的位置的速度,就可以用牛顿方程来计算任意时间时它们的位置,其实现实中没有那么简单的,虽然牛顿方程如此精妙,它们只在一种情况下能求得行星运动简单解,那就是只有两个星体相互围绕而不受宇宙引力作用时,只要再加上一个星体,大部分情况下其运动都将变得完全混乱,不再存在简单解,这就是三体问题。

19450d8af7453a6654650b9379178e9d.png

物理与科学界求解了300多年,找到三体运行轨道解意味着什么呢?

牛顿运动定律和万有引力定律给出了一组微分方程,可以用牛顿的另一伟大发明—“微分法”,求得一个简单方程,把具体数值代入即可求解,这些值就是引力质点的起始位置与速度,再代入时间值,就能通过这个方程得到任意时间下该系统的状态,这种简单、有确定解的方程叫做“解析表达式”,这意味着它能由有限数学运算及函数表达。

两个引力体时,牛顿定律的解就是星体行进路径的方程。路径有可能是抛物线,也可能是行星轨道的圆或椭圆,或是星际彗星所走的双曲线,总的来说都是圆锥截面,即切圆锥时得到的形状,这些方程很简单,开普勒在牛顿定律闻名前70年就得到了大部分行星运动的椭圆解,而牛顿定律发表之后,科学家就开始寻找更复杂系统的简单“解析解”。

三个引力体的系统,就是科学家一直想要解答的目标,但仅仅多一个星体的额外影响就将使得确切解成为不可能。三体的问题成为许多伟大的数学家的痴迷,但之后的三个世纪,也只有一些特殊情况下的解被求得。

19世纪00年代末,数学家“布伦斯”和“庞加莱”严词断言不存在一般“解析解”。

其实三体问题中几乎所有初始配置的演化都由混沌动力学决定

1、打开已有的一组三体配置文件(.tbc)并运行(点击播放按钮)。 "File"菜单下有导入(Import)、导出(Export)功能,在不能上传附件时方便以纯文本方式交流自己搜索出来的三体配置! 2、手工设定初始条件的全部数值(点击魔术棍按钮)。分别指定三个物体的初始条件(X、Y、Z坐标,质量,初始速度的幅度、在XY平面上的角度0~360、在XZ平面上的角度0~360)。四个圆形选项(Radio Button)是参照系选择:默认的"Normalize to Centroid"是按三体系统质心作为参照系进行速度平衡,相当于观察者总是跟随三体的质心运动。另外三个选项分别是以第一、二、三个天体作为参照系,即总是把这个天体放在中心位置从不移动--注意这是非惯性参照系!(一般应选取行星主要围绕的那个恒星,方便观察行星轨道) 如果XZ平面上的初始速度角度都是0,则退化为二维的三体。 不过手工设定的条件通常都很难稳定运行。 3、设定搜索条件,让软件自动搜索。搜索分为两步: 3.1、搜索稳定的三体(点击望远镜按钮) 第一部分是每个物体的约束条件:坐标最大值、最小质量、最大质量、最小速度幅度、最大速度幅度。 第二部分是是否要求三体在最初N步里超出一个边长为M的方框范围。这样看起来比较有趣,但搜索起来可能很慢。 第三部分是三体必须在N步里不超出一个边长为M的方框范围。否则它们很快发散就不好玩了。 然后那个复选框是:是否只进行二维搜索。 搜索结束后会出现一组初始条件值,点OK就开始运行了。 3.2、在三体的基础上,搜索稳定的行星(点击右下有小球的望远镜按钮) 手工设定或者自动搜索出来的,如果喜欢的话,可以存盘,也可以导出为纯文本贴在论坛上与大家共享。压缩包里的.tbc也是偶自己用这个软件搜出来的。 四个播放按钮: 第一个播放形状的,是开始或者继续运行; 第二个暂停形状的,是暂停; 第三个短箭头,是减速运行; 第三个双箭头,是加速运行。 速度有很多档次,从减速6倍到加速运行100倍,直到加速100倍跳3125帧(相当于加倍312500倍,但每隔3125帧才显示一帧,所以看起来很不连续),每5倍为一个档次。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值