MSysの优快云

本文参考百度经验如何用Mathematica计算三体问题？，进行了三体问题在三维空间的数值求解，给出了核心部分的计算代码。注意我下面这些初值设定来自一篇2013年的一篇论文（感兴趣的可以私信我），在这种设定下能构成闭合的三体周期轨道。以下这个部分是求解动力学方程的数值计算部分：{{r1xo,r1yo,r1zo,v1xo,v1yo,v1zo,r2xo,r2yo,r2zo,v2xo,v2yo,v2zo,r3xo,r3yo,r3zo,v3xo,v3yo,v3zo}}=NDSolve[{(*第一组*

微分代数（Differential Algebra）简单来说是一种利用计算机代数近似实际函数的手段，也就是说函数是以代数变量加减乘除的形式在计算机中保存的，而不是实际数值。因此，它可以作为主流符号运算软件的底层表达方式。它的核心方法是把函数进行Taylor展开，然后用nnn次多项式的系数表达出来，这其中会用到函数的微分，可能这是它称为“微分”代数的原因。在实际运算的时候，用这个“多项式近似”的结果来代替原来的函数关系，就可以有很多方便。典型例子比如，我们手头有一个函数f(x)=sinh⁡(arctan⁡

分析了受限轨迹优化问题数值求解的思路，按照直接转化法处理约束，给出了最优控制问题转化为非线性规划问题多重打靶法的描述，用MATLAB版本CasADi求解。然后又给出了一个直接配点法的计算程序。...

总结了数值求解轨迹优化常用的直接法、间接法，列举了我见到的最优控制求解软件。

以下我简单地给出解决思路，拒绝伸手党。（以后我有时间了会把整个代码发上来）主要用到的函数是plot(longitude,latitude)，其中long为经度，在±180°之间；lat为纬度，在±90°之间。我们所需要的星下点是连续的曲线，即使航天器星下点穿过了国际日界线，即±180°经线。如下图但是可惜的是，如果简单的用下面这句函数，就无法画出上面的曲线plot(Long*rad2deg, lat*rad2deg,'b-');axis([-180,180,-90,90]);grid on;可

角动量定理&刚体转动的姿态动力学，四元数形式、矩阵形式、姿态角形式的3种姿态运动学方程。