根据状态方程，求解系统的阶跃响应及其离散化后的阶跃响应

A=[0 1 0; 0 0 1; -6 -11 -6];
B=[1 0; 2 -1; 0 2];
C=[1 -1 0; 2 1 -1];
D=[0 0; 0 0];
sys=ss(A,B,C,D);
step(sys)

T=0.1;
sysD=c2d(sys,T);
step(sysD)

最后三行是系统离散化之后的阶跃响应。结果如图：

离散化后：

非常简单，不再赘述