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