duffing01.m

非线性动力学中lyapunov、duffing等方程的求解

function dx=ode45(t,x)
global m c k w f a


dx=[x(2);
    -c/m*x(2)+a/m*x(1)-k/m*x(1)^3+f/m*cos(w*t)];