che01.m

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

function dx=ode45(t,x)
global ee n alpha dx0 h gamma


dx=[x(2);
    -x(1)+ee*(-2*n*x(2)-alpha*sign(x(2)-dx0)+h*cos(gamma*t))];