Lorenz系统

Lorenz.m

function dx=Lorenz(t,x); % 洛仑滋方程[不显含时间t的自治系统] % dx=Lorenz(t,[x;y;z;a;b;r]) % t-可以取任何数字，x,y,z-为自变量，a,b,r-为如下方程所示的参数 % eg: dx=Lorenz(0,[1;1;1;16;4;49]) % % 方程如下: % dx/dt=a*(y-x) % dy/dt=x*(r-z)-y % dz/dt=x*y-b*z % % Example(函数图象): % [T,Y]=ode45('Lorenz',[0,50],[1;1;1;16;4;49]); % plot3(Y(:,1),Y(:,2),Y(:,3)); % % Example(分岔图): % % Z=[]; % for r=linspace(1,500,350); % [T,Y]=ode45('Lorenz',[0,50],[1;1;1;10;8/3;r]); % n=length(Y); % for k=round(n/2):n % if abs(Y(k,1))<1 % Z=[Z,r+abs(Y(k,2))*i]; % end % end % end % plot(Z,'.','markersize',1) % title('Lorenz映射分岔图') % xlabel('r'),ylabel('|y| where x=0') % % Author：yujunjie % Author's email: ustb03-07@yahoo.com.cn % dx(1,1)=x(4)*(x(2)-x(1)); dx(2,1)=x(1)*(x(6)-x(3))-x(2); dx(3,1)=x(1)*x(2)-x(5)*x(3); dx(4,1)=0; dx(5,1)=0; dx(6,1)=0;