ex754.m

来自「详细介绍了Mathlab的使用方法与技巧」· M 代码 · 共 10 行

syms z n k x1 x2
A=[0 -1;-0.4 0.3];B=[0;1];C=[0 1];                    %离散状态方程的系数矩阵
x=[x1;x2];X=A*x+B;                               %离散状态方程的状态方程
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                         %应用输出方程求出