状态空间方程离散化的matlab实现。

 

A=[0 1 ;1 -2];

B=[0;1];
a=0;
b=0;
T=0.01;
for i=1:2
    c=factorial(i);
    a=a+((A^(i-1))*(T^i))/c;
end
BB=a*B%零阶离散B

for i=1:2
    c=factorial(i);
    b=b+((A^i)*(T^i))/c;
end
I=eye(2);
AA=I+b%零阶离散A

以上实现为二阶矩阵ZOH离散方法。

如果要变阶次的话，以三阶举例，只需将i=1：3；I=eye(3)即可。

下图为离散后的矩阵A和B。