syms z n k x1 x2
A=[0 -1;-0.4 0.3];B=[0;1];C=[0 1];                    %離散狀態(tài)方程的系數(shù)矩陣
x=[x1;x2];X=A*x+B;                               %離散狀態(tài)方程的狀態(tài)方程
zIf=eye(2)*z-A;                                    % 
izIF=inv(zIf)                                      % 
zxGU=z*[1;-1]+B*(z/(z-1));                          % 
Z=izIF*zxGU;Z=factor(Z)             % [ ]，并化簡
xn=iztrans(Z);xkT=subs(xn,n,k)         %求反變換，并用 取代 求得 
ykT=C*xkT                         %應(yīng)用輸出方程求出