模拟洛仑滋系统

Lorenz3.m

function dx=Lorenz3(t,x); % 模拟洛仑滋系统[不显含时间t的自治系统] % dx=Lorenz(t,[x;y;z;a;b;c;m]) % t-可以取任何数字，x,y,z-为自变量，a,b,r-为如下方程所示的参数 % eg: dx=Lorenz(0,[1;1;1;.9;.1;2;0]) % % 方程如下: % dx/dt=a*(y-x) % dy/dt=sign(x)*(c-z)+m % dz/dt=|x|-b*z % % % Author：yujunjie % Author's email: ustb03-07@yahoo.com.cn % dx(1,1)=x(4)*(x(2)-x(1)); dx(2,1)=sign(x(1))*(x(6)-x(3))+x(7); dx(3,1)=abs(x(1))-x(5)*x(3); dx(4,1)=0; dx(5,1)=0; dx(6,1)=0; dx(7,1)=0;

>>    [T,Y]=ode45('Lorenz3',[0,10],[1;1;1;.9;.1;2;0]);
      x=Y(length(Y),:);
      [T,Y]=ode45('Lorenz3',[0,500],x);
      plot3(Y(:,1),Y(:,2),Y(:,3));
 

>>    [T,Y]=ode45('Lorenz3',[0,10],[1;1;1;.9;.1;2;1]);
      x=Y(length(Y),:);
      [T,Y]=ode45('Lorenz3',[0,500],x);
      plot3(Y(:,1),Y(:,2),Y(:,3));
 

>>    [T,Y]=ode45('Lorenz3',[0,10],[1;1;1;.9;.1;2;-1]);
      x=Y(length(Y),:);
      [T,Y]=ode45('Lorenz3',[0,500],x);
      plot3(Y(:,1),Y(:,2),Y(:,3));