Chen and Rossler

Chen.m

function dx=Chen(t,x); % Chen方程[不显含时间t的自治系统] % dx=Chen(t,[x;y;z;a;b;c]) % t-可以取任何数字，x,y,z-为自变量，a,b,c-为如下方程所示的参数 % eg: dx=Chen(0,[1;1;1;35;3;28]) % % 方程如下: % dx/dt=a*(y-x) % dy/dt=(c-a)*x+c*y-x*z % dz/dt=x*y-b*z % % Example(函数图象): % [T,Y]=ode45('Chen',[0,50],[1;1;1;35;3;28]); % plot3(Y(:,1),Y(:,2),Y(:,3)); % % Author's email: ustb03-07@yahoo.com.cn % dx(1,1)=x(4)*(x(2)-x(1)); dx(2,1)=(x(6)-x(4))*x(1)+x(6)*x(2)-x(1)*x(3); dx(3,1)=x(1)*x(2)-x(5)*x(3); dx(4,1)=0; dx(5,1)=0; dx(6,1)=0;

Rossler.m

function dx=Rossler(t,x); % Rossler方程[不显含时间t的自治系统] % dx=Rossler(t,[x;y;z;a;b;c]) % t-可以取任何数字，x,y,z-为自变量，a,b,c-为如下方程所示的参数 % eg: dx=Rossler(0,[1;1;1;.15;.2;10]) % % 方程如下: % dx/dt=-y-z % dy/dt=x+a*y % dz/dt=b+z*(x-c) % % Example(函数图象): % [T,Y]=ode45('Rossler',[0,400],[1;1;1;.15;.2;10]); % plot3(Y(:,1),Y(:,2),Y(:,3)); % % Author's email: ustb03-07@yahoo.com.cn % dx(1,1)=-x(2)-x(3); dx(2,1)=x(1)+x(4)*x(2); dx(3,1)=x(5)+x(3)*(x(1)-x(6)); dx(4,1)=0; dx(5,1)=0; dx(6,1)=0;