fdb.zip_fdb_模拟混沌_范德波方程

% Y. X. Li, G. R. Chen, K. S. Wallace, IEEE Trans. on Cir. Sys.��52 (2005) 204 % Author: Thomas Lee % E-mail: lixf1979@126.com % Corresponding: School of Mathematics, Physics and Software Engineering, Lanzhou Jiaotong University, Lanzhou 730070, China function xdot=lorrenz(t,x) for t=0:0.1:100; if(cos(5.3*t)>=0) p=1; else p=-1; end end xdot=[-8/3*x(1)+p*x(2)*x(3); (-x(2)+x(3))*(25-10*cos(5.3*t)); -p*x(1)*x(2)-(17.5*cos(5.3*t)+10.5)*x(2)+(13.3-14*cos(5.3*t))*x(3)]; tic; x0=[0;0;1e-10];tspan=[0,100]; [t,x]=ode45(@lorrenz,tspan,x0); plot(t,x(:,1),'k'); r=length(t)-5000:length(t); figure; plot(x(r,2),x(r,3),'k') toc;