tu1.asv

本程序画出了时间序列下的各变量变化图

clear;
t_space = [0 200];

[t,x] = ode45('fun_drive',t_space,[2,1,2]);

[tt,xx] = ode45('fun_drive',t_space,[1.5,0.5,1]);
[ttt,q] = ode45('fun_resp',[0,200],[2;1;2;1.5;0.5;0.5;1.5;0.05;0.05;0.05]);
figure(1);
subplot(3,1,1);
plot(t,x(:,1),tt,xx(:,1),'--');
ylabel('x1,x2');
xlabel('t');
subplot(3,1,2);
plot(t,x(:,2),tt,xx(:,2),'--');
ylabel('y1,y2');
xlabel('t');
subplot(3,1,3);
plot(t,x(:,3),tt,xx(:,3),'--');
ylabel('z1,z2');
xlabel('t');
%plot3(x(:,1),x(:,2),x(:,3));

figure(2);
plot(ttt,q(:,1),ttt,q(:,4));
ylabel('x1,x2');
xlabel('t');
figure(3);
plot(ttt,q(:,2),ttt,q(:,5));
ylabel('y1,y2');
xlabel('t');
figure(4);
plot(ttt,q(:,3),ttt,q(:,6));
ylabel('z1,z2');
xlabel('t');
figure(5);
plot(ttt,(q(:,1)-q(:,4)),ttt,(q(:,2)-q(:,5)),ttt,(q(:,3)-q(:,6)));
ylabel('ex,ey,ez');
xlabel('t');
figure(6);
plot(ttt,q(:,7));
ylabel('ex,ey,ez');
xlabel('t');