clear all;
options=odeset('RelTol',1e-7,'AbsTol',1e-7);
[T,Y]=ode45('chua',[0,300],[0.1,0.1,0.1],options);%解微分方程???
figure(1);
plot3(Y(:,1),Y(:,2),Y(:,3),'-');
xlabel('x');ylabel('y');zlabel('z');
title('x-y-z立体相图');
figure(2);
plot(T,Y(:,1),'-');
xlabel('t/s');ylabel('x');
title('x时域波形');
figure(3);
plot(T,Y(:,1),'-');
xlabel('t/s');ylabel('x');
title('y时域波形');
figure(4);
plot(T,Y(:,3),'-');
xlabel('t/s');ylabel('z');
title('z时域波形');
figure(5);
plot(Y(:,1),Y(:,2),'-');
xlabel('x');ylabel('y');
title('x-y平面相图');
figure(6);
plot(Y(:,1),Y(:,3),'-');
xlabel('x');ylabel('z');
title('x-z平面相图');
plot(Y(:,2),Y(:,3),'-');
xlabel('y');ylabel('z');
title('y-z平面相图');