weifen.m

来自「自己刚开始接触混沌方面」· M 代码 · 共 22 行

 function weifen
[t,x]=ode45(@weifen,[0:0.001:500],[1 1 1]);




figure,plot(x(:,1),x(:,2),'k');xlabel('x1');ylabel('x2');
figure,plot(x(:,1),x(:,3),'k');xlabel('x1');ylabel('x3');
figure,plot(x(:,2),x(:,3),'k');xlabel('x2');ylabel('x3');
figure,plot3(x(:,3),x(:,1),x(:,2),'k-'); xlabel('x3');ylabel('x1');zlabel('x2'); 

function dx=weifen(t,x)
dx=zeros(3,1);


a=0.5;b=-10;c=-4;
dx(1)=a*x(1)-x(2)*x(3);
dx(2)=b*x(2)+x(1)*x(3);
dx(3)=c*x(3)+x(1)*x(2);