庞加莱(Poincaré)截面

使用带有 WhenEventNDSolve 计算庞加莱截面.

    
对于周期驱动的平面例子,可以通过在规则时间间隔采样获得庞加莱截面. 可以使用 WhenEvent 进行采样,其中带有 Sow 行为和 Reap 封闭调用 NDSolve.
求解达芬方程 ,并在驱动力 的每个周期保存样本:
可以使用 ListPlot 可视化这些参数值的奇怪吸引子:
在高维系统中,庞加莱截面实际上在相空间的一个切片上采样. 这一般需要一个依赖于解的事件. 其中一个例子是 ArnoldBeltramiChildress (ABC) 流,它被用于模拟三维欧拉方程组的层流中的混沌:
庞加莱截面可以通过把样本作为跨过 平面的解进行计算. 为了获得空间的完整形状,有必要使用多个不同的初始化条件:
将函数运用到多个初始条件: