我大学时学的理论物理专业，直到现在还有很多问题不知道如何解决，比如限制性＊＊＊问题。虽然看过一些理论方面的资料，但是极其晦涩难懂，复杂的公式以及海量的计算根本就不是人力所能承受的。所以在学生时代非常的疑惑，这些理论是怎么被计算机实现的。

我知道利用计算机是可以进行求解的，我指的不是普通的有限元法或者有限差分法，而是利用计算机求解解析解的无穷级数的低频或低次项的系数，并且能估算被忽略掉的高频或高次在某一时刻的误差。

比如我先建立月球的运动模型，模型包含了月球的平动、进动、摄动、章动等运动参数，因为月球的运动是有周期的，可以把以上的运动都建模称三角级数的形式。接下来我就去根据力学定律来优化这个模型的参数，使其与力学定律的误差达到最小。在物理学中有最小作用量原理，可以通过优化最小作用量来实现对模型参数的调优。

思路是有了，接下来用PyTorch进行建模并训练模型参数应该就可以了。

我打算用Python 实现一下，加深一下自己对物理学的理解，解决我学生时代的疑惑。当然还有很多我至今没整明白的问题，比如高斯到底怎样用二乘法计算谷神星轨道的。亚当斯是怎样通过纸笔的计算就预测出天王星外面存在的海王星，并能给出海王星的轨道的。我在科普类的节目中刷到过这些，都是介绍这些人怎么怎么牛逼，但从来不对他们用到的方法展开具体介绍。我想真正搞明白还是要动手重新发明一遍吧，我不怕别人说我重复造轮子，我只怕自己对知识的理解停留在应用层而无法触及柏拉图洞穴之外。