liangfun.m

用matlab 计算李雅谱诺夫指数功率谱哦

function ydot=liangfun(t,y,flag,a,b,p,w0)
ydot=[y(2);-sin(y(1))-a*y(2)+b*cos(w0*t).*sin(y(1))+p];