untitled.m

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

clear;
[t,z] = ode45(@fun_resp,[0,20000],[0.82;0.29;0.48;0.9;0.1;0.1;0.05;0.05;0.05;0.05]);
figure(1);
plot(t,z(:,1),'--',t,z(:,4),'r');
ylabel('x1,x2');
xlabel('t');