www.gusucode.com > 直流电机拖动系统simulink仿真源码程序 > code/exam2_3.m

    clear
x0=[620 10 70];
tspan=[0,30];
[t,x]=ode45('fun2_3',tspan,x0);
plot(t,x(:,1),'k',t,x(:,2),'k--',t,x(:,3),'k:');
xlabel('time(天）  t0=0,tf=30');
ylabel('x(人)；x1(0)=620,x2(0)=10;x3(0)=70');
legend('x1','x2','x3');
grid;