根據狀態方程，求解系統的階躍響應及其離散化後的階躍響應

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)

最後三行是系統離散化之後的階躍響應。結果如圖：

離散化後：

非常簡單，不再贅述