che1.m

来自「非线性动力学中lyapunov、duffing等方程的求解」· M 代码 · 共 29 行

clear all
pack 

global ee n alpha dx0 h gamma

ee=0.01;
n=0.6;
alpha=0.32;
dx0=0.66;
h=3.6;
gamma=2/8.944;


step=2*pi/8.944/20;

t=0:step:1000;

x01=-0.02899486874725;
x02=0.00462025726431;
x0=[x01;x02];


[T,X]=ode45('che01',t,x0);

xx=X(:,1);

N=length(t);