状态观测器与加性扰动

需要运行下面的代码，得到矩阵E，然后才能进行仿真，更具体思路的可以参考 https://blog.****.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])';