woo.m

Poincare映射图像有不同的做法

function xdot=woo(t,x,flag,w)
m=10;
k=1000;
c=15;
p0=5;
p1=10;
xdot=[x(2);-k/m*x(1)-c/m*x(2)+p0*sin(1.0005*t)+p1*sin(w*t)];