u(k)=-10.0; end u_1=u(k); e_1=e(k); end figure(1);
plot(time,rin,'r',time,yout,'b'); xlabel('time(s)'),ylabel('rin,yout'); figure(2);
plot(time,rin-yout,'r');
xlabel('time(s)'),ylabel('error');
连续对象子函数:chap1_6f.m
function dy = PlantModel(t,y,flag,para) u=para; J=0.0067;B=0.1;
dy=zeros(2,1);