matlab中如何对状态空间方程的输出画图，以及lsim函数怎么用？

dx/du=ax+bu
y=cx+du
我设定的u是一个41的矩阵，数据均为恒定值。
a是1616的矩阵
c=[0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0]
d=[0 0 0 0]

sys = (a,b,c,d)
t=0:1e-7:0.001
y=lsim(sys,u,t')

怎样设定u矩阵可以满足lsim函数能够运行，能够输出我所设定的y