在[CNSNS 16,2845(2011)]中,我们提出了一种方法,该方法通过使用系统化的试验和误差收敛程序来扩展这种基于标准牛顿法的吸引力盆地,从而确定周期轨道。 为了做到这一点,我们将确定性算法与模拟退火算法相结合来近似周期轨道。 换句话说,这种随机方法的目标是能够以足够的精度确定初始猜测,以便将它们放入像 Newton-Raphson 算法这样的快速收敛算法的吸引力盆中。 由于该算法的潜在随机性质,它使人们能够通过多次启动该算法来确定所考虑的动力系统的几个不同的周期轨道。
本程序对应于参考文献中提出和使用的算法。 函数“SA_He_PO.m”作为一个测试用例来说明该方法的可行性和算法在寻找周期轨道方面的成功。 所考虑的示例描述了一个二自由度原子模型。