matlab动力学分析程序详解

1.微分方程的定义对于duffing方程，先将方程写作function dy=duffing(t,x)omega=1;%定义参数f1=x(2);f2=-omega^2*x(1)-x(1)^3;dy=[