我用了这个程序，将其中的lorenz_ext.m文件改成我所需要的Duffing模型，然后计算Lyapunov指数以判断系统状态，其值与相点曲线的状态不对应。也就是说，明明系统已经大周期运动了，max Lyapunov仍然为正。不对的地方请高手指点。