clc
clear
close all
a=[-21,19,-20;19 -21 20;40 -40 -40];
b=[0;1;2];
c=[1 0 2];
d=0;
%繪制閉環系統的階躍響應曲線
%g=ss(a,b,c,d);
%[y,t,x]=step(g);
[y,x,t]=step(a,b,c,d);
figure(1)
plot(t,y)
title('the step responce')
xlabel('time-sec')
figure(2)
%繪制狀態變量的軌跡
plot(t,x)