由给出的方程系数矩阵求出H(z),公式是H(z)=C(zI-A)^(-1)B+D,然后用H(z)=y(z)/f(z)写差分方程就行了。