日志
- 分享 duffing方程的时间历程，相图，频谱和庞加莱映射
-
function dy=duffing(t,x) omega=1;% 定义参数 f1=x(2); f2=-omega^2*x(1)-x(1)^3; dy= ; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% clear step=0.01 y0= ; tspan= ; =ode45('duffing',tspan,y0); %%%%%%%%%%%%%%%%%%%% ...
