需要运行下面的代码,得到矩阵E,然后才能进行仿真,更具体思路的可以参考 https://blog.csdn.net/XSTX1996/article/details/83117733
A = [1 1; 0 0];
C = [1 0];
syms e1 e2 lambda; E = [e1; e2];
T=0.01;
eigPloy = det(A-E*C - lambda*eye(2));
lambdaVal = roots([2*T^2, 2*T, 1]);
eigPloy = subs(eigPloy, lambda, lambdaVal);
[e1, e2] = solve(eigPloy, e1, e2);
E = double([e1 e2])';
1