smtu.m

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

[t,x]=ode45('smode' ,[0 50],[0.1 0.1 0.1 0.1])
 x1=x(:,1);
 x2=x(:,2);
 x3=x(:,3);
 x4=x(:,4);
 subplot(3,1,1),plot(t,x1);ylabel('x1');xlabel('t');
 subplot(3,1,2),plot(t,x2);ylabel('x2');xlabel('t');
 subplot(3,1,3),plot(t,x3);ylabel('x3');xlabel('t');
 subplot(3,1,4),plot(t,x4);ylabel('x4');xlabel('t');
 % subplot(3,1,5),plot(x1,x2)
end