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。