定理：给定系统（A,B,C,D),可以求出 α（s)=Sn+a(n-1)s(n-1)+`````+a1s+a0（不好意思在S后面的字母是幂，在a后面的字母是下标），En-1=CB
En-2=CAB+an-1CB
````````
E1=CA(n-2)B+a(n-1)CA(n-2)B+````+a2CB
E0=CA(n-1)B+a(n-1)CA(n-2)B+````+a1CB
(注意 字母A后面的是幂，a后面的是下标）
则 G(s)=[E(n-1)s(n-1）+E(n-2)s(n-2)+````+ E1s+E0]/α（s) （这里E后面的是下标，s后面是幂）

莱夫勒算法：特征多项式系数ai (i=1，2，3，````，n-1)
R(N-1) =I (单位矩阵） a(n-1)= - trR(n-1)A/1
R(N-2) =R(N-1)A+a(n-1)I a(n-2)= - trR(n-2)A/2
R(N-3) =R(N-2)A+a(n-2)I a(n-3)= - trR(n-3)A/3
```````
R1 =R2A+a2I a1= - trR2A/n-1
R0 =R1A+a1I a0= - trR1A/n (这里I均为单位矩阵）
青根据以上定理和莱夫勒算法用MATLAB编写从状态空间方程转化为传递函数阵的程序，并应用一个具体的系统进行验证。
谢谢大家的帮忙！