请教一个编程技巧的问题

gwei

现在我要求一个输入的相应，这个系统是个二阶的系统，系统为model_f=tf([1,5],[1,13,40]);输入U；问题可以简单的写成下面的代码。
u(1)=0;
time_step=0.001;
time_end=10;
t=0:time_step:time_end;
count=length(t);
ym(1)=0;
for i=2:count
   
    u(i)=ym(i-1)+1;
    [ym(i),t]=lsim(model_f,u(i),t);
   
end
??? Error using ==> rfinputs
Input data U and time vector T must have the same number of samples.
我想实现u随着ym变化，再作为输入得到下一个点的ym值。
请问我怎么实现这个功能啊？